Dependence of the energy released during earthquake on ambient shear stress

Starting with dislocation model, using the result of the fracture mechanics: the slip displacement at the crack tip is proportional to the length of the crack and the applied ambient shear stressτ ₀ ² , we consider the dislocation in the earthquake to be the slip displacement at the crack tip and have obtained the analysis expresses of displacement and velocity pulse for the circular crack and have calculated the seismic wave energy radiated by earthquake. The seismic wave energyE ∞M ₀ τ ₀ ² f(v) _r , i. e.E is proportional to the seismic momentM ₀ and the square of the ambient shear stressτ ₀ ² and increases with the rupture velocityv _r . In frequency domain, integrating the square of source velocity spectrum derived from our the scaling law model, we have also obtained the seismic wave energyE released by earthquake and earthquake radiated effficiencyη.E ∞M ₀ τ ₀ ² also. If takingτ ₀ = 10.0 MPa, E=4.79M ₀. This result is consistent with the estimate by Vassiliou and Kanamori (1982). Theη=5.26%. The distribution of the seismic wave energy is that most of the energy contains in the frequency range between the first corner frequencyf _c1 and thirdf _c3, amount to 92.3% the energy in the rangef<f _c1 is about 3.85% and 3.85% whenf>f _c3. Thef _c3 is about 8Hz forM ⩾ 6, thus most of radiated energy is below 2Hz. This phenomenon had been verified by Vassiliou Kanamori. Previous results show the energy radiated by earthquake to be strongly dependent on ambient shear stress. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 146–152, 1993. This work was supported by the Deutsche Forschungsgemeinschaft, Bonn, F. R. Germany. The support is grateful acknowledged. The authors are also grateful to Professor Klussmann and Mr. Hasthoff for their lots of help.     

Instability of flows near the onset of convection in a rotating layer with stress-free horizontal boundaries

We investigate instability of convective flows of simple structure (rolls, standing and travelling waves) in a rotating layer with stress-free horizontal boundaries near the onset of convection. We show that the flows are always unstable to perturbations, which are linear combinations of large-scale modes and short-scale modes, whose wave numbers are close to those of the perturbed flows. Depending on asymptotic relations of small parameters α (the difference between the wave number of perturbed flows and the critical wave number for the onset of convection) and ε (ε² being the overcriticality and the perturbed flow amplitude being O(ε)), either small-angle or Eckhaus instability is prevailing. In the case of small-angle instability for rolls the largest growth rate scales as ε^8/5, in agreement with results of Cox and Matthews (Cox, S.M. and Matthews, P.C., Instability of rotating convection. J. Fluid. Mech., 2000, 403, 153–172) obtained for rolls with k = k _c . For waves, the largest growth rate is of the order ε^4/3. In the case of Eckhaus instability the growth rate is of the order of α².     

Discussion: ‘Macroscale surface roughness and frictional resistance in overland flow’ by D. S. L. Lawrence

Lawrence argued that the inundation ratio Λ, defined as the mean flow depth d divided by the roughness height k, is the dominant control of flow resistance f and should be used as the primary variable when evaluating the hydraulics of overland flow on rough surfaces. Lawrence defined three flow regimes on the basis of Λ and developed an expression for f in terms of Λ for each regime. Common sense, however, suggests that f is independent of Λ where Λ < 1 because when roughness elements protrude through the flow, the value of f for the flow is the same regardless of the height of the elements. The error appears to have crept in as a result of Lawrence's representation of roughness elements by hemispheres. Lawrence found that f ∝ d/k, which she interpreted to mean f ∝ Λ. However, in her model the length dimension denoted by k is in fact half the breadth b/2 of the roughness elements. The distinction between k and b/2 is important, especially for roughness elements where k ≠ b/2. Thus, contrary to Lawrence's claim, f is not generally a function of Λ. Instead, f is a function of Λ only where Λ > 1. Where Λ < 1, f is a function of d/(b/2) or d/b. © 1998 John Wiley & Sons, Ltd.     

Energy of acoustic signals in a borehole

