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.     

The mean energetic level. theory

Summary One of the important atmospheric levels, the mean energetic level (MEL), which in a sense reflects the energetics of the whole atmosphere, is defined. Its fundamental properties are shown. In order to describe the MEL correctly a new vertical coordinate is introduced and discussed. The new coordinate, , is defined as the ratio of height and temperature. The MEL is shown to be a level with constant value of . Some incorrect conclusions concerning the MEL, derived in the past, have been corrected.List of symbols used c _p specific heat of air at constant pressure - c _v specific heat of air at constant volume - e base of natural logarithms - E total potential energy - f Coriolis parameter - g acceleration of gravity - i specific internal energy - I internal energy - J enthalpy - k unit vector pointing upwards - p pressure - Q diabatic heating rate - R gas constant of the air - t time - T temperature - v horizontal velocity - v ⁽³⁾ three-dimensional velocity - w vertical velocity in thez-system - z height - temperature growth rate (T/z) - Pechala's vertical coordinate (z/T) - generalized vertical velocity in the -system (d/dt) - specific potential energy - potential energy - density of the air - Ruppert function - T(1–)^–1 - ( ) _S quantity at the sea level - ( )* quantity at the MEL     

Determination of the geopotential scale factor from satellite altimetry

Summary The geopotential scale factor R ₀ =GM/W ₀ has been determined on the basis of satellite altimetry as R _{0=(6 363 672·5±0·3) m} and/or the geopotential value on the geoid W ₀ =(62 636 256·5±3) m ² s ^–2 . It has been stated that R ₀ and/or W ₀ is independent of the tidal distortion of surface W=W ₀ due to the zero frequency tide.

---

¶rt;a nmu amumuu u ama amnmuaa R ₀ =GM/W ₀ =(6 363 672,5±0,3) m u/uu aunmuaa a nmuu¶rt;a W ₀ =(62 636 256,5±3) m² s^–2. m, m R ₀ u/uu W ₀ auum m nm amu a a nuu ¶rt;au nmu W=W ₀ .

     

The role of wavegroup velocity in the blocking problem

Considering the blocking problem as a baroclinic instability problem in a dispersive wave system with diabatic heating effects, it is of great interest to investigate the role of wavegroup velocityv _gr in blocking processes, becausev _gr controls the energy transfer in the wave field. Using a Newtonian Cooling —type of forcing with a phase differencek to the main field and taking the linearized version of a two-level model, the phase speedc _r, the group velocityv _gr and the growth ratekc _i have been obtained as analytical functions of the mean zonal windU, the thermal windU _T, the coefficient of diabatic heating x, the phase differencek and the wavelengthL. Now the hypothesis is introduced, that a blocking should be expected, ifv _gr has a maximum value in the vicinity ofL _o, for whichc _r vanishes and thee-folding timet=1/kc _i (kc _i>0) is smaller than 6 days (see condition (20) in the text). One finds, that for special parameter combinations (U _T, U, ), where 15 m/secU _T25m/sec,U=10m/sec, 0.8·10^–51.5·10^–5 [sec^–1], certain valuesL _o with an appropriate phase differencek exist, which satisfy the conditions mentioned above (for values see Table 2). TherebyL _o varies within the range 8500 km <L _o<11000 km corresponding to the preferred planetary blocking wavenumber 2 in middle latitudes 50°<<70° N.     

The gross features of the lithosphere-asthenosphere system in Europe from seismic surface waves and body waves

Long-period recordings of dispersive Rayleigh waves along numerous station lines, or profiles, in Europe have for the first time permitted a uniform inversion of these observations based on a new method of phase velocity regionalization.Regional dispersion relations obtained by this method have then been subjected to a complete inversion procedure commonly known as the hedgehog method. The results are presented in a map outlining the thickness of the lower lithosphere (lid) and the shear (S) velocities in both the lid and the asthenosphere channel.A comparison of these results with the minimum compressional (P) wave velocities in the asthenosphere and their corresponding depths provides an estimate of theV _p /V _s ratio for the asthenosphere in the European area.Contribution No. 314, Institute of Geophysics, ETH-Zürich, Switzerland.     

