Age and kinematics of ductile deformation in the Cerro Durazno area, NW Argentina: Significance for orogenic processes operating at the western margin of Gondwana during Ordovician – Silurian times

The Cerro Durazno Pluton belongs to a suite of Paleozoic granitoid intrusions in NW-Argentina, that are central for understanding the tectonic setting of the western margin of Gondwana in Ordovician and Silurian times. The pluton and its host rocks were tectonically overprinted by metamorphic mineral shape fabrics formed under middle greenschist-facies metamorphic conditions and associated with the nearby Agua Rosada Shear Zone. Kinematic analysis of the shear zone based on the geometric relationship between individual segments of the shear plane and principal axes of mineral fabric ellipsoids indicates reverse-sense of shear with a minor component of left-lateral displacement. This is compatible with the kinematics of other ductile deformation zones in this area, collectively forming a network, which accomplished orogen-parallel extension in addition to vertical thickening. Using the Rb–Sr isotopic system, an undeformed pegmatite dike of the Cerro Durazno Pluton was dated at 455.8 ± 3.6 Ma and mineral fabrics of the Agua Rosada Shear Zone formed at middle greenschist-facies metamorphism gave deformation ages of 437.0 ± 3.8 Ma and 428.4 ± 4.5 Ma. Thus, tectonic overprint at low metamorphic grade occurred about 20–30 Ma after terminal magmatism in the Cerro Durazno area. Our data from the Cerro Durazno area and regional considerations suggest that the western margin of Gondwana was characterized by orogen-parallel extension in addition to crustal thickening as well as episodes of magmatism and ductile deformation that varied greatly in time and space.     

Phase transitions in the system CaCO<Subscript>3</Subscript> at high P and T determined by in situ vibrational spectroscopy in diamond anvil cells and first-principles simulations

The stability of the high-pressure CaCO₃ calcite (cc)-related polymorphs was studied in experiments that were performed in conventional diamond anvil cells (DAC) at room temperature as a function of pressure up to 30 GPa as well as in internally heated diamond anvil cells (DAC-HT) at pressures and temperatures up to 20 GPa and 800 K. To probe structural changes, we used Raman and FTIR spectroscopy. For the latter, we applied conventional and synchrotron mid-infrared as well as synchrotron far-infrared radiation. Within the cc-III stability field (2.2–15 GPa at room temperature, e.g., Catalli and Williams in Phys Chem Miner 32(5–6):412–417, 2005), we observed in the Raman spectra consistently three different spectral patterns: Two patterns at pressures below and above 3.3 GPa were already described in Pippinger et al. (Phys Chem Miner 42(1):29–43, 2015) and assigned to the phase transition of cc-IIIb to cc-III at 3.3 GPa. In addition, we observed a clear change between 5 and 6 GPa that is independent of the starting material and the pressure path and time path of the experiments. This apparent change in the spectral pattern is only visible in the low-frequency range of the Raman spectra—not in the infrared spectra. Complementary electronic structure calculations confirm the existence of three distinct stability regions of cc-III-type phases at pressures up to about 15 GPa. By combining experimental and simulation data, we interpret the transition at 5–6 GPa as a re-appearance of the cc-IIIb phase. In all types of experiments, we confirmed the transition from cc-IIIb to cc-VI at about 15 GPa at room temperature. We found that temperature stabilizes cc-VI to lower pressure. The reaction cc-IIIb to cc-VI has a negative slope of ?7.0 × 10^?3 GPa K^?1. Finally, we discuss the possibility of the dense cc-VI phase being more stable than aragonite at certain pressure and temperature conditions relevant to the Earth’s mantle.     

Dust deposition and climate in the Carpathian Basin over an independently dated last glacial–interglacial cycle

Loess in the Carpathian Basin is some of the thickest and most complete in Europe. Located in the Vojvodina region of the southern Carpathian Basin the Crvenka loess-palaeosol section appears to preserve a detailed climate proxy archive of the last glacial–interglacial cycle. Central to the interpretation of the site is a detailed and independent age model. Here, the results of detailed optically stimulated luminescence (OSL) dating and elevated temperature post-IR infrared stimulated luminescence (post-IR IRSL) dating are presented. Quartz OSL ages appear accurate to about 50–60 ka, where 2D₀ values are reached, while elevated temperature post-IR IRSL yields more accurate ages below this. In line with recent results, the latter signal appears to show negligible fading rates. Two age models are developed that combine (a) OSL and post-IR IRSL ages and (b) OSL ages and ‘expected’ ages from tying unit boundaries to the marine record. If the luminescence model is regarded as accurate, differences between this and the OSL/marine age model raise questions over the accuracy of the latter, as well as the processes controlling the zeroing of luminescence dates. The luminescence based age model is then used to derive the first fully independent reconstruction of climate proxies and accumulation rates from Carpathian loess. Such reconstructions can be used to compare to other independent records without assumptions inherent in correlation-based approaches. The findings demonstrate how variable accumulation rate is at the site, and compared to other independently dated Carpathian loess records. Average values vary north–south but are of similar order throughout the basin. Accumulation rate was highest during the later part of the last glacial, but variation on millennial timescales does not always match shifts in grain-size, suggesting diverse and complex influences. Environmental reconstructions using grain-size and magnetic susceptibility show that no one atmospheric system or air mass can explain the changes in the Carpathian Basin and that millennial-scale variability can only intermittently be tied to North Atlantic Heinrich events. Expanded ice sheets during the peak last glacial, combined with other atmospheric teleconnections, may have served to develop a strong anticyclone in the region. It was likely windier during earlier parts of the last glacial, but Atlantic and Mediterranean moisture was probably less abundant than during more humid interglacials.     

