Seismic analysis of lock–soil–fluid systems by hybrid BEM–FEM

A study of soil–structure–fluid interaction (SSFI) of a lock system subjected to harmonic seismic excitation is presented. The water contained lock is embedded in layered soils supported by a half-space bedrock. The ground excitation is prescribed at the soil–bedrock interface. The response is numerically obtained through a hybrid boundary element (BEM) finite element method (FEM) formulation. The semi-infinite soil and the fluid are modeled by the BEM and the lock is modeled by the FEM. The equilibrium equation for the lock system is obtained by enforcing compatibility and equilibrium conditions at the fluid–structure, soil–structure and soil–layer interfaces under conditions of plane strain. To the authors’ knowledge this is the first study of a lock system that considers the effects of dynamic soil–fluid–structure interaction through a BEM–FEM methodology. A numerical example and parametric studies are presented to examine the effects of the presence of water, lock stiffness, and lock embedment on the response.     

Nutrient-chlorophyll trajectories across trophic gradients

We estimate the response of chl-a (mg · m^–3) to changes in concentrations of total phosphorus (TP) by calculating the slopeS = chl-a/TP in chl-a =f(TP) graphs. Results show that in years where algae are P-limited oligotrophic lakes respond less (median slope 0.21) to changes in nutrient concentrations than eutrophic lakes, (median slope 0.31) and these again less than hypereutrophic lakes, (median slope 1.02). We find no saturation value for the slope within the TP range considered (6–480 mg · m^–3). Chl-a in eutrophic lakes responds more frequently to non-nutrient factors than oligotrophic and hypereutrophic lakes. Results obtained by replacing TP with a new nutrient parameter, TP = 0.056 · TP · IN^0.226, in which inorganic nitrogen, IN, is factored in, suggest that nitrogen has an influence on chl-a in oligotrophic lakes. Blue-green algae respond less to changes in TP than other algal species, e.g., diatoms.     

Two-dimensional scattering of SH waves in a soil layer underlain with a sloping bedrock

A two-dimensional analysis is applied to examine the effect that a sloping bedrock half-space has on the amplification of an anti-plane shear wave. The direct boundary integral equation method is used for the two-dimensional analysis. The particular soil–rock configuration investigated includes a homogeneous soil layer underlain by a sloping rock half-space. The rock half-space dips for a horizontal distance L and then becomes horizontal so that the overlying soil layer has a thickness H that remains constant from this point to infinity. The materials in the soil–rock configuration are considered viscoelastic except in the rock half-space below soil layer thickness H, which is considered elastic. This limitation in damping is due to the correction used for the truncation of the half-space boundary. Four cases are used to study the relationship between rock slope and surface displacement, vertical, 1:2, 1:4, 1:8. Surface displacements are determined for each of these cases for half-space incidence angles of 90, 75, and 60°. To allow for applicability to a wide range of problems, results are determined as a function of dimensionless parameters. In addition, solutions from a one-dimensional analysis are compared with the results of the two-dimensional analysis to establish limits outside of which a one-dimensional analysis suffices.     

Geomorphic imprint of landslides on alpine river systems,southwest New Zealand

An inventory of 846 mass movements, mainly landslides, in two alpine regions of southwest New Zealand was created to explore the geomorphic impacts of slope‐failure processes on river channels and valley floors. In total, 213 (i.e. 27 per cent) of the slope failures descended to valley floors, affecting the geomorphology of trunk channels (catchment area A_C > 10 km²) and valley floors in recurring patterns. A nominal classification system is introduced for characterizing (a) the physical contact nature between landslides and river channels, and (b) the resulting geomorphic consequences for drainage. Although landslide area A is useful for estimating the length of channel directly impacted by debris, it does not necessarily predict the direction of fluvial response or type of impact. Dominant persistent geomorphic imprints of bedrock landslides include channel occlusions and landslide dams in South Westland and Fiordland, respectively. Differences in size distribution and geomorphic effects on river systems between the two study regions are attributed to bedrock geology, tectonics and sediment flux. Although South Westland rivers are more frequently affected by landslides, disrupting long‐term effects such as blockage are more persistent in Fiordland. Copyright © 2005 John Wiley & Sons, Ltd.     

International standard (ISO) on seismic actions for designing geotechnical works–An overview