Flow and transport through two-dimensional isotropic media of binary conductivity distribution. Part 2: NUMERICAL simulations and comparison with theoretical results

In the present part the results of numerical simulations of flow and transport in media made up from circular inclusions of conductivity K that are submerged in a matrix of conductivity K ₀, subjected to uniform mean velocity, are presented. This is achieved for a few values of =K/K ₀ (0.01, 0.1 and 10) and of the volume fraction n (0.05, 0.1 and 0.2). The numerical simulations (NS) are compared with the analytical approximate models presented in Part 1: the composite elements (CEA), the effective medium (EMA), the dilute system (DSA) and the first-order in the logconductivity variance (FOA). The comparison is made for the longitudinal velocity variance and for the longitudinal macrodispersivity. This is carried out for n<0.2, for which the theoretical and simulation models represent the same structure of random and independent inclusions distribution. The main result is that transport is quite accurately modeled by the EMA and CEA for low , for which _L is large, whereas in the case of =10, the EMA matches the NS for n<0.1. The first-order approximation is quite far apart from the NS for the values of examined . This material is based upon work supported by the National Science Foundation under Grant No. 0218914. Authors also wish to thank the Center of Computational Research, University at Buffalo for assistance in running numerical simulations.     

Coulomb constitutive laws for friction: Contrasts in frictional behavior for distributed and localized shear

We describe slip-rate dependent friction laws based on the Coulomb failure criteria. Frictional rate dependence is attributed to a rate dependence of cohesionc and friction angle . We show that differences in the stress states developed during sliding result in different Coulomb friction laws for distributed shear within a thick gouge layer versus localized shear within a narrow shear band or between bare rock surfaces. For shear within gouge, shear strength is given by =c cos + _n sin, whereas for shear between bare rock surfaces the shear strength is =c cos + _n tan, where and _n are shear and normal stress, respectively. In the context of rate-dependent Coulomb friction laws, these differences mean that for a given material and rate dependence of the Coulomb parameters, pervasive shear may exhibit velocity strengthening frictional behavior while localized shear exhibits velocity weakening behavior. We derive from experimental data the slip-rate dependence and evolution ofc and for distributed and localized shear. The data show a positive rate dependence for distributed shear and a negative rate dependence for localized shear, indicating that the rate dependence ofc and are not the same for distributed and localized shear, even after accounting for differences in stress state. Our analysis is consistent with the well-known association of instability with shear localization in simulated fault gouge and the observation that bare rock surfaces exhibit predominantly velocity weakening frictional behavior whereas simulated fault gouge exhibits velocity strengthening followed by a transition to velocity weakening with increasing displacement. Natural faults also exhibit displacement dependent frictional behavior and thus the results may prove useful in understanding the seismic evolution of faulting.     

Hilbert transform in the interpretation of magnetic anomalies of various components due to a thin infinite dike

The Hilbert transformH(x) applicable to vertical (Z), horizontal (H), and total (T) magnetic anomalies due to a thin dike of infinite depth extent is derived from the generalised expression of magnetic effectF(x). The depth and dip of the dike is extracted by a simple procedure making use ofF(x) andH(x). A modified version of the amplitude of the analytic signal is given to locate the origin. The abscissa of the point of intersection ofF(x) and the discrete Hilbert transformH(1.x) directly yields the depth to the top. An example for each case is considered theoretically to illustrate the process. Applicability of the method is examined on the vertical component of the well-known magnetic anomaly at Kiirunavaara in northern Sweden, originally described by Von Carlheim Gyllenskjold, as well as on total magnetic anomaly of Bensons Mines, U.S.A.     

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      

The derivatives of the lunar disturbing potential harmonics with respect to euler's angles

