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.     

Dependence of the concentration of a passive substance on the duration of the measurement

Zusammenfassung Der Artikel hatte zur Aufgabe den Einfluss von Ver?nderung Beobachtungs-Intervalls auf die Konzentration der passiven, von erh?hter und steitiger Punktquelle ausgelassenen Substanz auszudrücken. Zu diesem Zwecke wurde ein Modell der sog, nichtstation?ren Rauchspur—einer durchschnittlichen Rauchspur im Messintervall T—konstruiert. Man setzt voraus, dass die Verteilung von Teilchen in der nichtstation?ren Rauchspur eine Gauss'sche ist mit der Dispersion in den Richtungen y, z. Für die Dispersion wurde die Beziehung (15) abgeleitet. Analoge Beziehung gilt auch in der Richtung z. Wenn uns die Auslenkung der Achse der nichtstation?ren Rauchspur von der Achse der station?ren Rauchspur, der Geraden (t,O, H), bekannt ist, so k?nnen wir die Konzentration q_T(t,y,z) explizit ausdrücken.

---

---

Address: Dúbravská cesta 4, Bratislava-Patr?nka.     

Magnetoconvection in Dependence on Roberts Number

Marginal instabilities in a horizontal unstably stratified rapidly rotating fluid layer permeated by an azimuthal magnetic field growing linearly with distance from the vertical rotation axis, are investigated in dependence on electromagnetic boundary conditions and Roberts number . Electrical conductivities of perfectly thermally conducting boundaries are either infinite or finite Earthlike ones. Following Soward's (1979) results the case of gives the instabilities of MAC-waves with much greater frequencies than for the case of q 1. Critical frequencies of MAC waves are most sensitive on Roberts number for q 1 and lose this sensitivity for q2. On the other hand the critical Rayleigh numbers are slightly dependent on q 1. Many qualitative differences in the dependences on q do exist for infinitely and finitely conducting boundaries. The former do not allow fast eastward modes for q 1 and the latter, e.g. allow eastward MC-waves with azimuthal wave number m = 1. The MC-waves are almost independent on Roberts number q.     

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

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

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.

     

Empirical ground-motion relations using moderate earthquakes recorded by Medellín–Aburrá Valley (Colombia) strong-motion networks

We tested attenuation relations obtained for different regions of the world to verify their suitability to predict strong-motion data recorded by Medellín and Aburrá Valley Accelerographic Networks. We used as comparison criteria, the average of the difference between the observed and the predicted data as a function of epicenter distance and its standard deviation. We also used the approach developed by Sherbaum et al. (Bull Seism Soc Am 94:2164–2185, 2004) that provides a method to evaluate the overall goodness-of-fit of ground-motion prediction equations. The predictive models selected use a generic focal depth. We found that this parameter has an important influence in the ground-motion predictions and must be taken into account as an independent variable. We also found important to characterize the local soil amplification to improve the attenuation relations. We found empirical relations for peak horizontal acceleration PGA and velocity PGV based on the Kamiyama and Yanagisawa (Soils Found 26:16–32, 1986) approach. $$\begin{aligned} \log _{10} (PGA)=0.5886M_L -1.0902\log _{10}(R)-0.0035H+C_{st}\pm 0.\text{29} \end{aligned}$$ $$\begin{aligned} \log _{10} (PGV)=0.7255M_L -1.8812\log _{10}(R)-0.0016H+C_{st}\pm 0.36 \end{aligned}$$ where PGA is measured in cm/s $^{2}$ and PGV in cm/s, $M_{L}$ is local magnitude in the range 2.8–6.5, $R$ is epicentral distance up to 290 km, $H$ is focal depth in km and $C_{st}$ is a coefficient that accounts for the site response due to soil conditions of each recording station. The introduction of focal depth and local site conditions as independent variables, minimize the residuals and the dispersion of the predicted data. We conclude that $H$ and $C_{st}$ are sensitive parameters, having a strong influence on the strong-motion predictions. Using the same functional form, we also propose an empirical relation for the root mean square acceleration a $_\mathrm{rms}$ : $$\begin{aligned} \log _{10} \left( {a_{rms} } \right)=0.4797M_L -1.1665\log _{10} (R)-0.00201H+C_{st}\pm 0.40 \end{aligned}$$ where a $_\mathrm{rms}$ is measured in cm/s $^{2}$ , from the S-wave arrival and using a window length equal to the rupture duration. The other variables are the same as those for PGA and PGV. The site correction coefficients $C_{st}$ found for PGA, PGV and a $_\mathrm{rms}$ show a similar trend indicating a good correlation with the soil conditions of the recording sites.     