The seismic performance of geotechnical works is significantly affected by ground displacement. In particular, soil–structure interaction and effects of liquefaction play major roles and pose difficult problems for engineers. An International Standard, ISO23469, is being developed for addressing these issues in a systematic manner within a consistent framework. The objective of this paper is to give an overview of this International Standard.In this International Standard, the seismic actions are determined through two stages. The first stage determines basic seismic action variables, including the earthquake ground motion at the site, the potential for earthquake-associated phenomena such as liquefaction and induced lateral ground displacement. These basic variables are used, in the second stage, for specifying the seismic actions for designing geotechnical works. In the second stage, the soil–structure interaction plays a major role. Types of analyses are classified based on a combination of static/dynamic analyses and the procedure for soil–structure interaction classified as follows:

– simplified: soil–structure interaction of a global system is modeled as an action on a substructure;

– detailed: soil–structure interaction of a global system is modeled as a coupled system.

Earthquake damage detection in the Imperial County Services Building III: Analysis of wave travel times via impulse response functions

The majority of structural health monitoring methods are based on detecting changes in the modal properties, which are global characteristics of the structure, and are not sensitive to local damage. Wave travel times between selected sections of a structure, on the other hand, are local characteristics, and are potentially more sensitive to local damage. In this paper, a structural health monitoring method based on changes in wave travel times is explored using strong motion data from the Imperial Valley Earthquake of 1979 recorded in the former Imperial County Services (ICS) Building, severely damaged by this earthquake. Wave travel times are measured from impulse response functions computed from the recorded horizontal seismic response in three time windows—before, during, and after the largest amplitude response, as determined from previous studies of this building, based on analysis of novelties in the recorded response. The results suggest initial spatial distribution of stiffness consistent with the design characteristics, and reduction of stiffness following the major damage consistent with the spatial distribution of the observed damage. The travel times were also used to estimate the fundamental fixed-base frequency of the structure f₁ (assuming the building deformed as a shear beam), and its changes during this earthquake. These estimates are consistent with previous estimates of the soil–structure system frequency, f_sys, during the earthquakes (f₁<f_sys as expected from soil–structure interaction studies), and with other estimates of frequency (f₁ from ETABS models, and f_sys from ambient vibration tests, and “instantaneous” f₁ from high-frequency pulse propagation).     

Experimental investigation of the elastic modulus of a fractal system—A model of fractured rocks

The elastic properties of a physical model representing a damaged rock matrix were studied using a square lattice deformed under tensile stress. The elastic modulusM of such a system varies in agreement with percolation theory as|x–x _c | ^f , wherex is the damage parameter andx _c the threshold value of the damage parameter,f3.6. Atxx _c the scale dependence ofM can be expressed asML ^–f/v , whereL is the size of the sample andv the correlation exponent in percolation theory.The experimental results are of interest in assessing elastic properties in earthquake focal zones and fault zones in general.     

Determination of characteristics of steam reservoirs by Radon-222 measurements in geothermal fluids

The authors conducted a Rn²²² survey in wells of the Larderello geothermal field (Italy) and observed considerable variations in concentrations. Simple models show that flow-rate plays an important part in the Rn²²² content of each well, as it directly affects the fluid transit time in the reservoirs. Rn²²² has been sampled from two wells of the Serrazzano area during flow-rate drawdown tests. The apparent volume of the steam reservoir of each of these two wells has been estimated from the Rn²²² concentration versus flow-rate curves.List of symbols Q Flow-rate (kg h^–1) - Decay constant of Rn²²² (=7.553×10^–3 h^–1) - Porosity of the reservoir (volume of fluid/volume of rock) - ₁ Density of the fluid in the reservoir (kg m^–3) - ₂ Density of the rock in the reservoir (kg m^–3) - M Stationary mass of fluid filling the reservoir (kg). - E Emanating power of the rock in the reservoir (nCi kg _rock ^–1 h^–1). - P Production rate of Rn²²² in the reservoir: number of atoms of Rn²²² (divided by 1.764×10⁷) transferred by the rock to the mass unit of fluid per unit time (nCi kg _fluid ^–1 h^–1). - N Specific concentration of Rn²²² in the fluid (nCi kg^–1) - Characteristic time of the steam reservoir at maximum flow-rate (=M/Q)     