In the present article, the dependencies of the acoustic signal total energy and the energies of the wave packets of different types of the waves on the elastic parameters and permeability of rocks have been studied. We have considered traditional logging tools containing acoustical monopole source. Calculations were performed in a frequency range of dozens of kilohertz, typical for acoustic well logging. It was shown that in a typical high-velocity formation (v_s > v_f, where v_s and v_f are the velocities of the shear wave in the rock and of the compressional wave in the borehole fluid, respectively), the pseudo-Rayleigh waves, whose elastic properties depend mainly on the shear modulus of the rock, contributed significant energy to the total signal energy in the borehole. The energies of different wave packets depend on the permeability in different ways. The greatest sensitivity to permeability change has been shown by the acoustic signal total energy and the energy of the low-velocity part of the pseudo-Rayleigh wave packet. The theoretical analysis was illustrated by real sonic log data.     

Accuracy of kinematic wave and diffusion wave for spatial‐varying rainfall excess over a plane

The kinematic‐wave and diffusive‐wave approximations were investigated for unsteady overland flow resulting from spatially varying rainfall excess. Three types of boundary conditions were adopted: zero flow at the upstream end, and critical flow and zero depth‐gradient at the downstream end. Errors were derived by comparing the dimensionless profiles of the flow depth over the plane with those computed from the dynamic‐wave solution. It was found that the mean errors for both the approximations were independent of the type of rainfall excess distribution for KF₀² > 5, where K is the kinematic‐wave number and F₀ is the Froude number. Therefore, the regions (KF₀², F₀) where the kinematic‐wave and diffusive‐wave solutions would be fairly accurate and for any distribution of spatially varying rainfall, were characterized. The kinematic‐wave approximation was reasonably accurate, with a mean error of less than 5% and for the critical depth at the downstream end, for KF₀² ≥ 20 with F₀ ≤ 1; if the rainfall excess was concentrated in a portion of the plane, the field where the kinematic‐wave solution was found accurate, it was more limited and characterized for KF₀² > 35 with F₀ ≤ 1. The diffusive‐wave solution was in good agreement with the dynamic‐wave solution with a mean error of less than 5%, in the flow depth, for KF₀² ≥ 15 with F₀ ≤ 1; for rainfall excess concentrated in a portion of the plane, the accuracy of the diffusion wave solution was in a region more restricted and defined for KF₀² ≥ 30 with F₀ ≤ 1. For zero‐depth gradient at the downstream end, the accuracy field of the kinematic‐wave was found to be greater and characterized for KF₀² > 10 with F₀ ≤ 1; for rainfall excess concentrated in a portion of the plane, the region was smaller and defined for KF₀² > 15 with F₀ ≤ 1. The diffusive‐wave solution was found accurate in the region defined for KF₀² > 7·5, whereas for rainfall excess concentrated in a portion of the plane, the field of accuracy was for KF₀² > 12·5 with F₀ ≤ 1. The lower limits of the regions, defined on KF₀², can be considered generally valid for both approximations, but for F₀ < 1 smaller lower limits were also characterized. Finally, the accuracy of these approximations was influenced significantly by the downstream boundary condition. Copyright © 2002 John Wiley & Sons, Ltd.     

A probabilistic model of seismicity: Kamchatka earthquakes

The catalog of Kamchatka earthquakes is represented as a probability space of three objects {Ω, $ \tilde F $ \tilde F P}. Each earthquake is treated as an outcome ω _i in the space of elementary events Ω whose cardinality for the period under consideration is given by the number of events. In turn, ω _i is characterized by a system of random variables, viz., energy class k_i, latitude φ _i , longitude λ _i , and depth h _i . The time of an outcome has been eliminated from this system in this study. The random variables make up subsets in the set $ \tilde F $ \tilde F and are defined by multivariate distributions, either by the distribution function $ \tilde F $ \tilde F (φ, λ, h, k) or by the probability density f(φ, λ, h, k) based on the earthquake catalog in hand. The probabilities P are treated in the frequency interpretation. Taking the example of a recurrence relation (RR) written down in the form of a power law for probability density f(k), where the initial value of the distribution function f(k ₀) is the basic data [Bogdanov, 2006] rather than the seismic activity A ₀, we proceed to show that for different intervals of coordinates and time the distribution f _elim(k) of an earthquake catalog with the aftershocks eliminated is identical to the distribution f _full(k), which corresponds to the full catalog. It follows from our calculations that f ₀(k) takes on nearly identical numeral values for different initial values of energy class k ₀ (8 ≤ k ₀ ≤ 12) f(k ₀). The difference decreases with an increasing number of events. We put forward the hypothesis that the values of f(k ₀) tend to cluster around the value 2/3 as the number of events increases. The Kolmogorov test is used to test the hypothesis that statistical recurrence laws are consistent with the analytical form of the probabilistic RR based on a distribution function with the initial value f(k ₀) = 2/3. We discuss statistical distributions of earthquake hypocenters over depth and the epicenters over various areas for several periods     

