On the heat budget of the Arabian Sea

Summary The rate of oceanic heat storage of the upper 200m of the Arabian Sea is explained in terms of net air-sea heat flux (Q _F), heat change due to horizontal divergence and vertical motion (Q _V) and heat change due to lateral advection (Q _A). The analysis revealed that the heat storage of the Arabian Sea is mainly controlled byQ _V while the effect ofQ _A is much larger than expected. Parameterisation of summer cooling revealed that the depletion of energy from the mixed layer is mainly due to upwelling and horizontal advection though large amount of heat is accumulated due to net air-sea heat flux. The annual heat balance of the upper 200m of the Arabian Sea suggested large heat gain by air-sea exchange processes. About two third of this heat gain is compensated by horizontal advection and one third by vertical advection.With 4 Figures     

Multiple episodes of tectonothermal processes in the Eastern Ghats granulite belt

A suite of rocks from Borra Carbonate Granulite Complex (BCGC) in the Eastern Ghats granulite belt displays superposed structures and overprinted mineral assemblages that reveal multiple episodes of tectonothermal reworking of the complex under granulite facies condition. Five distinct episodes of deformation (D₁, D₂, D₃, D₄ and D₅) and four phases of metamorphism (M₁, M₂, M₃ and M₄) are recorded. The signature of the earliest tectonothermal event, D₁ is a gneissic foliation (S₁) denned by segregation of peak granulite facies mineral assemblages corresponding to prograde M₁ metamorphism. M₂ metamorphic overprint represents an episode of near-isobaric cooling of the complex under a static condition. D₂ represents an episode of ductile deformation manifested by isoclinal folding (F₂) and associated extensional structures, within a broad framework of coaxial bulk deformation. The present study reveals that D₂ took place subsequent to M₂ - Subsequent deformation, D₃, produced F₃ folds and also deformations of boudins formed during D₂. M₃, which is synchronous with F₃, represents a near isothermal decompression of the BCGC. This was followed by a weak structural readjustment (D₄), producing E-W cross folds. The latter was not, however, associated with any recognizable petrological reworking. In the terminal events, deformation (D₅) and mineral reactions (M₄) were localized along narrow intersecting shear zones. The latter acted as channelways for carbonic and still later hydrous fluid infiltration. The available thermobarometric data from BCGC and other areas of the Eastern Ghats belt reveal that reworking during M₂ and M₃ ensued in a thermally perturbed regime. The high thermal regime might also have persisted during carbonic fluid infiltration related to terminal reworking (M₄).     

Electron cyclotron waves in the presence of parallel electric fields in the Earth’s auroral plasma

The electron cyclotron waves that originate at low altitudes (<0.5 R_E) and observed by ground facilities have been studied in the presence of a weak parallel electric field in auroral magnetoplasma consisting of trapped energetic auroral electrons and cold background electrons of ionospheric origin. The model distribution for auroral trapped electrons is taken as Maxwellian ring distribution. An expression for the growth rate has been obtained in the presence of parallel electric field assuming that the real frequency in the whistler mode is not affected by the presence of the electric field. The results show that waves grow (or damp) in amplitude for a parallel (or antiparallel) electric field. The influence of the electric field is more pronounced at a shorter wavelength spectrum. An increase in population of energetic electrons increases the growth rate and thus, plays a significant role in the wave excitation process in the auroral regions.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Seismic horizontal pullout capacity of shallow vertical anchors in sand

The concept of the equivalent free surface has been extended to determine the seismic horizontal pullout capacity of shallow vertical strip plate anchors buried in sand. The analysis has been done rigorously by using the method of stress characteristics. The results have been expressed in the form of non-dimensional charts. The pullout resistance has been found to reduce quite extensively with increase in the magnitude of horizontal earthquake acceleration. The results were compared with the previously published data, and it was seen that the computed pullout resistance with the proposed method was found to be lowest. This revised version was published online in July 2006 with corrections to the Cover Date.