Anelasticity of the crust and upper mantle beneath north and central India from the inversion of observed Love and Rayleigh wave attenuation data

The fundamental mode Love and Rayleigh waves generated by earthquakes occurring in Kashmir, Nepal Himalaya, northeast India and Burma and recorded at Hyderabad, New Delhi and Kodaikanal seismic stations are analysed. Love and Rayleigh wave attenuation coefficients are obtained at time periods of 15–100 seconds, using the spectral amplitude of these waves for 23 different paths along northern (across Burma to New Delhi) and central (across Kashmir, Nepal Himalaya and northeast India to Hyderabad and Kodaikanal) India. Love wave attenuation coefficients are found to vary from 0.0003 to 0.0022 km^–1 for northern India and 0.00003 km^–1 to 0.00016 km^–1 for central India. Similarly, Rayleigh wave attenuation coefficients vary from 0.0002 km^–1 to 0.0016 km^–1 for northern India and 0.00001 km^–1 to 0.0009 km^–1 for central India. Backus and Gilbert inversion theory is applied to these surface wave attenuation data to obtainQ ^–1 models for the crust and uppermost mantle beneath northern and central India. Inversion of Love and Rayleigh wave attenuation data shows a highly attenuating zone centred at a depth of 20–80 km with lowQ for northern India. Similarly, inversion of Love and Rayleigh wave attenuation data shows a high attenuation zone below a depth of 100 km. The inferred lowQ value at mid-crustal depth (high attenuating zone) in the model for northern India can be by underthrusting of the Indian plate beneath the Eurasian plate which has caused a low velocity zone at this shallow depth. The gradual increase ofQ ^–1 from shallow to deeper depth shows that the lithosphere-asthenosphere boundary is not sharply defined beneath central India, but rather it represents a gradual transformation, which starts beneath the uppermost mantle. The lithospheric thickness is 100 km beneath central India and below that the asthenosphere shows higher attenuation, a factor of about two greater than that in the lithosphere. The very lowQ can be explained by changes in the chemical constitution taking place in the uppermost mantle.     

Charge distribution on the atmospheric aerosol over the North Atlantic and the Indian Ocean

The distribution of electric charge on the marine aerosol was determined near the sea surface of the Indian Ocean and the North Atlantic during the final stage of the Snellius II-expedition. Mean values for small ion concentrationsn ⁺=455 cm^–3 andn ^–=340 cm^–3 were found over the Atlantic, whilen ⁺=310 andn ^–=250 cm^–3 were the mean values over the Indian Ocean. The ration ⁺/n^– increased from 1.2 to 2.0 with decreasing wind velocity. At wind velocities below 5 m/s 75% of the net space charge near the ocean surface was found to be carried by small ions.     

CodaQ as a combination of scattering and intrinsic attenuation: Numerical simulations with the boundary integral method

Numerical modelling ofSH wave seismograms in media whose material properties are prescribed by a random distribution of many perfectly elastic cavities and by intrinsic absorption of seismic energy (anelasticity) demonstrates that the main characteristics of the coda waves, namely amplitude decay and duration, are well described by singly scattered waves in anelastic media rather than by multiply scattered waves in either elastic or anelastic media. We use the Boundary Integral scheme developed byBenites et al. (1992) to compute the complete wave field and measure the values of the direct waveQ and coda wavesQ in a wide range of frequencies, determining the spatial decay of the direct wave log-amplitude relation and the temporal decay of the coda envelope, respectively. The effects of both intrinsic absorption and pure scattering on the overall attenuation can be quantified separately by computing theQ values for corresponding models with (anelastic) and without (elastic) absorption. For the models considered in this study, the values of codaQ ^–1 in anelastic media are in good agreement with the sum of the corresponding scatteringQ ^–1 and intrinsicQ ^–1 values, as established by the single-scattering model ofAki andChouet (1975). Also, for the same random model with intrinsic absorption it appears that the singly scattered waves propagate without significant loss of energy as compared with the multiply scattered waves, which are strongly affected by absorption, suggesting its dominant role in the attenuation of coda waves.     

Zero-frequency tidal distrotion due to the orbital elements of tide-forming bodies

