Some comparisons between mining-induced and laboratory earthquakes

Although laboratory stick-slip friction experiments have long been regarded as analogs to natural crustal earthquakes, the potential use of laboratory results for understanding the earthquake source mechanism has not been fully exploited because of essential difficulties in relating seismographic data to measurements made in the controlled laboratory environment. Mining-induced earthquakes, however, provide a means of calibrating the seismic data in terms of laboratory results because, in contrast to natural earthquakes, the causative forces as well as the hypocentral conditions are known. A comparison of stick-slip friction events in a large granite sample with mining-induced earthquakes in South Africa and Canada indicates both similarities and differences between the two phenomena. The physics of unstable fault slip appears to be largely the same for both types of events. For example, both laboratory and mining-induced earthquakes have very low seismic efficiencies where _a is the apparent stress and is the average stress acting on the fault plane to cause slip; nearly all of the energy released by faulting is consumed in overcoming friction. In more detail, the mining-induced earthquakes differ from the laboratory events in the behavior of as a function of seismic momentM ₀. Whereas for the laboratory events 0.06 independent ofM ₀, depends quite strongly onM ₀ for each set of induced earthquakes, with 0.06 serving, apparently, as an upper bound. It seems most likely that this observed scaling difference is due to variations in slip distribution over the fault plane. In the laboratory, a stick-slip event entails homogeneous slip over a fault of fixed area. For each set of induced earthquakes, the fault area appears to be approximately fixed but the slip is inhomogeneous due presumably to barriers (zones of no slip) distributed over the fault plane; at constant , larger events correspond to larger _a as a consequence of fewer barriers to slip. If the inequality _a / 0.06 has general validity, then measurements of _a =µE _a /M ₀, where is the modulus of rigidity andE _a is the seismically-radiated energy, can be used to infer the absolute level of deviatoric stress at the hypocenter.     

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

Three-dimensional inverse problem for inhomogeneous transversely isotropic media

Summary The basic formula used in the presented paper gives the relation between the P wave travel-time perturbation and the perturbation of an inhomogeneous transversely isotropic medium, expressed by four perturbations of elastic parameters and by two angles of orientation of the axis of symmetry of transverse isotropy in space. The travel time perturbation is computed along the ray in the unperturbed inhomogeneous isotropic medium. Four elastic parameters and two angles are parametrized in the model under study and a system of equations for many rays is constructed. The equations are linear in the sought elastic parameters and nonlinear in the sought angles, and the iterative Levenberg-Marquardt algorithm is thus used to solve them. The theoretical 3-D inverse problem was solved in the presented numerical example. The data, simulating teleseismic data, were computed in the direct problem and then inverted. The results indicate the applicability and limitation of the presented algorithm in real problems.

---

a a, unaa n¶rt;aa am, ¶rt;am mu ¶rt; uu u na u uu ¶rt;¶rt; nn umn ¶rt;, a m nuu naamau u ¶rt; au umauu u umuu nn umnuu nmam. u u na um ¶rt; a aa ¶rt;¶rt; umn ¶rt;. nu naam u ¶rt;a a naamuua ¶rt;u u nma uma au ¶rt; u . au u n um nu naama u u n um a umauu umuu, nm un m umamu aum a-aa¶rt;ma ¶rt; u u. am nu¶rt; ¶rt; m u nu. nu muu ¶rt;a aaa a na a¶rt;aa u am ¶rt;a a. mam naam auu u mu nuu nu¶rt;uma a a.

     

Changes in the shape of the circumpolar vortex due to heating from above

Summary Based on model considerations it is shown that, under certain assumptions, zonalization of tropospheric circulation may be expected in the region of the auroral oval as a result of heat released at the time energetic electrons penetrate from the Sun into the lower stratosphere.

---

a auu ¶rt; a¶rt;u naa, m nu m n¶rt;nu u¶rt;am auau mn uuu amu aa aa mam ¶rt;u mna nuu mumu m u a u mam.

     

An experimental study of triggered stick-slip