The role of climate in the spread of modern humans into Europe

The spread of anatomically modern humans (AMH) into Europe occurred when shifts in the North Atlantic meridional overturning circulation triggered a series of large and abrupt climate changes during the last glacial. However, the role of climate forcing in this process has remained unclear. Here we present a last glacial record that provides insight into climate-related environmental shifts in the eastern Mediterranean region, i.e. the gateway for the colonisation of Europe by AMH. We show that the environmental impact of the Heinrich Event H5 climatic deterioration c. 48 kyr ago was as extreme as that of the glacial maximum of Marine Isotope Stage (MIS) 4 when most of Europe was deserted by Neanderthals. We argue that Heinrich H5 resulted in a similar demographic vacuum so that invasive AMH populations had the opportunity to spread into Europe and occupy large parts before the Neanderthals were able to reoccupy this territory. This spread followed the resumption of the Atlantic meridional overturning circulation at the beginning of Greenland Interstadial (GIS) 12 c. 47 kyr ago that triggered an extreme and rapid shift from desert-steppe to open woodland biomes in the gateway to Europe. We conclude that the extreme environmental impact of Heinrich H5 within a situation of competitive exclusion between two closely related hominids species shifted the balance in favour of modern humans.     

Ion tracks in apatite at high pressures: the effect of crystallographic track orientation on the elastic properties of fluorapatite under hydrostatic compression

Static elasticity measurements at high pressures were carried out on oriented fluorapatite single crystals, some of which contained oriented amorphous ion tracks (ITs) implanted with relativistic Au ions (2.2 GeV) from the UNILAC linear accelerator at GSI, Darmstadt. High-pressure experiments on irradiated and non-irradiated crystal sections were carried out in diamond-anvil high-pressure cells under hydrostatic conditions. In situ single-crystal diffraction was performed to determine the high-precision lattice parameters, simultaneously monitoring the widths of X-ray diffraction Bragg peaks. High-pressure Raman spectra were analyzed with respect to the frequency shift and widths of bands, which correspond to the Raman-active vibrational modes of the phosphate tetrahedra. Swift heavy ion irradiation was found to induce anisotropic lattice expansion and tensile strain within the host lattice dependent on the ion-track orientation. The relatively low Grüneisen parameter for the ν _1b(A _g) mode, which has been assigned to originate from the volume fraction of the amorphous tracks, and the γ(ν _1a)/γ(ν _1b) ratio reveals compressive strain on the amorphous ITs. The comparative compressibilities for the host lattice reveal approximately equivalent bulk moduli, but significantly different pressure derivatives (K _T = 88.4 ± 0.7 GPa, ∂K/∂P = 6.3 ± 0.3 for non-irradiated, K _T = 90.0 ± 1.7 GPa, ∂K/∂P = 3.8 ± 0.5 for irradiated samples). The axial compressibility moduli β ⁻¹ reveal significant differences, which correlate with the ion-track orientation [b_a ^{- 1} \beta_{a}^{ - 1}  = 240 ± 5 GPa, b_c ^{- 1} \beta_{c}^{ - 1}  = 361 ± 14 GPa, ∂( b_a ^{- 1} ) \left( {\beta_{a}^{ - 1} } \right) /∂P = 11.3 ± 1.2, ∂( b_c ^{- 1} ) \left( {\beta_{c}^{ - 1} } \right) /∂P = 11.6 ± 3.4 for irradiation ⊥(100); 246 ± 9 GPa, 364 ± 57 GPa, 9.5 ± 2.9, 14.7 ± 14.1 for irradiation ⊥(001), 230.7 ± 3.6 GPa, 373.5 ± 5.1 GPa, 19.2 ± 1.4, 20.1 ± 1.8 for no irradiation]. Line widths of XRD Bragg peaks in irradiated apatites confirm the strain of the host lattice, which appears to decrease with increasing pressure. By contrast, the bandwidths of Raman modes increase with pressure, and this is attributed to increasing strain gradients on the length scale of the short-range order. The investigations reveal considerable deviatoric stress on the [100]-oriented tracks due to the anisotropic elasticity, while the compression is uniform for the directions perpendicular to the tracks, which are aligned parallel to the c-axis. This difference might be considered to control the diffusion properties related to the annealing kinetics and its observed anisotropy, and hence to cause potential pressure effects on track-fading rates.     

Noble gas record and cosmic‐ray exposure history of the new CO3 chondrite Moss––Comparison with Lancé and other CO chondrite falls