Summary It has been proved that orbital elements of perturbing bodies should be taken into account when 10^–10 accuracy is required for zero-frequency tidal distortion in the second zonal Stokes parameter of the geopotential. A solution at the 10^–15 level of magnitude has been presented. The zero-frequency tidal distortion in the fourth zonal Stokes parameter has been derived as 1·3×10^–10 if the secular Love number is unity. It should be reflected in the geopotential models respecting the 10^–10 level of magnitude.     

A determination of the wind-driven ocean circulation from the vorticity budget of the atmosphere

Summary The annual mean distribution of the surface stress curl over the Northern Hemisphere has been estimated from the horizontal vorticity advection in the atmosphere by using the upper-wind statistics as published byCrutcher [2]³). The results are used to compute the wind-driven mass transport (Sverdrup transport) in North Atlantic and North Pacific. The calculated intensity of the Gulf Stream is largest at the latitude 35°N, where a mass transport of 45×10⁶ tons sec^–1 is obtained; for the maximum intensity of the Kuroshio current a value of 60×10⁶ tons sec^–1 is obtained.Research supported in part by the Section of Atmospheric Sciences, National Science Foundation, Grant GP-2561.The research for this study was started by the author at the Department of Meteorology and Oceanography, The University of Michigan, Ann Arbor, Mich.     

Fine scale distributions of porosity and particulate excessPb, organic carbon and CaCO3 in surface sediments of the deep equatorial Pacific

Sediment samples were recovered from the central equatorial Pacific Ocean, sectioned at 1-mm intervals, and analyzed for porosity, organic carbon, excess²¹⁰Pb and CaCO₃. Steep porosity gradients were measured in the upper 1 cm of the sediment column with extremely high values observed near the sediment surface. Similarly, particulate organic carbon contents are highest at the sediment surface, decrease sharply in the upper 1 cm, and are relatively constant between 1 and 5 cm. CaCO₃ values, on the other hand, are lowest at the sediment surface and increase to a constant value below 5–10 mm depth. At the carbonate ooze sites, excess²¹⁰Pb is present throughout the upper 5 cm of the sediments suggesting relatively rapid particle mixing rates. However, extremely high excess²¹⁰Pb activities (> 100 dpm/g) are observed at the sediment surface with sharp gradients present in the upper 1 cm which would suggest slow rates of mixing. This apparent contradiction along with the major features of the CaCO₃ and particulate organic carbon profiles can be explained by a particle-selective feeding mechanism in which organic carbon, excess²¹⁰Pb-enriched particles are preferentially maintained at the sediment surface via ingestion and defecation by benthic organisms.     

2-D Crustal Velocity Structure Along Hirapur-Mandla Profile in Central India: An Update

The 2-D crustal velocity model along the Hirapur-Mandla DSS profile across the Narmada-Son lineament in central India (Murty et al., 1998) has been updated based on the analysis of some short and discontinuous seismic wide-angle reflection phases. Three layers, with seismic velocities of 6.5–6.7, 6.35–6.40 and 6.8 km s^–1, and upper boundaries located approximately at 8, 17 and 22 km depth respectively, have been identified between the basement (velocity 5.9 km s^–1) and the uppermost mantle (velocity 7.8 km s^–1). The layer with 6.5–6.7 km s^–1 velocity is thin and is less than 2-km deep between the Narmada north (at Katangi) and south (at Jabalpur) faults. The upper crust shows a horst feature between these faults, which indicates that the Narmada zone acts as a ridge between two pockets of mafic intrusion in the upper crust. The Moho boundary, at 40–44 km depth and the intra-crustal layers exhibit an upwarp suggesting that the Narmada faults have deep origins, involving deep-seated tectonics. A smaller intrusive thickness between the Narmada faults, as compared to those beyond these faults, suggests that the intrusive activities on the two sides are independent. This further suggests that the two Narmada faults may have been active at different geological times. The seismic model is constrained by 2-D gravity modeling. The gravity highs on either side of the Narmada zone are due to the effect of the high velocity/high density mafic intrusion at upper crustal level.     

Enhanced ion acoustic fluctuations and ion outflows

