Attenuation of P,S, and coda waves in Koyna region,India

The attenuation properties of the crust in the Koyna region of the Indian shield have been investigated using 164 seismograms from 37 local earthquakes that occurred in the region. The extended coda normalization method has been used to estimate the quality factors for P waves and S waves , and the single back-scattering model has been used to determine the quality factor for coda waves (Q _c). The earthquakes used in the present study have the focal depth in the range of 1–9 km, and the epicentral distance vary from 11 to 55 km. The values of and Q _c show a dependence on frequency in the Koyna region. The average frequency dependent relationships (Q = Q ₀ f ⁿ) estimated for the region are , and . The ratio is found to be greater than one for the frequency range considered here (1.5–18 Hz). This ratio, along with the frequency dependence of quality factors, indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of Q _c and in the present study shows that for frequencies below 4 Hz and for the frequencies greater than 4 Hz. This may be due to the multiple scattering effect of the medium. The outcome of this study is expected to be useful for the estimation of source parameters and near-source simulation of earthquake ground motion, which in turn are required in the seismic hazard assessment of a region.     

The viscosity of hydrous dacitic liquids: implications for the rheology of evolving silicic magmas

The viscosity of a series of six synthetic dacitic liquids, containing up to 5.04 wt% dissolved water, was measured above the glass transition range by parallel-plate viscometry. The temperature of the 10¹¹ Pa s isokom decreases from 1065 K for the anhydrous liquid, to 864 K and 680 K for water contents of 0.97 and 5.04 wt% H₂O. Including additional measurements at high temperatures by concentric-cylinder and falling-sphere viscometry, the viscosity (η) can be expressed as a function of temperature and water content w according to: where η is in Pa s, T is temperature in K, and w is in weight percent. Within the conditions of measurement, this parameterization reproduces the 76 viscosity data with a root-mean square deviation (RMSD) of 0.16 log units in viscosity, or 7.8 K in temperature. The measurements show that water decreases the viscosity of the dacitic liquids more than for andesitic liquids, but less than for rhyolites. At low temperatures and high water contents, andesitic liquids are more viscous than the dacitic liquids, which are in turn more viscous than rhyolitic liquids, reversing the trend seen for high temperatures and low water contents. This suggests that the relative viscosity of different melts depends on temperature and water content as much as on bulk melt composition and structure. At magmatic temperatures, rhyolites are orders of magnitude more viscous than dacites, which are slightly more viscous than andesites. During degassing, all three liquids undergo a rapid viscosity increase at low water contents, and both dacitic and andesitic liquids will degas more efficiently than rhyolitic liquids. During cooling and differentiation, changing melt chemistry, decreasing temperature and increasing crystal content all lead to increases in the viscosity of magma (melt plus crystals). Under closed system conditions, where melt water content can increase during crystallization, viscosity increases may be small. Conversely, viscosity increases are very abrupt during ascent and degassing-induced crystallization.     

State of stress in the lithosphere: Inferences from the flow laws of olivine

The experimental flow data for rocks and minerals are reviewed and found to fit a law of the form $$\dot \varepsilon = A'\left[ {sinh (\alpha \sigma )} \right]^n \exp \left[ {{{ - (E * + PV * )} \mathord{\left/ {\vphantom {{ - (E * + PV * )} {RT}}} \right. \kern-\nulldelimiterspace} {RT}}} \right]$$ where \(\dot \varepsilon \) This law reduces to the familiar power-law stress dependency at low stress and to an exponential stress dependency at high stress. Using the material flow law parameters for olivine, stress profiles with depth and strain rate are computed for a representative range of temperature distributions in the lithosphere. The results show that the upper 15 to 25 km of the oceanic lithosphere must behave elastically or fail by fracture and that the remainder deforms by exponential law flow at intermediate depths and by power-law flow in the rest. A model computation of the gravitational sliding of a lithospheric plate using olivine rheology exhibits a very sharp decoupling zone which is a consequence of the combined effects of increasing stress and temperature on the flow law, which is a very sensitive function of both.     

The TKE dissipation rate in the northern South China Sea