Abstract– The Moss meteorite is the first CO chondrite fall after a time period of 70 yr and the least terrestrially contaminated member of its group. Its cosmic‐ray exposure (CRE) age (T₃ ～ 13.5 Ma; T₂₁ ～ 14.6 Ma) is distinct among CO chondrites and, within witnessed falls is the shortest after Lancé, which we have reanalyzed. Gas retention ages are approximately 3.95 × 10⁹ yr (U/Th‐He) and approximately 4.43 × 10⁹ yr (K/Ar), respectively. Trapped Ar, Kr, and Xe are present in Moss in abundances typical for CO chondrites, with “planetary” elemental and isotopic compositions. Presence of HL‐xenon from presolar diamonds is observed in the stepwise release analysis of Lancé. It may also be present in Moss, but it is difficult to ascertain in single‐step bulk analyses. It follows from our new data combined with a survey of the literature that the abundance of trapped gases in CO chondrites is not a good indicator of their petrological subtype.     

Back analysis of landslide susceptibility zonation mapping for the 2005 Kashmir earthquake: an assessment of the reliability of susceptibility zoning maps

The October 2005 earthquake triggered several thousand landslides in the Lesser Himalaya of Kashmir in northern Pakistan and India. Analyses of ASTER satellite imagery from 2001 were compared with a study undertaken in 2005; the results show the extent and nature of pre- and co-/post-seismic landsliding. Within a designated study area of ~2,250 km², the number of landslides increased from 369 in 2001 to 2,252 in October 2005. Assuming a balanced baseline landsliding frequency over the 4 years, most of the new landslides were likely triggered by the 2005 earthquake and its aftershocks. These landslides mainly happened in specific geologic formations, along faults, rivers and roads, and in shrubland/grassland and agricultural land. Preliminary results from repeat photographs from 2005 and 2006 after the snowmelt season reveal that much of the ongoing landsliding occurred along rivers and roads, and the extensive earthquake-induced fissuring. Although the susceptibility zoning success rate for 2001 was low, many of the co-/post-seismic landsliding in 2005 occurred in areas that had been defined as being potentially dangerous on the 2001 map. While most of the area in 2001 was (very) highly susceptible to future landsliding, most of the area in 2005 was only moderate to low susceptible, that is, most of the landsliding in 2005 actually occurred in the potentially dangerous areas on the 2001 map. This study supports the view that although susceptibility zoning maps represent a powerful tool in natural hazard management, caution is needed when developing and using such maps.     

Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems—the MoMaS benchmark case

Although multicomponent reactive transport modeling is gaining wider application in various geoscience fields, it continues to present significant mathematical and computational challenges. There is a need to solve and compare the solutions to complex benchmark problems, using a variety of codes, because such intercomparisons can reveal promising numerical solution approaches and increase confidence in the application of reactive transport codes. In this contribution, the results and performance of five current reactive transport codes are compared for the 1D and 2D subproblems of the so-called easy test case of the MoMaS benchmark (Carrayrou et al., Comput Geosci, 2009, this issue). This benchmark presents a simple fictitious reactive transport problem that highlights the main numerical difficulties encountered in real reactive transport problems. As a group, the codes include iterative and noniterative operator splitting and global implicit solution approaches. The 1D easy advective and 1D easy diffusive scenarios were solved using all codes, and, in general, there was a good agreement, with solution discrepancies limited to regions with rapid concentration changes. Computational demands were typically consistent with what was expected for the various solution approaches. The differences between solutions given by the three codes solving the 2D problem are more important. The very high computing effort required by the 2D problem illustrates the importance of parallel computations. The most important outcome of the benchmark exercise is that all codes are able to generate comparable results for problems of significant complexity and computational difficulty.     

A hydrogeochemical survey of Kilimanjaro (Tanzania): implications for water sources and ages

Kilimanjaro, Tanzania, the highest mountain in Africa, has undergone extensive hydrologic changes over the past century in an area where water resources are critical. A hydrochemical and isotopic synoptic sampling program in January 2006 is used to characterize hydrogeology, hydrology, and water quality of the area. Samples were collected from the summit and southern side of Kilimanjaro and the Moshi region (Tanzania). Sample sources included four glaciers, seven groundwater wells, 12 rivers, 10 springs, precipitation, and a lake. Analyses included major ion chemistry, stable isotopes of water (¹⁸O and D); in addition, seven samples were analyzed for tritium. The samples generally have good water quality with the exception of three samples with elevated fluoride concentrations (>3 mg/L) and elevated nitrate concentrations (>2.5 mg/L NO₃ as N). There is a strong elevation control on stable isotopes, with an apparent elevation effect of – 0.1 ‰ δ¹⁸O per 100 m rise in elevation (R ² = 0.79). The results, including the tritium values, show that the hydrogeologic system is comprised of both local and regional flow systems, and that regional rivers are receiving significant inflow from shallow groundwater, and at very high elevations the hydrologic system is derived from groundwater, precipitation, and glacial melt water.