Kinetics of heterogeneous bubble nucleation in rhyolitic melts: implications for the number density of bubbles in volcanic conduits and for pumice textures

We performed decompression experiments to simulate the ascent of a phenocryst-bearing rhyolitic magma in a volcanic conduit. The starting materials were bubble-free rhyolites water-saturated at 200 MPa–800°C under oxidizing conditions: they contained 6.0 wt% dissolved H₂O and a dense population of hematite crystals (8.7 ± 2 × 10⁵ mm⁻³). Pressure was decreased from the saturation value to a final value ranging from 99 to 20 MPa, at constant temperature (800°C); the rate of decompression was either 1,000 or 27.8 kPa/s. In all experiments, we observed a single event of heterogeneous bubble nucleation beginning at a pressure P _N equal to 63 ± 3 MPa in the 1,000 kPa/s series, and to 69 ± 1 MPa in the 27.8 kPa/s series. Below P _N, the degree of water supersaturation in the liquid rapidly decreased to a few 0.1 wt%, the nucleation rate dropped, and the bubble number density (BND) stabilized to a value strongly sensitive to decompression rate: 80 mm⁻³ at 27.8 kPa/s vs. 5,900 mm⁻³ at 1,000 kPa/s. This behaviour is like the behavior formerly described in the case of homogeneous bubble nucleation in the rhyolite-H₂O system and in numerical simulations of vesiculation in ascending magmas. Similar degrees of water supersaturation were measured at 27.8 and 1,000 kPa/s, implying that a faster decompression rate does not result in a larger departure from equilibrium. Our experimental results imply that BNDs in acid to intermediate magmas ascending in volcanic conduits will depend on both the decompression rate and the number density of phenocrysts, especially the number density of magnetite microphenocrysts (1–100 mm⁻³), which is the only mineral species able to reduce significantly the degree of water supersaturation required for bubble nucleation. Very low BNDs (≈1 mm⁻³) are predicted in the case of effusive eruptions ( ≈ 0.1 kPa/s). High BNDs (up to 10⁷ mm⁻³) and bimodal bubble size distributions are expected in the case of explosive eruptions: (1) a relatively small number density of bubbles (1–100 mm⁻³) will first nucleate in the lower part of the conduit ( ≈ 10 kPa/s), either at high pressure on magnetite or at lower pressure on quartz and feldspar (or by homogeneous nucleation in the liquid) and (2) then, extreme decompression rates near the fragmentation level ( ≈ 10³ kPa/s) will trigger a major nucleation event leading to the multitude of small bubbles, typically a few micrometers to a few tens of micrometers in diameter, which characterizes most silicic pumices.     

Multiphase lattice Boltzmann simulations for porous media applications

Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this article, we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.     

Molybdenum isotope evidence for mild environmental oxygenation before the Great Oxidation Event

We use the molybdenum isotope paleoredox proxy to look for evidence of small amounts of O₂ in the environment ∼50 Ma before the Great Oxidation Event (GOE) in a high resolution profile from the ∼2.5 Ga Mt. McRae Shale. The molybdenum isotope compositions (δ^98/95Mo) from samples throughout the sequence span a range from 0.99‰ to 1.86‰. All samples have heavier δ^98/95Mo compared to average upper continental crust. In addition, the upper (S1) and lower (S2) black shale units within the Mt. McRae Shale exhibit systematic differences in average isotopic compositions and distinct patterns of δ^98/95Mo variation. Heavier δ^98/95Mo values occur in the S1 unit, where δ^98/95Mo correlates with Mo enrichments. In the S2 unit, δ^98/95Mo is not as heavy and is relatively invariant.Based on sedimentary Fe proxies we infer that S1 sediments were deposited under euxinic conditions, so that Mo removal from bottom waters was likely quantitative. Thus, δ^98/95Mo in this interval likely records coeval seawater. The lighter δ^98/95Mo values in the S2 unit may indicate a less fractionated ocean Mo inventory relative to the S1 unit. However, sedimentary Fe proxies suggest that S2 sediments accumulated under a water column that was ferruginous rather than euxinic, raising the possibility of non-quantitative Mo scavenging and hence an expressed δ^98/95Mo fractionation relative to coeval seawater. Because any associated fractionations during this process would have favored the light isotope in sediments, the lighter δ^98/95Mo values in the S2 unit represent a lower limit on the value in contemporaneous seawater.After evaluating a range of hypotheses, we conclude that the isotopically heavy δ^98/95Mo values seen throughout the Mt. McRae Shale likely reflect the effects of oxidative weathering and adsorption of Mo to oxide mineral surfaces on land or in surface oceans. The extent of environmental oxygenation in either unit is difficult to assess due to uncertainties in the global Mo isotope budget. Because of the small ocean Mo inventory in the Late Archean, documented by low concentrations of Mo and low Mo/TOC, the extent of oxygenation required to account for the observed fractionations is much smaller than in modern oceans. However, when juxtaposed against the record of δ^98/95Mo through time, our findings provide further evidence of the onset of environmental oxygenation before the GOE.     

