Shock‐darkening in ordinary chondrites: Determination of the pressure‐temperature conditions by shock physics mesoscale modeling

We determined the shock‐darkening pressure range in ordinary chondrites using the iSALE shock physics code. We simulated planar shock waves on a mesoscale in a sample layer at different nominal pressures. Iron and troilite grains were resolved in a porous olivine matrix in the sample layer. We used equations of state (Tillotson EoS and ANEOS) and basic strength and thermal properties to describe the material phases. We used Lagrangian tracers to record the peak shock pressures in each material unit. The post‐shock temperatures (and the fractions of the tracers experiencing temperatures above the melting point) for each material were estimated after the passage of the shock wave and after the reflections of the shock at grain boundaries in the heterogeneous materials. The results showed that shock‐darkening, associated with troilite melt and the onset of olivine melt, happened between 40 and 50 GPa with 52 GPa being the pressure at which all tracers in the troilite material reach the melting point. We demonstrate the difficulties of shock heating in iron and also the importance of porosity. Material impedances, grain shapes, and the porosity models available in the iSALE code are discussed. We also discuss possible not‐shock‐related triggers for iron melt.     

Estimates of Horizontal Groundwater Flow Velocities in Boreholes

Currently, monitoring tools can be deployed in observation boreholes to better assess groundwater flow, flux of dissolved contaminants and their mass discharge in an aquifer. The relationship between horizontal water velocity in observation boreholes and Darcy fluxes in the surrounding aquifer has been studied for natural flow conditions (i.e., no pumping). Interpretation of measurements taken with dilution tests, the colloidal borescope, the Heat Pulse Flowmeter, and other techniques require the conversion of observed borehole velocity u to aquifer Darcy flux q_∞ . This conversion is typically done through a proportionality factor α = u/q_∞ . In experimental studies as well as in theoretical developments, reported values of α vary almost three orders of magnitude (from 0.5 to 10). This large variability in reported values of α could be explained by: (1) unclear distinction between Darcy flux and water seepage velocity, (2) unclear definition of water velocity in the borehole, (3) effects of well screen and the presence of the measurement device itself on the observable velocities, and (4) hydraulic conditions in the borehole annulus. We address (1), (2) from a conceptual/theoretical perspective, and (3) by means of numerical simulations. We show that issue (1) in low porosity aquifers can yield to order-of-magnitude discrepancies in estimates of q_∞ ; (2) may result in discrepancies of up to 50%, and (3) can cause differences up to 20% of water velocity in the borehole void space compared to the theoretical case of an open borehole.     

Patterns in the spatial distribution of Peruvian anchovy (Engraulis ringens) revealed by spatially explicit fishing data

Peruvian anchovy (Engraulis ringens) stock abundance is tightly driven by the high and unpredictable variability of the Humboldt Current Ecosystem. Management of the fishery therefore cannot rely on mid- or long-term management policy alone but needs to be adaptive at relatively short time scales. Regular acoustic surveys are performed on the stock at intervals of 2 to 4 times a year, but there is a need for more time continuous monitoring indicators to ensure that management can respond at suitable time scales. Existing literature suggests that spatially explicit data on the location of fishing activities could be used as a proxy for target stock distribution. Spatially explicit commercial fishing data could therefore guide adaptive management decisions at shorter time scales than is possible through scientific stock surveys. In this study we therefore aim to (1) estimate the position of fishing operations for the entire fleet of Peruvian anchovy purse–seiners using the Peruvian satellite vessel monitoring system (VMS), and (2) quantify the extent to which the distribution of purse–seine sets describes anchovy distribution. To estimate fishing set positions from vessel tracks derived from VMS data we developed a methodology based on artificial neural networks (ANN) trained on a sample of fishing trips with known fishing set positions (exact fishing positions are known for approximately 1.5% of the fleet from an at-sea observer program). The ANN correctly identified 83% of the real fishing sets and largely outperformed comparative linear models. This network is then used to forecast fishing operations for those trips where no observers were onboard. To quantify the extent to which fishing set distribution was correlated to stock distribution we compared three metrics describing features of the distributions (the mean distance to the coast, the total area of distribution, and a clustering index) for concomitant acoustic survey observations and fishing set positions identified from VMS. For two of these metrics (mean distance to the coast and clustering index), fishing and survey data were significantly correlated. We conclude that the location of purse–seine fishing sets yields significant and valuable information on the distribution of the Peruvian anchovy stock and ultimately on its vulnerability to the fishery. For example, a high concentration of sets in the near coastal zone could potentially be used as a warning signal of high levels of stock vulnerability and trigger appropriate management measures aimed at reducing fishing effort.     

Trophic interactions in the coastal ecosystem of Sri Lanka: An ECOPATH preliminary approach

This study attempts to assemble and summarize existing information in order to build a general representation of the trophic interactions within the shallow coastal ecosystem of Sri Lanka. A multispecific ecosystem-based approach on trophic relationships and their possible variations was performed using ECOPATH. Thirty-nine functional groups were considered representing all trophic levels in the food web.     

