Holocene survival of the wild horse in Europe: a matter of open landscape?

The wild horse Equus ferus was one of the most frequent species of the Late Pleistocene large ungulate fauna in Eurasia and played an important role in the subsistence of human groups, especially at the end the Late Glacial. It is frequently assumed that E. ferus became extinct in Europe at the beginning of the Holocene because of the development of woodlands and loss of open habitats. Because of its preference for open habitats and in spite of its adaptability, the appearance or disappearance of the wild horse could therefore be a suitable palaeoecological indicator for the opening of the Holocene primeval woodlands. We revised the dating and reliability of the subfossil record and dated several bones by atomic mass spectrometry ¹⁴C dating. From the beginning of the Holocene (9600 cal a BC) to the end of the Atlantic Period (3750 cal a BC) there are 207 archaeological sites with wild horse records available in Europe. E. ferus survived the Pleistocene Holocene transition in Europe, but the spatiotemporal dynamics of populations fluctuated remarkably in the early and middle Holocene. Small and sparse populations increasingly became extinct during the early Holocene, until between 7100 and 5500 cal a BC the wild horse was almost absent in central parts of the European Lowlands. Particular conditions in natural open patches in the canopy forests, chalklands and floodplains may have maintained the local survival of the horse in some regions of the Lowlands, however. In the Late Atlantic, between 5500 and 3750 cal a BC the range of the wild horse was again extended. It re‐immigrated into central and western Europe, probably as a consequence of increasing landscape opening by Neolithic peoples. The data presented here may be a valuable part of the debate on the degree of openness of the early and middle Holocene landscape. Copyright © 2011 John Wiley & Sons, Ltd.     

Seasonal temperatures for the past ∼400 years reconstructed from diatom and chironomid assemblages in a high-altitude lake (Lej da la Tscheppa,Switzerland)

We analysed a 42 cm long sediment record from Lej da la Tscheppa, a high-altitude lake (2,616 m a.s.l.) in the Upper Engadine valley (Switzerland) for subfossil diatoms, chironomids and pollen. The chronology of the top 21 cm of the record was established using ²¹⁰Pb analysis using a constant-rate-of-supply model, and validated with ¹³⁷Cs measurements and the content of spheroidal carbonaceous particles. A tentative chronology for the lower part of the core was obtained through extrapolation of the sedimentation rates in the uppermost part of the record. Pollen assemblages in the record reflect regional changes in forestation and land-use patterns in the Upper Engadine valley and show no evidence of significant local human activity in the lake’s catchment. Diatom assemblages record a distinct increase in planktonic taxa since the early 19th century, suggesting a decrease in the duration of ice-cover. In contrast, chironomid assemblages remained stable during a large part of the record. We applied an established chironomid-based July air temperature transfer function and a newly developed diatom-based spring air temperature transfer function to reconstruct past seasonal air temperature changes at Lej da la Tscheppa. The reconstructions indicate a diatom-inferred warming trend in spring temperatures during the past ca. 400 years, whereas chironomid-inferred summer temperatures suggest a slight cooling trend. These biota-based reconstructions are in good agreement with the centennial-scale temperature trend in an independent reconstruction of regional temperatures in the Upper Engadine region based on instrumental records and documentary proxy evidence from the Alps. Our results suggest that, in high-altitude lakes, independent chironomid- and diatom-based seasonal temperature reconstruction is possible and can be successfully used to track seasonal temperature trends.     

Evaluation of magmatic processes for the products of the Nevado del Ruiz Volcano, Colombia from geochemical and petrological data

Trace-element concentrations of the eruptive products from the Nevado del Ruiz volcano indicate that at least for the last 1 Ma, the major controlling factor in the evolution of the magma has been a simple fractional crystallization.The Zr/Hf ratios suggest that the source material is mantle-like in origin and whilst subduction-related contamination of the source material has taken place, there is no evidence that assimilation of continental crust has affected the series.No convincing geochemical evidence for magma mixing was found, and data on banded pumices of Pleistocene, Historic and the 1985 eruptions strongly support a crystal fractionation model to explain the compositional variation of their glasses. Therefore, if mixing has taken place it can only be a mixing of layers within a stratified or partially stratified magma chamber.Modeling of fractional crystallization using a stepwise program has been undertaken and least-squares approximations with small residuals are commensurate with analyzed samples. Crystal fractionation models with plagioclase dominating two-pyroxenes and iron oxides give good fits for all calculations. We conclude that simple crystal fractionation is the process most consistent with major- and trace-element variations in the Ruiz series. This closed-system model requires at least 77% crystal fractionation of a basaltic parent to generate the observed compositions of dacites.     

