The planets capture model of V838 Monocerotis: conclusions for the penetration depth of the planet(s)

V838 Mon is the prototype of a new class of objects. Understanding the nature of its multistage outburst and similar systems is challenging. So far, several scenarios have been invoked to explain this group of stars. In this work, the planets-swallowing model for V838 Mon is further investigated, taking into account the findings that the progenitor is most likely a massive B-type star. We find that the super-Eddington luminosity during the eruption can explain the fast rising times of the three peaks in the optical light curve. We used two different methods to estimate the location where the planets were consumed. There is a nice agreement between the values obtained from the luminosities of the peaks and from their rising time-scale. We estimate that the planets were stopped at a typical distance of one solar radius from the centre of the host giant star. The planets-devouring model seems to give a satisfying explanation to the differences in the luminosities and rising times of the three peaks in the optical light curve of V838 Mon. The peaks may be explained by the consumption of three planets or alternatively by three steps in the terminal falling process of a single planet. We argue that only the binary merger and the planets-swallowing models are consistent with the observations of the new type of stars defined by V838 Mon.     

MANGALA VALLES, MARS: ASSESSMENT OF EARLY STAGES OF FLOODING AND DOWNSTREAM FLOOD EVOLUTION

The Mangala Valles system is an ∼ ∼900 km fluvially carved channel system located southwest of the Tharsis rise and is unique among the martian outflow channels in that it heads at a linear fracture within the crust as opposed to a collapsed region of chaos as is the case with the circum-Chryse channels. Mangala Valles is confined within a broad, north–south trending depression, and begins as a single valley measuring up to 350 km wide that extends northward from a Memnonia Fossae graben, across the southern highlands toward the northern lowlands. Approximately 600 km downstream, this single valley branches into multiple channels, which ultimately lose their expression at the dichotomy boundary. Previous investigations of Mangala Vallis suggested that many of the units mapped interior to the valley were depositional, related to flooding, and that a minimum of two distinct periods of flooding separated by tens to hundreds of millions of years were required to explain the observed geology. We use infrared and visible images from the THermal EMission Imaging System (THEMIS), and topographic data from the Mars Orbiting Laser Altimeter (MOLA), to investigate the nature of the units mapped within Mangala Vallis. We find that the geomorphology of the units, as well as their topographic and geographic distribution, are consistent with most of them originating from a single assemblage of volcanic flow deposits, once continuous with volcanic flows to the south of the Memnonia Fossae source graben. These flows resurfaced the broad, north–south trending depression into which Mangala Vallis formed prior to any fluvial activity. Later flooding scoured and eroded this volcanic assemblage north of the Mangala source graben, resulting in the present distribution of the units within Mangala Vallis. Additionally, our observations suggest that a single period of catastrophic flooding, rather than multiple periods separated by tens to hundreds of millions of years, is consistent with and can plausibly explain the interior geology of Mangala Vallis. Further, we present a new scenario for the source and delivery of water to the Mangala source graben that models flow of groundwater through a sub-cryosphere aquifer and up a fracture that cracks the cryosphere and taps this aquifer. The results of our model indicate that the source graben, locally enlarged to a trough near the head region of Mangala, would have required less than several days to fill up prior to any spill-over of water to the north. Through estimates of the volume of material missing from Mangala (13,000–20,000 km³), and calculation of mean discharge rates through the channel system (∼ ∼5 × 10⁶ m³ s⁻¹), we estimate that the total duration of fluvial activity through the Mangala Valles was 1–3 months.     

Investigation and sensitivity analysis of a mechanistic phytoplankton model implemented in a new modular numerical tool (Eco3M) dedicated to biogeochemical modelling

A new class of phytoplankton models with a mechanistic basis has been presented in a companion paper (Baklouti, M., Diaz, F., Pinazo, C., Faure, V., Queguiner, B., 2006. Investigation of mechanistic formulations depicting phytoplankton dynamics for models of marine pelagic ecosystems. Progress in Oceanography). It is the default class of models implemented in our new numerical tool Eco3M, which is dedicated to Ecological, Mechanistic and Modular Modelling. A brief overview of its main features is given in Section 2 of the present paper. In the next sections, a particular phytoplankton model among the aforementioned class has been tested with special emphasis on the mechanistic photosynthesis component relating the photosynthetic rate to the proportion of open photosystems II. The present study encompasses several essential steps that are inherent to any modelling, including model reduction, model sensitivity analysis and comparison of model outputs with experiments. The global sensitivity analysis of the plankton model for one-at-a-time parameter perturbations revealed a restricted set of parameters having major influence on the model outputs. Sensitivity tests involving simultaneous parameter perturbations within the range actually encountered in the literature provided a confidence interval for the outputs. Chemostat experiments performed on nitrate-limited diatoms grown under low (LL) and high-light (HL) conditions have been used for comparison with model outputs. The good fit between measured data and model outputs using the same parameter values in both the LL and HL cases demonstrates the ability of our model to represent the main features of phytoplankton dynamics including photoacclimation. Finally, Eco3M is ultimately intended to include explicit bacterial and zooplankton compartments, as well as to be coupled with ocean circulation models, but the intrinsic behavior of the phytoplankton model has been investigated first, independently of physical forcing.     