Analysis of the relations between slope failure distribution and seismic ground motion during the 2008 Wenchuan earthquake

The 2008 Wenchuan earthquake with M_w7.9 occurred at Yingxiu County in Longmenshan thrust fault belt, southwest China, having triggered a huge amount of slope failures. This paper applied a detailed inventory with more than 190,000 slope failures and strong ground motion records of 187 seismic stations to analyze the qualitative and quantitative relations between slope failure distribution and seismic parameters. The results revealed that slope failure distribution was exponentially decreased with the increment of epicentral distance and distance from surface fault rupture; peak ground acceleration (PGA) on the hanging wall side was apparently larger than that on footwall side. Linear correlation between concentration of slope failures (LNC) and the percentage of the area affected by slope failure (LAP) and PGA was demonstrated by statistical analysis, which revealed that 0.18~0.21 g horizontal PGA was the threshold value of the occurrence of slope failure. Furthermore, this paper presented an empirical model for the attenuation relation of slope failure distribution.     

Ionisations-Neutralisationsbilanz in derE-Schicht nach Angaben von der totalen Sonnenfinsternis am 15. Februar 1961

Zusammenfassung Es werden die Grundgleichungen der Ionisations-Neutralisationsbilanz in derE-Schicht präzisiert und auf die Notwendigkeit der Berücksichtigung des Einflusses der lokalen ionisierenden Strahlungsquellen auf der Sonnenscheibe, der dynamischen ionosphärischen Vorgängen und der Veränderungen des äquivalenten Rekombinationskoeffizienten hingewiesen. Ferner wird eine Methodik zur Bestimmung der charakteristischen Grössen der Ionisations-Rekombinationsbilanz (q _0m ,q _d, ) dargelegt. Zu dem Zweck werden die Messergebnisse aus der Periode der Sonnenfinsternis von zwei unweit voneinander gelegenen Ionosphärenstationen verwendet. Die Methodik wird auf die während der totalen Sonnenfinsternis am 15. Februar 1961 auf den Ionosphärenstationen in Sofia und Nesebar erhaltenen Ergebnisse angewandt. Für den Rekombinationskoeffizienten werden Werte zwischen 0,63·10^–7 cm³sec^–1 und 2,32·10^–7 cm³sec^–1 und für die Elektronenproduktion unter dem Einfluss der ionisierenden Strahlung von der homogenen Sonnenscheibeq _0m 1700 cm^–3sec^–1 erhalten.

---

Summary The basic equations for the ionisation-neutralisation balance in theE layer have been shown the necessity to take in account the influence of the local ionisation sources on the disk of the sun, the dynamic ionospheric processes and the variations in the equivalent recombination coefficient. The method for determining the characteristic quantities of the ionisation-recombination balance (q _0m ,q _d, ) has been exhibited for this purpose are used measurement data from two ionospheric stations located not far from each other, in the period of the solar eclips on 15 February, 1961. The quantity of the recombination coefficient lay between 0.63×10^–7 sec^–1 cm³ and 2.32×10^–7 sec^–1 cm³; the electron production under the action of the ionising radiation of the homogen disk showsq _0m 1700 sec^–1 cm^–3.

     

Ambient vibration tests of a seven-story reinforced concrete building in Van Nuys, California, damaged by the 1994 Northridge earthquake

Results of two detailed ambient vibration surveys of a 7-story reinforced concrete building in Van Nuys, California, are presented. Both surveys were conducted after the building was severely damaged by the 17 January 1994, Northridge earthquake (M_L=5.3, epicenter 1.5 km west from the building site) and its early aftershocks. The first survey was conducted on 4 and 5 February 1994, and the second one on 19 and 20 April 1994, about one month after the 20 March aftershock (M_L=5.3, epicenter 1.2 km north–west from the building site). The apparent frequencies and two- and three-dimensional mode shapes for longitudinal, transverse and vertical vibrations were calculated. The attempts to detect the highly localized damage by simple spectral analyses of the ambient noise data were not successful. It is suggested that very high spatial resolution of recording points is required to identify localized column and beam damage, due to the complex building behavior, with many interacting structural components. The loss of the axial capacity of the damaged columns could be seen in the vertical response of the columns, but similar moderate or weak damage typically would not be noticed in ambient vibration surveys. Previous analysis of the recorded response of this building to 12 earthquakes suggests that, during large response of the foundation and piles, the soil is pushed sideways and gaps form between the foundation and the soil. These gaps appear to be closing during “dynamic compaction” when the building site is shaken by many small aftershocks. The apparent frequencies of the soil–foundation–structure system appear to be influenced significantly by variations in the effective soil–foundation stiffness. These variations can be monitored by a sequence of specialized ambient vibration tests.     

