Balloon measurement of vertical and horizontal atmospheric electric fields

Summary Twenty-eight balloons have been flown in Canada to measure vertical and horizontal electric fields at balloon altitudes. A horizontal electric field of magnitude typically between 10 and 50 millivolts/meter has been measured at night at high latitudes in association with auroras and magnetic disturbances. Its origin is in the magnetosphere and it maps to balloon alitudes with small attenuation according to Maxwell's equations. The vertical electric field at the balloon displays variations as the horizontal ionospheric field changes in order to maintain . Thus, magnetospheric processes affect both vertical and horizontal atmospheric electric fields and the potential differences induced by these processes may be comparable to weather induced potential differences. Weather processes have also been observed to produce large horizontal electric fields at balloon altitudes. Methods of distinguishing horizontal fields of ionospheric origin from those of weather origin are investigated with the conclusion that a determination of the source of a given event can often if not generally be made.This paper was presented byU. Fahleson     

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.     

Strong Potential Magnetic Field Influence on Slightly Differentially Rotating Spherical Shell

An electrically conducting viscous fluid-filled spherical shell is permeated by an axisymmetric strong potential magnetic field with large Elssaser number ² 1. We describe analytically the steady flow driven by a slightly faster rotation of the conducting inner boundary of the shell. The main flow is controlled by Ekman-Hartmann boundary layers with a small thickness /, where ² is the Ekman number. Asymptotics based on small ^–1 1 reveal the nature of a free shear layer O((/)^1/2) and a super-rotation that allows a part of the fluid to rotate faster than the inner sphere. The free shear is following an imposed field line that is tangent to the inner or outer sphere. Meridional flux is concentrated in the shear and boundary layers. Fluid tends to rotate with the inner sphere and to expel azimuthal magnetic field from an -region restricted by the free shear in the spherical shell. For an imposed axial uniform magnetic field, this -region is outside the cylinder tangent to the inner sphere and rotates with the outer sphere. Weak differential rotation O(/) is inside the cylinder, while almost all difference in rotation rates between spheres is accommodated in the thin O((/)^1/2) free shear. For an imposed dipole magnet, the region has a shape of a lobe touching the outer equator. Inside a super-rotation exists; this is the common case for such when the source of the imposed field is inside.     

Secular variation of the geomagnetic virtual dipole field during the last 2000 years as inferred from the spherical harmonic analysis of the available archeomagnetic data

Summary Recently excellent archeomagnetic data sequences have been bublished from several parts of the world. Using these sequences, an attempt is made to trace the secular variation of the virtual geomagnetic dipole field characterized by the three first order spherical harmonic coefficients . The archeomagnetic data (declination, inclination and total intensity) are transformed into the first order coefficients mentioned by a simple mathematical method. The secular variations of these coefficients, however, contain both dipole and non-dipole components. The separation of these is also attempted.Paper presented at the IAGA Symposium, Madrid, September 1969.     

Sull'approssimazione numerica di una variabile casuale normale,con applicaztone agli errori d'osservazione

Riassunto Data una variabile casuale X che segue la legge normale di probabilitl con valor medio a ed error medio y 1'A. considera un'altra variabile casuale che prende il valore intero r quando r–1/2

Summary Given a random variable X following the normal probability law, with expectation a and standard error p, the author considers another random variable , that takes the entire value r when r–1/2     

The disturbance due to a point source in a homogeneous liquid layer over a heterogeneous liquid half-space

Summary The problem of a periodic point source in a homogeneous liquid layer overlying a heterogeneous liquid half-space is discussed. After obtaining the formal solution, the path of integration for the displacement potential of the layer is transformed from the positive real axis to the positive imaginary axis and the Sommerfeld contour and the latter is further distorted to the modified Sommerfeld contour. The residues of the integrand at the poles contained within the Sommerfeld loop constitute the normal mode solution to the problem. The integrands in the expressions for the integrals along the imaginary axis are expanded in a series of negative powers of exponentials and then some of the terms in these expansions are evaluated approximately. This gives various waves reflected from the interface and the integral along the Sommerfeld loop vanishes. The frequency equation is obtained, also by the principle of constructive interference. An expression for the reflection coefficient at an interface of two liquid media, the upper medium being homogeneous and the lower one inhomogeneous, is obtained.     