Deep submarine pyroclastic eruptions: theory and predicted landforms and deposits

Submarine pyroclastic eruptions at depths greater than a few hundred meters are generally considered to be rare or absent because the pressure of the overlying water column is sufficient to suppress juvenile gas exsolution so that magmatic disruption and pyroclastic activity do not occur. Consideration of detailed models of the ascent and eruption of magma in a range of sea floor environments shows, however, that significant pyroclastic activity can occur even at depths in excess of 3000 m. In order to document and illustrate the full range of submarine eruption styles, we model several possible scenarios for the ascent and eruption of magma feeding submarine eruptions: (1) no gas exsolution; (2) gas exsolution but no magma disruption; (3) gas exsolution, magma disruption, and hawaiian-style fountaining; (4) volatile content builds up in the magma reservoir leading to hawaiian eruptions resulting from foam collapse; (5) magma volatile content insufficient to cause fragmentation normally but low rise speed results in strombolian activity; and (6) volatile content builds up in the top of a dike leading to vulcanian eruptions. We also examine the role of bulk-interaction steam explosivity and contact-surface steam explosivity as processes contributing to volcaniclastic formation in these environments. We concur with most earlier workers that for magma compositions typical of spreading centers and their vicinities, the most likely circumstance is the quiet effusion of magma with minor gas exsolution, and the production of somewhat vesicular pillow lavas or sheet flows, depending on effusion rate. The amounts by which magma would overshoot the vent in these types of eruptions would be insufficient to cause any magma disruption. The most likely mechanism of production of pyroclastic deposits in this environment is strombolian activity, due to the localized concentration of volatiles in magma that has a low rise rate; magmatic gas collects by bubble coalescence, and ascends in large isolated bubbles which disrupt the magma surface in the vent, producing localized blocks, bombs, and pyroclastic deposits. Another possible mode of occurrence of pyroclastic deposits results from vulcanian eruptions; these deposits, being characterized by the dominance of angular blocks of country rocks deposited in the vicinity of a crater, should be easily distinguishable from strombolian and hawaiian eruptions. However, we stress that a special case of the hawaiian eruption style is likely to occur in the submarine environment if magmatic gas buildup occurs in a magma reservoir by the upward drift of gas bubbles. In this case, a layer of foam will build up at the top of the reservoir in a sufficient concentration to exceed the volatile content necessary for disruption and hawaiian-style activity; the deposits and landforms are predicted to be somewhat different from those of a typical primary magmatic volatile-induced hawaiian eruption. Specifically, typical pyroclast sizes might be smaller; fountain heights may exceed those expected for the purely magmatic hawaiian case; cooling of descending pyroclasts would be more efficient, leading to different types of proximal deposits; and runout distances for density flows would be greater, potentially leading to submarine pyroclastic deposits surrounding vents out to distances of tens of meters to a kilometer. In addition, flows emerging after the evacuation of the foam layer would tend to be very depleted in volatiles, and thus extremely poor in vesicles relative to typical flows associated with hawaiian-style eruptions in the primary magmatic gas case. We examine several cases of reported submarine volcaniclastic deposits found at depths as great as 3000 m and conclude that submarine hawaiian and strombolian eruptions are much more common than previously suspected at mid-ocean ridges. Furthermore, the latter stages of development of volcanic edifices (seamounts) formed in submarine environments are excellent candidates for a wide range of submarine pyroclastic activity due not just to the effects of decreasing water depth, but also to: (1) the presence of a summit magma reservoir, which favors the buildup of magmatic foams (enhancing hawaiian-style activity) and episodic dike emplacement (which favors strombolian-style eruptions); and (2) the common occurrence of alkalic basalts, the CO₂ contents of which favor submarine explosive eruptions at depths greater than tholeiitic basalts. These models and predictions can be tested with future sampling and analysis programs and we provide a checklist of key observations to help distinguish among the eruption styles.     

Scaling of the Vertical Spreading of a Plume of a Passive Tracer in a Rough-Wall Neutral Boundary Layer