Reoccurrence formulas for the largest earthquake magnitudes based on the modified cumulative frequency-magnitude relationship

A modified formula of the cumulative frequency-magnitude relation has been formulated and tested in a previous paper by the authors of this study. Based on the modified relationship, the following reoccurrence formulas have been obtained.

1. For the ‘T-years period’ larger earthquake magnitude,M _T $$M_T = \frac{1}{{A_3 }}ln\frac{{A_2 }}{{(1/T) + A_1 }}.$$
2. For the value of the maximum earthquake magnitude, which is exceeded with probabilityP inT-years period,M _PT $$M_{PT} = \frac{{ln(A_2 .T)}}{{A_3 }} - \frac{{ln[A_1 .T - ln(1 - P)]}}{{A_3 }}.$$
3. For the probability of occurrence of an earthquake of magnitudeM in aT-years period,P _MT $$P_{MT} = 1 - \exp [ - T[ - A_1 + A_2 \exp ( - A_3 M)]].$$

The above formulas provide estimates of the probability of reoccurrence of the largest earthquake events which are significantly more realistic than those based on the Gutenberg-Richter relationships; at least for numerous tested earthquake samples from the major area of Greece.     

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.     

Factors controlling permeability–porosity relationships in magma

A model has been developed to investigate the sensitivity of magma permeability, k, to various parameters. Power-law relationships between k and porosity J are revealed, in agreement with previous experimental and theoretical studies. These relationships take the form % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavTnhis1MBaeXatLxBI9gBae % rbd9wDYLwzYbWexLMBbXgBcf2CPn2qVrwzqf2zLnharyavP1wzZbIt % LDhis9wBH5garqqtubsr4rNCHbGeaGqiVCI8FfYJH8sipiYdHaVhbb % f9v8qqaqFr0xc9pk0xbba9q8WqFfeaY-biLkVcLq-JHqpepeea0-as % 0Fb9pgeaYRXxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaabau % aaaOqaaiqbdUgaRzaajaGaeyypa0Jaem4AaSMaei4la8IaemOCai3a % aWbaaSqabeaacqaIYaGmaaGccqGH9aqpcqWGHbqycqGGOaakcqaHgp % GzcqGHsislcqaHgpGzdaWgaaWcbaGaem4yamMaemOCaihabeaakiab % cMcaPmaaCaaaleqabaGaemOyaigaaaaa!4CE4! [^(k)] = k/r² = a(f- f_cr )^b\hat k = k/r^2 = a(\phi - \phi _{cr} )^b where r is the mean bubble radius, J_cr is the percolation threshold below which permeability is zero, and a and b are constants. It is discovered that % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavTnhis1MBaeXatLxBI9gBae % rbd9wDYLwzYbWexLMBbXgBcf2CPn2qVrwzqf2zLnharyavP1wzZbIt % LDhis9wBH5garqqtubsr4rNCHbGeaGqiVCI8FfYJH8sipiYdHaVhbb % f9v8qqaqFr0xc9pk0xbba9q8WqFfeaY-biLkVcLq-JHqpepeea0-as % 0Fb9pgeaYRXxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaabau % aaaOqaaiqbdUgaRzaajaGaeyOeI0IaeqOXdygaaa!3CDB! [^(k)] - f\hat k - \phi relationships are independent of bubble size. The percolation threshold was found to lie at around 30% porosity. Polydisperse bubble-size distributions (BSDs) give permeabilities around an order of magnitude greater than monodisperse distributions at the same porosity. If bubbles are elongated in a preferred direction then permeability in this direction is increased, but, perpendicular to this direction, permeability is unaffected. In crystal-free melts the greatest control on permeability is the ease of bubble coalescence. In viscous magmas, or when the cooling time-scale is short, bubble coalescence is impeded and permeability is much reduced. This last effect can cause variations in permeability of several orders of magnitude.     

