El Chichón volcano, April 4, 1982: volcanic cloud history and fine ash fallout

This retrospective study focuses on the fine silicate particles (<62 µm in diameter) produced in a large eruption that was otherwise well studied. Fine particles represent a potential hazard to aircraft, because as simple particles they have very low terminal velocities and could potentially stay aloft for weeks. New data were collected to describe the fine particle size distributions of distal fallout samples collected soon after eruption. Although, about half of the mass of silicate particles produced in this eruption of ~1 km³ dense rock equivalent magma were finer than 62 µm in diameter, and although these particles were in a stratospheric cloud after eruption, almost all of these fine particles fell to the ground near (<300 km) the volcano in a day or two. Particles falling out from 70 to 300 km from the volcano are mostly <62 µm in diameter. The most plausible explanation for rapid fallout is that the fine ash nucleates ice in the convective cloud and initiates a process of meteorological precipitation that efficiently removes fine silicates. These observations are similar to other eruptions and we conclude that ice formation in convective volcanic clouds is part of an effective fine ash removal process that affects all or most volcanic clouds. The existence of pyroclastic flows and surges in the El Chichón eruption increased the overall proportion of fine silicates, probably by milling larger glassy pyroclasts.     

Atmospheric circulations of terrestrial planets orbiting low-mass stars

Circulations and habitable zones of planets orbiting low-mass stars are investigated. Many of these planets are expected to rotate synchronously relative to their parent stars, thereby raising questions about their surface temperature distributions and habitability. We use a global circulation model to study idealized, synchronously rotating (tidally locked) planets of various rotation periods, with surfaces of all land or all water, but with an Earth-like atmosphere and solar insolation. The dry planets exhibit wide variations in surface temperature: >80 °C on the dayside to <−110 °C on the nightside for the 240-h rotator, for example. The water-covered aquaplanets are warmer and exhibit narrower ranges of surface temperatures, e.g., ∼40 °C to >−60 °C for the 240-h orbiter. They also have a larger habitable area, defined here as the region where average surface temperatures are between 0 °C and 50 °C. This concept has little relevance for either dry or aquaplanets, but might become relevant on a planet with both land area and oceans.The circulations on these tidally locked planets exhibit systematic changes as the rotation period is varied. However, they also reveal abrupt transitions between two different circulation regimes and multiple equilibria. For the dry planet, the transition occurs between a 4-day and a 5-day period, while for the aquaplanet, it occurs between a 3-day and a 4-day period. For both dry and aqua planets, this transition occurs when the Rossby deformation radius exceeds half the planetary radius. Further investigation on the dry planet reveals that multiple equilibria exist between 100- and 221-h periods. These multiple equilibria may be relevant for real planets within the habitable zones of late K and M stars, because these planets are expected to have rotation periods between 8 and 100 Earth days.     

Coupled convection and tidal dissipation in Europa’s ice shell using non-Newtonian grain-size-sensitive (GSS) creep rheology

We present self-consistent, fully coupled two-dimensional (2D) numerical models of thermal evolution and tidal heating to investigate how convection interacts with tidal dissipation under the influence of non-Newtonian grain-size-sensitive creep rheology (plausibly resulting from grain boundary sliding) in Europa’s ice shell. To determine the thermal evolution, we solved the convection equations (using finite-element code ConMan) with the tidal dissipation as a heat source. For a given heterogeneous temperature field at a given time, we determined the tidal dissipation rate throughout the ice shell by solving for the tidal stresses and strains subject to Maxwell viscoelastic rheology (using finite-element code Tekton). In this way, the convection and tidal heating are fully coupled and evolve together. Our simulations show that the tidal dissipation rate can have a strong impact on the onset of thermal convection in Europa’s ice shell under non-Newtonian GSS rheology. By varying the ice grain size (1-10 mm), ice-shell thickness (20-120 km), and tidal-strain amplitude (0-4 × 10⁻⁵), we study the interrelationship of convection and conduction regimes in Europa’s ice shell. Under non-Newtonian grain-size-sensitive creep rheology and ice grain size larger than 1 mm, no thermal convection can initiate in Europa’s ice shell (for thicknesses <100 km) without tidal dissipation. However, thermal convection can start in thinner ice shells under the influence of tidal dissipation. The required tidal-strain amplitude for convection to occur decreases as the ice-shell thickness increases. For grain sizes of 1-10 mm, convection can occur in ice shells as thin as 20-40 km with the estimated tidal-strain amplitude of 2 × 10⁻⁵ on Europa.     

