On the problem of standardization of the magnitude-frequency relation

au u uu mu a amu am a¶rt; mum ma¶rt;amuauau nmm m a, ¶rt;a m u ¶rt;u naam — aum¶rt; a, nu¶rt; a¶rt;u uu na¶rt;u am — nu muaua muam ¶rt; m ¶rt;a. ¶rt;aa n ma¶rt;amuauu nm naam ¶rt;u munaua nmmu ¶rt;.     

A note on the stress-dilatancy relation for simulated fault gouge

Theoretical constraints on the stress-dilation relation for a deforming Coulomb material requirev ifC=0 andv sin^-1( _m / _m ) always, wherev is the dilation angle, is the friction angle,C is cohesion, _m is the maximum shear stress, and _m is the mean effective stress. Recent laboratory measurements of friction and dilatancy of simulated fault gouge show that small amplitude shear-load cycling causes compaction and consolidation. Comparison of the data with theory indicates that such load cycling produces: (1) increased coefficient of friction (or friction angle), (2) increased cohesion, and (3) increased dilatancy rate (or dilation angle). Under certain conditions of load cycling without significant plastic shear strain accumulation ( ^p <0.005) we find thatv exceeds both and, in contrast to theory, sin^-1( _m / _m ). This result is interpreted in terms of enhanced cohesion and overconsolidation, which lead to residual stresses within the gouge. An analogy is drawn between these special loading conditions and those extant on natural faults. In particular, our results imply that jostling and minor stress variations associated with microearthquakes may produce strengthening of fault gouge and changes in the fault zone's stress-dilatancy relation. Hence, compaction associated with microseismicity may lead to subsequent dilation of fault gouge, even for faults with large displacement rates and large net offsets (e.g., San Andreas). In regions where such dilation persists over sufficient displacements (on the order of the critical slip distance for seismic faulting) it may tend to inhibit unstable slip.     

A model for the simulation of turbulent and eddy diffusion processes at heights of 0–65 km

A generalized turbulent diffusion model has been developed which evaluates the time rate of growth of a simulated cloud of particles released into a turbulent (i.e. diffusive) atmosphere. The general model, in the form of second-order differential equations, computes the three-dimensional size of the cloud as a function of time. Parameters which influence the cloud growth, and which are accounted for in the model equations, are: (1) length scales and velocity magnitudes of the diffusive field, (2) rate of viscous dissipation , (3) vertical stability as characterized by the relative adiabatic lapse rate (1/T)(g/C _p +T/z), and (4) vertical shear in the mean horizontal winds , and , for a given height and of spatial extent equal to that of the diffusing cloud. Sample results for near ground level and for upper stratospheric heights are given. For the atmospheric boundary layer case, the diffusive field is microscale turbulence. In the upper stratospheric case it is considered to be a field of highly interactive and dispersive gravity waves.     

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.     

On the propagation of wave from a cylindrical cavity

Summary This paper discusses the disturbance produced in an infinite layer of non-homogeneous elastic material characterised by =0 ⁿ and =0ⁿ(n>0) where and are the density and shear modulus respectively of the material, due to periodic torsional force applied on the wall of a cylindrical hole in the layer. The variation of the displacement component with the radius vector is shown graphically and compared with the corresponding homogeneous case.     

Molecular refraction in the Earth's mantle

The Drude law (molecular refraction) for the temperature radiation in a monoatomic model of the Earth's mantle is derived. The considerations are based on the Lorentz electron theory of solids. The characteristic frequency (or eigenfrequency) of independent electron oscillators (in energy units, ) is identified with the band gapE _G of a solid. The only assumption is that solid material related to the Earth's mantle has the mean atomic weight A21 g/mole, and its energy gap (E _G) is about 9 eV. In this case the value of molecular refraction (in cm³/g) is (n ²–1)/=0.5160.52, where andn are the density and the refractive index at wavelength _D=0.5893 m (sodium light), respectively. The average molecular refraction of important silicate and oxide minerals with A21, obtained byAnderson andSchreiber (1965) from laboratory data, is , where denotes the mean arithmetic value calculated from three principal refractive indices of crystal. For the rock-forming minerals with 19A<24 g/mole the new relation was found byAnderson (1975).     