Geodynamic implications of deep mantle upwelling in the source of Tertiary volcanics from the Veneto region (South-Eastern Alps)

Major and trace element and Sr–Nd–Hf–Pb isotopic data for the most primitive Tertiary lavas from the Veneto region (South-Eastern Alps, Italy) show the typical features of HIMU hotspot volcanism, variably diluted by a depleted asthenospheric mantle component (⁸⁷Sr/⁸⁶Sr_i=0.70306–0.70378; _Ndi=+3.9 to +6.8; _Hfi=+6.4 to +8.1, ²⁰⁶Pb/²⁰⁴Pb_i=18.786–19.574). P-wave seismic tomography of the mantle below the Veneto region shows the presence of low-velocity anomalies at depth, which is consistent with possible upwellings of plume material. Between the depths of 100–250 km the velocity anomalies are approximately 2–2.5% slower than average, implying a temperature excess of about 220–280 K, in agreement with estimates for other mantle plumes in the world. In this context, the Veneto volcanics may represent the shallow expression of a mantle upflow. The presence of a HIMU-DM component in a collision environment has significant geodynamic implications. Slab detachment and ensuing rise of deep mantle material into the lithospheric gap is proposed to be a viable mechanism of hotspot magmatism in a subduction zone setting.     

Triggered earthquakes as stress gauge: Implication for the uppercrust behavior in the Grenoble area,France

A sequence of moderate shallow earthquakes (3.5M _L5.3) was located within the Vercors massif (France) in the period 1961–1984. This subalpine massif has been a low seismic area for at least 5 centuries. During the period 1962–1963, 12 shallow earthquakes occurred in the neighborhood (10 km) of the Monteynard reservoir, 30 km south of the city of Grenoble. The latest fourM _L4.0 earthquakes occurred in 1979–1984 either at larger distance (35 km) or greater depth (10 km) from the reservoir. Two triggering mechanisms are suggested for this sequence: (i) the direct effect of elastic loading through either increased shear stress or strength reducing by increased pore pressure at depth; (ii) the pore pressure diffusion induced by poroelastic stress change due to the reservoir filling.The weekly water levels, local balanced geological cross sections, and focal mechanisms argue for two types of mechanical connection between the earthquake sequence and the filling cycles of the Monteynard reservoir. The seismic sequence started with the 1962–1963 shallow earthquakes that occurred during the first filling of the reservoir and are typical of the direct effect of elastic loading. The 1979 deeper earthquake is located at a 10 km depth below the reservoir. This event occurred 16 years after the initial reservoir impoundment, but one month after the previous 1963 maximum water level was exceeded. Moreover the yearly reservoir level increased gradually in the period 1962–1979 and has decreased since 1980. Accordingly we suggest that the gradual diffusion of water from reservoir to hypocentral depths decreases the strength of the rock matrices through increased pore pressure. The transition between the two types of seismic response is supported by the analysis ofM _L3.5 earthquakes which all occurred in the period 1964–1971, ranging between 10 and 30 km distance from the reservoir. The three other delayed earthquakes of the 1961–1984 seismic sequence (M _L4 during the 1979–1984 period) are all located 35 km away from the reservoir. Based on the seismic activity, the estimates for the hydraulic diffusivities range between 0.2–10 m²/s, except for the first event that occurred 30 km north of the reservoir, the filling just started. The lack ofin situ measurements of crustal hydrological properties in the area, shared by most of the Reservoir-Induced-Seismicity cases, prevents us from obtaining absolute evidence for the triggering processes. These observations and conceptual models attest that previous recurrence times for moderate natural shocks (4.5M _L5.5) estimated within this area using historical data, could be modified by 0.1–1 MPa stress changes. These small changes in deviatoric stress suggest that the upper crust is in this area nearly everywhere at a state of stress near failure. Although the paucity of both number and size of earthquakes in the French subalpine massif shows that aseismic displacements prevail, our study demonstrates that triggered earthquakes are important tools for assessing local seismic risk through mapping fault zones and identifying their possible seismic behavior.     