Stick-slips have been studied in the laboratory on granite, labrodorite and sandstone samples of two different sizes. Different roughness was achieved on the sawcut surfaces by finishing them with different grinding compounds ranging from grit 40 to grit 1000. Stick-slips occurred as a result of 1) slowly increasing the shear and normal stresses, 2) superimposing a sinusoidal stress modulation (0.1 and 10 Hz) on the slowly increasing stresses, 3) triggering by a stress impulse when the shear stress was well below the levels where stick-slips occurred without the impulse, and 4) foreshocks.Stick-slips triggered by impulses or foreshocks occurred long after the beginning of the triggering events, i.e., long in comparison with elastic wave travel times through the sample. All triggered events were very rich in high frequencies (corner frequency of 100 kHz). The untriggered stick-slips did not contain much energy at the high frequencies (corner frequency of 10 kHz). The dynamic friction coefficients for the triggered stick-slips were smaller than for the untriggered events.The long delay between the onset of the trigger and the stick-slip, and the high frequencies may be a consequence of corrosion of asperities. The ultimate triggering and the rate of corrosion are likely related to the interplay of the normal and shear stresses as they load and unload the fault surface. The consistent shape of the high frequency spectra is probably due to sample resonances which are excited rather than being characteristic of the details of the stick-slips. If these laboratory observations are directly applicable to earthquake seismology, the spectra of earthquakes which were triggered by other earthquakes should be anomalously rich in high frequencies.     

Optimum processing of gravimetric observations affected by the drift in a one-stage net,and the analysis of the accuracy of the results

Summary By using a general model of optimum elimination of systematic effects[5], the procedure has been developed for the optimum processing of gravimetric observations affected by the drift in a one-stage net, and statistically justified characteristics of the accuracy of the points of this net have been determined.

---

u uuu u ¶rt;u nmua uuauu um amuu uu [5] aamaa m¶rt;ua nmua amuaumuu uu ¶rt; aa ¶rt;-man mu; nu¶rt; mamumuu aamumuu mmu m¶rt; m m mu.

     

High heat flow measured in Ostrava-Karviná coal basin

aamuam mam uu mn nma ma-au a. a¶rt; uua¶rt;uma mnam, uuma mnn¶rt;mu u mn nma m ¶rt;uamau au uu uma, n¶rt; m n¶rt; mn nm nu¶rt; a. 2. a u au Q=1,99 u 2,06×10^–6 /² am ¶rt;au mu n ua amu, ma mu (2,5×10^–6 a/² ) aum nam mm ¶rt;uana. um mm am na m¶rt;, ¶rt;¶rt; nu¶rt;m am mam.     

Observed Response of the Ionospheric F2–Layer and Lower Thermosphere to Geomagnetic Storms during High Solar Activity

The response of the critical frequency of the ionosphere F2–layer, described by its main Fourier components (daily constant, diurnal and semidiurnal waves) and the lower thermosphere dynamics to the geomagnetic storms in July 1991 and February 1992 is studied. The daily constant displays a negative response, however, the magnitude of reaction depends on the season and latitude. The amplitudes of diurnal and semidiurnal waves increase during a geomagnetic storm, as this enhancement is very strong at high latitudes in winter. The prevailing neutral wind, especially the zonal wind, shows an inclination to decrease during the geomagnetic storm (the effect is more distinct in summer). The amplitudes of diurnal and semidiurnal tides also demonstrate a tendency toward reduction during high geomagnetic activity.     

Regional trends in European induction data

n ¶rt;a, n¶rt;mau 531 au ¶rt; u aum m u¶rt;uu n mumu ana¶rt;, ¶rt; u -m n, aauum ¶rt;um u u amuaa n¶rt;naa ma ua aama. uu nmam an¶rt;u ¶rt; u aum aam ¶rt;au cuu uP _n ^m , n¶rt;am mn n=1, 2, 3 u 5 (m n). u uua ¶rt;a¶rt;amu uu n¶rt;mauu uum au (a. 1) u u n aumam uu nmu, m n¶rt;mam u¶rt; am uuu ¶rt; u aum (u. 1–4). annuau 2 u 5 mn nm ma am mmmu m (u. 5, 6). ama uuu u m aam amu uu uma.     

A three-dimensional numerical model of advection and diffusion of urban source pollutants

ama mam maum n¶rt; n¶rt;um ma u ¶rt;, ¶rt;u ¶rt;unu nu nu au muu u. ¶rt;aa ¶rt; aa a uu au mm ¶rt;uuu n m¶rt;a am. u ¶rt;a n¶rt;aam m ¶rt; muu m nau am u u a um anu u m umu au, m ¶rt;aa ¶rt; a nm n¶rt;um nmam an¶rt;u mauu nu .     