The microstructure measurements taken during the summer seasons of 2009 and 2010 in the northern South China Sea (between 18°N and 22.5°N, and from the Luzon Strait to the eastern shelf of China) were used to estimate the averaged dissipation rate in the upper pycnocline 〈ε _p〉 of the deep basin and on the shelf. Linear correlation between 〈ε _p〉 and the estimates of available potential energy of internal waves, which was found for this data set, indicates an impact of energetic internal waves on spatial structure and temporal variability of 〈ε _p〉. On the shelf stations, the bottom boundary layer depth-integrated dissipation \( {\widehat{\varepsilon}}_{\mathrm{BBL}} \) reaches 17–19 mW/m², dominating the dissipation in the water column below the surface layer. In the pycnocline, the integrated dissipation \( {\widehat{\varepsilon}}_{\mathrm{p}} \) was mostly ∼10–30 % of \( {\widehat{\varepsilon}}_{\mathrm{BBL}} \). A weak dependence of bin-averaged dissipation \( \overline{\varepsilon} \) on the Richardson number was noted, according to \( \overline{\varepsilon}={\varepsilon}_0+\frac{\varepsilon_{\mathrm{m}}}{{\left(1+ Ri/R{i}_{\mathrm{cr}}\right)}^{1/2}} \), where ε ₀ + ε _m is the background value of \( \overline{\varepsilon} \) for weak stratification and Ri _cr = 0.25, pointing to the combined effects of shear instability of small-scale motions and the influence of larger-scale low frequency internal waves. The latter broadly agrees with the MacKinnon–Gregg scaling for internal-wave-induced turbulence dissipation.

     

Electromagnetic induction in the earth due to stationary and moving sources

A new approach to the theory of electromagnetic induction is developed that is applicable to moving as well as stationary sources. The source field is considered to be a standing wave generated by two waves travelling in opposite directions along the surface of the earth. For a stationary source the incident waves have velocities of the same magnitude, however for a moving source the velocities of the two incident waves are respectively increased and decreased by the velocity of the source. Electromagnetic induction in the earth is then considered as refraction of these waves and gives, for both stationary and moving sources, the magnetotelluric relation: $$\frac{{ - E_y }}{{H_x }} = \left( {\frac{{i\omega \mu }}{\sigma }} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} \left( {1 - i\frac{{v^2 }}{{\omega \mu \sigma }}} \right)^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} $$ where ν is the wavenumber of the source, μ is the permeability (4π·10^?7) and σ is the conductivity of the earth. ω is the angular frequency of the variation observed on the earth. For a stationary source the observed frequency is the same as the source frequency, however the effect of moving a time-varying source is to make the observed frequency different from the frequency of the source. Failure to recognise this in previous studies led to some erroneous conclusions. This study shows that a moving source isnot “electromagnetically broader” than a stationary source as had been suggested.     

“Pandora box” — Matrix in the calculus of observations

Summary If the condition R(A)=k(n), whereA is the design matrix of the type n × k and k the number of parameters to be determined, is not satisfied, or if the covariance matrixH is singular, it is possible to determine the adjusted value of the unbiased estimable function of the parameters f(), its dispersion D( (x)) and ² as the unbiased estimate of the value of ² by means of an arbitrary g-inversion of the matrix . The matrix , because of its remarkable properties, is called the Pandora Box matrix. The paper gives the proofs of these properties and the manner in which they can be employed in the calculus of observations.     

Generous statistical tests

A common statistical problem is deciding which of two possible sources, A and B, of a contaminant is most likely the actual source. The situation considered here, based on an actual problem of polychlorinated biphenyl contamination discussed below, is one in which the data strongly supports the hypothesis that source A is responsible. The problem approach here is twofold: One, accurately estimating this extreme probability. Two, since the statistics involved will be used in a legal setting, estimating the extreme probability in such a way as to be as generous as is possible toward the defendant’s claim that the other site B could be responsible; thereby leaving little room for argument when this assertion is shown to be highly unlikely. The statistical testing for this problem is modeled by random variables {X _i } and the corresponding sample mean the problem considered is providing a bound ɛ for which for a given number a ₀. Under the hypothesis that the random variables {X _i } satisfy E(X _i ) ≤ μ, for some 0  < μ < 1, statistical tests are given, described as “generous”, because ɛ is maximized. The intent is to be able to reject the hypothesis that a ₀ is a value of the sample mean while eliminating any possible objections to the model distributions chosen for the {X _i } by choosing those distributions which maximize the value of ɛ for the test used.     

