Introducing a hybrid radiative transfer method for smoothed particle hydrodynamics

A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods – the polytropic cooling approximation as devised by Stamatellos et al. and flux-limited diffusion. This hybrid method preserves the strengths of its individual components, while removing the need for atmosphere matching or other boundary conditions to marry optically thick and optically thin regions. The code uses a non-trivial equation of state to calculate temperatures and opacities of SPH particles, which captures the effects of H₂ dissociation, H⁰ ionization, He⁰ and He⁺ ionization, ice evaporation, dust sublimation, molecular absorption, bound-free and free–free transitions and electron scattering. The method is tested in several scenarios, including (i) the evolution of a  0.07 M_⊙  protoplanetary disc surrounding a  0.5 M_⊙  star; (ii) the collapse of a  1 M_⊙  protostellar cloud and (iii) the thermal relaxation of temperature fluctuations in a static homogeneous sphere.     

Capelinhos 1957–1958, Faial, Azores: deposits formed by an emergent surtseyan eruption

The 1957–1958 eruption of Capelinhos, Faial island, Azores, involved three periods of surtseyan, hydromagmatic activity: two in 1957 and one in 1958. Deposits from this eruption are exposed both in sea cliffs cut into the flanks of the tuff cone and more distally >1 km from the vent. Five lithofacies are identified: lithofacies I is composed of even thickness beds with laterally continuous internal stratigraphy and is interpreted to have been formed by fallout. Lithofacies II consists of beds with internally discontinuous lenses, and has sand-wave structures that increase in abundance toward the outer margins of the tuff cone. This lithofacies is interpreted as having been deposited from pyroclastic surges. Lithofacies III is composed of mantle-bedded deposits with laterally discontinuous internal stratigraphy. This lithofacies is interpreted to have been formed by hybrid processes where fallout of tephra occurred simultaneously with pyroclastic surges. In the outer flanks of the tuff cone, lithofacies III grades laterally into fallout beds of lithofacies I. Lithofacies IV consists of alternating beds of coarse ash aggregates and non-aggregated fine ash, and is particularly well developed in distal regions. Some of this facies was formed by fallout. Alternating beds also occur plastered against obstacles up to 2 km from the vent, indicating an origin from wet pyroclastic surges. The orientation of plastered tephra indicates that the surges were deflected by topography as they decelerated. The distinction between surge and fallout in distal regions is uncertain because wind-drifted fallout and decelerating surge clouds can generate similar deposits. Lithofacies V consists of scoria lapilli beds interpreted to be fallout from hawaiian-style fire-fountaining in the later stages of the eruption. Juvenile pyroclasts within hydromagmatic deposits are predominantly poorly vesicular (25–60% of clasts <30% vesicles). However, on both micro- and macroscopic scales, there is a wide range in clast vesicularity (up to 70% vesicles) indicating that, although fragmentation was predominantly hydromagmatic, vesiculation and magmatic-volatile-driven fragmentation operated simultaneously.     

The creep of potash salt rocks from Saskatchewan

Summary The results of creep tests on the Esterhazy-and the Patience-Lake-types of potash salt rocks from Saskatchewan, Canada are presented. The investigations involved over 6 years of time-dependent experiments in uniaxial compression using potash from the Rocanville and the Lanigan mines of the Potash Corporation of Saskatchewan. A creep test at a given load would last from 2 to 8 months, with most tests conducted over a 4-month period.Since the yield stress of both types of potash lies between 9 and 11 MPa, there is very little creep below 11 MPa. Between 11 and 13 MPa, creep strain production increases sharply through plastic deformation. Above about 13 MPa, however, plastic creep is dominated by brittle creep caused by microcracking. As a result, the lateral and volume creep strain curves may then display the transient and the steady-state, or all three stages of creep, while the axial strain, which is not affected by microcracking, usually attenuates for the whole duration.Two different interpretations of the results are offered. Identifying the last (the fourth) month of testing with the steady-state model, the stress dependence of the steady-state rate has been established for both rock types. A unimodal rate model for the axial strain and a bimodal model for the lateral strain are suggested. The alternate interpretation proceeds on the assumption that under 13 MPa, both the axial and the lateral strain can be modelled through the power function formulation of transient strain.     

Geochronology of basalts from the ninetyeast ridge and continental dispersion in the eastern Indian Ocean