Benchmarking spatial joins a la carte

Spatial joins are join operations that involve spatial data types and operators. Spatial access methods are often used to speed up the computation of spatial joins. This paper addresses the issue of benchmarking spatial join operations. For this purpose, we first present a WWW-based benchmark generator to produce sets of rectangles. Using a Web browser, experimenters can specify the number of rectangles in a sample, as well as the statistical distributions of their sizes, shapes, and locations. Second, using the generator and a well-defined set of statistical models we define several tests to compare the performance of three spatial join algorithms: nested loop, scan-and-index, and synchronized tree traversal. We also added two real-life data sets from the Sequoia 2000 storage benchmark. Our results show that the relative performance of the different techniques mainly depends on the selectivity factor of the join predicate. All of the statistical models and algorithms are available on the Web, which allows for easy verification and modification of our experiments.     

Mineralogical and chemical evolution of the Ernest Henry Fe oxide–Cu–Au ore system, Cloncurry district, northwest Queensland, Australia

The Ernest Henry Cu–Au deposit was formed within a zoned, post-peak metamorphic hydrothermal system that overprinted metamorphosed dacite, andesite and diorite (ca 1740–1660 Ma). The Ernest Henry hydrothermal system was formed by two cycles of sodic and potassic alteration where biotite–magnetite alteration produced in the first cycle formed ca 1514±24 Ma, whereas paragenetically later Na–Ca veining formed ca 1529 +11/−8 Ma. These new U–Pb_titanite age dates support textural evidence for incursion of hydrothermal fluids after the metamorphic peak, and overlap with earlier estimates for the timing of Cu–Au mineralization (ca 1540–1500 Ma). A distal to proximal potassic alteration zone correlates with a large (up to 1.5 km) K–Fe–Mn–Ba enriched alteration zone that overprints earlier sodic alteration. Mass balance analysis indicates that K–Fe–Mn–Ba alteration—largely produced during pre-ore biotite- and magnetite-rich alteration—is associated with K–Rb–Cl–Ba–Fe–Mn and As enrichment and Na, Ca and Sr depletion. The aforementioned chemical exchange almost precisely counterbalances the mass changes associated with regional Na–Ca alteration. This initial transition from sodic to potassic alteration may have been formed during the evolution of a single fluid that evolved via alkali exchange during progressive fluid-rock interaction. Cu–Au ore, dominated by co-precipitated magnetite, minor specular hematite, and chalcopyrite as breccia matrix, forms a pipe-like body at the core of a proximal alteration zone dominated by K-feldspar alteration. Both the core and K-feldspar alteration overprint Na–Ca alteration and biotite–magnetite (K–Fe) alteration. Ore was associated with the concentration of a diverse range of elements (e.g. Cu, Au, Fe, Mo, U, Sb, W, Sn, Bi, Ag, F, REE, K, S, As, Co, Ba and Ca). Mineralization also involved the deposition of significant barite, K(–Ba)–feldspar, calcite, fluorite and complexly zoned pyrite. The complexly zoned pyrite and variable K–(Ba)–feldspar versus barite associations are interpreted to indicate fluctuating sulphur and/or barium supply. Together with the alteration zonation geochemistry and overprinting criteria, these data are interpreted to indicate that Cu–Au mineralization occurred as a result of fluid mixing during dilation and brecciation, in the location of the most intense initial potassic alteration. A link between early alteration (Na–Ca and K–Fe) and the later K-feldspathization and the Cu–Au ore is possible. However, the ore-related enrichments in particular elements (especially Ba, Mn, As, Mo, Ag, U, Sb and Bi) are so extreme compared with earlier alteration that another fluid, possibly magmatic in origin, contributed the diverse element suite geochemically independently of the earlier stages. Structural focussing of successive stages produced the distinctive alteration zoning, providing a basis both for exploration for similar deposits, and for an understanding of ore genesis.     

High-T, low-P metamorphism in the Palaeoproterozoic Halls Creek Orogen, northern Australia: the middle crustal response to a mantle-related transient thermal pulse