A number of observations showing enhanced ion acoustic echoes observed by means of incoherent scatter radars have been reported in the literature. The Received power is extremely enhanced by up to 1 or 2 orders of magnitude above usual values, and it is mostly contained in one of the two ion acoustic lines. This spectral asymmetry and the intensity of the received signal cannot be resolved by the standard analysis procedure and often causes its failure. As a result, and in spite of a very clear spectral signature, the analysis is unable to fit the plasma parameters inside the regions of ion acoustic turbulence. We present European Incoherent Scatter radar (EISCAT) observations of large ion outflows associated with the simultaneous occurrence of enhanced ion acoustic echoes. The ion fluxes can reach 10¹⁴ m^–2 s^–1 at 800 km altitude. From the very clear spectral signatures of these echoes, a method is presented to extract estimates of the electron temperature and the ion drift within the turbulent regions. It is shown that the electron gas is strongly heated up to 11 000 K. Also electron temperature gradients of about 0.02 K/m exist. Finally, the estimates of the electron temperature and of the ion drift are used to study the possible implications for the plasma transport inside turbulent regions. It is shown that strong electron temperature gradients cause enhancement of the ambipolar electric field and can account for the observed ion outflows.     

Coupling of early diagenetic processes and sedimentary dynamics in tropical shelf environments: the Gulf of Papua deltaic complex

Tropical mobile mud belts represent a major class of biogeochemical and diagenetic systems characterized by extensive and frequent physical reworking of fine-grained, organic-rich deposits underlying oxygenated waters. Large regions of the Gulf of Papua, Papua New Guinea deltaic complex are dominated by such conditions. A reworked mud belt lies within the inner shelf between 10 and 20 m depth on a sedimentary clinoform derived from coalescing deltas. Deposits across the topset are typically suboxic, nonsulfidic over the upper 0.5–1 m, and have low to moderate maximum pore water concentrations of dissolved Fe(II) and Mn(II) (100–200, but up to 800 μM). Sediments are reactive, with surficial ΣCO₂ production 0.1–0.3 mM d⁻¹ and benthic O₂ fluxes 23±15 mmol m⁻² d⁻¹ (upper 20 cm). The highest rates occur within inner topset deposits (10–20 m) and near the high accumulation rollover region of the topset–foreset beds (40–50 m). Lower rates are found inshore along intertidal channels—mangrove fringe and within scoured or exposed consolidated deposits of the middle topset region. Remineralization rate patterns are independent of relative dominance by terrestrial or marine carbon in sediments. Dissolved O₂ usually penetrates 2–5 mm into surface sediments when macrofaunal burrows are absent. More than 75% of the highly reactive sedimentary Fe(III) pool (350–400 μmol g⁻¹) is typically diagenetically reduced in the upper 0.5 m. Pore water can be measureably depleted at depths >0.5 m, but dissolved H₂S generally remains below detection over the upper 1–2 m. As in other deltaic topset regions, concentration gradients often indicate that compared to many marine deposits of similar sediment accumulation rates, relatively refractory C_org is supplied to the SO₄ reducing zone. Sedimentary C/S ratios are 4–6 within the suboxic topset regions but decrease to <3 in offshore foreset beds where sulfidic diagenesis dominates. Only 15–20% of the diagenetically reduced Fe(II) is pyritic and a maximum of 10–25% is carbonate, implying that most Fe(II) is associated with authigenic or lithogenic silicates or oxides. The dominance of suboxic, nonsulfidic diagenetic processes reflect coupling between delivery of oxide-rich terrestrial debris, remobilization and reoxidation of deposits, and repetitive entrainment/remineralization of both labile and refractory organics. Distinct sedimentary indicators of reactive, suboxic mobile mud belts within tropical climatic zones are: abundant total highly reactive Fe (ΣFeR )>300 μmol g⁻¹; most reactive Fe is diagenetically reduced (ΣFe(II)/ΣFeR0.7–0.8); the proportion of diagenetically reduced Fe present as pyrite is low (Py–Fe(II)<0.2); C/S 4–8; and C_org/particle surface area <0.4 (mg C m⁻²). These depositional environments must be most common in tropical climates during high sea stand.     

Diffraction ofSH-waves originating from a moving point source by a rigid quarter space