Conventional K-Ar and ⁴⁰Ar/³⁹Ar total-fusion and incremental-heating ages from basalts recovered from the Ninetyeast Ridge confirm earlier indication from paleontology of basal sediments that basement crystallization ages increase systematically from south to north. A rate of migration of volcanism of 9.4 ± 0.3 cm/yr best fits the age-distance relationship determined from these basalts. The geometry, distribution of ages, paleomagnetism and geochemistry are most compatible with an origin of the Ninetyeast Ridge by northward movement of the Indian plate over a hotspot near Heard Island in the southern Indian Ocean, from mid-Cretaceous to Early Oligocene time. The Ninetyeast Ridge and nearly parallel Chagos-Laccadive Ridge provide an absolute frame of reference for the reconstruction of the eastern Indian Ocean.     

Morphologic and stratigraphic evolution of muddy ebb-tidal deltas along a subsiding coast: Barataria Bay, Mississippi River delta

The Barataria barrier coast formed between two major distributaries of the Mississippi River delta: the Plaquemines deltaic headland to the east and the Lafourche deltaic headland to the west. Rapid relative sea‐level rise (1·03 cm year^?1) and other erosional processes within Barataria Bay have led to substantial increases in the area of open water (> 775 km² since 1956) and the attendant bay tidal prism. Historically, the increase in tidal discharge at inlets has produced larger channel cross‐sections and prograding ebb‐tidal deltas. For example, the ebb delta at Barataria Pass has built seaward > 2·2 km since the 1880s. Shoreline erosion and an increasing bay tidal prism also facilitated the formation of new inlets. Four major lithofacies characterize the Barataria coast ebb‐tidal deltas and associated sedimentary environments. These include a proximal delta facies composed of massive to laminated, fine grey‐brown to pale yellow sand and a distal delta facies consisting of thinly laminated, grey to pale yellow sand and silty sand with mud layers. The higher energy proximal delta deposits contain a greater percentage of sand (75–100%) compared with the distal delta sediments (60–80%). Associated sedimentary units include a nearshore facies consisting of horizontally laminated, fine to very fine grey sand with mud layers and an offshore facies that is composed of grey to dark grey, laminated sandy silt to silty clay. All facies coarsen upwards except the offshore facies, which fines upwards. An evolutionary model is presented for the stratigraphic development of the ebb‐tidal deltas in a regime of increasing tidal energy resulting from coastal land loss and tidal prism growth. Ebb‐tidal delta facies prograde over nearshore sediments, which interfinger with offshore facies. The seaward decrease in tidal current velocity of the ebb discharge produces a gradational contact between proximal and distal tidal delta facies. As the tidal discharge increases and the inlet grows in dimensions, the proximal and distal tidal delta facies prograde seawards. Owing to the relatively low gradient of the inner continental shelf, the ebb‐tidal delta lithosome is presently no more than 5 m thick and is generally only 2–3 m in thickness. The ebb delta sediment is sourced from deepening of the inlet and the associated channels and from the longshore sediment transport system. The final stage in the model envisages erosion and segmentation of the barrier chain, leading to a decrease in tidal discharge through the former major inlets. This process ultimately results in fine‐grained sedimentation seaward of the inlets and the encasement of the ebb‐tidal delta lithosome in mud. The ebb‐tidal deltas along the Barataria coast are distinguished from most other ebb deltas along sand‐rich coasts by their muddy content and lack of large‐scale stratification produced by channel cut‐and‐fills and bar migration.     

Radar facies of paraglacial barrier systems: coastal New England, USA

Analysis of a large data base of ground-penetrating-radar (GPR) profiles from both natural and developed paraglacial barriers along the coast of New England has allowed identification of eight reflection configurations that characterize this type of mid- to high-latitude coastal environment. Bedrock anchor points yield primarily hyperbolic configurations, whereas glacial anchor points and sediment-source areas are characterized by chaotic, parallel, and tangential-oblique configurations. Beaches and dunes produce predominantly sigmoidal oblique, hummocky, reflection-free, and bounding-surface configurations. Back-barrier sediments may yield basin-fill configurations, but generally include abundant signal-attenuating units.
The GPR data, calibrated with information from cores, were collected across swash-aligned and drift-aligned barriers in a variety of wave- and tidal-energy settings. Application of a 120-MHz antenna, as used in this study, enables portrayal of a range of sedimentary units, from individual bedforms (on single records) to entire barrier elements (using large numbers of intersecting GPR sections), at maximum vertical resolutions that vary between 0·2 m and 0·7 m.
The most important drawback of GPR in the coastal environment is attenuation of the electromagnetic (EM) signal by layers of salt-marsh peat or by brackish or salty groundwater, primarily along barrier edges. This disadvantage is offset by many benefits. Data can be collected at rates of several km per day, making GPR an excellent reconnaissance tool. A core that is used in the calibration of GPR data can be matched with great accuracy to its position on the complementary GPR record, allowing detailed correlation between lithostratigraphy and reflection configuration.     