Galaxy groups in the 2dF Galaxy Redshift Survey: the catalogue

We construct a galaxy groups catalogue from the public 100-K data release of the 2dF Galaxy Redshift Survey. The group identification is carried out using a slightly modified version of the group-finding algorithm developed by Huchra & Geller. Several tests using mock catalogues allow us to find the optimal conditions to increase the reliability of the final group sample. A minimum number of four members, an outer number density enhancement of 80 and a linking radial cut-off of 200 km s⁻¹ are the best obtained values from the analysis. Using these parameters, approximately 90 per cent of groups identified in real space have a redshift space counterpart. On the other hand, the level of contamination in redshift space reaches 30 per cent, including ∼6 per cent of artificial groups and ∼24 per cent of groups associated with binaries or triplets in real space. The final sample comprises 2209 galaxy groups covering the sky region described by Colless et al.spanning over the redshift range of  0.003 ≤ z ≤ 0.25  with a mean redshift of 0.1.     

Double-difference Relocation of the Aftershocks of the Tecom��n, Colima, Mexico Earthquake of 22 January 2003

On 22 January 2003, the M _w?=?7.6 Tecomán earthquake struck offshore of the state of Colima, Mexico, near the diffuse triple junction between the Cocos, Rivera, and North American plates. Three-hundred and fifty aftershocks of the Tecomán earthquake with magnitudes between 2.6 and 5.8, each recorded by at least 7 stations, are relocated using the double difference method. Initial locations are determined using P and S readings from the Red Sismológica Telemétrica del Estado de Colima (RESCO) and a 1-D velocity model. Because only eight RESCO stations were operating immediately following the Tecomán earthquake, uncertainties in the initial locations and depths are fairly large, with average uncertainties of 8.0?km in depth and 1.4?km in the north?Csouth and east?Cwest directions. Events occurring between 24 January and 31 January were located using not only RESCO phase readings but also additional P and S readings from 11 temporary stations. Average uncertainties decrease to 0.8?km in depth, 0.3?km in the east?Cwest direction, and 0.7?km in the north?Csouth direction for events occurring while the temporary stations were deployed. While some preliminary studies of the early aftershocks suggested that they were dominated by shallow events above the plate interface, our results place the majority of aftershocks along the plate interface, for a slab dipping between approximately 20° and 30°. This is consistent with the slab positions inferred from geodetic studies. We do see some upper plate aftershocks that may correspond to forearc fault zones, and faults inland in the upper plate, particularly among events occurring more than 3?months after the mainshock.     

Magma ascent rate and initial water concentration inferred from diffusive water loss from olivine-hosted melt inclusions