Propagation of love waves in a non-homogeneous stratum of finite depth sandwitched between two semi-infinite isotropic media

Summary The aim of this paper is to study a problem in which the intermediate layer is non-homogeneous, the rigidity varying exponentially with depth i.e. ₂=₂ v ₀ ² e ^2pz , the density being constant, velocity varies also exponentially with depth according to the law =v ₀ e ^pz . The variability ofKH with the change of phase velocity is shown graphically.     

Analysis of aeromagnetic anomalies over the central part of the Narmada-Son lineament

The Narmada-Son lineament (NSL) forms a major tectonic feature on the Indian subcontinent. The importance of this lineament lies in its evolution as well as its tectonic history. The lineament seems to have been active since Precambrian times. In order to understand the history of its evolution, it is necessary to know what igenous activity has been taking place along this lineament, and how the Deccan trap volcanics, which cover large areas along this lineament, have erupted.For the study of this problem an analysis of the aeromagnetic anomaly map lying between 76°15 to 77°30E and 21°45 to 22°50N has been carried out. Four different profiles (B ₁ B ₁,B ₂ B ₂,B ₃ B ₃ andB ₄ B ₄) have been drawn in N-S direction over this area and interpreted in terms of the intrusive bodies present within or below the surface of Deccan trap exposures. Inversion and forward modelling techniques have been adopted for interpretation purposes. An analysis of frequency spectra along the profiles has also been carried out to estimate the average depth of the different magnetic bodies. These results have been correlated with the available geological information. It has been found that most of the small wavelength anomalies are caused by dyke-like bodies within or below the Deccan trap at a depth of less than 0.5 km.     

The dynamics of magma mixing in a rising magma batch

The conditions under which two magmas can become mixed within a rising magma batch are investigated by scaling analyses and fluid-dynamical experiments. The results of scaling analyses show that the fluid behaviours in a squeezed conduit are determined mainly by the dimensionless number where ₁ is the viscosity of the fluid, U is the velocity, g is the acceleration due to gravity, is the density difference between the two fluids, and R is the radius of the tube. The parameter I represents a balance between the viscous effects in the uppermost magma which prevent it from being moved off the conduit walls, and the buoyancy forces which tend to keep the interface horizontal. The experiments are carried out using fluid pairs of various density and viscosity contrasts in a squeezed vinyl tube. They show that overturning of the initial density stratification and mixing occur when I>order 10^-1; the two fluids remain stratified when I 10^-3. Transitional states are observed when 10^-3<I<10^-1. These results are nearly independent of Reynolds number and viscosity ratio in the range of and Re ₁<300. Applying these results to magmas shows that silicic to intermediate magmas overlying mafic magma will be prone to mixing in a rising magma batch. This mechanism can explain some occurrences of small-volume mixed lava flows.     

Propagation of Rayleigh waves in an axially symmetric heterogeneous layer lying between two homogeneous halfspaces

Summary The propagation of Rayleigh type waves in an axially symmetric inhomogeneous layer lying between two halfspaces is studied. The halfspaces are supposed to be identical in their elastic properties. The variation of the parameters in the layer is assumed to be of the form where is a constant andz is the distance measured from one interface into the layer. With this assumption, the vector wave equation for the layer is separable. The solution is obtained in terms ofWhittaker's functions and the frequency equation of Rayleigh type waves is derived.     

Uniform heat flow in a semi-infinite medium disturbed by a body of different thermal conductivity