Comments on the development of rifts and transform faults during continental breakup; examples from the Gulf of Aden and northern Red Sea

A number of basins are observed to extend inland from the coasts on both sides of the Gulf of Aden. The basins are orientated at approximately right angles to the spreading direction and intersect the coasts at the meeting of sheared and rifted continental margins. They appear to be grabens, one wall of which is continuous with the half graben of the neighbouring rifted margin. It is suggested that these were once parts of a number of discrete rifts arranged en-echelon along a zone of lithospheric weakness during the early opening of the Gulf of Aden, which became redundant when transform faults formed. The proposed development of rifts and transform faults is similar to that of a spreading centre, transform fault, spreading centre pattern developed in the freezing wax model of Oldenburg and Brune (1975). The Gulf of Suez at the northern end of the Red Sea is interpreted in a similar way since it has a number of features in common with the basins in the continents adjacent to the Gulf of Aden.     

Manganese chemistry in rivers and streams

The speciation of Mn has been determined in 15 rivers and streams representing a wide variety of physico-chemical conditions. Using the technique of anodic stripping voltammetry (asv), specific for reduced Mn(II) species, it is found that a major part of the <0.015 μm Mn size fraction is present in a reduced Mn(II), asv-labile, form. In some waters there is also a significant asv inactive Mn fraction considered to be present as a ‘small colloidal’ species. The soluble (<0.015 μm) Mn fraction represents 15–95% of total Mn and does not appear to be dependent upon pH, alkalinity, specific conductance or humic substance concentration in the water. It is argued that under the dynamic, short residence time, conditions that apply in most rivers the paniculate and soluble Mn fractions are decoupled, their respective presence being dependent principally upon the catchment hydrogeological conditions. This contrasts with a previously held view that the paniculate phase is coupled to the dissolved phase by the pH dependent oxidation of dissolved Mn(II) to highly insoluble Mn(IV) species (Graham et al., 1976). Consideration of manganese speciation in waters which were incubated for five months showed that pH becomes the controlling factor when equilibrium is approached.     

Reconstructing the ancestral Yellowstone plume from accreted seamounts and its relationship to flat-slab subduction

Recent geodynamic analyses have emphasized the relationship between modern flat-slab subduction zones and the overriding of buoyant oceanic crust. Although most models for the evolution of the Late Mesozoic–Cenozoic Laramide orogeny in the southwestern United States involve flat-slab subduction, the mechanisms proposed are controversial. An examination of the geological evolution of the 60–50-Ma Crescent terrane of the Coast Ranges indicates that it was formed in a shallowing-upward Loihi-type oceanic setting culminating in the eruption of subaerial lavas. Plate reconstructions indicate that the Crescent terrane was emplaced into ca. 20-Ma crust, and the presence of subaerial lavas implies an uplift due to the plume of ca. 4.2 km, which we use to calculate a minimum buoyancy flux of 1.1 Mg s⁻¹, similar to that of the modern Yellowstone plume.Published paleomagnetic data indicate that the Crescent terrane was formed at a paleolatitude similar to that of the Yellowstone plume. The Crescent seamount was accreted within 5 My of the cessation of plume magmatism. Plate reconstructions indicate that it would have originated about 750 km to the west of the North American plate margin if it developed above a fixed Yellowstone plume, and are therefore consistent with the recorded very short interval between its formation and tectonic emplacement.We interpret the Crescent terrane as due to the ancestral Yellowstone plume. Such a plume would have generated an elongate swell and related plateau that would have been overridden by the North American margin. Taken together, the relationship between flat-slab and overriding of oceanic plateau in Laramide times would have been analogous to the relationship between modern Andean flat-slab subduction zones and the Juan Fernandez and Nazca oceanic plateaus.