Energy mode distribution: An analysis of the ratio of anti-Stokes to Stokes amplitudes generated by a pair of counterpropagating Langmuir waves

Results from plasma wave experiments in spacecrafts give support to nonlinear interactions involving Langmuir, electromagnetic, and ion-acoustic waves in association with type III solar radio bursts. Starting from a general form of Zakharov equation (Zakharov, V.E., 1985. Collapse and self-focusing of Langmuir waves. Hand-book of Plasma Physics Cap.2, 81–121) the equations for electric fields and density fluctuations (density gratings) induced by a pair of counterpropagating Langmuir waves are obtained. We consider the coupling of four triplets. Each two triplets have in common the Langmuir pump wave (forward or backward wave) and a pair of independent density gratings. We numerically solve the dispersion relation for the system, extending the work of (Alves, M.V., Chian, A.C.L., Moraes, M.A.E., Abalde, J.R., Rizzato, F.B., 2002. A theory of the fundamental plasma emission of type- III solar radio bursts. Astronomy and Astrophysics 390, 351–357). The ratio of anti-Stokes (AS) (ω₀+ω) to Stokes (S) (ω₀-ω^*) electromagnetic mode amplitudes is obtained as a function of the pump wave frequency, wave number, and energy. We notice that the simultaneous excitation of AS and S distinguishable modes, i.e., with Re{ω}=ω_r≠0, only occurs when the ratio between the pump wave amplitudes, r is ≠1 and the pump wave vector k₀ is , W₀ being the forward pump wave energy. We also observe that the S mode always receives more energy.     

Migration velocity analysis for tilted transversely isotropic media

Tilted transversely isotropic formations cause serious imaging distortions in active tectonic areas (e.g., fold‐and‐thrust belts) and in subsalt exploration. Here, we introduce a methodology for P‐wave prestack depth imaging in tilted transversely isotropic media that properly accounts for the tilt of the symmetry axis as well as for spatial velocity variations. For purposes of migration velocity analysis, the model is divided into blocks with constant values of the anisotropy parameters ε and δ and linearly varying symmetry‐direction velocity V_P0 controlled by the vertical (k_z) and lateral (k_x) gradients. Since determination of tilt from P‐wave data is generally unstable, the symmetry axis is kept orthogonal to the reflectors in all trial velocity models. It is also assumed that the velocity V_P0 is either known at the top of each block or remains continuous in the vertical direction. The velocity analysis algorithm estimates the velocity gradients k_z and k_x and the anisotropy parameters ε and δ in the layer‐stripping mode using a generalized version of the method introduced by Sarkar and Tsvankin for factorized transverse isotropy with a vertical symmetry axis. Synthetic tests for several models typical in exploration (a syncline, uptilted shale layers near a salt dome and a bending shale layer) confirm that if the symmetry‐axis direction is fixed and V_P0 is known, the parameters k_z, k_x, ε and δ can be resolved from reflection data. It should be emphasized that estimation of ε in tilted transversely isotropic media requires using nonhyperbolic moveout for long offsets reaching at least twice the reflector depth. We also demonstrate that application of processing algorithms designed for a vertical symmetry axis to data from tilted transversely isotropic media may lead to significant misfocusing of reflectors and errors in parameter estimation, even when the tilt is moderate (30°). The ability of our velocity analysis algorithm to separate the anisotropy parameters from the velocity gradients can be also used in lithology discrimination and geologic interpretation of seismic data in complex areas.     