Summary The derivatives of the harmonicsP _n ^(k) (sin _O)cos kT_O andP _n ^(k) (sin _O)sin kT_O, occurring in the development of the lunar disturbing potential, are derived upto n=4 and for k== 0, 1, ..., n. The equatorial co-ordinates _OT_O are referred to the Moon's mass centre; the procedure for the solar disturbing potential is formally identical.     

The effect of the velocity of the centre of a cyclone on the generation of microseisms

Summary The influence of the velocity of the movement of the centre of the cycloneV _c.c. on the rate of amplitudes' change A/t and periods' change T/t of storm microseisms is investigated. The dependence A/t=k V _c.c. and T/t=k ¹ V _c.c. is obtained. Unmovable depression (V _c.c. =0) does not stipulate the change of A/t and T/t.

---

u V _c.c. A/t T/t . A/t=k V _c.c. T/t=¹ V _c.c. . (V _c.c. =0) A/t T/t.

---

---

Presented as a scientific communication to the IASPEI Assembly in Madrid, Sept. 1969.     

Hydroxyl (6/2) airglow emission intensity ratios for rotational temperature determination

OH(6/2) Q1/P1 and R₁/P₁ airglow emission intensity ratios, for rotational states up to j = 4.5, are measured to be lower than implied by transition probabilities published by various authors including Mies, Langhoff et al. and Turnbull and Lowe. Experimentally determined relative values of j transitions yield OH(6/2) rotational temperatures 2 K lower than Langhoff et al., 7 K lower than Mies and 13 K lower than Turnbull and Lowe.     

Simple model for long-range forecasts of AT500 hPa level fluctuations

Summary A forecast model to predict the fluctuations of level AT 500 hPa in a selected grid of points is derived. The solution of the compensation equation is sought in the form of a trigonometric polynomial in three idependent variables. It constitutes the basis of a numerical solution of the prediction problem with the use of a high-speed computer. A three-month forecast of the altitude fluctuation of level 500 hPa is evaluated by means of the daily values of the correlation coefficient. The results are satisfactory and the general evaluation shows the model to be prospective.

---

¶rt;um nmua ¶rt; ¶rt; na au m AT 500a a uum m m. u au nauu um u¶rt; mu mu nua m m auu n. m u u nmu a¶rt;au a . u¶rt;um a m na au m 500a nu nu m au uuma uu. mam m anumu, u nu ¶rt; n¶rt;mam nnmu.

---

List of symbols used * coupling coefficient between two conjugate atmospheric elements - * ageostrophicity coefficient of the atmospheric system - Coriolis' parameter (=2 sin) - , * geographic latitude, geographic longitude - *, geopotential reference and pressure levels - , * compensation and coupling frequencies - integration field over the entire atmospheric system - A** constant (A*=2 ²(*+*)) - A _r ,B _r ,C _r ,D _r constants related to subsystem -r- - A _s ,B _s ,C _s ,D _s constants related to subsystem -s- - B** constant (B**= ⁴(*²+2**)) - C _r constant (i=1, 2, 3, 4) - E _k ,E _v ,E _p energies of the atmospheric system: kinetic, internal and potential - K transformation constant - m total number of generalized frequencies - R(T) frequency characteristics of the numerical band filter - r0(t),r1(t) daily values of the correlation coefficient - Q heat - x, y coordinates in the reference plane - t time - _p ² Laplace's differential operator in thep-system     

Water uptake and equilibrium sizes of aerosol particles at high relative humidities: Their dependence on the composition of the water-soluble material

