Mineral equilibria constraints on open‐system melting in metamafic compositions

The recent development of activity–composition relations for mineral and melt phases in high‐grade metamafic rocks allows mineral equilibria tools to be used to further aid our understanding of partial melting and the mineralogical consequences of melt segregation in these rocks. We show that bulk compositional data from natural amphibolites cover a wide compositional range, with particular variability in the content and ratios of Ca, Na and K indicating that low‐grade metasomatic alteration can substantially alter the igneous protolith chemistry and potentially affect the volume and composition of melt generated. Mineral equilibria calculations for five samples that span the compositional variability in our data set indicate that melting occurs primarily via the fluid‐absent breakdown of amphibole+quartz to produce a pressure‐sensitive peritectic assemblage of augite, orthopyroxene and/or garnet. The introduction of orthopyroxene at the onset of the amphibolite‐to‐granulite‐facies transition at lower pressure results in an increased rate of melt production until quartz is typically exhausted, and this is similarly seen for the introduction of garnet at higher pressure. Calculated melt compositions are dependent on the protolith composition, but initial solidus melting and biotite breakdown produce 1–3 mol.% of K‐rich granitic melts. As hornblende melting proceeds, 15–20 vol.% of either more granodioritic‐to‐tonalitic or granodioritic‐to‐trondhjemitic melt is produced. Once quartz is exhausted, intermediate to mafic melt compositions are produced at ultrahigh‐temperature conditions. Quartz‐rich lithologies with high Ca coupled to low Na and K are the most fertile under orogenic conditions, yielding up to 25 mol.% of sub‐alkalic granitic melt by 850°C. Such rocks did not experience significant subsolidus alteration. Altered compositions with low Ca and elevated Na and K are not as fertile, yielding less than 15 mol.% of alkalic granitic melt by 850°C. These melt volumes are enough to be segregated, and can make a contribution to granite magmatism and intracrustal differentiation that should not be overlooked.     

Internal vein texture and vein evolution of the epithermal Shila-Paula district,southern Peru

The epithermal Shila-Paula Au–Ag district is characterized by numerous veins hosted in Tertiary volcanic rocks of the Western Cordillera (southern Peru). Field studies of the ore bodies reveal a systematic association of a main E–W vein with secondary N55–60°W veins—two directions that are also reflected by the orientation of fluid-inclusion planes in quartz crystals of the host rock. In areas where this pattern is not recognized, such as the Apacheta sector, vein emplacement seems to have been guided by regional N40°E and N40°W fractures. Two main vein-filling stages are identified. stage 1 is a quartz–adularia–pyrite–galena–sphalerite–chalcopyrite–electrum–Mn silicate–carbonate assemblage that fills the main E–W veins. stage 2, which contains most of the precious-metal mineralization, is divided into pre-bonanza and bonanza substages. The pre-bonanza substage consists of a quartz–adularia–carbonate assemblage that is observed within the secondary N45–60°W veins, in veinlets that cut the stage 1 assemblage, and in final open-space fillings. The two latter structures are finally filled by the bonanza substage characterized by a Fe-poor sphalerite–chalcopyrite–pyrite–galena–tennantite–tetrahedrite–polybasite–pearceite–electrum assemblage. The ore in the main veins is systematically brecciated, whereas the ore in the secondary veins and geodes is characteristic of open-space crystallization. Microthermometric measurements on sphalerite from both stages and on quartz and calcite from stage 2 indicate a salinity range of 0 to 15.5 wt% NaCl equivalent and homogenization temperatures bracketed between 200 and 330°C. Secondary CO₂-, N₂- and H₂S-bearing fluid inclusions are also identified. The age of vein emplacement, based on ⁴⁰Ar/³⁹Ar ages obtained on adularia of different veins, is estimated at around 11 Ma, with some overlap between adularia of stage 1 (11.4±0.4 Ma) and of stage 2 (10.8±0.3 Ma). A three-phase tectonic model has been constructed to explain the vein formation. Phase 1 corresponds to the assumed development of E–W sinistral shear zones and associated N60°W cleavages under the effects of a NE–SW shortening direction that is recognized at Andean scale. These structures contain the stage 1 ore assemblage that was brecciated during ongoing deformation. Phase 2 is a reactivation of earlier structures under a NW–SE shortening direction that allowed the reopening of the preexisting schistosity and the formation of scarce N50°E-striking S2-cleavage planes filled by the stage 2 pre-bonanza minerals. Phase 3 coincides with the bonanza ore emplacement in the secondary N45–60°W veins and also in open-space in the core of the main E–W veins. Our combined tectonic, textural, mineralogical, fluid-inclusion, and geochronological study presents a complete model of vein formation in which the reactivation of previously formed tectonic structures plays a significant role in ore formation.     

Stochastic representation for anisotropic permeability tensor random fields