Measurement of the radius of Mercury by radio occultation during the MESSENGER flybys

The MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft completed three flybys of Mercury in 2008–2009. During the first and third of those flybys, MESSENGER passed behind the planet from the perspective of Earth, occulting the radio-frequency (RF) transmissions. The occultation start and end times, recovered with 0.1 s accuracy or better by fitting edge-diffraction patterns to the RF power history, are used to estimate Mercury's radius at the tangent point of the RF path. To relate the measured radius to the planet shape, we evaluate local topography using images to identify the high-elevation feature that defines the RF path or using altimeter data to quantify surface roughness. Radius measurements are accurate to 150 m, and uncertainty in the average radius of the surrounding terrain, after adjustments are made from the local high at the tangent point of the RF path, is 350 m. The results are consistent with Mercury's equatorial shape as inferred from observations by the Mercury Laser Altimeter and ground-based radar. The three independent estimates of radius from occultation events collectively yield a mean radius for Mercury of 2439.2±0.5 km.     

Saturn’s tropospheric composition and clouds from Cassini/VIMS 4.6-5.1 μm nightside spectroscopy

The latitudinal variation of Saturn’s tropospheric composition (NH₃, PH₃ and AsH₃) and aerosol properties (cloud altitudes and opacities) are derived from Cassini/VIMS 4.6-5.1 μm thermal emission spectroscopy on the planet’s nightside (April 22, 2006). The gaseous and aerosol distributions are used to trace atmospheric circulation and chemistry within and below Saturn’s cloud decks (in the 1- to 4-bar region). Extensive testing of VIMS spectral models is used to assess and minimise the effects of degeneracies between retrieved variables and sensitivity to the choice of aerosol properties. Best fits indicate cloud opacity in two regimes: (a) a compact cloud deck centred in the 2.5-2.8 bar region, symmetric between the northern and southern hemispheres, with small-scale opacity variations responsible for numerous narrow light/dark axisymmetric lanes; and (b) a hemispherically asymmetric population of aerosols at pressures less than 1.4 bar (whose exact altitude and vertical structure is not constrained by nightside spectra) which is 1.5-2.0× more opaque in the summer hemisphere than in the north and shows an equatorial maximum between ±10° (planetocentric).Saturn’s NH₃ spatial variability shows significant enhancement by vertical advection within ±5° of the equator and in axisymmetric bands at 23-25°S and 42-47°N. The latter is consistent with extratropical upwelling in a dark band on the poleward side of the prograde jet at 41°N (planetocentric). PH₃ dominates the morphology of the VIMS spectrum, and high-altitude PH₃ at p < 1.3 bar has an equatorial maximum and a mid-latitude asymmetry (elevated in the summer hemisphere), whereas deep PH₃ is latitudinally-uniform with off-equatorial maxima near ±10°. The spatial distribution of AsH₃ shows similar off-equatorial maxima at ±7° with a global abundance of 2-3 ppb. VIMS appears to be sensitive to both (i) an upper tropospheric circulation (sensed by NH₃ and upper-tropospheric PH₃ and hazes) and (ii) a lower tropospheric circulation (sensed by deep PH₃, AsH₃ and the lower cloud deck).     

Recycling of Coastal Dredged Sediments from the Northern Nile Delta,Egypt, for Heavy Minerals Exploitation

Mineralogical and textural analyses of 45 undisturbed short cores and 80 grab sediment samples, collected from five frequently dredged navigational areas within harbors and water pathways of the Nile delta littoral system, were utilized for evaluation of these sediments as potential source of economic heavy minerals (EHMs). Results of mineralogical characterization indicate that the average total heavy mineral (HM) concentrations are as follows: Abu Qir Bay (1.7%), Rosetta estuary (3.1%), Burullus fishing port (4.5%), Damietta Harbor (2.9%), and El Gamil lagoon inlet (1.9%). Assessment of HM grades indicates predominance of magnetite, ilmenite, hematite, leucoxene, garnet, zircon, and rutile. Results of the feasibility analysis indicate that dredged sediments at these study areas are considered as a potential source of EHMs and economically promising to be mined for HMs. The present study suggests a practical operative plan of two successive phases for HMs recycling: (1) in-situ initial separation of HMs on the dredger deck using wet-gravity spirals, then (2) transportation of the recovered HM concentrates to an onshore processing plant to selectively separate individual HMs via wet and dry magnetic and electrostatic separators. Alternatively, dredged sediments can be directly pumped to a nearby onshore area as stockpile to be recycled afterward in inland processing plant.     

