Contact metamorphism/hydrothermal alteration of Tertiary basalts from the Isle of Skye,northwest Scotland

Alkali olivine basalts from Skye were simultaneously contact metamorphosed by Tertiary gabbro and granite intrusions and altered by the hydrothermal convection system that the plutons induced. Four metamorphic zones were mapped around the plutons. Furthest from the intrusions, in the primary olivine zone, metabasalts are composed of combinations of igneous olivine, augite, plagioclase, titaniferous magnetite, ilmenite, zeolites, gyrolite, sulfides, and chlorite-smectite intergrowths. Closer to the plutons, in the smectite zone, saponite and carbonate appear, primary olivine and gyrolite disappear, and zeolites decrease dramatically in abundance. Still closer to the plutons, in the amphibole zone, actinolite, edenite, chlorite, sphene, epidote, andradite, and quartz appear and saponite and chlorite-smectite intergrowths disappear. Along parts of the contact between gabbro and basalt, in the orthopyroxeneolivine zone, orthopyroxene, metamorphic olivine, and biotite appear and amphibole, chlorite, sphene, epidote, andradite, carbonate, and quartz disappear. Whole-rock chemical data indicate only minor change in the major-element chemical composition of the metabasalts during progressive metamorphism/hydrothermal alteration. Two-pyroxene eothermometry and various mineral-fluid equilibria suggest the range of peak temperatures attained in the metamorphic zones: orthopyroxene-olivine zone, 900°1, 030° C; amphibole zone, 400°–900° C; smectite and primary olivine zones, < 400° C. Mineralogical and oxygen isotopic alteration of the metabasalts were closely coupled: Basalts from the primary olivine zone with nearly unaltered igneous mineralogies have normal or near-normal wholerock ¹⁸O>+5 (SMOW); mineralogically more altered basalts from the smectite zone have whole-rock ¹⁸O=+2 to +5; still more mineralogically altered basalts from the amphibole zone (with one exception) have ¹⁸O<+ 2; completely recrystallized hornfelses from the orthopyroxene-olivine zone have ¹⁸O<0. The principal mechanism of isotope exchange between basalt and metamorphic/ hydrothermal fluid probably was heterogeneous mineralfluid reaction.Metabasalts from the orthopyroxene-olivine zone are mineralogically fresh pyroxene hornfelses that record crystallization temperatures > 1,000° C yet have highly altered whole-rock oxygen isotope compositions, ¹⁸O<0%. The hornfelses chemically interacted with metamorphic/hydrothermal fluids either at very high temperatures or while they were heated to > 1,000° C or both. Their mineralogy, however, rules out significant water-rock interaction after they cooled below 900° C. Hydrothermal convection on Skye was a two-stage process: (a) fluid flow through wall rocks initially was pervasive while they are heated; (b) fluid flow after the thermal peak in the wall rocks was sufficiently channelized that rocks such as those in the orthopyroxeneolivine zone were isolated from further fluid-rock interaction during all or almost all of the cooling history of the hydrothermal system.     

The direction of fluid flow during contact metamorphism of siliceous carbonate rocks: new data for the Monzoni and Predazzo aureoles,northern Italy,and a global review