The influence of surface roughness on the dispersion of a passive scalar in a rough wall turbulent boundary layer has been studied using wind-tunnel experiments. The surface roughness was varied using different sizes of roughness elements, and different spacings between the elements. Vertical profiles of average concentration were measured at different distances downwind of the source, and the vertical spread of the plume was computed by fitting a double Gaussian profile to the data. An estimate of the integral length scale is derived from the turbulence characteristics of the boundary layer and is then used to scale the measured values of plume spread. This scaling reduces the variability in the data, confirming the validity of the model for the Lagrangian integral time scale, but does not remove it entirely. The scaled plume spreading shows significant differences from predictions of theoretical models both in the near and in the far field. In the region immediately downwind of the source this is due to the influence of the wake of the injector for which we have developed a simple model. In the far field we explain that the differences are mainly due to the absence of large-scale motions. Finally, further downwind of the source the scaled values of plume spread fall into two distinct groups. It is suggested that the difference between the two groups may be related to the lack of dynamical similarity between the boundary-layer flows for varying surface roughness or to biased estimates of the plume spread.     

Turbulent Transfer Between Street Canyons and the Overlying Atmospheric Boundary Layer

The turbulent exchange of momentum between a two-dimensional cavity and the overlying boundary layer has been studied experimentally, using hot-wire anemometry and particle image velocimetry (PIV). Conditions within the boundary layer were varied by changing the width of the canyons upstream of the test canyon, whilst maintaining the square geometry of the test canyon. The results show that turbulent transfer is due to the coupling between the instabilities generated in the shear layer above the canyons and the turbulent structures in the oncoming boundary layer. As a result, there is no single, unique velocity scale that correctly characterizes all the processes involved in the turbulent exchange of momentum across the boundary layer. Similarly, there is no single velocity scale that can characterize the different properties of the turbulent flow within the canyon, which depends strongly on the way in which turbulence from the outer flow is entrained into the cavity and carried round by the mean flow. The results from this study will be useful in developing simple parametrizations for momentum exchange in the urban canopy, in situations where the street geometry consists principally of relatively long, uniform streets arranged in grid-like patterns; they are unlikely to be applicable to sparse geometries composed of isolated three-dimensional obstacles.     

Consequences of explosive eruptions on small Solar System bodies: the case of the missing basalts on the aubrite parent body

Aubrite meteorites are composed of constituents which are almost certainly of igneous origin. If they were generated by the melting and fractionation of enstatite chondrite-like parental material, as seems very likely, then plagioclase-rich, basaltic complements to the aubrites should have formed. However, such materials are not known as individual meteorites, and the compositions of two plagioclase-silica clasts and one albite-silica-(diopside-anorthite) clast (probably an impact melt) in the Norton County aubrite breccia suggest that they are not the putative enstatite-plagioclase basalts. We propose a new mechanism that explains the absence of these materials, showing that the expansion of even very small amounts of volatiles present in a melt approaching the surface of a small, low-gravity body will be enough to disrupt the melt into a spray of droplets moving faster than the local escape velocity. This explosive volcanic process of melt removal requires larger melt volatile contents on larger bodies, and data on the solubility of volatiles in basaltic melts suggest that the process was limited to bodies smaller than about 100 km in radius.     

Volumetric characteristics of lava flows from interferometric radar and multispectral satellite data: the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos

We have used a suite of remotely sensed data, numerical lava flow modeling, and field observations to determine quantitative characteristics of the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos Islands. Flank lava flow areas, volumes, instantaneous effusion rates, and average effusion rates were all determined for these two eruptions, for which only limited syn-eruptive field observations are available. Using data from SPOT, TOPSAR, ERS-1, and ERS-2, we determined that the 1995 Fernandina flow covers a subaerial area of 6.5×10⁶ m² and has a subaerial dense rock equivalent (DRE) volume of 42×10⁶ m³. Field observations, ATSR satellite data, and the FLOWGO numerical model allow us to determine that the effusion rate declined exponentially from a high of ~60–200 m³ s^-1 during the first few hours to <5 m³ s^-1 prior to ceasing after 73 days, with a mean effusion rate of 4–16 m³ s^-1. Integrating the ATSR-derived, exponentially declining effusion rate over the eruption duration produces a total (subaerial + submarine) DRE volume of between 27 and 100×10⁶ m³, the range in values being due to differing assumptions about heat loss characteristics; only values in the higher part of this range are consistent with the independently derived subaerial volume. Using SPOT, TOPSAR, ERS-1, and ERS-2 data, we determine that the 1998 Cerro Azul flow is 16 km long, covers 16 km², and has a DRE volume of 54×10⁶ m³. FLOWGO produces at-vent velocity and effusion rate values of 11 m s^-1 and ~600 m³ s^-1, respectively. The velocity value agrees well with the 12 m s^-1 estimated in the field. The mean effusion rate (total DRE volume/duration) was 7–47 m³ s^-1. Dike dimensions, fissure lengths, and pressure gradients along the conduit based on magma chamber depth estimates of 3–5 km produce mean effusion rates for the two eruptions that range over nearly four orders of magnitude, the range being due to uncertainty in the magma viscosity, dike dimensions, and pressure gradient between magma chamber and vent. Although somewhat consistent with mean effusion rates from other techniques, their wide range makes them less useful. The exponentially declining effusion rates during both eruptions are consistent with release of elastic strain being the driving mechanism of the eruptions. Our results provide independent input parameters for previously published theoretical relationships between magma chamber pressurization and eruption rates that constrain chamber volumes and increases in volume prior to eruption, as well as time constants of exponential decay during the eruption. The results and theoretical relationships combine to indicate that at both volcanoes probably 25–30% of the volumetric increase in the magma chamber erupted as lava onto the surface. In both eruptions the lava flow volumes are less than 1% of the magma chamber volume.     

