Across-arc variation of lava chemistry in the Izu-Bonin Arc: identification of subduction components

In order to understand the role of the subducted lithosphere in producing the geochemical characteristics of arc magmas, major- and trace-element along with Sr- and Nd-isotope compositions have been determined for Quaternary volcanic rocks from the Izu-Bonin intra-oceanic arc. ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd ratios decrease away from the volcanic front of this arc and lie on mixing lines between the assumed isotopic compositions of fluid phases mainly derived from the basalt layer of the subducted lithosphere and upper-mantle materials in the sub-arc wedge. This across-arc variation can be explained through a simple sequence of processes involving initial release of fluid phases from the subducted oceanic crust to produce hydrous peridotite at the base of the mantle wedge. This hydrous peridotite is dragged downward with the slab and releases a second-stage metasomatizing fluid beneath the volcanic arc. The higher concentrations of both Sr and Nd in the fluid beneath the volcanic front than those beneath the back-arc side may be a possible cause of the observed across-arc variation in Sr-Nd isotopic ratios. The difference in compositions of fluid phases is attributed to the different hydrous phases which decompose in the hydrous peridotite layer; amphibole beneath the volcanic front and phlogopite beneath the back-arc side of the volcanic arc. The mineralogically controlled fluid addition may also be responsible for the across-arc variation in Rb/K and Rb/Zr ratios, increasing away from the volcanic front.     

A posteriori local error estimation and adaptive time-stepping for newmark integration in dynamic analysis

A simple a posteriori local error estimate for Newmark time integration schemes in dynamic analysis is presented, based on the concept of a so called ‘post-processing’ technique. In conjunction with the error estimate, an adaptive time-stepping algorithm is described, which adjusts the time step size so that the local error of each time step is within a prescribed error tolerance. Numerical examples given in the paper indicate that the error estimate is asymptotically convergent, computationally efficient and convenient, and the adaptive time-stepping scheme can predict a nearly optimal step size from time to time, thus making the numerical solution reliable in an efficient manner.     

Dynamic interaction analysis of a circular cylindrical shell of finite length in a half-space

A study on the transient response of a circular cylindrical shell of finite length embedded in a homogeneous, isotropic and linear elastic half-space is presented. The soil-structure system is subjected to suddenly applied explosion waves. The numerical method employed is a combination of the time domain semi-analytical boundary element method used for the semi-infinite soil medium and the finite strip method used for the circular cylindrical shell. The two methods are combined through equilibrium and compatibility conditions at the soil-structure interface. The dynamic responses at the interface between the soil medium and the structure for every time step are obtained. Numerical examples are presented in detail to demonstrate the use and versatility of the proposed method. The following parameters are found to affect the response: (1) the slenderness ratio of the length over the diameter of the shell, L/D; (2) the relative wall thickness, h/a; (3) the relative stiffness ratio between the shell and the medium, E_s/E_m; and (4) the incidence angle of the explosion wave, α.     

Lateral migration,thermoluminescence chronology and colour variation of longitudinal dunes near Birdsville in the Simpson Desert,central Australia

The longitudinal dunes of the Simpson Desert, in the vicinity of Birdsville, have been reworked largely during the Holocene from dunes deposited up to 80000 years ago or earlier. The widespread asymmetry of these roughly northward-trending dunes, with steeper eastern faces and more gentle western faces, supports wind-rose data showing sand-transporting winds from the southwest obliquely intersecting the dunes. While this suggests a change in the wind pattern since the dune field was oriented, it does not indicate that the dunes are necessarily shifting leeward (eastward) as a consequence. It is hypothesized that the direction of migration is controlled by the extent to which the dunes are vegetated. Relatively well-vegetated dunes can accrete sand on their gentle stoss slopes and erode on their lee slopes causing them to shift westward and hence obliquely into the wind, a condition that probably prevails in wetter regions and during episodes of relatively humid climate. In contrast, in very dry areas or during arid phases, sand can move unimpeded up a sparsely vegetated stoss face and over the crest to form an avalanche or slip face on the lee side, thereby causing the dunes to shift eastward. Despite evidence that longitudinal dune crests can shift laterally to some extent, the dunes in the western part of the Simpson Desert have not migrated, either westward or eastward, more than 100m or so from their Pleistocene cores. Aeolian transport and partial or complete removal of iron cutans from around quartz grains results in dunes of widely varying colour yet of similar age.     

Hydrography and circulation in the northwestern Bay of Bengal during the retreat of southwest monsoon

The distribution of temperature and salinity in the upper 500 m of the northwestern Bay of Bengal, adjoining the east coast of India, during the retreat of southwest monsoon (September) of 1983 is presented. This study reveals coastal upwelling (limited to the upper 40 m) induced by the local winds. Waters of higher surface salinity near the coast characterize the upwelling. The freshwater influx near the head of the Bay diluted the surface salinity to as low as 26.0 × 10⁻³. The surface circulation was weak and led to a net transport of 2.0 × 10⁶m³.s⁻¹ directed towards northeast.     

Estimation of the transfer function by autoregressive deconvolution techniques—An application to time series analysis