Summary The approximate solution for the disturbance of a uniform heat flow in a homogeneous semi-infinite medium by an object of different thermal conductivity buried in it is generally used in the interpretation of heat flow anomalies on the ocean bottom. In order to know the accuracy of the approximate solution, a comparison between the approximate solution and the exact solution is given in the case of a very long horizontal cylinder in a semi-infinite medium. The computed results indicate that the two solutions agree to within 10% whend>1.3 and 0.5<<2, whered is ratio of the depth to the radiusR ₀ of the cylinder and is the factor of the contrast of the thermal conductivities between the medium and the body. As for the cases when 1 and 1, the same accuracy can be obtained only whend>2. A similar approach is also applied to the case of a spherical conductor in a semi-infinite medium by using a bispherical harmonic solution. The results of both the bipolar solution and the bispherical solution show that when 1 andd1, the vertical thermal gradient at the surface of the semi-infinite medium is always positive and tends to zero, but a negative vertical gradient may be obtained for the approximate solutions.     

The intensity of the main phase of geomagnetic storms in relation to the parameters of the solar wind

Summary On the basis of investigating 10 storms (1965–1967) good correlation was found between the density of the solar wind energy (₂=1/2mNv²) and the intensity of the main phase of the geomagnetic storms, expressed in terms of the maximum decrease of the horizontal intensity (B=H/cos). The relation between ₂, or Nv², and B could then be used to determine the quantities and ₀ ( is the factor expressing the increase in energy density in the magnetosphere, ₀ is the energy density of the particles in a quiet magnetosphere). A comparison with the directly observed distribution of the energy density of the particles in the magnetosphere indicates that the computed value of ₀ seems to be realistic. The magnitude of the factor will have to be checked again.     

Nonlinear Magnetoconvection and the Geostrophic Flow - Subcritical Solutions

The magnetoconvection problem under the magnetostrophic approximation is investigated as the nonlinear regime is entered. The model consists of a fluid filled sphere, internally heated, and rapidly rotating in the presence of a prescribed, axisymmetric, toroidal magnetic field. For simplicity only a dipole parity and a single azimuthal wavenumber (m = 2) is considered here. The leading order nonlinearity at small amplitude is the geostrophic flow U _g which is introduced to the previously linear model (Walker and Barenghi, 1997a, b). Walker and Barenghi (1997c) considered parameter space above critical and found that U _g acts as an equilibration mechanism for moderately supercritical solutions. However, for solutions well above critical a Taylor state is approached and the system can no longer equilibrate. More importantly though, in the context of this paper, is that subcritical solutions were found. Here subcritical solutions are considered in more detail. It was found that, at is strongly dependent on . ( is the critical value of the modified Rayleigh number is a measure of the maximum amplitude of the generated geostrophic flow while , the Elsasser number, defines the strength of the prescribed toroidal field.) R_m at proves to be the key measure in determining how far into the subcritical regime the system can advance.     

Integral Equations for Kinematic Dynamo Models

A new technique for the treatment of the kinematic dynamo problem is presented. The method is applicable when the dynamo is surrounded by a medium of finite conductivity and is based on a reformulation of the induction equation and boundary conditions at infinity into an integral equation. We show that the integral operator involved here is compact in the case of homogeneous conductivity, which is important for both mathematical and numerical treatment. A lower bound for the norm of then yields a necessary condition for the generation of magnetic fields by kinematic dynamos. Numerical results are presented for some simple ²-dynamo models. The far-field asymptotics for stationary and time-dependent field modes are discussed.     

Variation of continental mantle heat flow with age. Possibility of discriminating between thermal models of the lithosphere

Continental mantle heat flow values are obtained by subtracting the radiogenic heat produced in the lower crust and lithosphere beneath the crust from reduced heat flow values reported for various heat flow provinces. The significance of continental mantle heat flow values thus obtained is that they can be considered essentially as representing the residual heat of cooling of the continental lithosphere. A plot of these mantle heat flow values against 1/t where t is the geologic age of the last thermal event suggests a linear trend. It is also found that the recently proposed relationQ=500 (1/t) for the variation of oceanic heat flowQ (in mW/M²) with aget (in million years) provides a reasonably good fit to the mantle heat flow data. The constant thickness plate model however, is found to be unsatisfactory in explaining the variation of continental mantle heat flow with age.     