In this paper, we introduce a novel stochastic model for the permeability tensor associated with stationary random porous media. In the light of recent works on mesoscale modeling of permeability, we first discuss the physical interpretation of the permeability tensor randomness. Subsequently, we propose a nonparametric prior probabilistic model for non‐Gaussian permeability tensor random fields, making use of the information theory and a maximum entropy procedure, and provide a physical interpretation of the model parameters. Finally, we demonstrate the capability of the considered class of random fields to generate higher levels of statistical fluctuations for selected stochastic principal permeabilities. This unique flexibility offered by the parameterization of the model opens up many new possibilities for both forward simulations (e.g. for uncertainty propagation in predictive simulations) and stochastic inverse problem solving. Copyright © 2011 John Wiley & Sons, Ltd.     

Progressive formation of modern drumlins at Múlajökull,Iceland: stratigraphical and morphological evidence

The drumlin field at Múlajökull, Iceland, is considered to be an active field in that partly and fully ice‐covered drumlins are being shaped by the current glacier regime. We test the hypothesis that the drumlins form by a combination of erosion and deposition during successive surge cycles. We mapped and measured 143 drumlins and studied their stratigraphy in four exposures. All exposures reveal several till units where the youngest till commonly truncates older tills on the drumlin flanks and proximal slope. Drumlins inside a 1992 moraine are relatively long and narrow whereas drumlins outside the moraine are wider and shorter. A conceptual model suggests that radial crevasses create spatial heterogeneity in normal stress on the bed so that deposition is favoured beneath crevasses and erosion in adjacent areas. Consequently, the crevasse pattern of the glacier controls the location of proto‐drumlins. A feedback mechanism leads to continued crevassing and increased sedimentation at the location of the proto‐drumlins. The drumlin relief and elongation ratio increases as the glacier erodes the sides and drapes a new till over the landform through successive surges. Our observations of this only known active drumlin field may have implications for the formation and morphological evolution of Pleistocene drumlin fields with similar composition, and our model may be tested on modern drumlins that may become exposed upon future ice retreat.     

Polyphase movement on the Lavanttal Fault Zone (Eastern Alps): reconciling the evidence from different geochronological indicators

The Lavanttal Fault Zone (LFZ) is generally considered to be related to Miocene orogen-parallel escape tectonics in the Eastern Alps. By applying thermochronological methods with retention temperatures ranging from ~450 to ~40°C we have investigated the thermochronological evolution of the LFZ and the adjacent Koralm Complex (Eastern Alps). ⁴⁰Ar/³⁹Ar dating on white mica and zircon fission track (ZFT) thermochronology were carried out on host rocks (HRs) and fault-related rocks (cataclasites and fault gouges) directly adjacent to the unfaulted protolith. These data are interpreted together with recently published apatite fission track (AFT) and apatite (U-Th)/He ages. Sample material was taken from three drill cores transecting the LFZ. Ar release spectra in cataclastic shear zones partly show strongly rejuvenated incremental ages, indicating lattice distortion during cataclastic shearing or hydrothermal alteration. Integrated plateau ages from fault rocks (~76 Ma) are in parts slightly younger than plateau ages from HRs (>80 Ma). Incremental ages from fault rock samples are in part highly reduced (~43 Ma). ZFT ages within fault gouges (~65 Ma) are slightly reduced compared to the ages from HRs, and fission tracks show reduced lengths. Combining these results with AFT and apatite (U-Th)/He ages from fault rocks of the same fault zone allows the recognition of distinct faulting events along the LFZ from Miocene to Pliocene times. Contemporaneous with this faulting, the Koralm Complex experienced accelerated cooling in Late Miocene times. Late-Cretaceous to Palaeogene movement on the LFZ cannot be clearly proven. ⁴⁰Ar/³⁹Ar muscovite and ZFT ages were probably partly thermally affected along the LFZ during Miocene times.     

Microdiamond grade as a regionalised variable – some basic requirements for successful local microdiamond resource estimation of kimberlites

Microdiamonds offer several advantages as a resource estimation tool, such as access to deeper parts of a deposit which may be beyond the reach of large diameter drilling (LDD) techniques, the recovery of the total diamond content in the kimberlite, and a cost benefit due to the cheaper treatment cost compared to large diameter samples. In this paper we take the first step towards local estimation by showing that micro-diamond samples can be treated as a regionalised variable suitable for use in geostatistical applications and we show examples of such output. Examples of microdiamond variograms are presented, the variance-support relationship for microdiamonds is demonstrated and consistency of the diamond size frequency distribution (SFD) is shown with the aid of real datasets. The focus therefore is on why local microdiamond estimation should be possible, not how to generate such estimates. Data from our case studies and examples demonstrate a positive correlation between micro- and macrodiamond sample grades as well as block estimates. This relationship can be demonstrated repeatedly across multiple mining operations. The smaller sample support size for microdiamond samples is a key difference between micro- and macrodiamond estimates and this aspect must be taken into account during the estimation process. We discuss three methods which can be used to validate or reconcile the estimates against macrodiamond data, either as estimates or in the form of production grades: (i) reconcilliation using production data, (ii) by comparing LDD-based grade estimates against microdiamond-based estimates and (iii) using simulation techniques.

     

