Arctic freshwater ostracods from modern periglacial environments in the Lena River Delta (Siberian Arctic,Russia): geochemical applications for palaeoenvironmental reconstructions

The aim of this study is to describe ostracods from freshwater habitats in the Siberian Arctic in order to estimate the present-day relationships between the environmental setting and the geochemical properties of ostracod calcite. A special focus is on the element ratios (Mg/Ca, Sr/Ca), and the stable isotope composition (δ¹⁸O, δ¹³C), in both ambient waters and ostracod calcite. The most common species are Fabaeformiscandona pedata and F. harmsworthi with the highest frequency in all studied waters. Average partition coefficients D(Sr) of F. pedata are 0.33 ± 0.06 (1σ) in females, and 0.32 ± 0.06 (1σ) in males. A near 1:1 relationship of δ¹⁸O was found, with a mean shift of Δ_mean = 2.2‰ ± 0.5 (1σ) to heavier values in ostracod calcite of F. pedata as compared to ambient waters. The shift is not dependent on δ¹⁸O_water, and is caused by metabolic (vital) and temperature effects. Temperature-dependence is reflected in the variations of this shift. For ostracod calcite of F. pedata a vital effect as compared to inorganic calcite in equilibrium was quantified with 1.4‰. Results of this study are valuable for the palaeoenvironmental interpretation of geochemical data of fossil ostracods from permafrost deposits.     

Retrieving earthquake signature in grace gravity solutions

The GRACE satellites have been orbiting the Earth since 2002, monitoring the time variable gravity field. Some of the observed fluctuations are due to geodynamic causes, but they are often hidden in the complex signal, composed of hydrology, ocean, atmosphere, and geodynamics, the signal of geodynamic origin being usually the smallest. In addition, dealiasing residuals and noise make the separation of the signal from the different causes more difficult. We proposed a method based on the Empirical Orthogonal Function decomposition to extract the signal of physical origin, under the hypothesis that the physical signal is spatially more consistent than the noise and aliasing incomplete correction. We used synthetic geoid variations associated with earthquakes located at nearly 2000 positions at the Earth surface, based on several examples of large actual subduction events. We show that, with the present day accuracy, we can retrieve the geoid variations associated with more than 98 per cent of the earthquakes of magnitude 9 or above, around 60 per cent for magnitude 8.8, 40 per cent for magnitude 8.6 and 33 per cent for magnitude 8.3. Some events, with the right properties and location, can be detected with magnitude as low as 8. We then applied the method to the GRACE solutions, and retrieved the Hokkaido event (2003) and the Sumatra event (2004), which is in agreement with the retrieval rates mentioned here above.     

On the estimation of continuous 24-h precipitation maxima

Extreme value analysis of precipitation is of great importance for several types of engineering studies and policy decisions. For return level estimation of extreme 24-h precipitation, practitioners often use daily measurements (usually 08:00–08:00 local time) since high-frequency measurements are scarce. Annual maxima of daily series are smaller or equal to continuous 24-h precipitation maxima such that the resulting return levels may be systematically underestimated. In this paper we use a rule, derived earlier, on the conversion of the generalized extreme value (GEV) distribution of daily to 24-h maxima. We develop an estimator for the conversion exponent by combining daily maxima and high-frequency sampled 24-h maxima in one joint log-likelihood. Once the conversion exponent has been estimated, GEV-parameters of 24-h maxima can be obtained at sites where only daily data is available. The new methodology has been extended to spatial regression models.     

Seismic characterization of switching platform geometries and dominant carbonate producers (Miocene,Las Negras,Spain)

