Direct dating of thick‐ and thin‐skin thrusts in the Peruvian Subandean zone through apatite (U–Th)/He and fission track thermochronometry

Although the structure of the central Peruvian Subandean zone is well defined, the timing of thrust‐related exhumation and Cenozoic sedimentation remain poorly constrained. In this study, we report new apatite (U–Th)/He (AHe) and fission track (AFT) ages from thrust‐belt and foreland strata along three published balanced cross sections. AHe data from the northern, thick‐skinned domain (i.e. Shira Mountain, Otishi Cordillera and Ucayali Basin) show young AHe ages (ranging from 2.6 ± 0.2 to 13.1 ± 0.8 Ma) compared with AFT ages (ranging from 101 ± 5 to 133 ± 11 Ma). In the southern Camisea Basin, where deformation is mainly thin‐skinned, AHe and AFT ages have been both reset and show young cooling ages (3.7 ± 0.8 Ma and 8 ± 2 Ma respectively). Using low‐temperature thermochronology data and the latest fission track annealing and He diffusion codes, the thermal history of the study area has been reconstructed using inverse modelling. This history includes two steps of erosion: Early Cretaceous and late Neogene, but only Neogene sedimentation and exhumation varies in the different sectors of the study area. From a methodological point of view, large AHe data dispersion point to the need for refinement of AHe damage and annealing models. The influence of grain chemistry on damage annealing, multiple age components and the possibility of fission tracks as traps for He need further consideration. For the central Peruvian Subandes, AHe and AFT ages combined with balanced cross sections emphasize the dominant control of Paleozoic inheritance rather than climate on Cenozoic infilling and exhumation histories. Finally, our data provide the first field example of how thick‐skinned thrust‐related deformation and exhumation in the Subandes can be directly dated through AHe thermochronology.     

Structure and nuclear density distribution in the cheralite—CaTh(PO4)2: studies of its behaviour under high pressure (36?GPa)

The crystal structure of the cheralite—CaTh(PO₄)₂—has been revisited by neutron diffraction and its behaviour under high pressure investigated by X-ray diffraction up to 36?GPa. The neutron diffraction data at ambient pressure gave a more accurate determination of the Ca/Th cation position than previous XRD data, taking advantage that the neutron scattering lengths of calcium and thorium are of same order of magnitude. The nuclear density distribution was also determined using the maximum entropy method (MEM) confirming that the two cations are not located at the same position in the unit cell but are slightly displaced from one another along a specific direction in order to minimize the electrostatic repulsion with the surrounding phosphorus atoms. At high pressure, the compound did not show any phase transition or amorphization. From the evolution of the unit-cell volume as a function of the pressure, the zero-pressure bulk modulus B₀ and its pressure derivative B₀ ^′ have been determined by fitting the experimental compressibility curve to the Birch–Murnaghan equation of state. The results are B₀?=?140(2) GPa and B ₀ ^′ ?=?4.4(4) GPa.     

Micro- and nano-textural evidence of Ti(–Ca–Fe) mobility during fluid–rock interactions in carbonaceous lawsonite-bearing rocks from New Zealand

Understanding the mobility of chemical elements during fluid–rock interactions is critical to assess the geochemical evolution of a rock undergoing burial and metamorphism and, more generally, to constrain the geochemical budget of the subduction factory. In particular, determining the behavior and mobility of Ti in aqueous fluids constitutes a great challenge that is still under scrutiny. Here, we study plant fossils preserved in blueschist metasedimentary rocks from the Marybank Formation (New Zealand). Using scanning and transmission electron microscopies (SEM and TEM), we show that the carbonaceous material (CM) composing the fossils contains abundant nano-inclusions of Ti- and Fe-oxides. These nanocrystals are mainly anatase, rutile, and Fe–Ti oxides. The mineral composition observed within the fossils is significantly different from that detected in the surrounding rock matrix. We propose that Ti and Fe might have been mobilized by the alteration of a detrital Ti–Fe-rich protolith during an early diagenetic event under acidic and reducing conditions. Aqueous fluids rich in organic ligands released by the degradation of organic matter may have been involved. Moreover, using mass balance and petrological observations, we show that the contrasted mineralogy between the rock matrix and the fossil CM might be the consequence of the chemical isolation of fossil CM during the prograde path of the rock. Such an isolation results from the early formation of quartz and Fe-rich phyllosilicate layers enclosing the fossil as characterized by SEM and TEM investigations. Overall, this study shows that investigating minerals associated with CM down to the nanometer scale in metamorphic rocks can provide a precious record of early prograde geochemical conditions.     

Links between NAO fluctuations and inter-annual variability of winter-months precipitation in the Seine River watershed (north-western France)