As the water concentration in magma decreases during magma ascent, olivine-hosted melt inclusions will reequilibrate with the host magma through hydrogen diffusion in olivine. Previous models showed that for a single spherical melt inclusion in the center of a spherical olivine, the rate of diffusive reequilibration depends on the partition coefficient and diffusivity of hydrogen in olivine, the radius of the melt inclusion, and the radius of the olivine. This process occurs within a few hours and must be considered when interpreting water concentration in olivine-hosted melt inclusions. A correlation is expected between water concentration and melt inclusion radius, because small melt inclusions are more rapidly reequilibrated than large ones when the other conditions are the same. This study investigates the effect of diffusive water loss in natural samples by exploring such a correlation between water concentration and melt inclusion radius, and shows that the correlation can be used to infer the initial water concentration and magma ascent rate. Raman and Fourier transform infrared spectroscopy measurements show that 31 melt inclusions (3.6–63.9 μm in radius) in six olivines from la Sommata, Vulcano Island, Aeolian Islands, have 0.93–5.28 wt% water, and the host glass has 0.17 wt% water. The water concentration in the melt inclusions shows larger variation than the data in previous studies (1.8–4.52 wt%). It correlates positively with the melt inclusion radius, but does not correlate with the major element concentrations in the melt inclusions, which is consistent with the hypothesis that the water concentration has been affected by diffusive water loss. In a simplified hypothetical scenario of magma ascent, the initial water concentration and magma ascent rate are inferred by numerical modeling of the diffusive water loss process. The melt inclusions in each olivine are assumed to have the same initial water concentration and magma ascent rate. The melt inclusions are assumed to be quenched after eruption (i.e., the diffusive water loss after eruption is not considered). The model results show that the melt inclusions initially had 3.9–5.9 wt% water and ascended at 0.002–0.021 MPa/s before eruption. The overall range of ascent rate is close to the lower limit of previous estimates on the ascent rate of basalts.     

http://www.sciencedirect.com/science/article/pii/S1674987113000029

The Preandean geological configuration of the eastern North Patagonian Massif is established through the use of geological and geophysical analysis.The positive gravity anomalies located near the Atlantic coast are due to 535 and 540 Ma old rocks belonging to the Pampean Orogeny （Precambrian-middle Cambrian）,which are widely recognized in central and northern Argentina.The Famatinian Cycle （Ordovician-Devonian） is represented by a Silurian-Devonian marine basin equivalent to those of eastern-central Argentina and South Africa,and which was deformed at the end of the Devonian byan～E-W to WNW-ESE compressional event,part of the Famatinian Orogeny.Containing strong gravity gradients,the NW-SE belt is coincident with fault zones which were originated during the Gondwanide Orogeny.This event also produced NW-SE overthrusting of the Silurian-Devonian sequences and strike-slip faults that displaced blocks in the same direction.This deformation event belongs to the Gondwanide Orogeny that includes movements related to a counterclockwise rotation of blocks in northern Patagonia.The strong negative anomalies located in the western part of the area stem from the presence of rocks of the Jurassic Ca（n）adón Asfalto basin interbedded in the Marifil Complex.These volcaniclastic sequences show mild deformation of accommodation zones in a pre-Jurassic paleorelief.     

Insights from a Qualitative Study of Rural Communes: Physical and Social Dimensions of Place

Existing research on people–place dynamics often focuses on either physical environment or social dimensions, but many theoretical models posit that both factors are important. I argue that rural communes provide unique settings to explore both social and environmental aspects of place. Based on qualitative interview data with members of two communes, I address emotional connection to the land, changing social relationships over time, and the interplay of physical and social aspects of place. Results demonstrate how the connections to land and social ties are both prevalent. However, social interaction has decreased while commitment to sharing and protecting the physical environment has endured.     

Overcoming unreported violence using place-based ambulance data: The case for mapping hotspots based on health data for crime prevention initiatives

A key concern in crime analysis is the “hidden crime” problem. Crime events unaccounted for in police records limit the external validity of official statistics and, more importantly, hinder the ability of the police to manage crime and utilize their resources effectively. The problem is exacerbated in proactive initiatives aimed at curbing violence through hotspot policing, where inaccuracies and imprecision, or, worse, no data at all, diminish prevention efforts. Previous studies have sought to overcome the data problem by juxtaposing police records with ambulance data on assault callouts and have found profound disparities. Specifically, researchers matched “crime hotspots” with “ambulance hotspots” (rather than individual events) because patient confidentiality considerations have prevented health professionals from sharing subject-level data with the police. However, health services can safely share spatial data on wider areas that do not disclose personal information. We build on this line of inquiry by analyzing data from the Thames Valley, United Kingdom, and observing spatial hotspots of different sizes. The results demonstrate that while the police and ambulance services attend to the same communities and similar types of facilities, the police are “blinded” to the location of nearly 8 out of 10 assaults. The incongruency is shown even with severe assaults, but to a lesser extent. We then simulate the reduction in injuries if the police had access to health data at different spatial levels and show that even under the most conservative set of assumptions, such an approach can prevent between 113 and 116 violent injuries each year that might otherwise require hospitalization.     