The primary goals of seismic interpretation and quantification are to understand and define reservoir architecture and the distribution of petrophysical properties. Since seismic interpretation is associated with major uncertainties, outcrop analogues are used to support and improve the resulting conceptual models. In this study, the Miocene carbonates of Cerro de la Molata (Las Negras, south‐east Spain) have been selected as an outcrop analogue. The heterogeneous carbonate rocks of the Cerro de la Molata Platform were formed by a variety of carbonate‐producing factories, resulting in various platform morphologies and a wide range of physical properties. Based on textural (thin sections) and petrophysical (porosity, density, carbonate content and acoustic properties) analyses of the sediments, eleven individual facies types were determined. The data were used to produce synthetic seismic profiles of the outcrop. The profiles demonstrate that the spatial distribution of the facies and the linked petrophysical properties are of key importance in the appearance of the synthetic seismic sections. They reveal that carbonate factory and facies‐specific reflection patterns are determined by porosity contrasts, diagenetic modifications and the input of non‐carbonate sediment. The reflectors of the seismograms created with high‐frequency wavelets are coherent with the spatial distribution of the predefined facies within the depositional sequences. The synthetic seismograms resulting from convolution with lower frequency wavelets do not show these details – the major reflectors coincide with: (i) the boundary between the volcanic basement and the overlying carbonates; (ii) the platform geometries related to changes in carbonate factories, thus sequence boundaries; and (iii) diagenetic zones. Changes in seismic response related to diagenesis, switching carbonate producers and linked platform geometries are important findings that need to be considered when interpreting seismic data sets.     

CrystalMom: a new model for the evolution of crystal size distributions in magmas with the quadrature-based method of moments

Nucleation and growth of crystals, and the resulting crystal size distribution, play a fundamental role in controlling the physical properties of magmas and consequently the dynamics of the eruptions. In the past decades, laboratory experiments demonstrated that size and shape of crystals strongly control the physical properties of magma and lava. Additionally, natural and experimental samples are usually characterized in terms of their crystal size distribution to link it with physical processes that are not directly observable, such as cooling or decompression mechanisms. In this paper, we present CrystalMoM, a new predictive model, based on the quadrature-based method of moments, developed for studying the kinetic of crystallization in volcanic systems. The quadrature-based method of moments, well established in the field of chemical engineering, represents a mesoscale modelling approach that rigorously simulates the space–time evolution of a distribution of particles, by considering its moments. The method is applied here, for the first time, for studying the equilibrium/disequilibrium crystallization in magma, modelling the temporal evolution of the moments of a crystal size distribution. The model, verified against numerical and experimental data, represents a valuable tool to infer the cooling and decompression rates from the crystal size distribution observed in natural samples.     

Stability of phlogopite in ultrapotassic kimberlite-like systems at 5.5–7.5 GPa

Hydrous K-rich kimberlite-like systems are studied experimentally at 5.5–7.5 GPa and 1200–1450?°C in terms of phase relations and conditions for formation and stability of phlogopite. The starting samples are phlogopite–carbonatite–phlogopite sandwiches and harzburgite–carbonatite mixtures consisting of Ol?+?Grt?+?Cpx?+?L (±Opx), according to the previous experimental results obtained at the same P–T parameters but in water-free systems. Carbonatite is represented by a K- and Ca-rich composition that may form at the top of a slab. In the presence of carbonatitic melt, phlogopite can partly melt in a peritectic reaction at 5.5 GPa and 1200–1350?°C, as well as at 6.3–7.0 GPa and 1200?°C: 2Phl?+?CaCO₃ (L)?Cpx?+?Ol?+?Grt?+?K₂CO₃ (L)?+?2H₂O (L). Synthesis of phlogopite at 5.5 GPa and 1200–1350?°C, with an initial mixture of H₂O-bearing harzburgite and carbonatite, demonstrates experimentally that equilibrium in this reaction can be shifted from right to left. Therefore, phlogopite can equilibrate with ultrapotassic carbonate–silicate melts in a?≥?150?°C region between 1200 and 1350?°C at 5.5 GPa. On the other hand, it can exist but cannot nucleate spontaneously and crystallize in the presence of such melts in quite a large pressure range in experiments at 6.3–7.0 GPa and 1200?°C. Thus, phlogopite can result from metasomatism of peridotite at the base of continental lithospheric mantle (CLM) by ultrapotassic carbonatite agents at depths shallower than 180–195 km, which creates a mechanism of water retaining in CLM. Kimberlite formation can begin at 5.5 GPa and 1350?°C in a phlogopite-bearing peridotite source generating a hydrous carbonate–silicate melt with 10–15 wt% SiO₂, Ca# from 45 to 60, and high K enrichment. Upon further heating to 1450?°C due to the effect of a mantle plume at the CLM base, phlogopite disappears and a kimberlite-like melt forms with SiO₂ to 20 wt% and Ca#?=?35–40.