Apparent and real bias in numerical transfer functions in palaeolimnology

We present an evaluation of the procedure by which model prediction bias is examined in palaeolimnological transfer function inference models. We argue that most of the prediction biases commonly reported in the literature are, in fact, fallacious, and are the artificial consequence of the inappropriate manner in which residuals are traditionally examined. We show that the extent of the specious model bias is entirely predictable from first principles and is essentially determined by the strength of the predictive model. We suggest that the analysis of residuals should always be examined as a function of the model's predictions and we discuss the implications of the old and new approaches.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

Determination of Total Organic Carbon Content and δ13C in Carbonate-Rich Detrital Sediments

The determination of total organic carbon content and composition in detrital sediments requires careful removal of their carbonate minerals. In detrital sediments containing large amounts of carbonates, including dolomite, this can only be achieved by liquid acid leaching that may solubilise a significant proportion of the organic carbon. For a set of detrital sediments from the Himalayan system and the Amazon River as well as five geological reference materials, we determined the proportion of organic carbon (C_org) solubilised during acid leaching. This proportion is significant for all analysed sediments and generally tends to increase with the organic carbon content. Compared to other types of sediments analysed, clay fractions extracted from river sediments and bed sediments with very low organic carbon content have high and low proportions of acid soluble C_org respectively. In Himalayan and Amazon river sediments, the proportion of C_org solubilised during acid leaching was relatively constant with average values of 14 and 19 % respectively. Thus, it is possible to correct the C_org content for the dissolved organic carbon content measured after decarbonation. Data presented here show that C_org dissolved during liquid acid leaching must be taken into account. After careful calibration, the method presented here should, therefore, be applied to any carbonate-rich detrital sediment.     

Formation by silicate–fluoride + phosphate melt immiscibility of REE-rich globular segregations within aplite dikes

Aplite dikes intruding the Proterozoic 1.42(±?3) Ga Longs Peak-St. Vrain Silver Plume-type peraluminous granite near Jamestown, Colorado, contain F, P, and rare earth element (REE)-rich globular segregations, with 40–46% REE, 3.7–4.8 wt% P₂O₅, and 5–8 wt% F. A combination of textural features and geochemical data suggest that the aplite and REE-rich globular segregations co-existed as two co-genetic liquids prior to their crystallization, and we propose that they are formed by silicate–fluoride?+?phosphate (+?S?+?CO₂) melt immiscibility following ascent, cooling, and decompression of what was initially a single homogeneous magma that intruded the granite. The REE distribution coefficients between the silica-rich aplites and REE-rich segregations are in good agreement with experimentally determined distribution coefficients for immiscible silicate–fluoride?+?phosphate melts. Although monazite-(Ce) and uraninite U–Th–Pb microprobe ages for the segregations yield 1.420(±?25) and 1.442(±?8) Ga, respectively, thus suggesting a co-genetic relationship with their host granite, ε_Nd1.42Ga values for the granites and related granitic pegmatites range from ??3.3 to ??4.7 (average ??3.9), and differ from the values for both the aplites and REE-rich segregations, which range from ??1.0 to ??2.2 (average ??1.6). Furthermore, the granites and pegmatites have (La/Yb)_N <50 with significant negative Eu anomalies, which contrast with higher (La/Yb)_N >100 and absence of an Eu anomaly in both the aplites and segregations. These data are consistent with the aplite dikes and the REE-rich segregations they contain being co-genetic, but derived from a source different from that of the granite. The higher ε_Nd1.42Ga values for the aplites and REE-rich segregations suggest that the magma from which they separated had a more mafic and deeper, dryer and hotter source in the lower crust or upper mantle compared to the quartzo-feldspathic upper crustal source proposed for the Longs Peak-St. Vrain granite.     

The mesoscale variability in the Caribbean Sea. Part I: Simulations and characteristics with an embedded model

The variability in the Caribbean Sea is investigated using high resolution (1/15°) general circulation model experiments. For the first time in this region, simulations were carried out with a 2-way nested configuration of the NEMO primitive equation model. A coarse North Atlantic grid (1/3°) reproduces the main features of the North Atlantic and Equatorial circulation capable of influencing ocean dynamics in the Caribbean Sea. This numerical study highlights strong dynamical differences among basins and modifies the view that dynamics are homogeneous over the whole Caribbean Basin. The Caribbean mean flow is shown to organize in two intense jets flowing westward along the northern and southern boundaries of the Venezuela Basin, which merge in the center of the Colombia Basin. Diagnostics of model outputs show that width, depth and strength of baroclinic eddies increase westward from the Lesser Antilles to the Colombia Basin. The widening and strengthening to the west is consistent with altimetry data and drifter observations. Although influenced by the circulation in the Colombia Basin, the variability in the Cayman Basin (which also presents a westward growth from the Chibcha Channel) is deeper and less energetic than the variability in the Colombia/Venezuela Basins. Main frequency peaks for the mesoscale variability present a westward shift, from roughly 50 days near the Lesser Antilles to 100 days in the Cayman Basin, which is associated with growth and merging of eddies.     

Zooplankton community structure in a highly turbid environment (Charente estuary,France): Spatio-temporal patterns and environmental control

