Numerical simulation of percussive drilling

This paper presents a novel dynamical model to analyze the long‐term response of a percussive drilling system. This departs from existing approaches that usually consider a single activation and bit/rock interaction cycle for the analysis of the process performance. The proposed model integrates the axial dynamics of an elastic piston and an elastic drill bit, a motion‐dependent pressure law to drive the piston, and a generalized bit/rock interaction law representative of the dynamic indentation taking place at the bit/rock interface. It applies to down‐the‐hole percussive drilling as well as top‐hole, with minor modifications. The model does not account for the angular motion or the hole cleaning, however. The model is first formulated mathematically; then, a finite‐dimensional approximation is proposed for computations. Numerical analyses of the model response, for a low‐size down‐the‐hole percussive system, follow. The period‐1 stationary response for the reference configuration is studied in detail, and parametric analyses assessing the influence on the rate of penetration of the bit/rock interaction parameters, the feed force, and the percussive activation parameters are conducted. These analyses reveal that the multiscale nature of the process is well captured by the model and recover expected trends for the influence of the parameters. They also suggest that a significant increase of the penetration rate can be achieved by increasing the percussive frequency. Copyright © 2015 John Wiley & Sons, Ltd.     

Kimberlite petrogenesis: Insights from clinopyroxene-melt partitioning experiments at 6 GPa in the CaO-MgO-Al2O3-SiO2-CO2 system

In this experimental study, we examine the mineral-melt partitioning of major and trace elements between clinopyroxene and CO₂-rich kimberlitic melts at a pressure of 6 GPa and temperatures of 1410°C and 1430°C. The melts produced contain ∼ 28 wt% dissolved CO₂, and are saturated with olivine and clinopyroxene. To assess the effects of temperature, crystal and melt compositions on trace element partitioning, experiments were performed in the model CaO-MgO-Al₂O₃-SiO₂-CO₂ system. Our results reveal that all the elements studied, except Al, Mg, Si, and Ga, are incompatible in clinopyroxene. Partition coefficients show a considerable range in magnitude, from ∼ 10⁻³ for D_U and D_Ba to ∼ 2.5 for D_Si. The two experimental runs show similar overall partitioning patterns with the D values being lower at 1430°C. Rare earth elements display a wide range of partition coefficients, D_La (0.012-0.026) being approximately one order of magnitude lower than D_Lu (0.18-0.23). Partition coefficients for the 2+ and 3+ cations entering the M2-site exhibit a near-parabolic dependence on radius of the incorporated cations as predicted from the lattice strain model. This underlines the contribution made by the crystal structure toward controlling the distribution of trace elements. Using data obtained in this study combined with that in the published literature, we also discuss the effects that other important parameters, namely, melt composition, pressure, and temperature, could have on partitioning.Our partition coefficients have been used to model the generation of the Group I (GI) kimberlites from South Africa. The numerical modeling shows that kimberlitic melts can be produced by ∼0.5% melting of a MORB-type depleted source that has been enriched by small-degree melts originating from a similar depleted source. This result suggests that the source of GI kimberlites may be located at the lithosphere-asthenosphere transition. Percolation of small degree melts from the asthenosphere would essentially create a metasomatic horizon near the bottom of the non-convecting sublithospheric mantle. Accumulation of such small degree melts together with the presence of volatiles and conductive heating would trigger melting of the ambient mantle and subsequently lead to eruption of kimberlitic melts. Additionally, our model shows that the GI source can be generated by metasomatism of a 2 Ga old MORB source ca. 1 Ga ago. Assuming that MORB-type mantle is the most depleted source of magmas on earth, then this is the oldest age at which the GI source could have existed. However, this age most likely reflects the average age of a series of metasomatic events than that of a single event.     

Distribution and sources of aliphatic hydrocarbons in surface sediments of Sergipe River estuarine system

The assessment of aliphatic hydrocarbons was performed in the Sergipe River estuarine system, northeastern Brazil. Aliphatic hydrocarbons concentration ranged from 9.9 ug g?1 up to 30.8 ug g?1 of dry sediment. The carbon preference index (CPI, based on nC?? to nC?? range), indicated predominance of petrogenic input in two of the sites analyzed (P4 and P5). The unresolved complex mixture (UCM) was found to be present in seven of the nine sites sampled (except for P4 and P5). Overall, the results of this work suggest that there is a mix of organic matter sources to the sediment. Although the coast of Sergipe has an intense off shore petroleum exploration and the Sergipe River crosses the entire city of Aracaju, the capital city of Sergipe, non-significant anthropogenic fingerprint was assessed.     

Magnetic fabrics and rock magnetism of Proterozoic dike swarm from the southern São Francisco Craton, Minas Gerais State, Brazil

