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.

     

The role of intracontinental deformation in supercontinent assembly: insights from the Ribeira Belt,Southeastern Brazil (Neoproterozoic West Gondwana)

In this study, we challenge the multiple collision model for the tectonic evolution of the Neoproterozoic Ribeira Belt in Southeastern Brazil. New U–Pb SHRIMP data reveal Palaeoproterozoic (2153 ± 15 Ma) and Cryogenian (783 ± 6 and 768 ± 8 Ma) granitic rocks in the Embu Domain, and detrital zircon data of metasedimentary units from the Embu and Costeiro domains suggest a coherent tectonic evolution for the whole Ribeira Belt. Rather than by multiple collisions, these data are best explained by a simpler tectonic model involving continent (craton)‐volcanic arc collisions in the Dom Feliciano and Brasilia belts that led to intracontinental crustal thickening of the adjacent thinned hinterland (Ribeira Belt) at ~640–610 Ma, followed by widespread post‐collisional magmatism and rift‐related sedimentation at ~600–540 Ma. We suggest that intracontinental orogeny is a relevant process during supercontinent assembly, as illustrated here by the evolution of significant parts of the Brasiliano orogen.     

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.     

Coordinated Ocean-ice Reference Experiments (COREs)

Coordinated Ocean-ice Reference Experiments (COREs) are presented as a tool to explore the behaviour of global ocean-ice models under forcing from a common atmospheric dataset. We highlight issues arising when designing coupled global ocean and sea ice experiments, such as difficulties formulating a consistent forcing methodology and experimental protocol. Particular focus is given to the hydrological forcing, the details of which are key to realizing simulations with stable meridional overturning circulations.The atmospheric forcing from [Large, W., Yeager, S., 2004. Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies. NCAR Technical Note: NCAR/TN-460+STR. CGD Division of the National Center for Atmospheric Research] was developed for coupled-ocean and sea ice models. We found it to be suitable for our purposes, even though its evaluation originally focussed more on the ocean than on the sea-ice. Simulations with this atmospheric forcing are presented from seven global ocean-ice models using the CORE-I design (repeating annual cycle of atmospheric forcing for 500 years). These simulations test the hypothesis that global ocean-ice models run under the same atmospheric state produce qualitatively similar simulations. The validity of this hypothesis is shown to depend on the chosen diagnostic. The CORE simulations provide feedback to the fidelity of the atmospheric forcing and model configuration, with identification of biases promoting avenues for forcing dataset and/or model development.     

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.