Ueber Herde elastischer Wellen in isotropen,homogenen Medien

Zusammenfassung 1) Es werden Multipollösungen der skalaren Wellengleichung ² f/t ² – c² div gradf=0 betrachtet. Einerseits kann man solche Lösungen direkt durch Kugelfunktionenn-ter Ordnung ausdrücken, anderseits aus der Einpollösungf=1/p F(t–p/c) durch Differentiation nachn Richtungen erhalten. Es wird der Zusammenhang zwischen den Ergebnissen der beiden Verfahren gezeigt. — 2) Für die Energiedichte und den Energiefluss durch Kugelflächen bei kleinen elastischen Verschiebungen werden Ausdrücke in Kugelkoordinaten angegeben. — 3) Für die Wellengleichung grad div –b ² rot rot werden rotationsfreie Multipollösungen angegeben und Ausdrücke für Energiedichte und Energiefluss hergeleitet. — 4) Das gleiche wird für divergenzfreie Multipollösungen durchgeführt. — 5) Es werden Multipole betrachtet, die weder rotationsfrei noch divergenzfrei sind. Als Spezialfälle werden Multipole mit zeitlich begrenzter und solche mit periodischer Erregung gezeigt, ferner Lösungen der Wellengleichung, die sowohl rotationsfrei wie divergenzfrei sind. — 6) Es wird gezeigt, wie man die elastischen Wellen, die im Sinne vonStokes von einem Herdgebiet endlicher Ausdehnung ausgehen, näherungsweise durch elastische Multipole darstellen kann. — 7) Es wird angedeutet, wie man durch Messung von Komponenten von oder u.s.w. in Punkten im Innern des Mediums die Erregung und Energie von elastischen Multipolen bestimmen kann. Ferner wird auf den Fall hingewiesen, wo ein rotationsfreier Einpol sich im Innern eines Halbraumes befindet und die Messungen an seiner Oberfläche ausgeführt werden.

---

Summary (On foci of elastic waves in isotropic homogeneous media) — 1) Multiplets as solutions of the scalar wave equation ² f/t ² – c² div gradf=0 are considered. Such solutions can be obtained either directly by aid of spherical harmonics of ordern, or by differentiating the single polef=1/p F(t–p/c) with respect ton directions. The relations between the results of those two procedures are shown. — 2) In the case of small elastic displacements , the density of energy and the flow of energy through spherical surfaces are expressed by spherical coordinates. — 3) Multiplets which satisfy the equation of motion =a ² grad div b ² curl curl and the equation curl = 0 are given, and expressions for the density and flow of energy are found. — 4) The same is done with multiplets satisfying the equation of motion and the equation div = 0. — 5) General multiplets which satisfy the equation of motion are treated. As special cases, multiplets with excitation of finite length and multiplets with periodic excitation are considered, furthermore solutions of the equation of motion and of the equations curl = 0 and div = 0 are given. — 6) It is shown how elastic waves whose origin is a region of finite extension in the sense given byStokes, can be approximated by elastic multiplets. — 7) Some indications are given on the problem of how to find the functions of excitation and the energy of an elastic multiplet by measuring components of or etc., at points in the interior of the medium. The same problem is considered in the case of the single elastic pole. = grad 1/p F (t–p/a), if the measurements are made at the surface of an elastic half space.

     

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.     

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.     

Frequency-time analysis of geomagnetic beating-type pulsationsPc3