Late orogenic rebound and oblique Alpine convergence: New constraints from subsidence analysis of the Austrian Molasse basin

Subsidence analysis of 16 wells in the Austrian Molasse basin documents major spatial and temporal changes in tectonic subsidence as well as a late-stage surface uplift. The timing of the main phase of tectonic subsidence shifted from early Oligocene in the western part of the peripheral foreland to the early Miocene in the eastern part. These temporal and spatial changes in tectonic subsidence reflect a change from oblique dextral to sinistral convergence between the Alpine nappe stack and its foreland. The main phase of sediment accumulation was delayed to the early Miocene and led to the infill of the basin and a major second, sediment-load driven phase of basement subsidence. Sediment accumulation rates in the basin reflect the build-up of topography in the Alpine mountain chain. Since approximately 6 Ma a pronounced regional uplift of the entire Molasse basin has taken place, marking the transition from lateral extrusion to orthogonal contraction within the Alpine system and deep-seated changes in geodynamic boundary conditions, possibly due to delamination of previously thickened lithosphere. Surface uplift is contemporaneous with similar processes in extra-Alpine Central Europe, where it is interpreted to reflect intra-plate stress changes.     

Multiproxy Analyses of Stratigraphy and Palaeoenvironment of the Late Palaeolithic Grabow Floodplain Site,Northern Germany

Changing river courses and fluctuations of the water table were some of the most fundamental environmental changes that humans faced during the Late Glacial, particularly as these changes affected areas intensively used for settlement and resource exploitation. Unfortunately, only a few stratigraphies have been documented in the North European plain that show the interaction between river development, vegetation history, and occupation by Late Palaeolithic humans. Here, we present the results of detailed stratigraphical studies (pedology, archaeology, chrono‐, tephra‐, and palynostratigraphy) at the Federmesser site Grabow 15 located in the broad Elbe River valley. The research aimed to produce a model of site formation based on a multiproxy approach, relating the local evidence to the palaeoenvironmental and settlement history of the wider region. After deposition of fluvial sands during the Late Pleniglacial in a braided setting, the river course developed locally toward a meandering system at the transition from the Older Dryas to the Allerød, while periodic flooding led to the deposition of floodplain sediments during the early Allerød. The floodplain was settled by people of the earliest “Federmessergruppen,” who are believed to have chosen this open floodplain area along the river for collecting and processing amber of local origin. Their artifacts became embedded in the aggrading floodplain sediments. In the late Allerød, floodplain sedimentation ceased and a Fluvisol‐type soil developed, indicating a trend toward geomorphic stability. The Fluvisol was then covered by silty floodplain sediments due to a rising water level during the late Younger Dryas resulting in the cessation of human occupation in the area. Subsequent organic‐rich Late Glacial/Holocene sediments preserved the settlement remains to the present.     

Patterns of decadal-scale Arctic warming events in simulated climate

Pronounced positive decadal-scale temperature anomalies occurred in the Arctic region in the first half of the twentieth century, an episode known as the early twentieth century warming (ETCW). Analyzing a 3,000-year unperturbed climate simulation performed with the Max Planck Institute Earth System Model, we demonstrate that internal variability of the Northern Hemisphere climate system is sufficient to reproduce warm events matching the observed ETCW. We perform a superposed epoch analysis on simulated data and identify 26 Arctic warming episodes compatible with the ETCW. The simulated events reproduce, in their ensemble average, magnitude as well as spatial and temporal extent of the observed ETCW. In individual realizations, the ETCW-like events indicate that different patterns of internally generated decadal Arctic warming are possible, including pan-Arctic warming events. We investigate the dynamics that typically lead to the simulated warming events: positive oceanic heat transport anomalies that form in the North Atlantic initialize the warming events and trigger an ocean-ice-albedo feedback in the Barents Sea region. The consequent reduction in sea-ice extent leads to enhanced multi-year surface warming through strengthened ocean heat release to the atmosphere. The oceanic heat transport anomalies reduce to pre-event levels around the year of the maximum warming. However, the warming events typically lasts for another 5–7 years until the sea-ice extent recovers to pre-event conditions.     

Analytical 2D GNSS PVT solutions from a hyperbolic positioning approach

GNSS receivers estimate 3D antenna position and receiver clock bias when at least four satellites are tracked. If only three satellites are available, a 2D antenna position solution is still possible. We derive an almost exact algorithm for the determination of two possible antenna positions and the corresponding receiver clock biases based on pseudorange measurements to three GNSS satellites and a height measurement. The two ambiguous solutions exactly reflect the same height measurement. One of the solutions can be eliminated if some prior knowledge of the user position, for example, near the Earth, is available. In general, a less accurate height measurement gives a less accurate 2D GNSS solution, and vice versa. The determination of the receiver antenna position is based upon the intersection of two confocal hyperboloid sheets and the ellipsoid, resulting in a hyperbola along which the user is located. The algorithm is verified by numerical computations.