Fractal approach of the temporal earthquake distribution in the Hellenic arc-trench system

The time clustering of earthquakes occurring in the Hellenic arc-trench system is quantitatively analyzed by means of the fractal dimension,D, of their time distribution in the time intervals of 1950–1985 (M _s >-4.5) and 1964–1985 (M _s 4.0). The results obtained imply that scale-invariant clustering holds over very large scale lengths of time,T, with 2²–2⁸T (in min) 2²⁰–2²², depending on the seismotectonic segment considered. In all segments a common feature is the relation between theD ₁,D ₂ andD ₃-values found for shallow, intermediate-depth and all-depth shocks, respectively:D ₃>D₁>D₂. TheD-values found for shallow shocks range between 0.137 and 0.191 with the exception of the Ionian Islands and Cretan segments where anomalously high values (D=0.221–0.251) have been determined. We discuss possible seismotectonic interpretations of the results.     

Unified approach to amplitude attenuation and coda excitation in the randomly inhomogeneous lithosphere

A unified model is proposed for explaining the frequency dependent amplitude attenuation and the coda wave excitation on the basis of the single scattering process in the randomly inhomogeneous lithosphere. Adopting Birch's law and a direct proportion between density and wave velocity, we statistically describe the inhomogeneous medium by one random function characterized by the von Karman autocorrelation function. We calculate the amplitude attenuation from the solid angle integral of scattered wave energy on the basis of the Born approxiimation after subtracting the travel-time fluctuation effect caused by slowly varying velocity inhomogeneities. This subtraction is equivalent to neglect energy loss by scattering within a cone around the forward direction. The random inhomogeneity of the von Karman autocorrelation function of order 0.35 with the mean square fractional fluctuation of 7.2×10^–3 1.3×10^–2 and the correlation distance of 2.15.1 km well explains observed backward scattering coefficientg and the ratioQ _P ^–1 /Q _S ^–1 , and observed and partially conjecturedQ _S ^–1 for frequencies between 0.5 Hz and 30 Hz.     

The Polar Cap (PC) indices: Relations to solar wind parameters and global magnetic activity

The polar geomagnetic activity resulting from solar wind–magnetosphere interactions can be characterized the Polar Cap (PC) indices, PCN and PCS. PC index values are derived from polar magnetic variations calibrated on a statistical basis such that the index approximate values in units of mV/m of the interplanetary “geo-effective” (or “merging”) electric field (E_M) conveyed by the solar wind. The timing and amplitude relations of the PC index to solar wind plasma and magnetic field parameters are reported. The solar wind effects are parameterized in terms of the geo-effective electric field (E_M) and the dynamical pressure (P_DYN). The PC index has a delayed and damped response to E_M variations and display saturation-like effects for E_M values exceeding 10 mV/m. Steady or slowly varying levels of solar wind dynamical pressure have little or no impact on the PC index above the effects related to E_M for which the solar wind velocity is also a factor. Sharp increases in the dynamical pressure generate impulsive variations in the PC index comprising a initial negative impulse of 5–10 min duration followed by a positive impulse lasting 10–20 min. Typical amplitudes of both the negative and the positive impulses are 0.2–0.5 units. A sharp decrease in the pressure produces the inverse sequence of pulses in the PC index. Auroral substorm activity represented by the AL index level has a marked influence on the average PC/E_M level at the transition from very quiet (AL0 nT) to disturbed conditions while more or less disturbed conditions (AL<−100 nT) have no systematic effect on the average PC/E_M values. At distinct substorm events the PC/E_M ratio has a minimum (0.8) in the pre-onset phase at around 20 min before substorm onset. The average ratio gradually increases in the expansion phase to reach a maximum value (1.1) at around 40 min after substorm onset (or 20 min after the largest (negative) peak in AL). At substorm recovery during the next 2 h the PC/E_M ratio decreases. Finally, we report on the application of polar magnetic variations to model the disturbance storm time (Dst) index development during magnetic storms by using the PC index as a source function to quantify the energy input to the ring current representing accumulated storm energy and characterized by the Dst index.     