Magnetic fabric and rock magnetism studies were performed on 32 mafic dikes of a Proterozoic dike swarm from the southern São Francisco Craton (SFC; Minas Gerais State, SE Brazil). Magnetic anisotropies were determined by applying anisotropy of low-field magnetic susceptibility (AMS) and anisotropy of remanent magnetization (ARM). The latter was performed imposing both anhysteretic (total (AAR) and partial pAAR)) and isothermal remanence magnetizations (AIRM). Partial anhysteretic remanence anisotropy was performed based on remanent coercivity spectra from a pilot specimen of each site. In most sites, AMS is dominantly carried by ferromagnetic minerals, however, in some sites, the paramagnetic contribution exceeds 70% of bulk susceptibility. Rock magnetism and thin section analysis allow classifying the dikes as non-hydrothermalized and hydrothermalized. Magnetic measurement shows that the mean magnetic susceptibility is usually lower than 5×10⁻³ (SI). Ti-poor titanomagnetites up to pure magnetite pseudo-single-domain (PSD) grain sizes carry the majority of magnetic fabrics for non-hydrothermalized dikes whereas coarse to fine grained Ti-poor titanomagnetites carry the majority of magnetic fabrics for hydrothermalized dikes.Three primary AMS fabrics are recognized which are coaxial with ARM fabric, except for two dikes, from both non-hydrothermalized and hydrothermalized dikes. Normal AMS fabric surprisingly is not dominant (31%). The parallelism between AMS, pAAR_0–30, pAAR_30–60 and pAAR_60–90 fabrics in the hydrothermalized dikes indicates that magnetic grains formed due to late-stage crystallization or to remobilization of iron oxides due to hydrothermal alteration after dike emplacement have acquired a mimetic fabric coaxial with the primary fabric given by coarse-grained early crystallized Ti-poor titanomagnetites. This fabric is interpreted as magma flow in which the analysis of K_max inclination permitted the inference that the dikes were fed by horizontal or subhorizontal fluxes (K_max<30°). Intermediate AMS fabric is the most important (41%) in the investigated swarm. It is interpreted as due to vertical compaction of a static magma column with the minimum stress along the dike strike. ARM determinations for these sites also remained intermediate except for two dikes. In one of them, AIRM fabric resulted in normal AMS fabric while for the other AAR fabric resulted in inverse AMS fabric. A combination of AMS and ARM fabrics suggest that magmatic fabric for both dikes were overprinted by some late local event, probably related to Brasiliano orogenic processes after dike emplacement. InverseInverse AMS fabric is a minority (four dikes). ARM determinations also remained inverse suggesting a primary origin for inverse AMS fabric.     

Biogenic carbon flows through the planktonic food web of the Amundsen Gulf (Arctic Ocean): A synthesis of field measurements and inverse modeling analyses

Major pathways of biogenic carbon (C) flow are resolved for the planktonic food web of the flaw lead polynya system of the Amundsen Gulf (southeast Beaufort Sea, Arctic Ocean) in spring-summer 2008. This period was relevant to study the effect of climate change on Arctic marine ecosystems as it was characterized by unusually low ice cover and warm sea surface temperature. Our synthesis relied on a mass balance estimate of gross primary production (GPP) of 52.5 ± 12.5 g C m⁻² calculated using the drawdown of nitrate and dissolved inorganic C, and a seasonal f-ratio of 0.64. Based on chlorophyll a biomass, we estimated that GPP was dominated by phytoplankton (93.6%) over ice algae (6.4%) and by large cells (>5 μm, 67.6%) over small cells (<5 μm, 32.4%). Ancillary in situ data on bacterial production, zooplankton biomass and respiration, herbivory, bacterivory, vertical particle fluxes, pools of particulate and dissolved organic carbon (POC, DOC), net community production (NCP), as well as selected variables from the literature were used to evaluate the fate of size-fractionated GPP in the ecosystem. The structure and functioning of the planktonic food web was elucidated through inverse analysis using the mean GPP and the 95% confidence limits of every other field measurement as lower and upper constraints. The model computed a net primary production of 49.2 g C m⁻², which was directly channeled toward dominant calanoid copepods (i.e. Calanus hyperboreus 20%, Calanus glacialis 10%, and Metridia longa 10%), other mesozooplankton (12%), microzooplankton (14%), detrital POC (18%), and DOC (16%). Bacteria required 29.9 g C m⁻², a demand met entirely by the DOC derived from local biological activities. The ultimate C outflow comprised respiration fluxes (82% of the initial GPP), a small sedimentation (3%), and a modest residual C flow (15%) resulting from NCP, dilution and accumulation. The sinking C flux at the model limit depth (395 m) supplied 60% of the estimated benthic C demand (2.8 g C m⁻²), suggesting that the benthos relied partly on other C sources within the bottom boundary layer to fuel its activity. In summary, our results illustrate that the ongoing decline in Arctic sea ice promotes the growth of pelagic communities in the Amundsen Gulf, which benefited from a ∼80% increase in GPP in spring-summer 2008 when compared to 2004 – a year of average ice conditions and relatively low GPP. However, 53% of the secondary production was generated within the microbial food web, the net ecological efficiency of zooplankton populations was not particularly high (13.4%), and the quantity of biogenic C available for trophic export remained low (6.6 g C m⁻²). Hence it is unlikely that the increase in lower food web productivity, such as the one observed in our study, could support new harvestable fishery resources in the offshore Beaufort Sea domain.     