Influence of the IMF sector boundaries on cosmic rays and tropospheric vorticity

Summary The effect of the IMF sector boundary crossing (IMF SBC) in the vorticity area index (VAI) — the well-known dip in the VAI after IMF SBC — is found to be independent of the IMF SBC effect in the cosmic ray flux. This finding refutes a recent suggestion by Lundstedt [1] that the IMF SBC effect in VAI is caused by a decrease in cosmic ray flux, but supports the concept of the IMF SBC effects in the ionosphere and atmosphere developed by Latovika [2–4]. Cosmic rays seem to affect the troposphere in another way.

---

¶rt;mu nu mau nam aum n ( ) a u¶rt; na¶rt;u aumu () — um uu n — a¶rt; auu m ma nm uu . mm mam nam ¶rt;a n¶rt;u ¶rt;m¶rt;a [1], m m a nuu nma uu , n¶rt;¶rt;uam nu m u u am, aum amu [2–4]. am m uu u m um a mn ¶rt;u a.

     

A steady-state model of a cumulus cloud

nua u mam n ma aa ¶rt;mn ¶rt;u aa, ma a n¶rt;a [1]. mma m ¶rt;u aa a ¶rt;auu nauMi (uuuu n) uMa (¶rt;uauu anm), m amm u nu [2–5]. u¶rt;m mam u nu ¶rt;u ¶rt; u¶rt; mu ¶rt; u a ¶rt; ¶rt;a. mu a aa, n¶rt;ua m a aa, nua a¶rt; u u _i uu a a¶rt; z_i nu nu aa ¶rt; ¶rt;a, m ¶rt;um mua m, a¶rt;u ¶rt; a, ¶rt;m, mnam, anu ¶rt; naa, nm a an u nu. a ¶rt;mn ¶rt;u m n a nu u ¶rt;mn ¶rt;u, ma a ma n¶rt;a [1].     

Micromechanics of the velocity and normal stress dependence of rock friction

Among the second-order effects on friction the most important are those of variable normal stress and of slip velocity. Velocity weakening, which is usually considered the source of the stick-slip instability in rock friction, has been observed in velocity stepping experiments with Westerly granite. The friction change, , was –0.01 to –0.008 for a tenfold velocity increase. Using normal closure measurements, we observed dilation upon each increase in sliding rate. We also observed, for the first time, time-dependent closure between surfaces during static loading. The dilation that occurred during the velocity stepping experiment was found to be that expected from the static time-dependent closure phenomenon. This change in closure was used to predict friction change with an elastic contact model. The calculated friction change which results from a change in contact area and asperity interlocking, is in good agreement with the observed velocity dependence of steady-state friction. Variable normal stress during sliding has two effects, first in creating new partial slip contacts and locking some existing fully sliding contacts and second in increasing interlocking, for instance when normal load is suddenly increased. As a result, a transient change in friction occurs upon a sudden change in normal load.     

The influence of a flood event on phytoplankton succession

The effect of the August 1987 River Reuss flood event on the phytoplankton biocoenosis in Lake Uri (Urnersee, part of Lake Lucerne, Switzerland) was investigated firstly by comparing biological, chemical and physical data sampled before the event with equivalent data sampled after the event; and secondly by comparing the phytoplankton succession in 1987 with that occurring in the floodfree year 1989. As a consequence of the flood, the physical and chemical environment of the phytoplankton was found to have undergone a change which resulted in an alteration in the composition of the phytoplankton community. The phytoplankton community existing previous to the flood event, which had been dominated byTabellaria fenestrata sensu Husted 1930 (K-strategist), was replaced by a biocoenosis characterized mainly by various species of flagellates, which represent a typical spring successional stage (r-strategists). After the externally-imposed perturbation, the return to stable physical and chemical conditions was followed by the re-establishment of the successional stage which had existed before the flood (termed reversion by Reynolds, 1980).     

Nonequilibrium thermodynamics of pressure solution