Re-Os and Mo isotope systematics of black shales from the Middle Proterozoic Velkerri and Wollogorang Formations, McArthur Basin, northern Australia

Joint application of the Mo isotope paleoredox proxy and Re-Os deposition-age geochronometer to euxinic black shales has potential for tracing the evolution of ocean redox chemistry over geological time. Here, we report new Re-Os and Mo isotope data for the Mesoproterozoic Velkerri Formation (Roper Group) and Paleoproterozoic Wollogorang Formation (Tawallah Group), McArthur Basin, northern Australia. New Re-Os ages of 1361 ± 21 Ma (2σ, n = 14, mean square of weighted deviates [MSWD] = 1.3, Model 1) and 1417 ± 29 Ma (2σ, n = 12, MSWD = 1.3, Model 1) constrain the depositional age of the Velkerri Formation and its contained biomarkers, as well as acritarchs and microfossils from the Roper Group. Black shales from the upper Velkerri Formation have high Mo abundances (105-119 ppm) and degree of pyritization [DOP] values (0.90-0.92) implying quantitative conversion of molybdate (MoO₄²⁻) to thiomolybdate (MoS₄²⁻) in overlying bottom waters. The average δ^97/95Mo (0.72 ± 0.10‰, 2σ, n = 6) of these shales is consistent with previous data, but represents a significantly more precise determination for global seawater δ^97/95Mo at 1.4 Ga. This value is lighter than present-day seawater by ∼0.85‰ and reflects expanded strongly euxinic deep ocean conditions ([H₂S]_aq > 11 μM) relative to oxic, suboxic, and weakly/intermittently euxinic ([H₂S]_aq < 11 μM) marine deposition in the 1.4 Ga oceans. Mass-balance modelling suggests Mo removal into strongly euxinic and oxic sediments may have comprised 30-70% and less than 15%, respectively, of the oceanic Mo sink at 1.4 Ga as opposed to 5% and 35% today, respectively.The Re-Os radioisotope system in organic-rich shales serves as a test for post-depositional alteration that could affect the usefulness of paleoredox tracers such as Mo stable isotopes. Re-Os isotope data for the Wollogorang Formation black shales are scattered and yield a highly imprecise date of 1359 ± 150 Ma (2σ, n = 21, MSWD = 85, Model 3). This age is younger than U-Pb zircon ages from interbedded tuffs that constrain the age of deposition at ca. 1730 Ma. In conjunction with previous petrological, geochemical, and paleomagnetic data, the Re-Os isotope data suggest hydrothermal fluid flow through the Wollogorang Formation, possibly associated with formation of the ca. 1640 Ma McArthur River Pb-Zn-Ag sedimentary exhalative deposit, resulted in post-depositional mobilization of Re and Os. Based on the degree of deviation of the Re-Os data from a 1730 Ma reference line, open-system behavior of Re and Os was greatest near the base of the black shale unit. Wollogorang Formation black shales are enriched in Mo (41-58 ppm), but are characterized by variable δ^97/95Mo (0.3-0.8‰) and DOP (0.57-0.92). The lightest δ^97/95Mo values occur near the base of the black shale unit. Based on the Re-Os systematics, hydrothermal fluids have probably overprinted the authigenic δ^97/95Mo signature in those shales. However, the heaviest δ^97/95Mo values in the Wollogorang Formation come from stratigraphically higher shales, and are similar to those observed for the Velkerri Formation, and thus may reflect seawater δ^97/95Mo at 1.73 Ga.