Two-dimensional kinematic and rheological modeling of the 1912 pyroclastic flow,Katmai, Alaska

Pyroclastic flow emplacement is strongly influenced by eruption column height. A surface along which kinetic energy is zero theoretically connects the loci of eruption column collapse with all coeval ignimbrite termini. This surface is reconstructed as a two-dimensional energy line for the 1912 Katmai pyroclastic flow in the Valley of Ten Thousand Smokes from mapped flow termini and the runup of the ignimbrite onto obstructions and through passes. Extrapolation of the energy line to the vicinity of the source vent at Novarupta suggests the eruption column which generated the ignimbrite eruption was approximately 425 m high. The 1912 pyroclastic flow travelled about 25 km downvalley. Empirical velocity data calculated from runup elevations and surveyed centrifugal superelevations indicate initial velocities near Novarupta were greater than 79–88 m s^–1. The flow progressively decelerated and was travelling only 2–8 m s^–1 when it crossed a moraine 16 km downvalley. The constant slope of the energy line away from Novarupta suggests the flow was systematically slowed by internal and basal friction. Using a simple physical model to calculate flow velocities and a constant kinetic friction coefficient (Heim coefficient) of 0.04 derived from the reconstructed energy line, the flow is estimated to have decelerated at an average rate of –0.16 m s^–2 and to have taken approximately 9.5 minutes to travel 25 km down the Valley of Ten Thousand Smokes. The shear strength of the flowing ignimbrite at the moraine was approximately 0.5 kPa, and its Bingham viscosity when it crossed the moraine was 3.5 × 10³ P. If the flow was Newtonian, its viscosity was 4.2 × 10³ P. Reynolds and Froude numbers at the moraine were only 41–62 and 0.84–1.04, respectively, indicating laminar, subcritical flow.     

Erosional power in the Swiss Alps: characterization of slope failure in the Illgraben

Landslides and rockfalls are key geomorphic processes in mountain basins. Their quantification and characterization are critical for understanding the processes of slope failure and their contributions to erosion and landscape evolution. We used digital photogrammetry to produce a multi‐temporal record of erosion (1963–2005) of a rock slope at the head of the Illgraben, a very active catchment prone to debris flows in Switzerland. Slope failures affect 70% of the study slope and erode the slope at an average rate of 0.39 ± 0.03 m yr¯¹. The analysis of individual slope failures yielded an inventory of ~2500 failures ranging over 6 orders of magnitude in volume, despite the small slope area and short study period. The slope failures form a characteristic magnitude–frequency distribution with a rollover and a power‐law tail between ~200 m³ and 1.6 × 10⁶ m³ with an exponent of 1.65. Slope failure volume scales with area as a power law with an exponent of 1.1. Both values are low for studies of bedrock landslides and rockfall and result from the highly fractured and weathered state of the quartzitic bedrock. Our data suggest that the magnitude–frequency distribution is the result of two separate slope failure processes. Type (1) failures are frequent, small slides and slumps within the weathered layer of highly fractured rock and loose sediment, and make up the rollover. Type (2) failures are less frequent and larger rockslides and rockfalls within the internal bedded and fractured slope along pre‐determined potential failure surfaces, and make up the power‐law tail. Rockslides and rockfalls of high magnitude and relatively low frequency make up 99% of the total failure volume and are thus responsible for the high erosion rate. They are also significant in the context of landscape evolution as they occur on slopes above 45° and limit the relief of the slope. Copyright © 2012 John Wiley & Sons, Ltd.     

Q measurements for PhaseX overtones

Linear stacking procedures are used to retrieve the attenuation of 91 modes belonging to the 3rd, 4th and 5th Rayleigh overtones branches in the 80–160 s period range, and contributing to the so-called PhaseX wave group. Our data show in general slightly less attenuation than expected from available models. Data space inversion shows that, when combined with previously measured fundamental modeQ's, this new dataset improves resolution significantly in the 1000–2000 km depth range. Based on this remark, we carry out a number of parameter space inversions. Our results suggest a narrow (80–200 km) zone of high attenuation (Q =75–90), low attenuation in the intermediate mantle (670–1500 km); (Q 350), and lower values in the deeper mantle (Q 200).     

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.