A co-accretional model of satellite formation

Numerical calculation of a simple accretion model including the effects of tidal friction indicate that coformation is tenable only if the planet's Q is less than about 10³. The parameter which most strongly affects the final mass ratio of the pair is the time at which the secondary embryo is introduced. Our model yields the proper Moon-Earth mass ratio if the Moon embryo is introduced when the Earth is only about $\text{1}\text{10}$ of its final mass. The lunar orbit remains at about 10 Earth radii throughout most of the growth.This model of satellite formation overcomes two difficulties of the “circumterrestrial cloud” model of Ruskol (1960, 1963, 1972): (1) The difficulty of accumulating a mass as great as the entire Moon before gravitational instability reduces the cloud to a small number of moonlets is removed. (2) The differences between terrestrial and outer planet satellite systems is easily understood in terms of the differences in Q between these planets. The high Q of the outer planets does not allow a satellite embryo to survive a significant portion of the accretion process, thus only small bodies which formed very late in the accumulation of the planet remain as satellites. The low Q of the terrestrial planets allows satellite embryos of these planets to survive during accretion, thus massive satellites such as the Earth's Moon are expected. The present lack of such satellites of the other terrestrial planets may be the result of tidal evolution, either infall following primary despinning (Burns, 1973) or escape due to increase in orbit eccentricity.     

30 keV (n, γ) cross-sections from the nuclear statistical model

Global input parameters for a Hauser-Feshbach calculation of astrophysical nuclear reaction rates are evaluated; wherever possible they are expressed as systematic functions of nucleon numbersN andZ rather than mass numberA. Average (n, ) cross-sections at 30 keV are calculated for all stable nuclei and nuclei with lifetimes 1 yr, taking into account the effects of width fluctuations. Agreement with experimental results is generally obtained to within 70%, but systematic discrepancies arise for even-even target nuclei between neutron shell closuresN=50 andN=82. A N curve for thes-process nuclei is in good agreement with theoretical predictions, except for a possible anomalous peak aroundA=180.     

Carbon flux in the northwest Mediterranean estimated from microbial production

Spring profiles of microbial production derived from the dark incorporation of tritiated leucine and tritiated thymidine in the northwest Mediterranean show an exponential decline with depth. Assuming this to represent a steady-state balance between microbial respiration and the downward flux of carbon, the downward flux is estimated as (1−/)p/b, where p is the microbial production, their gross growth efficiency and b the coefficient of exponential decline with depth. Summer profiles, ranging over about 3° of latitude and 4° of longitude, were well fitted by a two-component exponential decline, suggesting two distinct microbial substrates. Values of b for the more rapidly declining component varied between 0.01 and 0.06 m⁻¹ according to location. In the case of the slower component, b was estimated as 0.002 m⁻¹, and did not vary significantly over the region. Estimated fluxes of carbon at the surface are 123–335 mg m⁻² d⁻¹ for the fast and 95 mg m⁻² d⁻¹ for the slow component. Below about 200 m, carbon flux is dominated by the slow component. Flux estimates are compatible with flux estimates from sediment traps in the same region. The observed changes between the spring and summer profiles, combined with the horizontal homogeneity of the summer profiles below 200 m, are consistent with a downward transport of about 5–10 m d^–1, implying a significant dispersive component to the observed fluxes.     

Delineation of intrinsic geological units

Conventional methods for the estimation of mineral resources potential are broadly subdivided into two groups: subjective models based upon ore genetic theories and objective models based upon multiple geological observations. This paper examines some objective means that brings the two types of model together and jointly utilizes the information from both ore genesis and geological observations in delineation of promising areas for mineral exploration. These promising areas are here referred to as intrinsic geological units. They are so-called because they are not identified directly with respect to mineral deposits, but delineated in terms of some critical genetic factor, a necessary condition for formation of the mineral deposits. Four major steps are involved in this analysis. First, establish the structure of critical genetic factor-recognition criteria-multiple geodata for the deposit type of interest. Next, estimate probabilities of occurrence of each of the recognition criteria based upon geological observations. Third, a synthesized probability measure for occurrence of critical genetic factor is determined as an optimum linear combination of the probabilities estimated in the second stage. Finally, intrinsic geological units are delineated by optimally discretizing the probability measure for the critical genetic factor. This estimation procedure is formulated into a three-stage optimization model: logistic regression-quadratic programming-optimum discretization, which is demonstrated with a case study on epithermal Au-Ag mineral deposits in the Walker Lake quadrangle of California and Nevada. In this case study, heat source is selected as the critical genetic factor and it is identified by three selected recognition criteria: Au-Ag mineral occurrence, Tertiary intrusive, and hydrothermal alteration. The intrinsic geological units for epithermal Au-Ag deposits are delineated on the basis of six different preprocessed and integrated geological, geochemical, and geophysical fields.     