An algorithm for the optimum distribution of a regional seismic network — II. An analysis of the accuracy of location of local earthquakes depending on the number of seismic stations

Summary Seven optimal networks consisting of 4 to 10 stations are compared for a given region, where velocity-depth profiles and the distribution of seismic intensity are known. Assuming that the standard error of arrival time is _t =0.05 s and the standard errors of the parameters of velocity-depth profiles are equal to 5% of their values, the average standard errors of the origin time and focus coordinates are estimated. The application of optimum methods to the planning of seismic networks in the Lublin Coal Basin is presented, and maps of standard errors of origin time , depth and epicenter ( _xy ) for the case of an optimum network of 6 seismic stations are given.     

Empirical ground-motion models for point- and extended-source crustal earthquake scenarios in Europe and the Middle East

This article presents the latest generation of ground-motion models for the prediction of elastic response (pseudo-) spectral accelerations, as well as peak ground acceleration and velocity, derived using pan-European databases. The models present a number of novelties with respect to previous generations of models (Ambraseys et al. in Earthq Eng Struct Dyn 25:371–400, 1996, Bull Earthq Eng 3:1–53, 2005; Bommer et al. in Bull Earthq Eng 1:171–203, 2003; Akkar and Bommer in Seismol Res Lett 81:195–206, 2010), namely: inclusion of a nonlinear site amplification function that is a function of $\text{ V }_\mathrm{S30}$ and reference peak ground acceleration on rock; extension of the magnitude range of applicability of the model down to $\text{ M }_\mathrm{w}$ 4; extension of the distance range of applicability out to 200 km; extension to shorter and longer periods (down to 0.01 s and up to 4 s); and consistent models for both point-source (epicentral, $\text{ R }_\mathrm{epi}$ , and hypocentral distance, $\text{ R }_\mathrm{hyp}$ ) and finite-fault (distance to the surface projection of the rupture, $\text{ R }_\mathrm{JB}$ ) distance metrics. In addition, data from more than 1.5 times as many earthquakes, compared to previous pan-European models, have been used, leading to regressions based on approximately twice as many records in total. The metadata of these records have been carefully compiled and reappraised in recent European projects. These improvements lead to more robust ground-motion prediction equations than have previously been published for shallow (focal depths less than 30 km) crustal earthquakes in Europe and the Middle East. We conclude with suggestions for the application of the equations to seismic hazard assessments in Europe and the Middle East within a logic-tree framework to capture epistemic uncertainty.     

Aspherical correction for free vibration of elastically anisotropic solid by means of theXYZ algorithm

By applying the perturbation theory to theXYZ algorithm (a kind of variational method), the difference f in free vibration frequencies between sphere and ellipsoid was approximated as , where i and _i (i = x,y andz) (i=x, y andz) are aspherical coefficients and asphericities of the ellipsoid, respectively. We developed an analytic method to compute the aspherical coefficients719-4 by using theXYZ algorithm. A numerical example was given for an ellipsoidal olivine, and an attempt was made to estimate the asphericities of the specimen by a least-squares method, based on the relationship between frequency shift and asphericity.     

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

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.     

Estimation of coda wave attenuation in East Central Iran