The inter-annual to multi-decadal winter variability (DJFM) of precipitation on the Seine River watershed (France) was analysed using continuous wavelet transform analysis and compared to the winter North Atlantic Oscillation Index (NAOI). Nine weather stations were used over the 1951 to 2004 period and confirmed the homogeneity of inter-annual fluctuations for all stations but one. Wavelet coherence between SLP over the Icelandic and Azores regions and precipitation highlighted coherence for different scales of variability according to the centre of action considered. Segmentation and wavelet analysis and coherence between precipitation and NAOI over a long period of time (1873–2004) showed: i) increasing variability across the last century at most time scales, especially for NAOI; ii) the existence of change points for the mean and variance of both signals; iii) overall discontinuity of the coherence whatever the scale considered, especially between ∼1910 and ∼1955 for inter-decennial to pluri-decennial scales.     

Improving the assessment of ICESat water altimetry accuracy accounting for autocorrelation

Given that water resources are scarce and are strained by competing demands, it has become crucial to develop and improve techniques to observe the temporal and spatial variations in the inland water volume. Due to the lack of data and the heterogeneity of water level stations, remote sensing, and especially altimetry from space, appear as complementary techniques for water level monitoring. In addition to spatial resolution and sampling rates in space or time, one of the most relevant criteria for satellite altimetry on inland water is the accuracy of the elevation data. Here, the accuracy of ICESat LIDAR altimetry product is assessed over the Great Lakes in North America. The accuracy assessment method used in this paper emphasizes on autocorrelation in high temporal frequency ICESat measurements. It also considers uncertainties resulting from both in situ lake level reference data. A probabilistic upscaling process was developed. This process is based on several successive ICESat shots averaged in a spatial transect accounting for autocorrelation between successive shots. The method also applies pre-processing of the ICESat data with saturation correction of ICESat waveforms, spatial filtering to avoid measurement disturbance from the land–water transition effects on waveform saturation and data selection to avoid trends in water elevations across space. Initially this paper analyzes 237 collected ICESat transects, consistent with the available hydrometric ground stations for four of the Great Lakes. By adapting a geostatistical framework, a high frequency autocorrelation between successive shot elevation values was observed and then modeled for 45% of the 237 transects. The modeled autocorrelation was therefore used to estimate water elevations at the transect scale and the resulting uncertainty for the 117 transects without trend. This uncertainty was 8 times greater than the usual computed uncertainty, when no temporal correlation is taken into account. This temporal correlation, corresponding to approximately 11 consecutive ICESat shots, could be linked to low transmitted ICESat GLAS energy and to poor weather conditions. Assuming Gaussian uncertainties for both reference data and ICESat data upscaled at the transect scale, we derived GLAS deviations statistics by averaging the results at station and lake scales. An overall bias of −4.6 cm (underestimation) and an overall standard deviation of 11.6 cm were computed for all lakes. Results demonstrated the relevance of taking autocorrelation into account in satellite data uncertainty assesment.     

Influence of the orographic roughness of glacier valleys across the Transantarctic Mountains in an atmospheric regional model

Glacier valleys across the Transantarctic Mountains are not properly taken into account in climate models, because of their coarse resolution. Nonetheless, glacier valleys control katabatic winds in this region, and the latter are thought to affect the climate of the Ross Sea sector, frsater formation to snow mass balance. The purpose of this paper is to investigate the role of the production of turbulent kinetic energy by the subgrid-scale orography in the Transantarctic Mountains using a 20-km atmospheric regional model. A classical orographic roughness length parametrization is modified to produce either smooth or rough valleys. A one-year simulation shows that katabatic winds in the Transantarctic Mountains are strongly improved using smooth valleys rather than rough valleys. Pressure and temperature fields are affected by the representation of the orographic roughness, specifically in the Transantarctic Mountains and over the Ross Ice Shelf. A smooth representation of escarpment regions shows better agreement with automatic weather station observations than a rough representation. This work stresses the need to improve the representation of subgrid-scale orography to simulate realistic katabatic flows. This paper also provides a way of improving surface winds in an atmospheric model without increasing its resolution.     

Influence of coupling on atmosphere, sea ice and ocean regional models in the Ross Sea sector, Antarctica

Air–sea ice–ocean interactions in the Ross Sea sector form dense waters that feed the global thermohaline circulation. In this paper, we develop the new limited-area ocean–sea ice–atmosphere coupled model TANGO to simulate the Ross Sea sector. TANGO is built up by coupling the atmospheric limited-area model MAR to a regional configuration of the ocean–sea ice model NEMO. A method is then developed to identify the mechanisms by which local coupling affects the simulations. TANGO is shown to simulate realistic sea ice properties and atmospheric surface temperatures. These skills are mostly related to the skills of the stand alone atmospheric and oceanic models used to build TANGO. Nonetheless, air temperatures over ocean and winter sea ice thickness are found to be slightly improved in coupled simulations as compared to standard stand alone ones. Local atmosphere ocean feedbacks over the open ocean are found to significantly influence ocean temperature and salinity. In a stand alone ocean configuration, the dry and cold air produces an ocean cooling through sensible and latent heat loss. In a coupled configuration, the atmosphere is in turn moistened and warmed by the ocean; sensible and latent heat loss is therefore reduced as compared to the stand alone simulations. The atmosphere is found to be less sensitive to local feedbacks than the ocean. Effects of local feedbacks are increased in the coastal area because of the presence of sea ice. It is suggested that slow heat conduction within sea ice could amplify the feedbacks. These local feedbacks result in less sea ice production in polynyas in coupled mode, with a subsequent reduction in deep water formation.     