Sensitivity equations for hyperbolic conservation law-based flow models

The present paper focuses on the governing equations for the sensitivity of the variables to the parameters in flow models that can be described by one-dimensional scalar, hyperbolic conservation laws. The sensitivity is shown to obey a hyperbolic, scalar conservation law. The sensitivity is a conserved scalar except in the case of discontinuous flow solutions, where an extra, point source term must be added to the equations in order to enforce conservation. The propagation speed of the sensitivity waves being identical to that of the conserved variable in the original conservation law, the system of conservation laws formed by the original hyperbolic equation and the equation satisfied by the sensitivity is linearly degenerate. A consequence on the solution of the Riemann problem is that rarefaction waves for the variable of the original equation result in vacuum regions for the sensitivity. The numerical solution of the hyperbolic conservation law for the sensitivity by finite volume methods requires the implementation of a specific shock detection procedure. A set of necessary conditions is defined for the discretisation of the source term in the sensitivity equation. An application to the one-dimensional kinematic wave equation shows that the proposed numerical technique allows analytical solutions to be reproduced correctly. The computational examples show that first-order numerical schemes do not yield satisfactory numerical solutions in the neighbourhood of moving shocks and that higher-order schemes, such as the MUSCL scheme, should be used for sharp transients.     

Influence of the oceanic biology on the tropical Pacific climate in a coupled general circulation model

The influence of chlorophyll spatial patterns and variability on the tropical Pacific climate is investigated by using a fully coupled general circulation model (HadOPA) coupled to a state-of-the-art biogeochemical model (PISCES). The simulated chlorophyll concentrations can feedback onto the ocean by modifying the vertical distribution of radiant heating. This fully interactive biological-ocean-atmosphere experiment is compared to a reference experiment that uses a constant chlorophyll concentration (0.06 mg m⁻³). It is shown that introducing an interactive biology acts to warm the surface eastern equatorial Pacific by about 0.5°C. Two competing processes are involved in generating this warming: (a) a direct 1-D biological warming process in the top layers (0–30 m) resulting from strong chlorophyll concentrations in the upwelling region and enhanced by positive dynamical feedbacks (weaker trade winds, surface currents and upwelling) and (b) a 2-D meridional cooling process which brings cold off-equatorial anomalies from the subsurface into the equatorial mixed layer through the meridional cells. Sensitivity experiments show that the climatological horizontal structure of the chlorophyll field in the upper layers is crucial to maintain the eastern Pacific warming. Concerning the variability, introducing an interactive biology slightly reduces the strength of the seasonal cycle, with stronger SST warming and chlorophyll concentrations during the upwelling season. In addition, ENSO amplitude is slightly increased. Similar experiments performed with another coupled general circulation model (IPSL-CM4) exhibit the same behaviour as in HadOPA, hence showing the robustness of the results.     

Interannual variability of the South Pacific Convergence Zone and implications for tropical cyclone genesis

The interannual variability of the South Pacific Convergence Zone (SPCZ) and its influence on tropical cyclone (TC) genesis in the South Pacific are investigated using observations and ERA40 reanalysis over the 1979?C2002 period. In austral summer, the SPCZ displays four typical structures at interannual timescales. The first three are characterized by a diagonal orientation of the SPCZ and account for 85% of the summer seasons. One is close to climatology and the other two exhibit a 3° northward or southward departure from the SPCZ climatological position. In contrast, the fourth one, that only encompasses three austral summer seasons (the extreme 1982/1983 and 1997/1998 El Ni?o events and the moderate 1991/1992 El Ni?o event), displays very peculiar behaviour where the SPCZ largely departs from its climatological position and is zonally oriented. Variability of the western/central Pacific equatorial sea surface temperature (SST) is shown to modulate moisture transport south of the equator, thereby strongly constraining the location of the SPCZ. The SPCZ location is also shown to strongly modulate the atmospheric circulation variability in the South Pacific with specific patterns for each class. However, independently of its wide year-to-year excursions, the SPCZ is always collocated with the zero relative vorticity at low levels while the maximum vorticity axis lies 6° to the south of the SPCZ position. This coherent atmospheric organisation in the SPCZ region is shown to constrain tropical cyclogenesis to occur preferentially within 10° south of the SPCZ location as this region combines all the large-scale atmospheric conditions that favour the breeding of TCs. This analysis also reveals that cyclogenesis in the central Pacific (in the vicinity of French Polynesia) only occurs when the SPCZ displays a zonal orientation while this observation was previously attributed to El Ni?o years in general. Different characteristics of El Ni?o Southern Oscillation (ENSO)-related Pacific equatorial warming are shown to impact differently on the SPCZ position, suggesting that for regional scales, such as the South Pacific, the SPCZ classification is more appropriate than a simple ENSO index to characterize TC interannual variability. These findings suggest that forecasting the strength of El Ni?o through SST variations in the eastern Pacific may not be sufficient to accurately predict cyclogenesis in the South Pacific, especially east of the dateline.     

The Paris CM chondrite: Secondary minerals and asteroidal processing