The transfer function of time-dependent models is classically inferred by the ordinary least squares (OLS) techniques. This OLS technique assumes independence of the residuals with time. However, in practical cases, this hypothesis is often not justified producing inefficient estimation of the transfer function. When the residuals constitute an autoregressive process, we propose to apply the Box-Jenkins' method to model the residuals, and to modify in a simple manner the primary convolution equation. Then, a multivariate regression technique is used to infer the transfer function of the new equation producing time-independent residuals. This three-step autoregressive deconvolution technique is particularly efficient for time series analysis. The reconstitution and the forecasting of real data are improved efficiently. Theoretically, the proposed method can be extended to the convolution equations for which the residuals follow a moving average or an autoregressive-moving average process, but the mathematical formulation is no longer direct and explicit. For this general case, we propose to approximate the moving average or the autoregressive-moving average process by an autoregressive process of sufficient order, and then the transfer function. Two case studies in hydrogeology will be used to illustrate the procedure.     

Mössbauer studies of iron silicate spinel at high pressure

We have measured in situ Mössbauer transmission spectra of iron silicate spinel (γ-Fe₂SiO₄) in a diamond anvil cell at room temperature and pressures up to 16 GPa. The observed spectra show a doublet characteristic of the paramagnetic state. The isomer shift and quadrupole splitting at atmospheric pressure are 1.10 and 2.63 mm/s, respectively, which are smaller than those of fayalite (α-Fe₂SiO₄). Both the isomer shift and quadrupole splitting decrease linearly with pressure with slope of ?0.003(1) and ?0.020(1) mm/sec · GPa, respectively. This simple linear trend suggests that no electronic or polymorphic transitions occur under 16 GPa except for those due to the small and continuous changes of volume and local symmetry under pressure. On the basis of a crystalline field calculation, the negative pressure derivative of the quadrupole splitting is associated with a trend towards an ideal cubic symmetry of the oxygen sublattice.     

A method for quantifying living roots of Spartina (cordgrass) and Juncus (needlerush)

Quantifying living roots in a marsh is a necessary but difficult task in wetland research. The two main difficulties usually encountered are distinguishing living from dead roots and processing a dense mat of fine roots. We found that living roots of salt-marsh plants release much more dissolved organic carbon (DOC) in boiling water than dead roots. Based on the finding, we developed a DOC procedure to quantify living roots of Spartina alterniflora and Juncus roemerianus. The DOC released in boiling water is a function of root activity, and the amount released can be used to calculate the living root biomass of a sample. The amount of living roots determined by the DOC method correlated well with the amount of living roots determined by the manual, sorting method (r² = 0.78, p<0.01). The DOC method is more objective, precise, and much less tedious than the manual sorting method.     

The value ofC efor the Arabian Sea during summer monsoon

We estimate, from the moisture budget the bulk aerodynamic coefficient for latent heat flux (C _e) during the monsoon season over the central Arabian Sea. The average value ofC _eunder active monsoon conditions was found to be 2.25 × 10⁻³ which is nearly 60% higher than those previously used.     

Reassessment of the garnet-clinopyroxene Fe−Mg exchange thermometer: I. Evaluation of the Pattison and Newton (1989) experiments

This is Part 1 of a two-part study on the garnet(Grt)-clinopyroxene(Cpx) Fe−Mg exchange equilibrium widely used in geothermometry of amphibolites, granulites and eclogites. The experimental data set previously published by Pattison and Newton (PN) on this equilibrium comprises potentially the largest source of data used for extraction of thermodynamic properties for Grt and Cpx and for formulation of a new thermometer expression (Part 2 of this study). However, the experiments suggested unusual Grt and Cpx solution properties and resulted in a geothermometer expression that generally returns low temperatures for granulites. To verify the reliability of the experiments, and to test the experimental technique of Grt buffering used in the PN study, 69 of the 125 PN runs were reanalysed using wavelength dispersive spectroscopic analysis guided by backscattered electron images. Special attention was given to analysis of compositional changes in Grt in addition to Cpx. With minor exceptions, the reanalysis supports the analytical results of the original PN study. The mechanism of equilibration in the runs was a multicomponent solution-precipitation process involving dissolution of metastable Grt and Cpx starting materials and precipitation of newly formed Cpx and Grt, the latter assumed to represent equilibrim compositions. The experiments cannot be regarded either as reversed or as strictly bracketing. Nevertheless, Cpx compositional parameters, including Mg#, Na, Ca and Al, show the same smooth variations with temperature, pressure and composition as found by PN, suggestive of a close approach to equilibrium. Grt in run products retained the moderate heterogeneity of the starting Grt, although there are subtle compositional trends indicating small changes in Grt Mg# consistent with mass balance constraints. These results uphold the essential validity of the technique of Grt buffering used in the experiments. Nevertheless, it is possible that the final Grt compositional may not coincide with the nominal starting Grt compositions, requiring assignment of compositional uncertainties to Grt for purposes of thermodynamic modelling of the data. Cpx analyses also require assignment of non-zero compositional uncertainties.