Summary The mechanism of beating of Pc3 type pulsations is studied. Using the method of numerical computation of a sonagram (the method of frequency-time analysis) a set of samples of pulsations from the Budkov Observatory is treated (1968–1969) mostly at K-indices equal to 2–3. By comparing f–t diagrams with the spectra of the samples an attempt has been made at interpreting the beating as a superposition of the frequency components, contained in the pulsation signal. In most observed cases it is possible to determine two close frequencies, the difference of which is on the average =5.4 mHz. The average carrier frequency of the samples was =37.6 mHz, and the average frequency of the beating =2.7 mHz. The interval of observed values of f_B amounted to 1–5 mHz. A tendency was observed for f_B to increase with increasing degree of disturbance of the geomagnetic field.     

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.     

Integral Equations for Kinematic Dynamo Models

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

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.     

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.     

Empirische Untersuchung zur Frage der Beziehung zwischen durchschnittlicher und kennzeichnender Wellenperiode im Seegang

Zusammenfassung Als Zusammenhang zwischen der kennzeichnenden Wellenperiode und der durchschnittlichen Periode im Seegang wird die Formel angesetzt. Mit Hilfe empirischer Unterlagen wird nachgewiesen, daßc eine Funktion des von D. E. Cartwright und M. S. Longuet-Higgins [1956] eingeführten Spektralparameters ist. Es wird eine vorläufige quantitative Beziehung zwischenc und abgeleitet.

---

Empirical investigations of the relation between the mean and the significant wave period in the sea

Summary It is supposed that the formula represents the relation between the significant wave period and the mean period in the sea. With the aid of empirical data it is demonstrated thatc is a function of the spectral parameter introduced by D. E. Cartwright and M. S. Longuet-Higgins [1956]. A preliminary quantitative relation betweenc and is derived.

---

Etudes empiriques de la relation entre la période moyenne et la période significative des vagues dans la houle

Résumé On suppose que la formule représente la relation entre la période significative des vagues et la période moyenne dans la houle. A l'aide des données empiriques on montre quec est une fonction du paramètre spectral , introduit par D. E. Cartwright et M. S. Longuet-Higgins [1956]. Une relation quantitative préliminaire entrec et est dérivée.

---

     

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

Estimate of earthquake hazard in the Vrancea (Romania) region

A maximum likelihood method is used to estimate the earthquake hazard parameters maximum magnitudeM _max, annual activity rate , and theb value of the Gutenberg-Richter equation in the Vrancea (Romania) region. The applied procedure permits the use of mixed catalogs with incomplete historical as well as complete instrumental parts, the consideration of variable detection thresholds, and the incorporation of earthquake magnitude uncertainty.Our imput data, comprises 105 historical earthquakes which occurred between 984 and 1934, and a complete data file containing 1067 earthquakes which occurred during the period 1935–30 August, 1986. The complete part was divided into four subcatalogs according to different thresholds of completeness. Only subcrustal events were considered, and dependent events were removed.The obtained value (=0.65) is at the lower range of the previously reported results, but it appears concurrent with conceptual and observational facts. The same concerns inferred value of _max = 7.8 and activity rate _4.0 = 5.34.     

The dynamics of magma mixing in a rising magma batch

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

Molecular refraction in the Earth's mantle

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

An analysis of heat-, momentum-transports,and spectra for clear air turbulence in mid-stratosphere

Summary An analysis of the clear air turbulence in the mid-stratosphere indicates that the turbulence is characterized by an anisotropic field of turbulence with an intense lateral component of the turbulence, associated with strong thermally stable stratification. The Richardson number in the region of turbulence is generally smaller than , the Taylor's criterion for a stratified shear flow.The cospectra for the momentum transport by the streamwise, lateral, and vertical components of the velocity show similarity in their distribution, decreasing with increasing wavenumber. The cospectra for the heat transport by the streamwise, lateral, and vertical components of the velocity show similarity in the high and medium wavenumber ranges, but not in the low wavenumber range of the spectra.The power spectra of the temperature and wind speed are very similar, and are approximately proportional to the power of the wavenumber. The power spectra of the streamwise and lateral components of the velocity are approximately proportional to the –2 power of the wavenumber, whereas the spectrum of the vertical component of the velocity is approximately proportional to the –1 power of the wavenumber.sponsored by the National Science Foundation.