The morphology and evolution of the Stromboli 2002–2003 lava flow field: an example of a basaltic flow field emplaced on a steep slope

The use of a hand-held thermal camera during the 2002–2003 Stromboli effusive eruption proved essential in tracking the development of flow field structures and in measuring related eruption parameters, such as the number of active vents and flow lengths. The steep underlying slope on which the flow field was emplaced resulted in a characteristic flow field morphology. This comprised a proximal shield, where flow stacking and inflation caused piling up of lava on the relatively flat ground of the vent zone, that fed a medial–distal lava flow field. This zone was characterized by the formation of lava tubes and tumuli forming a complex network of tumuli and flows linked by tubes. Most of the flow field was emplaced on extremely steep slopes and this had two effects. It caused flows to slide, as well as flow, and flow fronts to fail frequently, persistent flow front crumbling resulted in the production of an extensive debris field. Channel-fed flows were also characterized by development of excavated debris levees in this zone (Calvari et al. 2005). Collapse of lava flow fronts and inflation of the upper proximal lava shield made volume calculation very difficult. Comparison of the final field volume with that expecta by integrating the lava effusion rates through time suggests a loss of ~70% erupted lava by flow front crumbling and accumulation as debris flows below sea level. Derived relationships between effusion rate, flow length, and number of active vents showed systematic and correlated variations with time where spreading of volume between numerous flows caused an otherwise good correlation between effusion rate, flow length to break down. Observations collected during this eruption are useful in helping to understand lava flow processes on steep slopes, as well as in interpreting old lava–debris sequences found in other steep-sided volcanoes subject to effusive activity.     

Towards quantifying uncertainty in transient climate change

Ensembles of coupled atmosphere–ocean global circulation model simulations are required to make probabilistic predictions of future climate change. “Perturbed physics” ensembles provide a new approach in which modelling uncertainties are sampled systematically by perturbing uncertain parameters. The aim is to provide a basis for probabilistic predictions in which the impact of prior assumptions and observational constraints can be clearly distinguished. Here we report on the first perturbed physics coupled atmosphere–ocean model ensemble in which poorly constrained atmosphere, land and sea-ice component parameters are varied in the third version of the Hadley Centre model (the variation of ocean parameters will be the subject of future study). Flux adjustments are employed, both to reduce regional sea surface temperature (SST) and salinity biases and also to admit the use of combinations of model parameter values which give non-zero values for the global radiation balance. This improves the extent to which the ensemble provides a credible basis for the quantification of uncertainties in climate change, especially at a regional level. However, this particular implementation of flux-adjustments leads to a weakening of the Atlantic overturning circulation, resulting in the development of biases in SST and sea ice in the North Atlantic and Arctic Oceans. Nevertheless, model versions are produced which are of similar quality to the unperturbed and un-flux-adjusted version. The ensemble is used to simulate pre-industrial conditions and a simple scenario of a 1% per year compounded increase in CO₂. The range of transient climate response (the 20 year averaged global warming at the time of CO₂ doubling) is 1.5–2.6°C, similar to that found in multi-model studies. Measures of global and large scale climate change from the coupled models show simple relationships with associated measures computed from atmosphere-mixed-layer-ocean climate change experiments, suggesting that recent advances in computing the probability density function of climate change under equilibrium conditions using the perturbed physics approach may be extended to the transient case.     

The changing morphology of an open lava channel on Mt. Etna