Flow velocity profiles in the Lower Scheldt estuary

Recent acoustic Doppler current profiler (ADCP)-measurements in the Scheldt estuary near Antwerp, Belgium, revealed anomalous, i.e. anti-clockwise circulations in a left bend during the major part of the flood period; these circulations were established shortly after the turn of the tide. During ebb, anti-clockwise circulations persisted, as predicted by classical theory. These data were analysed with a 3D and a 1DV-model. The 3D simulations reveal that the anomalous circulations are found when salinity is included in the computations—without salinity “normal” circulations were found. From analytical and 1DV simulations, it is concluded that a longitudinal salinity gradient ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-0em} {\partial x}$ may induce a near-bed maximum in flow velocity reversing the direction of the secondary currents. The 1DV-model was then used to assess the contribution of various processes one by one. It was found that because of a reduction in vertical mixing, the vertical velocity profile is not at equilibrium during the first phase of accelerating tide, further enhancing the effects of ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-0em} {\partial x}$ . A small vertical salinity gradient ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial z}}} \right. \kern-0em} {\partial z}$ appeared to have a very large effect as the crosscurrents of the secondary circulations induced by ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-0em} {\partial x}$ became an order of larger magnitude. However, at the site under consideration, the effects of transverse salinity gradients, generated by differential advection in the river bend, were dominant: adverse directions of the secondary circulations were found even when the vertical velocity profile became more regular with a more or less logarithmic shape, i.e. when the effects of ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-0em} {\partial x}$ and ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial z}}} \right. \kern-0em} {\partial z}$ did not play a dominant role anymore. It is argued that data on the secondary velocity structure, which can be measured easily owing to today’s developments in ADCP equipment, may serve as an indicator for the accuracy at which the salinity field is computed with 3D numerical models. Moreover, the large effect of the salinity structure on the velocity field must have a large impact on the morphological development of estuaries, which should therefore be accounted for in morphological modelling studies.     

Coseismic slip in the 1964 Prince William Sound earthquake: A new geodetic inversion

The 1964 Prince William Sound earthquake (March 28, 1964;M _w =9.2) caused crustal deformation over an area of approximately 140,000 km² in south central Alaska. In this study geodetic and geologic measurements of this surface deformation were inverted for the slip distribution on the 1964 rupture surface. Previous seismologic, geologic, and geodetic studies of this region were used to constrain the geometry of the fault surface. In the Kodiak Island region, 28 rectangular planes (50 by 50 km each) oriented 218°N, with a dip varying from 8^o nearest the Aleutian trench to 9^o below Kodiak Island, define the rupture surface. In the Prince William Sound region 39 planes with variable dimensions (40 by 50 km near the trench, 64 by 50 km inland) and orientation (218°N in the west and 270°N in the east) were used to approximate the complex faulting. Prior information was introduced to constrain offshore dip-slip values, the strike-slip component, and slip variation between adjacent planes. Our results suggest a variable dip-slip component with local slip maximums occurring near Montague Island (up to 30 m), further to the east near Kayak Island (up to 14 m), and trenchward of the northeast segment of Kodiak Island (up to 17m). A single fault plane dipping 30°NW, corresponding to the Patton Bay fault, with a slip value of 8 m modeled the localized but large uplift on Montague Island. The moment calculated on the basis of our geodetically derived slip model of 5.0×10²⁹ dyne cm is 30% less than the seismic moment of 7.5×10²⁹ dyne cm calculated from long-period surface waves (Kanamori, 1970) but is close to the seismic moment of 5.9×10²⁹ dyne cm obtained byKikuchi andFukao (1987).     

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).     

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).