The attenuation of coda waves, Q _c , has been estimated in Zarand, Jiroft, and Bam regions of east central Iran using a single back-scattering model of S-coda envelopes. For this purpose, the recordings of 97 earthquakes by three seismic networks and a local strong ground motion network have been used. In this research, the frequency-dependent Q _c values are estimated at central frequencies of 1.5, 3, 6, 8, 12, 16, and 24 Hz using different lapse time windows from 20 to 60 s. The frequency-dependent relationships obtained are for Zarand, for Jiroft, and for Bam region. From the strong ground motion data, we obtain the relation . The Q _c frequency-dependent relationship for the entire region of east central Iran from all data (both seismograms and accelerograms) is . The average Q _c values estimated and their frequency dependent relationships correlate well with a highly heterogeneous and highly tectonically active region. Results also show that the attenuation is higher in Bam region compared to Zarand and Jiroft regions.     

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.     

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.     

Magnetic interpretation of dyke anomalies using derivatives

The horizontal and vertical derivative profiles of magnetic anomalies of dykes show some interesting properties. The points of zero derivatives and the points where the derivatives are equal are conjugate point pairs. A method of interpretation of dyke anomalies is suggested, which utilizes the distances between these points.Notation F Magnetic anomaly in total intensity - Z Depth to top of the dyke - 2T Width of the dyke - Geological dip of the dyke - I Effective intensity of magnetisation in the plane of profile - Dip of effective magnetisation vector in the plane of profile - Strike angle of the dyke - i Magnetic dip - Q – - Q _f –+arctan (sin coti) - I _f      

Zur Interpretation von seismischen Refraktionsmessungen

Zusammenfassung Unter der Voraussetzung, dass die Frontgeschwindigkeitc eine stetige und monoton wachsende Funktion von der Tiefez ist, wird dargelegt, wie man aus einer gemessenen Laufzeitkurve () die zuc inverse Funktionz=z (c) auf einfache Weise berechnen kann. Weiter wird die Eindringtiefez _m in Funktion von ermittelt und abschliessend ein Beispiel gegeben.

---

Summary Based on a recorded travel-time curve (), a simple direct method is developed for calculating the functionz=z (c), under the asumption that the wave velocityc is a regularly monotone increasing function of the depthz. Finally a rumerical example is given.

     

On the seismic response of under-designed caisson foundations

The seismic behaviour of caisson foundations supporting typical bridge piers is analysed with 3D finite elements, with due consideration to soil and interface nonlinearities. Single-degree-of freedom oscillators of varying mass and height, simulating heavily and lightly loaded bridge piers, founded on similar caissons are studied. Four different combinations of the static ( $\text{ FS }_\mathrm{V}$ FS V ) and seismic ( $\text{ FS }_\mathrm{E}$ FS E ) factors of safety are examined: (1) a lightly loaded ( $\text{ FS }_\mathrm{V}= 5$ FS V = 5 ) seismically under-designed ( $\text{ FS }_\mathrm{E} < 1$ FS E < 1 ) caisson, (2) a lightly loaded seismically over-designed ( $\text{ FS }_\mathrm{E} >1$ FS E > 1 ) caisson, (3) a heavily loaded ( $\text{ FS }_\mathrm{V} = 2.5$ FS V = 2.5 ) seismically under-designed ( $\text{ FS }_\mathrm{E} < 1$ FS E < 1 ) caisson and (4) a heavily loaded seismically over-designed caisson. The analysis is performed with use of seismic records appropriately modified so that the effective response periods (due to soil-structure-interaction effects) of the studied systems correspond to the same spectral acceleration, thus allowing their inelastic seismic performance to be compared on a fair basis. Key performance measures of the systems are then contrasted, such as: accelerations, displacements, rotations and settlements. It is shown that the performance of the lightly loaded seismically under-designed caisson is advantageous: not only does it reduce significantly the seismic load to the superstructure, but it also produces minimal residual displacements of the foundation. For heavily loaded foundations, however ( $\text{ FS }_{V} = 2.5$ FS V = 2.5 ), the performance of the two systems (over and under designed) is similar.