Basaltic thermals and Subplinian plumes: Constraints from acoustic measurements at Shishaldin volcano, Alaska

The 1999 basaltic eruption of Shishaldin volcano (Alaska, USA) included both Strombolian and Subplinian activity, as well as a “pre-Subplinian” phase interpreted as the local coalescence within a long foam in the conduit. Although few visual observations were made of the eruption, a great deal of information regarding gas velocity, gas flux at the vent and plume height may be inferred by using acoustic recordings of the eruption. By relating acoustic power to gas velocity, a time series of gas velocity is calculated for the Subplinian and pre-Subplinian phases. These time series show trends in gas velocity that are interpreted as plumes or, for those signals lasting only a short time, thermals. The Subplinian phase is shown to be composed of a thermal followed by five plumes with a total expelled gas volume of .The initiation of the Subplinian activity is probably related to the arrival of a large overpressurised bubble close to the top of the magma column. A gradual increase in low-frequency (0.01–0.5 Hz) signal prior to this “trigger bubble” may be due to the rise of the bubble in the conduit. This delay corresponds to a reservoir located at ≈3.9 km below the surface, in good agreement with studies on other volcanoes.The presence of two thermal phases is also identified in the middle of the pre-Subplinian phase with a total gas release of and . Gas velocity at the vent is found to be and for the Subplinian plumes and the pre-Subplinian thermals respectively.The agreement is very good between estimates of the gas flux from modelling the plume height and those obtained from acoustic measurements, leading to a new method by which eruption physical parameters may be quantified. Furthermore, direct measurements of gas velocity can be used for better estimates of the flux released during the eruption.     

A semiempirical mathematical model of the secondary pollution of water bodies by soluble iron and manganese forms