Summary ApplyingWiener-Hopf techniques the diffraction ofSH-waves originating from a moving point source due to the presence of a step change in the elevation of the surface of a homogeneous half-space have been obtained. An interesting feature of resonance associated with the motion of a load at a velocityv ₀ ^–1 approachingv ₀ ^–1 has been observed.     

An investigation of the size frequency resolution of an aerosol by the static method of diffusion

Summary Theoretical considerations are put forward for identifying the Dynamic and Static Methods of Diffusion measurements of aerosol particles. On the basic of this equality the Exhaustion Method of Diffusion developed byPollak andMetnieks for the size-frequency resolution of a heterogeneous aerosol by Dynamic Diffusion Methods should be applicable to Static Methods. This hypothesis is tested using a Pollak-Nolan nucleus counter as decay vessel. It is shown that turbulence in the counter is the most serious obstacle against this type of approach to finding the size components of an aerosol. The time involved, the high humidity in the counter and a hypothetical nucleus-free zone are other major difficulties for relatively large Aitken nuclei. It is established, however, that for particles whose diffusion coefficient is of the order of 10^–3 cm² sec^–1, the proposed method should be quite workable.     

A fluid-filled fracture as possible mechanism of ground deformation at Phlegraean Fields,Italy

Two uplift episodes have recently been recorded at Phlegraean Fields, a volcanic region near Naples (south-central Italy). The first episode occurred in 1970 and lasted up to 1972; the second has begun at the end of 1982 and is still in progress.An attempt to model the ground deformations which occurred during the 1970–1972 event is made in this paper. The source has been assumed to be a two-dimensional fluid-filled fracture, similar to a sill. A good fit with experimental data has been obtained for a short (1–2 km long) shallow (3 km deep) source and a driving pressure ranging from 60 to 210 bars. Rigidity values have been fixed at 3.5–4.0 × 10¹⁰ c.g.s., with Poisson ratio equal to 0.3.This solution which differs from previous approaches byMogi (1958) andWalsh andDecker (1971) is non-unique, but the present results are in good agreement with observed horizontal and vertical displacements. Notwithstanding the oversimplification made in assuming a homogeneous elastic semi-infinite medium, the predicted stress field seems to be in agreement with the fault-plane solutions and the pattern of epicenters determined for the uplift-seismic swarm episode that is still in progress.     

A new measure of seismicity

Summary In this paper an attempt has been made to analyse the World shallow focus earth-quakes by the method ofGutenberg andRichter. Frequency-magnitude analysis of various earthquakes indicate that in the relation logN=a–b M, the ratio (b/a)^–1 satisfies fairly well the seismicity of a region and it is felt that this ratio may be used as a measure of seismicity for a given region.     

Global geodynamic long-term variations and expanding earth hypothesis

Summary The expanding Earth hypothesis has been discussed on the basis of recent LLR data and the long-term decrease observed in the angular velocity of the Earth's rotation. It has been proved that the observed data do not support the hypothesis. Realistic limits for hypothetical long-term variations in the Earth's dimension consistent with observations have been derived as ±3 × 10^–3 m cy^–1.     

Global vorticity budget over the tropics and subtropics at 200-mb during northern hemisphere summer

The large-scale terms in the vorticity equation are evaluated usingKrishnamurti's (1971a, b) summer mean winds at 200 mb for a global belt from 25°S to 45°N. The production of vorticity by the divergent wind field is found to be imbalanced over all of the tropical and subtropical belt. As a result there is a requirement for a sub-grid scale (space or time) mechanism which removes negative vorticity from the regions of strong divergence (Tibetan and Mexican highlands) and removes positive vorticity from the regions of strong convergence (mid-oceanic troughs) at 200 mb during northern summer at a rate of approximately 4×10^–10 sec^–2. As suggested byHolton andColton (1972), in regions of strong and persistent convection, such as the Tibetan Plateau, deep cumulus clouds can account for this transport. However, the mechanism for removing positive vorticity in the vicinity of the upper tropospheric mid-oceanic troughs is still an intriguing and open question.On leave-of-absence at the National Science Foundation, Climate Dynamics Research Section.     

The evolution of andesite volcano structures: new evidence from gravity studies in Costa Rica