Equilibrium water uptake and the sizes of atmospheric aerosol particles have for the first time been determined for high relative humidities, i.e., for humidities above 95 percent, as a function of the particles chemical composition. For that purpose a new treatment of the osmotic coefficient has been developed and experimentally confirmed. It is shown that the equilibrium water uptake and the equilibrium sizes of atmospheric aerosol particles at large relative humidities are significantly dependent on their chemical composition.List of symbols A proportionality factor - a _w activity of water in a solution - c _{p v} specific heat of water vapour at constant pressure - c _w specific heat of liquid water - f relative humidity - l _w specific heat of evaporation of water - M _i molar mass of solute speciesi - M _s mean molar mass of all the solute species in a solution - M _w molar mass of water - m ₀ mass of an aerosol particle in dry state - m _i mass of solute speciesi - m _s mass of solute - m _w mass of water taken up by an aerosol particle in equilibrium state - m total molality=number of mols of solute species in 1000 g of water - m _i molality of solute speciesi - m _k total molality of a pure electrolytek - O(m ²) remaining terms being of the second and of higher powers ofm - p ⁺ standard pressure - p total pressure of the gas phase - p pressure within a droplet - p ₁,p ₂,p ₃ coefficients in the expansion of M - p _1i, p_2i, p_3i specific parameters of ioni - p _s saturation vapour pressure - p _w water vapour pressure - R _w individual gas constant of water - r radius of a droplet - r ₀ equivalent volume radius of an aerosol particle in dry state - T temperature - T ₀ standard temperature - T ₁ temperature of the pure water drop in the osmometer - v _w specific volume of pure water - z _i valence of ioni - _i relativenumber concentration of ioni in a solution - correction term due to the adsorption of ions at liquid-solid interfaces - activity coefficient of solute speciesi in a solution, related to molalities - I bridge current - T temperature difference between solution and pure water drop in the osmometer - exponential mass increase coefficient - _w specific chemical potential of water vapour - _w specific chemical potential of water - ⁰ _w specific chemical potential of pure water vapour - ⁰ _w specific chemical potential of pure water - ₀ density of an aerosol particle in dry state - _w density of pure water - surface tension of a droplet - ₀ surface tension of pure water, i.e., at infinite dilution of the solute - osmotic coefficient - _k osmotic coefficient of a solution of a pure electrolytek - _k osmotic coefficient of a solution of a mixed solute - _M fugacity coefficient of water vapour - ^s _{i=1 i} z ² _i This work is part of a Ph.D. thesis carried out at the Meteorological Institute of the Johannes Gutenberg-Universität, Mainz.     

Deformation of a homogeneous gravitating sphere by internal dislocations

Summary An explicit solution is obtained for the system of equations describing the spheroidal motion in a homogeneous, isotropic, gravitating, elastic medium possessing spherical symmetry. This solution is used to derive the Green's dyad for a homogeneous gravitating sphere. The Green's dyad is then employed to obtain the displacement field induced by tangential and tensile dislocations of arbitrary orientation and depth within the sphere.Notation G Gravitational constant - a Radius of the earth - A _o =4/3 G - Perturbation of the gravitational potential - Circular frequency - V _p ,V _s Compressional and shear wave velocities - k _p =/V _p - k _s =/V _s - k _p [(2.8)] - , [(2.17)] - f _l ⁺ Spherical Bessel function of the first kind - f _l ^– Spherical Hankel function of the second kind - x =r - y =r - x _o =r _o - y _o =r_o - x =r k _s - y =r k _p - x _o =r _o k _s - y _o =r _o k _p - =a - =a - [(5.17)] - _{m, l}      

Accuracy in measurements and the temperature and volume dependence of thermoelastic parameters

To obtain the temperatureT and volumeV (or pressureP) dependence of the Anderson-Grüneisen parameter _T , measurements with high sensitivity are required. We show two examples:P, V, T measurements of NaCl done with the piston cylinder and elasticity measurements of MgO using a resonance method. In both cases, the sensitivity of the measurements leads to results that provide information about _T (,T), where V/V ₀ andV ₀ is the volume at zero pressure. We demonstrate that determination of _T leads to understanding of the volume and temperature dependence ofq=( ln / lnV) _T over a broadV, T range, where is the Grüneisen ratio.     

Quasi-static thermal deformations of a sphere by radiation from a point source