A composite Fe,Ni‐FeS and enstatite‐forsterite‐diopside‐glass vitrophyre clast in the Larkman Nunatak 04316 aubrite: Origin by pyroclastic volcanism

Abstract– We studied the mineralogy, petrology, and bulk, trace element, oxygen, and noble gas isotopic compositions of a composite clast approximately 20 mm in diameter discovered in the Larkman Nunatak (LAR) 04316 aubrite regolith breccia. The clast consists of two lithologies: One is a quench‐textured intergrowth of troilite with spottily zoned metallic Fe,Ni which forms a dendritic or cellular structure. The approximately 30 μm spacings between the Fe,Ni arms yield an estimated cooling rate of this lithology of approximately 25–30 °C s^?1. The other is a quench‐textured enstatite‐forsterite‐diopside‐glass vitrophyre lithology. The composition of the clast suggests that it formed at an exceptionally high degree of partial melting, perhaps approaching complete melting, and that the melts from which the composite clast crystallized were quenched from a temperature of approximately 1380–1400 °C at a rate of approximately 25–30 °C s^?1. The association of the two lithologies in a composite clast allows, for the first time, an estimation of the cooling rate of a silicate vitrophyre in an aubrite of approximately 25–30 °C s^?1. While we cannot completely rule out an impact origin of the clast, we present what we consider is very strong evidence that this composite clast is one of the elusive pyroclasts produced during pyroclastic volcanism on the aubrite parent body ( Wilson and Keil 1991 ). We further suggest that this clast was not ejected into space but retained on the aubrite parent body by virtue of the relatively large size of the clast of approximately 20 mm. Our modeling, taking into account the size of the clast, suggests that the aubrite parent body must have been between approximately 40 and 100 km in diameter, and that the melt from which the clast crystallized must have contained an estimated maximum range of allowed volatile mass fractions between approximately 500 and approximately 4500 ppm.     

The initial responses of hot liquid water released under low atmospheric pressures: Experimental insights

Experiments have been performed to simulate the shallow ascent and surface release of water and brines under low atmospheric pressure. Atmospheric pressure was treated as an independent variable and water temperature and vapor pressure were examined as a function of total pressure variation down to low pressures. The physical and thermal responses of water to reducing pressure were monitored with pressure transducers, temperature sensors and visible imaging. Data were obtained for pure water and for solutions with dissolved NaCl or CO₂. The experiments showed the pressure conditions under which the water remained liquid, underwent a rapid phase change to the gas state by boiling, and then solidified because of removal of latent heat. Liquid water is removed from phase equilibrium by decompression. Solid, liquid and gaseous water are present simultaneously, and not at the 611 Pa triple point, because dynamic interactions between the phases maintain unstable temperature gradients. After phase changes stop, the system reverts to equilibrium with its surroundings. Surface and shallow subsurface pressure conditions were simulated for Mars and the icy satellites of the outer Solar System. Freezing by evaporation in the absence of wind on Mars is shown to be unlikely for pure water at pressures greater than c. 670 Pa, and for saline solutions at pressures greater than c. 610 Pa. The physical nature of ice that forms depends on the salt content. Ice formed from saline water at pressures less than c. 610 Pa could be similar to terrestrial sea ice. Ice formed from pure water at pressures less than c. 100 Pa develops a low thermal conductivity and a ‘honeycomb’ structure created by sublimation. This ice could have a density as low as c. 450 kg m⁻³ and a thermal conductivity as low as 1.6 W m⁻¹ K⁻¹, and is highly reflective, more akin to snow than the clear ice from which it grew. The physical properties of ice formed from either pure or saline water at low pressures will act to reduce the surface temperature, and hence rate of sublimation, thereby prolonging the lifespan of any liquid water beneath.