Published gravity data on active volcanoes generally reflecteither the low density scoriaceous/pumiceous deposits that are localized within ring-fracture collapse depressions, such as the calderas of mature silicic volcanoes,or the high density frozen magma conduits that occur beneath basaltic shields and cones. The intensive gravity surveys reported here over three complex andesite volcanoes reveal features of both types. Their multi-component gravity fields have crater-centred positive anomalies (1–2 km diameter) surrounded by broader zones of negative gravity with similar amplitudes but greater width (5–10 km). The former are thought to reflect sub-crater magma pipes ofnormal density (ca. 2.5–2.6 Mg m⁻³) surrounded by pyroclastic scoria, ashes and occasional lava flows of muchlower net density (1.8–2.4 Mg m⁻³) which, in turn, account for the negative anomalous zones because the deeper, more consolidated and older parts of these andesite volcano edifices have more normal densities (2.3–2.6 Mg m⁻³).The low density materials are particularly interesting because they appear to have filled topographic depressions to depths of several hundred metres, especially where old caldera-like structures have been postulated from the steep gravity gradients over perimeter ring faults. A model is developed whereby short periods of caldera collapse, associated with intermittent, large high level magma bodies, are interspersed by normal crater-like activity with narrow sub-surface magma pipes. Dominantly pyroclastic activity from summit craters generates the materials that gradually fill earlier-formed topographic depressions. This study demonstrates the unique value of detailed gravity surveys, combined with surface geological information, for modelling and understanding the evolution of active volcano summit regions.     

Heavy metals in surface sediment and epibenthic macroinvertebrates from the coastal zone and continental slope of Kenya

In surface sediments, a statistically significant increase was measured in the concentration of copper and cadmium along all transects, and of zinc along only the most southern (Gazi) transect, radiating into the Indian Ocean perpendicular to the Kenyan coast. Mean copper and cadmium increased from 5 to 30 μg g⁻¹ dry wt and from 0.01 to 0.34 μg g⁻¹ dry wt in shallow coastal (±20 m depth) to deep-sea stations (±2000 m depth), respectively. These gradients were found both during the south-east monsoon and north-east monsoon period. The shallow estuarine zone of the Sabaki river mouth showed significantly enhanced levels of total organic carbon and nitrogen and all metals analysed, except cadmium.In crustaceans, the concentration of copper and particularly cadmium was significantly above baseline levels, varying from 45 to 90 μg g⁻¹ dry wt and 1.0 to 8.5 μg g⁻¹ dry wt, respectively. Zinc levels (49–102 μg g⁻¹ dry wt) were at about baseline levels or a little elevated. On the contrary, lead showed very low concentrations, varying from 0.1 to 0.6 μg g⁻¹ dry wt. Other species generally showed the same pattern.     

Travel-time residuals and three-dimensional velocity structure of Italy

Jeffreys-Bullen P and PKP travel-time residuals observed at more than 50 seismic stations distributed along Italy and surrounding areas in the time interval 1962–1979, indicate the complex velocity pattern of this region. Strong lateral velocity inhomogeneities and low velocity zones are required to explain the observed pattern of residuals. In particular, late arrivals of about 1 sec are observed in the Apenninic mountain range, requiring both greater crustal thickness and low velocity layers, coherent with seismic refraction data and surface wave dispersion measurements. The seismic stations located in the Western and Eastern Alps indicate the presence of high velocities. In the Western Alps the strong azimuthal variation of residuals and the high values of early arrivals have a close relationship to the Ivrea body, an intrusive crustal complex characterized by a velocity as high as 7–7.2 km/sec.A travel-time inversion performed with theAki et al. (1977) block model, confirms the peculiar characteristics and the sharp variations in the lithosphere of the whole Italian region, with values of velocity perturbations between many adjacent blocks, ranging in size from 50 to 100 km, and independent from the earth parametrization chosen, reaching values up to 10% in the lithospheric part and 5% in the asthenosphere. 3-D inversion requires also high velocity along the Tyrrhenian coastal margin, equivalent to an uprise of major crustal and lithospheric discontinuities along this part of the Italian peninsula. Moreover low velocity material must be present in the northern part of the Adriatic foreland, in the lithosphere-asthenosphere system, closely related to the stress and seismicity pattern, and the lateral bending of the lithosphere in the same region.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.