Periclase formed in siliceous dolomitic marbles during contact metamorphism in the Monzoni and Predazzo aureoles, the Dolomites, northern Italy, by infiltration of the carbonate rocks by chemically reactive, H₂O-rich fluids at 500 bar and 565-710 °C. The spatial distribution of periclase and oxygen isotope compositions is consistent with reactive fluid flow that was primarily vertical and upward in both aureoles with time-integrated flux ~5,000 and ~300 mol fluid/cm² rock in the Monzoni and Predazzo aureoles, respectively. The new results for Monzoni and Predazzo are considered along with published studies of 13 other aureoles to draw general conclusions about the direction, amount, and controls on the geometry of reactive fluid flow during contact metamorphism of siliceous carbonate rocks. Flow in 12 aureoles was primarily vertically upward with and without a horizontal component directed away from the pluton. Fluid flow in two of the other three was primarily horizontal, directed from the pluton into the aureole. The direction of flow in the remaining aureole is uncertain. Earlier suggestions that fluid flow is often horizontal, directed toward the pluton, are likely explained by an erroneous assumption that widespread coexisting mineral reactants and products represent arrested prograde decarbonation reactions. With the exception of three samples from one aureole, time-integrated fluid flux was in the range 10²-10⁴ mol/cm². Both the amount and direction of fluid flow are consistent with hydrodynamic models of contact metamorphism. The orientation of bedding and lithologic contacts appears to be the principal control over whether fluid flow was either primarily vertical or horizontal. Other pre-metamorphic structures, including dikes, faults, fold hinges, and fracture zones, served to channel fluid flow as well.     

Testing the genetic relationship between fluid alteration and brecciation in CM chondrites

Boriskino is a poorly studied CM chondrite with numerous millimeter‐ to centimeter‐scale clasts exhibiting sharp boundaries. Clast textures and mineralogies attest to diverse geological histories with various degrees of aqueous alteration. We conducted a petrographic, chemical, and isotopic study on each clast type of the breccia to investigate if there exists a genetic link between brecciation and aqueous alteration, and to determine the controlling parameter of the extent of alteration. Boriskino is dominated by CM2 clasts for which no specific petrographic type could be assigned based on the chemical compositions and modal abundances of constituents. One clast stands out and is identified as a CM1 lithology, owing to its lack of anhydrous silicates and its overall abundance of dolomite‐like carbonates and acicular iron sulfides. We observe that alteration phases near clast boundaries exhibit foliation features, suggesting that brecciation postdated aqueous alteration. We measured the O‐isotopic composition of Ca‐carbonates and dolomite‐like carbonates to determine their precipitation temperatures following the methodology of Verdier‐Paoletti et al. (2017). Both types of carbonates yield similar ranges of precipitation temperatures independent of clast lithology, ranging from ?13.9 ± 22.4 (2σ) to 166.5 ± 47.3 °C, precluding that temperature alone accounts for the differences between the CM1 and CM2 lithologies. Instead, we suggest that initial water/rock ratios of 0.75 and 0.61 for the CM1 and CM2 clasts, respectively, might control the extent of aqueous alteration. Based on these estimates, we suggest that Boriskino clasts originated from a single parent body with heterogeneous distribution of water either due to local differences in the material permeability or in the initial content of ice available. These conditions would have produced microenvironments with differing geochemical conditions thus leading to a range of degrees of aqueous alteration.     

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.     

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.     

Heterogeneities in the Microstructure and Composition of Aircraft Engine Combustor Soot: Impact on the Water Uptake

Size, morphology, microstructure, chemical composition and hygroscopic properties of aircraft engine combustor (AEC) soot particles are studied by using a combination of several methods, namely atomic force microscopy, transmission electron microscopy, gravimetry, ionic chromatography analysis and wetting observations. From the microstructure and the composition of soot agglomerates, we find that we can separate AEC soot in two fractions having quite different physico-chemical properties: a main fraction of particles containing essentially amorphous carbon with small amounts of oxygen, sulfur and iron and a fraction of impurities characterized by various structures and a large amount of impurities. These properties of aircraft engine combustor soot are compared to those of soot obtained by burning TC1 aviation kerosene in a laboratory burner. It is shown that TC1 soot can be a good surrogate of the AEC main fraction. Such a finding allows us to perform water uptake measurements and to conclude that the AEC main fraction is rather hydrophobic whereas the AEC fraction of impurities is highly hydrophilic The ability of the two fractions of aircraft engine combustor soot to act as cloud condensation nuclei is discussed with respect to their implication in contrail and cirrus formation.Associated to Université de la Méditerranée (Aix-Marseille II) and to Université Paul Cézanne (Aix-Marseille III)     

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.