This paper is concerned with the thermodynamic theory of solution and precipitation processes in wet crustal rocks and with the mechanism of steady pressure-solution slip in contact zones, such as grain-to-grain contacts, fracture surfaces, and permeable gouge layers, that are infiltrated by a mobile aqueous solution phase. A local dissipation jump condition at the phase boundary is fundamental to identifying the thermodynamic force driving the solution and precipitation process and is used here in setting up linear phenomenological relations to model near-equilibrium phase transformation kinetics. The local thermodynamic equilibrium of a stressed pure solid in contact with its melt or solution phase is governed by Gibbs's relation, which is rederived here, in a manner emphasizing its independence of constitutive assumptions for the solid while neglecting surface tension and diffusion in the solid. Fluid-infiltrated contact zones, such as those formed by rough surfaces, cannot generally be in thermodynamic equilibrium, especially during an ongoing process of pressure-solution slip, and the existing equilibrium formulations are incorrect in overlooking dissipative processes tending to eliminate fluctuations in superficial free energies due to stress concentrations near asperities, defects, or impurities. Steady pressure-solution slip is likely to exhibit a nonlinear dependence of slip rate on shear stress and effective normal stress, due to a dependence of the contact-zone state on the latter. Given that this dependence is negligible within some range, linear relations for pressure-solution slip can be derived for the limiting cases of diffusion-controlled and interface-reaction-controlled rates. A criterion for rate control by one of these mechanisms is set by the magnitude of the dimensionless quantityk/2C _pD, wherek is the interfacial transfer coefficient, is the mean diffusion path length,C _p is the solubility at pressurep, andD is the mass diffusivity.     

Spectral theory of electromagnetic induction in a radially and laterally inhomogeneous earth

Summary The partial differential equations of electromagnetic induction in a 3-D Earth of inhomogeneous conductivity are reduced to a system of ordinary differential equations of the 2nd order for the spectral coefficients of the field.

---

au am nu¶rt; ¶rt; maum u¶rt;uu u m ¶rt;¶rt; n n¶rt;umu n¶rt; um ¶rt;uua au m n¶rt;a ¶rt; nma uum n.

     

Source radiation pattern andPKP amplitudes from the Tonga region

Summary The relation of short period PKP1 and PKP2 amplitudes from deep earthquakes in the Tonga region observed at a seismological station in Central Europe, is compared with the theoretical relations computed for different source radiation parameters and depths. The Preliminary Earth Model was used in the computations. Qualitatively good agreement of observed and calculated values was found for the radiation pattern corresponding to the likely geometry of the Tonga subduction zone. The influence of the variation of some source and model parameters on the calculated PKP1 and PKP2 relation is discussed.

---

¶rt;a au mu anum¶rt; mnu¶rt; PKP1 u PKP2 mu mu a, anua ¶rt; u mau ma n, mmmuu mmuuu uuau, aumau ¶rt; a naam ua umua uu aa. ama unaa ¶rt; PREM. a¶rt; u auma uu ¶rt;am am au ¶rt; a uu, mmm mmuu ¶rt;uu a a. ¶rt;am uu auau naam umu u ¶rt;u a uu am mu anum¶rt; PKP1 u PKP2.

     

Seismic velocity-ratio in the crust and uppermost mantle in southeast Wellington province,New Zealand

The regional variation of the seismic velocity-ratio () over a 200 km long traverse has been studied by means of microearthquake surveys. The Wadati-plot method is used with a minimum of four P and S arrivals for each of 49 earthquakes. The area as a whole is found to be characterized by a value of 1.74–1.76 for earthquakes of depth 12–40 km, except in a 50 km long section near Wellington, where is low at 1.60. This low has been attributed to the fault zones in the region. A small change of is observed between the upper crust (5 km) and lower crust (12 km), but there is no change of between the lower crust and uppermost mantle.     

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.     

Magnetic field at the core boundary in the nearly symmetric hydromagnetic dynamo z

Summary The behaviour of the poloidal and toroidal magnetic field at the core-mantle boundary is analysed in more detail, assuming that the conductive layer in the lowest mantle is thin. We can conclude that, in the case of the Z-model of the nearly symmetric hydromagnetic dynamo, the poloidal field may be considered potential everywhere in the mantle and that the azimuthal field depends on the geostrophic azimuthal velocity in the same manner as derived in[1] and[3].

---

aau ¶rt;-amu n¶rt;nuu m n¶rt; amuu aauum n¶rt;u nu¶rt;a u mu¶rt;a n. am ¶rt;, m Z-¶rt;u nmu umuu¶rt;aum ¶rt;ua aum nu¶rt;a n umam nmua ¶rt; amuu a n¶rt;u . ¶rt;m¶rt;am na [1] u [3] auum auma aum n m auma mu.