Talc mineralisation associated with soft hematite ore, Gongo Soco deposit, Minas Gerais, Brazil: petrography, mineral chemistry and boron-isotope composition of tourmaline

Talc mineralisation occurs as hematite–talc schist between soft hematite ore and dolomitic itabirite at Gongo Soco, Quadrilátero Ferrífero of Minas Gerais, Brazil. The hematite–talc schist and soft hematite have a prominent tectonic foliation of tabular hematite. Tabular hematite without preferential orientation is superimposed on the tectonic foliation. The talcose schist is enriched in F and has a constant Fe/S ratio. Electron-microprobe analyses indicate trace amounts of S in different generations of hematite. The whole-rock Fe/S ratio possibly represents sulfate S from hematite-hosted fluid inclusions. Fluid inclusions in foliation-overprinting hematite and chlorite geothermometry from talcose rocks suggest, respectively, temperatures from <200°C to ~300°C. Tourmaline, a rarely observed mineral in the hematite–talc schist, belongs to the alkali group and falls in the dravite compositional field. Boron-isotope determinations of tourmaline crystals, using secondary ion mass spectrometry, vary from −20‰ to −12‰ δ¹¹B. This compositional isotopic range and the tourmaline chemical composition suggest a meta-evaporitic origin. A non-marine evaporitic setting is the most likely source of acidic, highly oxidising fluids, which resulted in the abundant F-bearing talc and the presence of otherwise immobile Ti in hematite. Oxidising brines were channelled along shear zones and converted dolomitic itabirite into the Gongo Soco soft hematite and the talc mineralisation. The latter is envisaged as the hydrothermal wall-rock alteration of dolomitic itabirite, which gave rise to the soft hematite ore.     

Physicochemical properties of fine aerosols at Plan d'Aups during ESCOMPTE

The physical and chemical properties of aerosol particles were investigated at Plan d'Aups, one of the ESCOMPTE sites located in the St. Baume mountain area (700 m a.s.l.), 50 km east of Marseilles (France). The site is ideally located for assessing the vertical and horizontal extent of the pollution plume from the Marseilles–Berre area.Our study showed that polluted air masses from the Marseilles–Berre area are advected to Plan d'Aups in the early afternoon. Average daily concentration of particles reaches up to 40 μg m⁻³ while 1-h average particle number concentration is greater than 30,000 cm⁻³. Most of the particle mass is composed of SO₄²⁻ and organic carbon (OC). The chemical properties of the particles revealed that an additional source, possibly from the industrial area of Gardanne, contributes to the aerosol mass. This last source is characterised by significant emissions of elements, such as Zn, V, Al and Si.In addition to transport, we found that gas-to-particle conversion takes place at the interface between the free troposphere and the boundary layer. We estimated that on average, 30% of the particle number is accounted for by direct nucleation. This is potentially a major aerosol source to the free troposphere.     

Long-term evolution of the spin of Venus: II. numerical simulations

We present here the numerical application of the theoretical results derived in Correia et al. (2003, Icarus 163, 1-23) for the spin evolution of Venus since its formation. We explore a large variety of initial conditions to cover the possible formation and evolutionary scenarios. In particular, we pay special attention to the evolutions which cross the chaotic zone resulting from secular planetary perturbations (Laskar and Robutel, 1993, Nature 361, 608-612). We demonstrate that Venus’ axis can be temporarily trapped in a secular resonance with the node of Neptune’s orbit, which can prevent it from being tilted to 180° and will drive it toward 0°. We test several dissipation models and parameters to evaluate their contribution to the planet’s spin history. We confirm that despite the variations in the models, only three of the four final spin states of Venus are possible (Correia and Laskar, 2001, Nature 411, 767-770) and that the present observed retrograde spin state of Venus can be attained by two different processes. In the first scenario (π−), the axis is tilted toward 180° while its rotation rate slows down, while in the second one, the axis is driven toward 0° obliquity and the rotation rate decreases, stops, and increases again in the reverse direction to a final equilibrium value (0−).