Mo isotope composition in Mo-rich high- and low-T hydrothermal systems from the Swiss Alps

We analyzed the molybdenum (Mo) isotope compositions (IC) of 59 samples from two molybdenite mineralizations (Alpjahorn and Grimsel) and from a Mo-rich hydrothermal breccia (Grimsel) from the Aar Massif, Switzerland. The formation temperature of the Late Paleozoic Mo mineralizations (300-600 °C) is much higher than that of the Pliocene breccia (100-160 °C). The Mo IC of the molybdenites varies over 1.35‰. Even in a single hand specimen it spans 0.45‰, indicating that fractionation processes during molybdenite precipitation can vary on a cm scale. The Mo IC of most molybdenites analyzed here are significantly heavier than that of the host rock (δ^98/95Mo = (0.05 ± 0.1)‰) and show a bimodal distribution centered around δ^98/95Mo ≈ 1.1‰ and 0.2‰. This result rules out single stage Rayleigh fractionation as the relevant formation mechanism and instead, redox variations are suggested to be a main factor controlling the Mo IC of the studied high-temperature Mo deposits. The range of the Mo IC in one single deposit, the Alpjahorn, overlaps with the variation range of almost all other published values for Mo IC in Mo deposits. Compared to the molybdenites, the breccia shows an even wider variation of 3.0‰ (δ^98/95Mo between −1.6‰ and +1.4‰). In contrast to the high-T molybdenite deposits, here the Mo was transported via oxidized surface waters into the breccia system, where it was reduced and precipitated. This indicates that oxidation and reduction of Mo complexes may lead to highly variable Mo IC in hydrothermal systems.     

HP–UHP metamorphism as an indicator of slab dip variations in the Alpine arc

HP/UHP and LT metamorphic units that commonly occur in the inner parts of mountain belts result from the subduction of continental and oceanic material, most often exhumed prior to continental collision. The prograde pressure–temperature history of HP–UHP rocks strongly depends on the convergence rate and on the subduction zone geometry. The maximum pressure recorded provides a proxy for the depth of shearing off and stacking of HP metamorphic nappes. A 2-D thermal model of continental subduction at lithospheric scale is used to compute the length and pressure peak of detached HP metamorphic units as a function of the slab dip angle and the convergence rate. Model results are applied to the metamorphic nappe pile of the inner Alps. A mean convergence rate of 1 cm/year during the subduction of the Briançonnais terrane is indicated by the paleogeographic reconstructions between 46 and 38 Ma. On this basis, the available petrological data and lengths of metamorphic units are used to compute the variations of the slab dip angle. The slab dip angle is shown to increase, from the northeast to the southwest, along the Alpine arc with estimated values of 20° for Suretta, 30–45° for Monte Rosa and Gran Paradiso, and 60° for Dora Maira. From Eocene to Oligocene times, the increase in slab dip angle is controlled by changes of buoyancy, due to the spatial configuration of the Valaisan trough and the incoming of crustal material within the subduction zone.     

Molluskan species richness and endemism on New Caledonian seamounts: Are they enhanced compared to adjacent slopes?

Seamounts were often considered as ‘hotspots of diversity’ and ‘centers of endemism’, but recently this opinion has been challenged. After 25 years of exploration and the work of numerous taxonomists, the Norfolk Ridge (Southwest Pacific) is probably one of the best-studied seamount chains worldwide. However, even in this intensively explored area, the richness and the geographic patterns of diversity are still poorly characterized. Among the benthic organisms, the post-mortem remains of mollusks can supplement live records to comprehensively document geographical distributions. Moreover, the accretionary growth of mollusk shells informs us about the life span of the pelagic larva. To compare diversity and level of endemism between the Norfolk Ridge seamounts and the continental slopes of New Caledonia we used species occurrence data drawn from (i) the taxonomic literature on mollusks and (ii) a raw dataset of mainly undescribed deep-sea species of the hyperdiverse Turridae. Patterns of endemism and species richness were analyzed through quantitative indices of endemism and species richness estimator metrics. To date, 403 gastropods and bivalves species have been recorded on the Norfolk Ridge seamounts. Of these, at least 38 species (∼10%) are potentially endemic to the seamounts and nearly all of 38 species have protoconchs indicating lecithotrophic larval development. Overall, our results suggest that estimates of species richness and endemism, when sampling effort is taken into account, were not significantly different between slopes and seamounts. By including in our analyses 347 undescribed morphospecies from the Norfolk Ridge, our results also demonstrate the influence of taxonomic bias on our estimates of species richness and endemism.