A semiempirical mathematical model of iron and manganese migration from bottom sediments into the water mass of water bodies has been proposed based on some basic regularities in the geochemistry of those elements. The entry of dissolved forms of iron and manganese under aeration conditions is assumed negligible. When dissolved-oxygen concentration is <0.5 mg/L, the elements start releasing from bottom sediments, their release rate reaching its maximum under anoxic conditions. The fluxes of dissolved iron and manganese (Me) from bottom sediments into the water mass (J _Me) are governed by the gradients of their concentrations in diffusion water sublayer adjacent to sediment surface and having an average thickness of h = 0.025 cm: \({J_{Me}} = - {D_{Me}}\frac{{{C_{Me\left( {ss} \right)}} - {C_{Me\left( w \right)}}}}{h}\) (D _Me ≈ 1 × 10^–9 m²/s is molecular diffusion coefficient of component Me in solution; C _Me(ss) and C _Me(w) ≈ 0 are Me concentrations on sediment surface, i.e., on the bottom boundary of the diffusion water sublayer, and in the water mass, i.e., on the upper boundary of the diffusion water sublayer). The value of depends on water saturation with dissolved oxygen (\({\eta _{{O_2}}}\)) in accordance with the empiric relationship \({C_{Me\left( {ss} \right)}} = \frac{{C_{_{Me\left( {ss} \right)}}^{\max }}}{{1 + k{\eta _{{O_2}}}}}\) (k is a constant factor equal to 300 for iron and 100 for manganese; C _Me(ss) ^max is the maximal concentration of Me on the bottom boundary of the diffusion water sublayer with C _Fe(ss) ^max ≈ 200 μM (11 mg/L), and C _Mn(ss) ^max ≈ 100 μM (5.5 mg/L).     

Ground-motion Attenuation Relation from Strong-motion Records of the 2001 Mw 7.7 Bhuj Earthquake Sequence (2001–2006), Gujarat, India

Predictive relations are developed for peak ground acceleration (PGA) from the engineering seismoscope (SRR) records of the 2001 M_w 7.7 Bhuj earthquake and 239 strong-motion records of 32 significant aftershocks of 3.1 ≤ M_w ≤ 5.6 at epicentral distances of 1 ≤ R ≤ 288 km. We have taken advantage of the recent increase in strong-motion data at close distances to derive new attenuation relation for peak horizontal acceleration in the Kachchh seismic zone, Gujarat. This new analysis uses the Joyner-Boore’s method for a magnitude-independent shape, based on geometrical spreading and anelastic attenuation, for the attenuation curve. The resulting attenuation equation is,

Steady flows driven by sources of random strength in heterogeneous aquifers with application to partially penetrating wells

Average steady source flow in heterogeneous porous formations is modelled by regarding the hydraulic conductivity K(x) as a stationary random space function (RSF). As a consequence, the flow variables become RSFs as well, and we are interested into calculating their moments. This problem has been intensively studied in the case of a Neumann type boundary condition at the source. However, there are many applications (such as well-type flows) for which the required boundary condition is that of Dirichlet. In order to fulfill such a requirement the strength of the source must be proportional to K(x), and therefore the source itself results a RSF. To solve flows driven by sources whose strength is spatially variable, we have used a perturbation procedure similar to that developed by Indelman and Abramovich (Water Resour Res 30:3385–3393, 1994) to analyze flows generated by sources of deterministic strength. Due to the linearity of the mathematical problem, we have focused on the explicit derivation of the mean head distribution G _d (x) generated by a unit pulse. Such a distribution represents the fundamental solution to the average flow equations, and it is termed as mean Green function. The function G _d (x) is derived here at the second order of approximation in the variance σ² of the fluctuation (where K _A is the mean value of K(x)), for arbitrary correlation function ρ(x), and any dimensionality d of the flow domain. We represent G _d (x) as product between the homogeneous Green function G _d ⁽⁰⁾(x) valid in a domain with constant K _A , and a distortion term Ψ _d (x) = 1 + σ²ψ _d (x) which modifies G _d ⁽⁰⁾(x) to account for the medium heterogeneity. In the case of isotropic formations ψ _d (x) is expressed via one quadrature. This quadrature can be analytically calculated after adopting specific (e.g.. exponential and Gaussian) shape for ρ(x). These general results are subsequently used to investigate flow toward a partially-penetrating well in a semi-infinite domain. Indeed, we construct a σ²-order approximation to the mean as well as variance of the head by replacing the well with a singular segment. It is shown how the well-length combined with the medium heterogeneity affects the head distribution. We have introduced the concept of equivalent conductivity K ^eq(r,z). The main result is the relationship where the characteristic function ψ^(w)(r,z) adjusts the homogeneous conductivity K _A to account for the impact of the heterogeneity. In this way, a procedure can be developed to identify the aquifer hydraulic properties by means of field-scale head measurements. Finally, in the case of a fully penetrating well we have expressed the equivalent conductivity in analytical form, and we have shown that (being the effective conductivity for mean uniform flow), in agreement with the numerical simulations of Firmani et al. (Water Resour Res 42:W03422, 2006).     

The sunspot areas and the relative sunspot numbers

Summary Within each sunspot cycle the yearly means (A) of the daily sunspot areas increase faster than the corresponding sunspot numbers (R) from the minimum to the maximum of solar activity and then decrease also faster than these numbers till the next minimum. Relation (A)=16.7 (R), frequently used so far, is approximately valid only for the years in the vicinity of the sunspot maximum. Instead of that, author gives the relations: for the years preceding the sunspot maximum, for the years following the sunspot maximum, wherea andb are constants,T _R is the time of rise of the corresponding sunspot cycle expressed in years, andk takes the valuek=0 for the year of maximum solar activity andk=1, 2, 3, ... for the first, second, third ... year preceding or following that of maximum solar activity. The monthly means of the daily sunspot areas show a similar variation, but in this case the ratioq=AR varies with a greater amplitude both within each sunspot cycle and from cycle to cycle. The values ofq corresponding to all months of a given year in the sunspot cycle are contained between two limits depending on the time of rise.

---

Résumé Les valeurs moyennes (A) des aires diurnes des taches solaires à chaque année depuis 1878 augmentent plus rapidement du minimum vers le maximum de l'activité solaire que les nombres de Wolf correspondants (R). Elles diminuent aussi plus rapidement que les nombres de Wolf du maximum vers le minimum de l'activité solaire. La relation adoptée (A)=16.7 (R) ne s'applique pas avec une approximation satisfaisante que seulement pour les années voisines celle du maximum de l'activité solaire. L'auteur propose les relations: pour les années qui précédent le maximum, pour les années qui suivent le maximum, oùa, b sont des constantes,T _R le temps d'ascension du cycle correspondant exprimé en années et la parametrek prend la valeurk=0 à l'année du maximum de l'activité solaire etk=1, 2, 3 ... pour les années qui précédent et qui suivent celle du maximum. Les valeurs moyennes des aires diurnes des taches à chaque mois, suivent la même marche mais dans ce cas le rapportq=AR present des larges variations. Il oscille pourtant extre deux limites qui dependent du temps d'ascension.

     

Stochastic finite-fault simulation of ground motion from the August 11, 2012, <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 6.4 Ahar earthquake,northwestern Iran

In this study, the 11 August 2012 M _w 6.4 Ahar earthquake is investigated using the ground motion simulation based on the stochastic finite-fault model. The earthquake occurred in northwestern Iran and causing extensive damage in the city of Ahar and surrounding areas. A network consisting of 58 acceleration stations recorded the earthquake within 8–217 km of the epicenter. Strong ground motion records from six significant well-recorded stations close to the epicenter have been simulated. These stations are installed in areas which experienced significant structural damage and humanity loss during the earthquake. The simulation is carried out using the dynamic corner frequency model of rupture propagation by extended fault simulation program (EXSIM). For this purpose, the propagation features of shear-wave including \( {Q}_s \) value, kappa value \( {k}_0 \), and soil amplification coefficients at each site are required. The kappa values are obtained from the slope of smoothed amplitude of Fourier spectra of acceleration at higher frequencies. The determined kappa values for vertical and horizontal components are 0.02 and 0.05 s, respectively. Furthermore, an anelastic attenuation parameter is derived from energy decay of a seismic wave by using continuous wavelet transform (CWT) for each station. The average frequency-dependent relation estimated for the region is \( Q=\left(122\pm 38\right){f}^{\left(1.40\pm 0.16\right)}. \) Moreover, the horizontal to vertical spectral ratio \( H/V \) is applied to estimate the site effects at stations. Spectral analysis of the data indicates that the best match between the observed and simulated spectra occurs for an average stress drop of 70 bars. Finally, the simulated and observed results are compared with pseudo acceleration spectra and peak ground motions. The comparison of time series spectra shows good agreement between the observed and the simulated waveforms at frequencies of engineering interest.     

Estimates of fracture parameters of earthquakes

A new estimate of the fracture parameters of earthquakes is provided in this paper. By theMuskhelishvili method (1953) a number of basic relations among fracture-mechanics parameters are derived. A scheme is proposed to evaluate the slip weakening parameters in terms of fault dimension, average slip, and rise time, and the new results are applied to 49 events compiled in the earthquake catalogue ofPurcaru andBerckhemer (1982). The following empirical relations are found in the paper: $$\begin{gathered} \frac{{\tau _B - \tau _f }}{{\tau _\infty - \tau _f }} = 2.339 \hfill \\ {{\omega _c } \mathord{\left/ {\vphantom {{\omega _c } {W = 0.113}}} \right. \kern-\nulldelimiterspace} {W = 0.113}} \hfill \\ \log G_c \left( {{{dyne} \mathord{\left/ {\vphantom {{dyne} {cm}}} \right. \kern-\nulldelimiterspace} {cm}}} \right) = 2 \log L (km) + 6.167 \hfill \\ \log \delta _c (cm) = 2 \log L (km) - 1.652 \hfill \\ \end{gathered} $$ whereG _c is the specific fracture energy,ω _c the size of the slip weakening zone,δ _c the slip weakening displacement,τ _B ?τ _f the drop in strength in the slip weakening zone,τ _∞?τ _f the stress drop,L the fault length, andW the fault width. The investigation of 49 shocks shows that the range of strength dropτ _B ?τ _f is from several doze to several hundred bars at depthh<400 km, but it can be more than 10³ bars ath>500 km; besides, the range of the sizeω _c of the strength degradation zone is from a few tenths of a kilometer to several dozen kilometers, and the range of the slip weakening displacementδ _c is from several to several hundred centimeters. The specific fracture energyG _c is of the order of 10⁸ to 10¹¹ erg cm^?2 when the momentM ₀ is of the order of 10²³ to 10²⁹ dyne cm.     

Lyapunov exponent and dimension of the strange attractor of elastic frictional system

ＬｙａｐｕｎｏｖｅｘｐｏｎｅｎｔａｎｄｄｉｍｅｎｓｉｏｎｏｆｔｈｅｓｔｒａｎｅａｔｔｒａｃｔｏｒｏｆｅｌａｓｔｉｃｆｒｉｃｔｉｏｎａｌｓｙｓｔｅｍＺｈｉ－ＲｅｎＮＩＵ（牛志仁）ａｎｄＤａｎｇ－ＭｉｎＣＨＥＮ（陈党民）（ＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒ...     

Variation of Seismic Coda Wave Attenuation in the Garhwal Region, Northwestern Himalaya

Seismic coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) characteristics in the Garhwal region, northwestern Himalaya is studied using 113 short-period, vertical component seismic observations from local events with hypocentral distance less than 250?km and magnitude range between 1.0 to 4.0. They are located mainly in the vicinity of the Main Boundary Thrust (MBT) and the Main Central Thrust (MCT), which are well-defined tectonic discontinuities in the Himalayas. Coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) is estimated using the single isotropic scattering method at central frequencies 1.5, 3, 5, 7, 9, 12, 16, 20, 24 and 28?Hz using several starting lapse times and coda window lengths for the analysis. Results show that the ( $ Q_{\text{c}}^{ - 1} $ ) values are frequency dependent in the considered frequency range, and they fit the frequency power law ( $ Q_{\text{c}}^{ - 1} \left( f \right) = Q_{0}^{ - 1} f^{ - n} $ ). The Q ₀ (Q _c at 1?Hz) estimates vary from about 50 for a 10?s lapse time and 10?s window length, to about 350 for a 60?s lapse time and 60?s window length combination. The exponent of the frequency dependence law, n ranges from 1.2 to 0.7; however, it is greater than 0.8, in general, which correlates well with the values obtained in other seismically and tectonically active and highly heterogeneous regions. The attenuation in the Garhwal region is found to be lower than the Q _c ^?1 values obtained for other seismically active regions of the world; however, it is comparable to other regions of India. The spatial variation of coda attenuation indicates that the level of heterogeneity decreases with increasing depth. The variation of coda attenuation has been estimated for different lapse time and window length combinations to observe the effect with depth and it indicates that the upper lithosphere is more active seismically as compared to the lower lithosphere and the heterogeneity decreases with increasing depth.     