Summary In this paper the quasi-static temperature and stress distributions set up in an elastic sphere by radiation from a point source at a finite distance from the centre of the sphere and out-side it, have been discussed. The temperature boundary condition has been taken in the general form involving an arbitrary function of time. The final solutions have been obtained in terms of series involving Legendre polynomials. Numerical calculations have been done on IBM 1620 Computer and a desk calculator. The results have been represented in graphs.Notation the del operator - u the displacement vector - T the excess of temperature over that at state of zero stress and strain - , Lamé's constants - /2(+) Poisson's ratio - coefficient of linear expansion - 2(1+) - a radius of the sphere - d distance of the point source from the centre of the sphere - d _o a/d - K coefficient of thermal conductivity - h heat transfer coefficient of the surface     

A method of searching for radial velocity anomalies in the lower mantle using interfering signals

Radial velocity anomalies in the lower mantle that give rise to triplications in the travel-time curve for short-periodP waves will produce arrivals havingdT/d values that differe by roughly 0.2–0.5 s/deg. The first two arrivals associated with such triplications will be separated by less than one second over a distance range of 4°–10° they may not, therefore, be separable visually on single seismograms, so that their presence can only be inferred from some measurable property that depends on their mutual interference. If there are lateral variations in the regions of anomalous velocity gradients, the interfering signals will also have different azimuths of arrival. Using two synthetic wavelets we have investigated the effect of interference on bothdT/d and azimuth measurements at the Yellowknife Array. We found that if the interfering pulses have a dominant frequencyv, there is a range of time separations (0.30/v0.55/v) over which the measureddT/d and azimuth values may fluctuate by much more than the differences indT/d and azimuth between the interfering signals. We have evaluated the following empirically defined functions for three different primary signals, and for three different relative amplitudes of the interfering signals:f (t), the drift function, which expresses how the measured slownesses,p, and azimuths, , differ from the slownesses and azimuths of the primary wavelets; f(), the range function, which describes the behaviour of the upper and lower bounds ofp and as a function of the difference in arrival times of the signals, andf , studied the properties of these functions, and have outlined how these properties provide criteria based on the numerical and statistical characteristics of the arrival vectors, and on the waveform of the signal that will enable small radial velocity anomalies to be more clearly delineated.Contribution No. 863 from the Earth Physics Branch.     

Scattering attenuation and the fractal geometry of fracture systems

Scattering of seismic waves can be shown to have a frequency dependenceQ ^{–1 3–v} if scattering is produced by arrays of inhomogeneities with a 3D power spectrumW _3D(k) k ^–v. In the earth's crust and upper mantle the total attenuation is often dominated by scattering rather than intrinsic absorption, and is found to be frequency dependent according toQ ^–1 , where –1<–0.5. IfD ₁ is the fractal dimension of the surface of the 3D inhomogeneities measured on a 2D section, then this corresponds respectively to 1.5<D ₁1.75, since it can be shown that =2(D ₁–2). Laboratory results show that such a distribution of inhomogeneities, if due to microcracking, can be produced only at low stress intensities and slow crack velocities controlled by stress corrosion reactions. Thus it is likely that the earth's brittle crust is pervaded by tensile microcracks, at least partially filled by a chemically active fluid, and preferentially aligned parallel to the maximum principal compressive stress. The possibility of stress corrosion implies that microcracks may grow under conditions which are very sensitive to pre-existing heterogeneities in material constants, and hence it may be difficult in practice to separate the relative contribution of crack-induced heterogeneity from more permanent geological heterogeneities.By constrast, shear faults formed by dynamic rupture at critical stress intensities produceD ₁=1, consistent with a dynamic rupture criterion for a power law distribution of fault lengths with negative exponentD. The results presented here suggest empirically thatD ₁-1/2(D+1), thereby providing the basis for a possible framework to unify the interpretation of temporal variations in seismicb-value (b-D/2) and the frequency dependence of scattering attenuation ().This is PRIS contribution 046.     

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.