We report a petrographic and mineralogical survey of Paris, a new CM chondrite considered to be the least‐altered CM identified so far (Hewins et al. 2014 ). Compared to other CMs, Paris exhibits (1) a higher concentration of Fe‐Ni metal beads, with nickel contents in the range 4.1–8.1 wt%; (2) the systematic presence of thin lamellae and tiny blebs of pentlandite in pyrrhotite grains; and (3) ubiquitous tochilinite/cronstedtite associations with higher FeO/SiO₂ and S/SiO₂ ratios. In addition, Paris shows the highest concentration of trapped ³⁶Ar reported so far for a CM chondrite (Hewins et al. 2014 ). In combination with the findings of previous studies, our data confirm the reliability of (1) the alteration sequence based on the chemical composition of tochilinite/cronstedtite associations to quantify the fluid alteration processes and (2) the use of Cr content variability in type II ferroan chondrule olivine as a proxy of thermal metamorphism. In contrast, the scales based on (1) the Fe³⁺ content of serpentine in the matrix to estimate the degree of aqueous alteration and (2) the chemical composition of Fe‐Ni metal beads for quantifying the intensity of the thermal metamorphism are not supported by the characteristics of Paris. It also appears that the amount of trapped ³⁶Ar is a sensitive indicator of the secondary alteration modifications experienced by chondrites, for both aqueous alteration and thermal metamorphism. Considering Paris, our data suggest that this chondrite should be classified as type 2.7 as it suffered limited but significant fluid alteration and only mild thermal metamorphism. These results point out that two separated scales should be used to quantify the degree of the respective role of aqueous alteration and thermal metamorphism in establishing the characteristics of CM chondrites.     

The Orgueil meteorite: 150 years of history

The goal of this paper is to summarize 150 yr of history of a very special meteorite. The Orgueil meteorite fell near Montauban in southwestern France on May 14, 1864. The bolide, which was the size of the full Moon, was seen across Western France, and almost immediately made the news in local and Parisian newspapers. Within a few weeks of the fall, a great diversity of analyses were performed under the authority of Gabriel Auguste Daubrée, geology professor at the Paris Museum, and published in the Comptes Rendus de l'Académie des Sciences. The skilled scientists reported the presence of iron sulfides, hydrated silicates, and carbonates in Orgueil. They also characterized ammonium salts which are now gone, and observed sulfates being remobilized at the surface of the stone. They identified the high water and carbon contents, and noted similarities with the Alais meteorite, which had fallen in 1806, 300 km away. While Daubrée and his colleagues noted the similarity of the Orgueil organic matter with some terrestrial humus, they were cautious not to make a direct link with living organisms. One century later, Nagy and Claus were less prudent and announced the discovery of “organized” elements in some samples of Orgueil. Their observations were quickly discredited by Edward Anders and others who also discovered that some pollen grains were intentionally placed into the rock back in the 1860s. Orgueil is now one of the most studied meteorites, indeed one of the most studied rocks of any kind. Not only does it contain a large diversity of carbon‐rich compounds, which help address the question of organo‐synthesis in the early solar system but its chemical composition is also close to that of the Sun's photosphere and serves as a cosmic reference. Secondary minerals, which make up 99% of the volume of Orgueil, were probably formed during hydrothermal alteration on the parent‐body within the first few million years of the solar system; their study is essential to our understanding of fluid–rock interaction in asteroids and comets. Finally, the Orgueil meteorite probably originated from a volatile‐rich “cometary” outer solar system body as indicated by its orbit. Because it bears strong similarities to other carbonaceous chondrites that originated on dark asteroids, this cometary connection supports the idea of a continuum between dark asteroids and comets.     

Mineralogy and petrography of C asteroid regolith: The Sutter's Mill CM meteorite

Based upon our characterization of three separate stones by electron and X‐ray beam analyses, computed X‐ray microtomography, Raman microspectrometry, and visible‐IR spectrometry, Sutter's Mill is a unique regolith breccia consisting mainly of various CM lithologies. Most samples resemble existing available CM2 chondrites, consisting of chondrules and calcium‐aluminum‐rich inclusion (CAI) set within phyllosilicate‐dominated matrix (mainly serpentine), pyrrhotite, pentlandite, tochilinite, and variable amounts of Ca‐Mg‐Fe carbonates. Some lithologies have witnessed sufficient thermal metamorphism to transform phyllosilicates into fine‐grained olivine, tochilinite into troilite, and destroy carbonates. One finely comminuted lithology contains xenolithic materials (enstatite, Fe‐Cr phosphides) suggesting impact of a reduced asteroid (E or M class) onto the main Sutter's Mill parent asteroid, which was probably a C class asteroid. One can use Sutter's Mill to help predict what will be found on the surfaces of C class asteroids such as Ceres and the target asteroids of the OSIRIS‐REx and Hayabusa 2 sample return missions (which will visit predominantly primitive asteroids). C class asteroid regolith may well contain a mixture of hydrated and thermally dehydrated indigenous materials as well as a significant admixture of exogenous material would be essential to the successful interpretation of mineralogical and bulk compositional data.