Attenuation of coda waves in the Northeastern Region of India

Coda wave attenuation quality factor Qc is estimated in the northeastern region of India using 45 local earthquakes recorded by regional seismic network. The quality factor Qc was estimated using the single backscattering model modified by Sato (J Phys Earth 25:27–41, 1977), in the frequency range 1–18 Hz. The attenuation and frequency dependence for different paths and the correlation of the results with geotectonics of the region are described in this paper. A total of 3,890 Qc measurements covering 187 varying paths are made for different lapse time window of 20, 30, 40, 50, 60, 70, 80, and 90 s in coda wave. The magnitudes of the analyzed events range from 1.2 to 3.9 and focal depths range between 7 and 38 km. The source–receiver distances of the selected events range between 16 and 270 km. For 30-s duration, the mean values of the estimated Qc vary from 50 ± 12 (at 1 Hz) to 2,078 ± 211(at 18 Hz) for the Arunachal Himalaya, 49 ± 14 (at 1 Hz) to 2,466 ± 197 (at 18 Hz) for the Indo-Burman, and 45 ± 13 (at 1 Hz) to 2,069 ± 198 (at 18 Hz) for Shillong group of earthquakes. It is observed that Qc increases with frequency portraying an average attenuation relation for the region. Moreover, the pattern of Qc^− 1 with frequency is analogous to the estimates obtained in other tectonic areas in the world, except with the observation that the Qc^− 1 is much higher at 1 Hz for the northeastern region. The Qc^− 1 is about 10^− 1.8 at 1 Hz and decreases to about 10^− 3.6 at 18 Hz indicating clear frequency dependence. Pertaining to the spatial distribution of Qc values, Mikir Hills and western part of Shillong Plateau are characterized by lower attenuation.