Attenuation of High Frequency P and S Waves in the Gujarat Region,India

The local earthquake waveforms recorded on broadband seismograph network of Institute of Seismological Research in Gujarat, India have been analyzed to understand the attenuation of high frequency (2–25 Hz) P and S waves in the region. The frequency dependent relationships for quality factors for P (Q _P) and S (Q _S) waves have been obtained using the spectral ratio method for three regions namely, Kachchh, Saurashtra and Mainland Gujarat. The earthquakes recorded at nine stations of Kachchh, five stations of Saurashtra and one station in mainland Gujarat have been used for this analysis. The estimated relations for average Q _P and Q _S are: Q _P = (105 ± 2) f ^{0.82 ± 0.01}, Q _S = (74 ± 2) f ^{1.06 ± 0.01} for Kachchh region; Q _P = (148 ± 2) f ^{0.92 ± 0.01}, Q _S = (149 ± 14) f ^{1.43 ± 0.05} for Saurashtra region and Q _P = (163 ± 7) f ^{0.77 ± 0.03}, Q _S = (118 ± 34) f ^{0.65 ± 0.14} for mainland Gujarat region. The low Q (<200) and high exponent of f (>0.5) as obtained from present analysis indicate the predominant seismic activities in the region. The lowest Q values obtained for the Kachchh region implies that the area is relatively more attenuative and heterogeneous than other two regions. A comparison between Q _S estimated in this study and coda Q (Qc) previously reported by others for Kachchh region shows that Q _C > Q _S for the frequency range of interest showing the enrichment of coda waves and the importance of scattering attenuation to the attenuation of S waves in the Kachchh region infested with faults and fractures. The Q _S/Q _P ratio is found to be less than 1 for Kachchh and Mainland Gujarat regions and close to unity for Saurashtra region. This reflects the difference in the geological composition of rocks in the regions. The frequency dependent relations developed in this study could be used for the estimation of earthquake source parameters as well as for simulating the strong earthquake ground motions in the region.     

Distribution of magnetic energy in αΩ-dynamos,I: The method

Abstract

In this paper a method for solving the equation for the mean magnetic energy <BB> of a solar type dynamo with an axisymmetric convection zone geometry is developed and the main features of the method are described. This method is referred to as the finite magnetic energy method since it is based on the idea that the real magnetic field B of the dynamo remains finite only if <BB> remains finite. Ensemble averaging is used, which implies that fields of all spatial scales are included, small-scale as well as large-scale fields. The method yields an energy balance for the mean energy density ε ≡ B ²/8π of the dynamo, from which the relative energy production rates by the different dynamo processes can be inferred. An estimate for the r.m.s. field strength at the surface and at the base of the convection zone can be found by comparing the magnetic energy density and the outgoing flux at the surface with the observed values. We neglect resistive effects and present arguments indicating that this is a fair assumption for the solar convection zone. The model considerations and examples presented indicate that (1) the energy loss at the solar surface is almost instantaneous; (2) the convection in the convection zone takes place in the form of giant cells; (3) the r.m.s. field strength at the base of the solar convection zone is no more than a few hundred gauss; (4) the turbulent diffusion coefficient within the bulk of the convection zone is about 10¹⁴cm²s^?1, which is an order of magnitude larger than usually adopted in solar mean field models.     

Scattering and anelastic attenuation of seismic energy in northern Greece