An open channel lava flow on Mt. Etna (Sicily) was observed during May 30–31, 2001. Data collected using a forward looking infrared (FLIR) thermal camera and a Minolta-Land Cyclops 300 thermal infrared thermometer showed that the bulk volume flux of lava flowing in the channel varied greatly over time. Cyclic changes in the channel's volumetric flow rate occurred over several hours, with cycle durations of 113–190 min, and discharges peaking at 0.7 m³ s⁻¹ and waning to 0.1 m³ s⁻¹. Each cycle was characterized by a relatively short, high-volume flux phase during which a pulse of lava, with a well-defined flow front, would propagate down-channel, followed by a period of waning flow during which volume flux lowered. Pulses involved lava moving at relatively high velocities (up to 0.29 m s⁻¹) and were related to some change in the flow conditions occurring up-channel, possibly at the vent. They implied either a change in the dense rock effusion rate at the source vent and/or cyclic-variation in the vesicle content of the lava changing its bulk volume flux. Pulses would generally overspill the channel to emplace pāhoehoe overflows. During periods of waning flow, velocities fell to 0.05 m s^–1. Blockages forming during such phases caused lava to back up. Occasionally backup resulted in overflows of slow moving ‘a‘ā that would advance a few tens of meters down the levee flank. Compound levees were thus a symptom of unsteady flow, where overflow levees were emplaced as relatively fast moving pāhoehoe sheets during pulses, and as slow-moving ‘a‘ā units during backup. Small, localized fluctuations in channel volume flux also occurred on timescales of minutes. Volumes of lava backed up behind blockages that formed at constrictions in the channel. Blockage collapse and/or enhanced flow under/around the blockage would then feed short-lived, wave-like, down-channel surges. Real fluctuations in channel volume flux, due to pulses and surges, can lead to significant errors in effusion rate calculations. Editorial responsibility: A. Woods     

An analysis of mechanisms of ice-wedge casting based on geotechnical centrifuge simulations

Here we interpret the outcomes of scaled geotechnical centrifuge simulation of ice-wedge casting in terms of the likely significance of Quaternary ice-wedge pseudomorphs observed within different host sediments. Six experiments were completed in which 1/30th scale models of an ice-wedge embedded within frozen host sediments beneath a 25 mm thick unfrozen active layer were thawed from the surface downwards in the geotechnical centrifuge under a 30 times gravity (30×g) acceleration. Host sediment granulometry and/or ice contents were varied in each model, with host materials comprising medium sands, fine sandy silts and silty clays. The model ice-wedge was 50 mm at the top, 150 mm deep, and extended across the full width of the 450 mm wide test box. Centrifuge scaling laws indicate that under an acceleration of 30×g, stress distribution was equivalent to a 13.5 m long section of a 4.5 m high and 1.5 m wide full-scale prototype ice-wedge, covered in an active layer of thickness equivalent to 0.75 m. Thermal regimes, measured pore pressures during thaw, observed thaw consolidation and measured host sediment geotechnical properties are utilised in the interpretation of casting mechanisms. During a single uniform thaw event it is shown that arching of infilling sediment and the formation of a void is likely if negative pore pressures are developed in the host sediment. In fine silt and clays high ice contents are more likely than in sands, thaw consolidation is greater, positive pore pressures encourage complete filling of the ice-wedge void, and soft sediment deformation is likely to cause deformation of the cast and reduce its width and depth. Though natural casting mechanisms are likely to be more complex than those simulated here, modelling experiments highlight the need for care when inferring original ice-wedge geometry from observed shape and size of Quaternary ice-wedge casts.     

New insights into the occurrence of C-depleted carbon in the mantle from two closely associated kimberlites: Letlhakane and Orapa, Botswana

Carbon isotope measurements on diamonds from the Letlhakane kimberlite, and the analyses of their inclusions, permit the examination of km-scale mantle-composition variations by comparing the results with those for the nearby Orapa kimberlite. Diamonds from Letlhakane have a wide range in carbon isotopic composition (−3‰ to −21‰); however, the relative abundance of diamonds depleted in ¹³C is significantly lower than in the Orapa kimberlite. Most of the ¹³C-depleted diamonds belong to the eclogictic or websteritic paragenesis. The relative abundance of inclusions in diamonds and their composition indicate that there are significant differences in petrology in the mantle below the two locations. At Letlhakane, peridotitic compositions are more prevalent than at Orapa and the protolith of P-Type inclusions in diamonds may have experienced a higher degree of partial melting at Letlhakane compared to Orapa. P/T estimates for both W- and E-Type diamonds indicate that a region of ¹³C-depletion may exist beneath the two kimberlites. The relationships between carbon isotopic composition of the host diamond and the Al₂O₃/Cr₂O₃ ratios of their websteritic and eclogitic garnet inclusions indicate that the low δ¹³C regions may represent a primary mantle feature, unrelated to a crustal component.