Log‐jams and avulsions in the San Antonio River Delta,Texas

The San Antonio River Delta (SARD), Texas, has experienced two major avulsions in the past 80 years, and a number of other historical and Holocene channel shifts. The causes and consequences of these avulsions – one of which is ongoing – were examined using a combination of fieldwork, geographic information system (GIS) analysis, and historical information to identify active, semi‐active, and paleochannels and the sequence of shifting flow paths through the delta. The role of deposition patterns and antecedent morphology, large woody debris jams, and tectonic influences were given special attention. Sedimentation in the SARD is exacerbated by tectonic effects. Channel aggradation is ubiquitous, and superelevation of the channel bed above the level of backswamp areas on the floodplain is common. This creates ideal setup conditions for avulsions, and stable, cohesive fine‐grained banks favor avulsions rather than lateral migration. Flood basins between the alluvial ridges associated with the aggraded channels exist, but avulsions occur by re‐occupation of former channels found within or connected to the flood basins. Large woody debris and channel‐blocking log‐jams are common, and sometimes displace flow from the channel, triggering crevasses. However, a large, recurring log‐jam at the site of the ongoing avulsion from the San Antonio River into Elm Bayou is not responsible for the channel shift. Rather, narrow, laterally stable channels resulting from flow splits lead to accumulation of wood. Some aspects of the SARD avulsion regime are typical of other deltas, while others are more novel. These includes avulsions involving tributaries and subchannels within the delta as well as from the dominant channel; tectonic influences on delta backstepping and on channel changes within the delta; avulsions as an indirect trigger for log‐jam formation (as well as vice‐versa); and maintenance of a multi‐channel flow pattern distinct from classic anastamosing or distributary systems. Copyright © 2012 John Wiley & Sons, Ltd.     

Geographical information system approaches for hazard mapping of dilute lahars on Montserrat, West Indies

Many research tools for lahar hazard assessment have proved wholly unsuitable for practical application to an active volcanic system where field measurements are challenging to obtain. Two simple routing models, with minimal data demands and implemented in a geographical information system (GIS), were applied to dilute lahars originating from Soufrière Hills Volcano, Montserrat. Single-direction flow routing by path of steepest descent, commonly used for simulating normal stream-flow, was tested against LAHARZ, an established lahar model calibrated for debris flows, for ability to replicate the main flow routes. Comparing the ways in which these models capture observed changes, and how the different modelled paths deviate can also provide an indication of where dilute lahars, do not follow behaviour expected from single-phase flow models. Data were collected over two field seasons and provide (1) an overview of gross morphological change after one rainy season, (2) details of dominant channels at the time of measurement, and (3) order of magnitude estimates of individual flow volumes. Modelling results suggested both GIS-based predictive tools had associated benefits. Dominant flow routes observed in the field were generally well-predicted using the hydrological approach with a consideration of elevation error, while LAHARZ was comparatively more successful at mapping lahar dispersion and was better suited to long-term hazard assessment. This research suggests that end-member models can have utility for first-order dilute lahar hazard mapping.     

Biofilm Detection in a Model Well‐Bore Environment Using Low‐Field NMR

This research addresses the challenges of the lack of non‐invasive methods and poor spatiotemporal resolution associated with monitoring biogeochemical activity central to bioremediation of subsurface contaminants. Remediation efforts often include growth of biofilm to contain or degrade chemical contaminants, such as nitrates, hydrocarbons, heavy metals, and some chlorinated solvents. Previous research indicates that nuclear magnetic resonance (NMR) is sensitive to the biogeochemical processes of biofilm accumulation. The current research focuses on developing methods to use low‐cost NMR technology to support in situ monitoring of biofilm growth and geochemical remediation processes in the subsurface. Biofilm was grown in a lab‐scale radial flow bioreactor designed to model the near wellbore subsurface environment. The Vista Clara Javelin NMR logging device, a slim down‐the‐borehole probe, collected NMR measurements over the course of eight days while biofilm was cultivated in the sand‐packed reactor. Measured NMR mean log T₂ relaxation times decreased from approximately 710 to 389 ms, indicating that the pore environment and bulk fluid properties were changing due to biofilm growth. Destructive sampling employing drop plate microbial population analysis and scanning electron and stereoscopic microscopy confirmed biofilm formation. Our findings demonstrate that the NMR logging tool can detect small to moderate changes in T₂ distribution associated with environmentally relevant quantities of biofilm in quartz sand.