The relative contribution of scattering (Q _s ^–1 ) and intrinsic (Q _i ^–1 ) attenuation to the totalS-wave attenuation for the frequencies of 1.5, 3.0, 6.0 and 12.0 Hz has been studied by applying the radiative energy transfer theory, Data of local earthquakes which occurred in northern Greece and were recorded by the permanent telementered network of the Geophysical Laboratory of the University of Thessaloniki have been used. The results show that in this area the scattering attenuation is dominant over all frequencies while intrinsic attenuation is significantly lower. The estimatedQ _s ^–1 andQ _i ^–1 values have frequency dependences off ^–0.72 andf ^–0.45, respectively. The frequency dependence ofQ _s ^–1 is the same as that of the codaQ _c ^–1 , obtained by applying the single scattering model, which probably implies that the frequency dependence of the coda wave attenuation is attributed to the frequency dependence of the scattering attenuation.Q _c ^–1 values are very close to scattering attenuation for short lapse times, (10–20 sec), and intermediate between scattering and intrinsic attenuation for the longer lapse times, (50–100 sec). This difference is explained as the result of the depth-dependent attenuation properties and the multiple scattering effects.     

Coda wave attenuation in the Parecis Basin,Amazon Craton,Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos Gaúchos, Amazon craton—Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q _c) of coda waves modeled using Q_c = Q₀ f^hQ_{\rm c} =Q_{\rm 0} f^\eta , where Q ₀ is the coda quality factor at frequency of 1 Hz and η is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239–252, 2008), based in the temporal attenuation coefficient, χ(f) instead of Q(f), given by the equation c(f)=g+\fracpfQ_e \chi (f)\!=\!\gamma \!+\!\frac{\pi f}{Q_{\rm e} }, for the calculation of the geometrical attenuation (γ) and effective attenuation (Q_e^-1 )(Q_{\rm e}^{-1} ). Q _c values have been computed at central frequencies (and band) of 1.5 (1–2), 3.0 (2–4), 6.0 (4–8), 9.0 (6–12), 12 (8–16), and 18 (12–24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin Q_c = (98±12)f^(1.14±0.08)Q_{\rm c} =(98\pm 12)f^{(1.14\pm 0.08)}, for the surrounding shield Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)}, and for the whole region of Porto dos Gaúchos Q_c = (99±19)f^(1.17±0.02)Q_{\rm c} =(99\pm 19)f^{(1.17\pm 0.02)}. Using the independent frequency model, we found: for the cratonic zone, γ = 0.014 s^− 1, Q_e^-1 = 0.0001Q_{\rm e}^{-1} =0.0001, ν ≈ 1.12; for the basin zone with sediments of ~500 m, γ = 0.031 s^− 1, Q_e^-1 = 0.0003Q_{\rm e}^{-1} =0.0003, ν ≈ 1.27; and for the Parecis basin with sediments of ~1,000 m, γ = 0.047 s^− 1, Q_e^-1 = 0.0005Q_{\rm e}^{-1} =0.0005, ν ≈ 1.42. Analysis of the attenuation factor (Q _c) for different values of the geometrical spreading parameter (ν) indicated that an increase of ν generally causes an increase in Q _c, both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, Q_c = (78±23)f^(1.17±0.14)Q_{\rm c} =(78\pm 23)f^{(1.17\pm 0.14)} (in the deepest part of the basin), than in the basement, Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)} (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.     

An inversion of site response and Lg attenuation using Lg waveform

Based on spectral ratio method, a joint inversion method was used to obtain parameters of Lg wave attenuation and site response. The inversion method allows simple and direct (two-parameter) determination of Lg wave attenua- tion and site response from sparse spectral data, which are not affected by radiation pattern factor and different response of same instrument after geometrical spreading. The method was used successfully for estimating site re- sponse of stations of Zhejiang Seismic Network and measuring Lg wave attenuation. The study is based on 20 earth- quakes occurred in northeast of Taiwan with magnitude MS5.0~6.7 and 960 seismic wave records from 16 stations in Zhejiang area from 2002 to 2005. The parameters of site response and Lg attenuation were calculated with a fre- quency interval of 0.2 Hz in the range of 0.5 Hz to 10 Hz. Lg wave attenuation coefficient corresponding to U-D, E-W and N-S components are γ ( f )=0.00175 f 0.43485, γ ( f )=0.00145f 0.48467 and γ ( f )=0.0021f 0.41241, respectively. It is found that the site response is component-independent. It is also found that the site response of QIY station is significant above the frequency of 1.5 Hz, and that the site response of NIB station is low for most frequency     

Sediment transport in high‐speed flows over a fixed bed: 2. Particle impacts and abrasion prediction

Single bed load particle impacts were experimentally investigated in supercritical open channel flow over a fixed planar bed of low relative roughness height simulating high‐gradient non‐alluvial mountain streams as well as hydraulic structures. Particle impact characteristics (impact velocity, impact angle, Stokes number, restitution and dynamic friction coefficients) were determined for a wide range of hydraulic parameters and particle properties. Particle impact velocity scaled with the particle velocity, and the vertical particle impact velocity increased with excess transport stage. Particle impact and rebound angles were low and decreased with transport stage. Analysis of the particle impacts with the bed revealed almost no viscous damping effects with high normal restitution coefficients exceeding unity. The normal and resultant Stokes numbers were high and above critical thresholds for viscous damping. These results are attributed to the coherent turbulent structures near the wall region, i.e. bursting motion with ejection and sweep events responsible for turbulence generation and particle transport. The tangential restitution coefficients were slightly below unity and the dynamic friction coefficients were lower than for alluvial bed data, revealing that only a small amount of horizontal energy was transferred to the bed. The abrasion prediction model formed by Sklar and Dietrich in 2004 was revised based on the new equations on vertical impact velocity and hop length covering various bed configurations. The abrasion coefficient k_v was found to be vary around k_v ~ 10⁵ for hard materials (tensile strength f_t > 1 MPa), one order of magnitude lower than the value assumed so far for Sklar and Dietrich's model. Copyright © 2017 John Wiley & Sons, Ltd.     

Nearshore circulation over transverse bar and rip morphology with oblique wave forcing

There is a paucity of field data to describe the transition in nearshore circulation between alongshore, meandering and rip current systems. A combination of in‐situ current meters and surf zone drifters are used to characterize the nearshore circulation over a transverse bar and rip morphology at Pensacola Beach, Florida in the presence of relatively low energy oblique waves. Current speeds vary in response to the relative wave height ratio (H_s/h), which defines the degree and extent of breaking over the shoal. In the absence of wave breaking the nearshore circulation was dominated by an alongshore current driven by the oblique waves. As waves begin to break across the shoal (0.2<H_s/ h<0.5) the nearshore circulation is characterized by a meandering alongshore current. As conditions became more dissipative (H_s/h>0.5), the meandering current is replaced by an unsteady rip circulation that moves offshore between the shoals before turning alongshore in the direction of wave advance outside the surf zone. The increase in wave dissipation is associated with an increase in very low frequency (VLF) variations in the current speed across the shoal and in the rip channel that caused the circulation to oscillate between an offshore and an alongshore flow. The unsteady nature of the nearshore circulation is responsible for 55% of all surf zone exits under these more dissipative conditions. In contrast, only 29% of the drifters released from the shoal exited the surf zone and bypassed the adjacent shoal with the alongshore‐meandering current. While the currents had a low velocity (maximum of ~0.4 m s^‐1) and would not pose a significant hazard to the average swimmer, the results of this study suggest that the transverse bar and rip morphology is sufficient to create an alongshore variation in wave dissipation that forces alongshore meandering and low‐energy rip circulation systems under oblique wave forcing. Copyright © 2013 John Wiley & Sons, Ltd.     

Reliability of rayleigh wave dispersion curve obtained from f–k spectral analysis of microtremor array measurement

The paper utilizes previously developed microtremor simulation technique to evaluate the reliability of Rayleigh wave dispersion curve estimated by f–k spectral analysis of microtremor array measurement. The simulated microtremors are obtained for a fictitious (TEST) site. Attempt is also made to obtain the dispersion curve for two real sites (OHDATE and SKC) by inverse analysis of the microtremor array measurement using f–k spectra method. The estimated dispersion curve from simulated microtremors (TEST site) compares well with the theoretical dispersion curve, demonstrating the reliability of f–k spectra method and indicating that the estimated dispersion curve from microtremor measurements could be adequately used as the target for inverse analysis purposes. It is also demonstrated that the dispersion curve from microtremor measurements can be utilized to estimate the soil profile at OHDATE and SKC sites by inverse analysis. Results show that the theoretical dispersion curve of the fundamental mode of Rayleigh wave after the end of inverse analysis are in good agreement with the dispersion curve obtained by f–k spectral analysis of microtremor array measurement.     

Heat transfer and thermal boundary layers in high amplitude annulus waves

Abstract

The heat transfer by a rotating, differentially-heated annulus of fluid is measured throughout the high amplitude wave regime. Only Δr_w T was varied (although v(T₁₅ ).K(T₁₅ ) varied by 46%), and it is found that Nu = C₁(λ)Ra^? away from the symmetry and low amplitude to wave transition curves and this is independent of ω. (λ is the wavelength.) On the wave side of these transition curves a region exists in which Nu (symmetry) λ Nu λ C₁(λ)Ra^?. The local heat transfer rate also varies strongly with wave phase.

Using a selection of measured internal thermal fields in the steady, high amplitude wave regime, the side-wall thermal boundary layer structure is examined. It is found that Nu, = C₂·Gr₂ ^A2; both C ₂ and A ₂ are independent of ω and λ to first order. For the time mean profiles, A ₂ ≈ 0.25; in the high heat transfer portion of the wave A ₂ < ¼ and in the low heat transfer portion of the wave A ₂ > ?. These relations hold over most of the vertical extent of the side walls. The deviations of the boundary layers from the above behavior which occur on the remainder of the walls is illustrated. The average thicknesses of the wall boundary layers ∞ Ra^?¼ except in that phase of the wave in which the wall to mid-gap temperature difference is the largest.     

Effect of consolidation ratios on maximum dynamic shear modulus of sands

The dynamic shear modulus (DSM) is the most basic soil parameter in earthquake or other dynamic loading conditions and can be obtained through testing in the field or in the laboratory. The effect of consolidation ratios on the maximum DSM for two types of sand is investigated by using resonant column tests. And, an increment formula to obtain the maximum DSM for cases of consolidation ratio κc>1 is presented. The results indicate that the maximum DSM rises rapidly when κc is near 1 and then slows down, which means that the power function of the consolidation ratio increment κc-1 can be used to describe the variation of the maximum DSM due to κc>1. The results also indicate that the increase in the maximum DSM due to κc>1 is significantly larger than that predicted by Hardin and Black's formula.     

On the propagation of topographically trapped waves

Abstract

In the context of ageostrophic theory in a homogeneous ocean, a nondimensional number is determined which corresponds to the Ursell number for long gravity waves. It is defined as Q = NL ²/h ³, where N is the amplitude of the wave travelling along the long length-scale direction, L is its length and h (which for gravity waves is the water depth) is given by h=(l ⁴ f ²/g)^1/3. where l is the short length-scale, f the Coriolis parameter and g the acceleration due to gravity. The physical meaning of Q is as follows: if Q? O(1) the free evolution of the wave is linear and weakly dispersive, if Q = O(1) nonlinear and dispersive effects balance out and finally if Q ?O(1) the evolution is nonlinear and non-dispersive. Expressions for the time scales for the development of dispersive and nonlinear effects are also determined. These results apply to topographically trapped waves, namely barotropic continental shelf and double Kelvin waves travelling along a rectilinear topographic variation.     

A note on edge waves on a flat shelf

Abstract

A variational approximation to the dispersion relation for trapped waves on a flat shelf of depth h ₁, bounded internally by a vertical coast and externally by a semi-infinite ocean of depth h ₂>h ₁, is obtained through an integral-equation formulation that accounts for all of the non-propagated modes that are excited at the discontinuity in depth (the conventional formulation of the edge-wave problem allows only for the propagated mode on the shelf and the dominant, non-propagated mode in the deep water). Coriolis effects are neglected. The exact result in the limit ω² h ₂/g↓0 (ω = angular frequency) is obtained by conformal mapping and compared with the variational approximation, which proves to be quite accurate over the entire range 1>h ₂/h ₁>x. The effects of the higher-order, non-propagated modes are found to be small for the long waves observed over the Southern California shelf by Snodgrass, Munk and Miller (1962).