Zooplankton assemblages were studied from January 2007 to January 2008 along the salinity gradient of the Charente estuary (France). A Lagrangian survey was performed monthly at five sampling stations defined by salinity (freshwater, 0.5, 5, 15 and 25) in order to collect zooplankton and measure the main environmental parameters (concentrations of suspended particulate matter, particulate organic carbon, chlorophyll a and phaeopigments). A combination of multivariate cluster analysis, species indicator index and canonical correspondence analysis was used to relate the spatio-temporal patterns of the zooplankton assemblages with environmental drivers. The estuary was divided into three different zones by means of environmental parameters while four zooplankton assemblages were identified along the salinity gradient. The Charente estuary appeared as one of the most turbid systems in Europe, with suspended particulate matter (SPM) concentration reaching 3.5 g l⁻¹ in the Maximum Turbidity Zone (MTZ). Algal heterotrophy and microphytobenthos resuspension from the wide mudflats could be responsible for the relatively high chlorophyll a concentrations measured within this MTZ. Salinity and SPM affected significantly the spatial distribution of zooplankton species while temperature and river flow seemed to control their temporal variations. From a zooplanktonic viewpoint, the highly turbid Charente estuary seemed to match an “ecotone–ecocline” model: the succession of species assemblages along the salinity gradient matched the concept of ecocline while the MTZ, which is a stressful narrow area, could be considered as an ecotone. Although such ecoclinal characteristics seemed to be a general feature of estuarine biocenoses, the ecotone could be more system-specific and biological compartment-specific.     

Petroleum accumulation and leakage in a deeply buried carbonate reservoir,Níspero field (Mexico)

The filling history of the Níspero deeply buried Lower Cretaceous carbonate reservoir (below 4000 m) from the south part of Gulf of Mexico was studied using a combination of data from petrography, stable isotopes and fluid inclusions and compared with a one-dimensional burial model to derive timing.     

Microstructure modifications of silicates induced by the collection in aerogel: Experimental approach and comparison with Stardust results

Abstract– Particles from comet 81P/Wild 2 were captured with silica aerogel during the flyby Stardust mission. A significant part of the collection was damaged during the impact at hypervelocity in the aerogel. In this study, we conducted impact experiments into aerogel of olivine and pyroxene powder using a light‐gas gun in similar conditions as that of the comet Wild 2 particles collection. The shot samples were investigated using transmission electron microscopy to characterize their microstructure. Both olivine and pyroxene samples show evidence of thermal alteration due to friction with the aerogel. All the grains have rounded edges after collection, whereas their shape was angular in the initial shot powder set. This is probably associated with mass loss of particles. The rims of the grains are clearly melted and mixed with aerogel. The core of olivine grains is fairly well preserved, but some grains contain dislocations in glide configuration. We interpret these dislocations as generated by the thermal stresses that have emerged due to the high temperature gradients between the core and the rim of the grains. Most of the pyroxene grains have been fully melted. Their high silica concentration reflects a strong impregnation with melted aerogel. The preferential melting of pyroxene compared with olivine is due to a difference in melting temperatures of 300°. This melting point difference probably induces a bias in the measurements of the ratio olivine/pyroxene in the Wild 2 comet. The proportion of pyroxene was probably higher on Wild 2 than expected from the samples collected into aerogel.     

Planetary perturbations on Mercury’s libration in longitude

Two space missions dedicated to Mercury (MESSENGER and BepiColombo) aim at understanding its rotation and confirming the existence of a liquid core. This double challenge requires much more accurate models for the spin-orbit resonant rotation of Mercury. The purpose of this paper is to introduce planetary perturbations on Mercury’s rotation using an analytical method and to analyse the influence of the perturbations on the libration in longitude. Applying a perturbation theory based on the Lie triangle, we were able to re-introduce short periodic terms into the averaged Hamiltonian and to compute the evolution of the rotational variables. The perturbations on Mercury’s forced libration in longitude mainly come from the orbital motion of Mercury (with an amplitude around 41 arcsec that depends on the momenta of inertia). It is completed by various effects from Jupiter (11.86 and 5.93 year-periods), Venus (with a 5.66 year-period), Saturn (14.73 year-period), and the Earth (6.58 year-period). The amplitudes of the oscillations due to Jupiter and Venus are approximately 33% and 10% of those from the orbital motion of Mercury and the amplitudes of the oscillations due to Saturn and the Earth are approximately 3% and 2%. We compare the analytical results with the solution obtained from the spin-orbit numerical model SONYR.     

Cumulating processes at the crust-mantle transition zone inferred from Permian mafic-ultramafic xenoliths (Puy Beaunit,France)

Ultramafic and mafic xenoliths of magmatic origin, sampled in the Beaunit vent (northern French Massif Central), derive from the Permian (257 Ma) Beaunit layered complex (BLC) that was emplaced at the crust-mantle transition zone (∼1 GPa). These plutonic xenoliths are linked to a single fractional crystallisation process in four steps: peridotitic cumulates; websteritic cumulates; Al-rich mafic cumulates (plagioclase, pyroxenes, garnet, amphibole and spinel) and finally low-Al mafic cumulates. This sequence of cumulates can be related to the compositional evolution of hydrous Mg basaltic magma that evolved to high-Al basalt and finally to andesitic basalt. Sr and Nd isotopic compositions confirm the co-genetic character of the various magmatic xenoliths and argue for an enriched upper mantle source comparable to present mantle wedges above subduction zones. LILE, LREE and Pb enrichment are a common feature of all xenoliths and argue for an enriched sub-alkaline transitional parental magma. The existence of a Permian magma chamber at 30 km depth suggests that the low-velocity zone observed locally beneath the Moho probably does not represent an anomalous mantle but rather a sequence of mafic/ultramafic cumulates with densities close to those of mantle rocks.