High‐T, low‐P metamorphic rocks of the Palaeoproterozoic central Halls Creek Orogen in northern Australia are characterised by low radiogenic heat production, high upper crustal thermal gradients (locally exceeding 40 °C km^?1) sustained for over 30 Myr, and a large number of layered mafic‐ultramafic intrusions with mantle‐related geochemical signatures. In order to account for this combination of geological and thermal characteristics, we model the middle crustal response to a transient mantle‐related heat pulse resulting from a temporary reduction in the thickness of the mantle lithosphere. This mechanism has the potential to raise mid‐crustal temperatures by 150–400 °C within 10–20 Myr following initiation of the mantle temperature anomaly, via conductive dissipation through the crust. The magnitude and timing of maximum temperatures attained depend strongly on the proximity, duration and lateral extent of the thermal anomaly in the mantle lithosphere, and decrease sharply in response to anomalies that are seated deeper than 50–60 km, maintained for <5 Myr in duration and/or have half‐widths <100 km. Maximum temperatures are also intimately linked to the thermal properties of the model crust, primarily due to their influence on the steady‐state (background) thermal gradient. The amplitudes of temperature increases in the crust are principally a function of depth, and are broadly independent of crustal thermal parameters. Mid‐crustal felsic and mafic plutonism is a predictable consequence of perturbed thermal regimes in the mantle and the lowermost crust, and the advection of voluminous magmas has the potential to raise temperatures in the middle crust very quickly. Although pluton‐related thermal signatures significantly dissipate within <10 Myr (even for very large, high‐temperature intrusive bodies), the interaction of pluton‐ and mantle‐related thermal effects has the potential to maintain host rock temperatures in excess of 400–450 °C for up to 30 Myr in some parts of the mid‐crust. The numerical models presented here support the notion that transient mantle‐related heat sources have the capacity to contribute significantly to the thermal budget of metamorphism in high‐T, low‐P metamorphic belts, especially in those characterised by low surface heat flow, very high peak metamorphic geothermal gradients and abundant mafic intrusions.     

Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

Total nitrogen content of deep sea basalts

Various models have been suggested concerning the origin and evolution of the earth's atmosphere. An estimate of the nitrogen content of the mantle could further constrain atmospheric models. Total nitrogen content was determined by thermal neutron activation analysis via ¹⁴N(n,p)¹⁴C. The ¹⁴C was converted to carbon dioxide and counted in miniature low level proportional counters. The total nitrogen content of U.S.G.S. standards BCR-1 and G-2 as determined by different laboratories is variable, probably due to atmospheric adsorption by the finely ground samples. Total nitrogen content was determined in deep sea basalt glasses from three regions: East Pacific Rise ($\text{15 ± 4, 18 ± 4}$, and $\text{7 ±}\text{ppm}\text{2 N}$), Mid-Atlantic Rift (FAMOUS Region:$\text{22 ± 5, 18 ± 3}$, and $\text{10 ± 2}\text{ppm}\text{N}$) and the Juan de Fuca Ridge ($\text{17 ± 4}\text{ppm}\text{N}$). Matrix material from the same samples as the glasses was available from the East Pacific Rise ($\text{37 ± 6, 26 ± 4}$, and $\text{34 ± 6}\text{ppm}\text{N}$) and the Mid-Atlantic Rift ($\text{39 ± 4}\text{ppm}\text{N}$) which are about 50 to 100% greater than the associated glasses. The increased matrix abundance may be due to incorporation of chemically bound nitrogen from sea water rather than dissolved molecular nitrogen. The nitrogen content of the FAMOUS samples are inconsistent with the model of Langmuir et al. (1977) for petrogenesis based on trace element data. Factors which can affect the observed nitrogen content in the basalts and the interpretation in terms of the mantle nitrogen abundance are discussed (e.g. partial melting and degassing of the basalts). A lower limit of about 2 ppm N in the mantle can be estimated.     

Rapid orbit determination of LEO satellites using IGS clock and ephemeris products

Different types of GPS clock and orbit data provided by the International GPS Service (IGS) have been used to assess the accuracy of rapid orbit determination for satellites in low Earth orbit (LEO) using spaceborne GPS measurements. To avoid the need for reference measurements from ground-based reference receivers, the analysis is based on an undifferenced processing of GPS code and carrier-phase measurements. Special attention is therefore given to the quality of GPS clock data that directly affects the resulting orbit determination accuracy. Interpolation of clock data from the available 15 min grid points is identified as a limiting factor in the use of IGS ultra-rapid ephemerides. Despite this restriction, a 10-cm orbit determination accuracy can be obtained with these products data as demonstrated for the GRACE-B spacecraft during selected data arcs between 2002 and 2004. This performance may be compared with a 5-cm orbit determination accuracy achievable with IGS rapid and final products using 5 min clock samples. For improved accuracy, high-rate (30 s) clock solutions are recommended that are presently only available from individual IGS centers. Likewise, a reduced latency and more frequent updates of IGS ultra-rapid ephemerides are desirable to meet the requirements of upcoming satellite missions for near real-time and precise orbit determination.