Recherches complémentaires sur le bases du calculateur météorologique «Temp»

Résumé La formule de base, traduisant une propriété analytique d'une classe très générale de fonctions, est un corollaire du théorème fondamental démontré dans un mémoire précédent, d'après lequel, étant donnés une fonction continue,p(, ,t) des points (, ) d'une surface régulière fermée et du temps et le champ d'un vecteur vitesse de transfert ou d'advection tangent à et ayant des lignes de flux fermées et régulières, il existe un opérateur spatial, linéaire, non singulierA tel que la fonctionA(p+Const.) soit purement advective par rapport a (sans creusement ni comblement). Ce théorème peut être exprimé par l'équation , où est un opérateur spatial, linéaire et non singulier, fonction deA.La détermination de peut être faite, soit en comparant deux formes différentes de la solution générale de l'équation en , soit en utilisant un raisonnement a priori très simple. On arrive ainsi au résultat pour un certain scalaireu(, ).Dans le cas oùp(, ,t) est la perturbation de la pression sur la surface du géoïde l'équation résulte aussi, comme nous l'avons montré dans le mémoire précédent, de notre théorie hydrodynamique des perturbations. On montre ici que la même équation peut encore être déduite de l'équation de continuité associée à la condition d'équilibre quasi statique selon la verticale.Comme applications de la formule de base (solution générale de l'équation enM), on étudie les problèmes suivants: 1^o creusement et comblement en général; 2^o creusement et comblement des centres et des cols; 3^o mouvement des centres et des cols; 4^o instabilité d'un champ moyen; 5^o propriétés spatiales des champsp(, ,t) et des vecteurs d'advection analytiques.Après une discussion des erreurs de la prévision d'un champp(, ,t) par la formule de base, du fait des erreurs des observations et du fonctionnement du calculateur, on examine quelques particularités du transfert ou advection d'un champf ₀(, ) par le vecteur . Enfin, le dernier chapitre du mémoire donne des éclaircissements complémentaires sur la structure du calculateur électronique «Temp» (qui effectue automatiquement les opérations mathématiques de la formule de base) et expose l'état actuel de sa construction.

---

Summary The basic formula, expressing an analytical property of a very general class of functions, is a corollary of the fundamental theorem, proved in a previous paper, according to which, given a functionp(, ,t) of the points (, ) of a closed regular surface and of the time, and a transfer or advection velocity vector tangent to and having regular closed streamlines, there is a spatial, linear, non singular operatorA such thatA(p+const.) is a purely advective function in respect to (no deepening). This theorem can be expressed by the equation where is a spatial, linear, non singular operator depending onA.The determination of can be attained, either by the comparison of two different forms of the general solution of the -equation, or by a simple a priori reasonning. The conclusion is thus reached that for a certain scalaru(, ).Whenp(, ,t) is the pressure perturbation at sea level, it was shown, in the preceding paper, that the equation can also be derived from our hydrodynamical perturbation theory. We now show that for this particular case, the same equation is also a consequence of the equation of continuity together with the condition of quasi statical vertical equilibrium.The following problems are then analysed by means of the basic formula: 1^o deepening and filling in general; 2^o deepening and filling of the centres and cols; 3^o motion of the centres and cols; 4^o instability of a mean field; 5^o spatial properties of the analytical fields and advection vectors .The errors in the forecast of a field,p(, ,t) by means of the basic formula, due to the observational and computational errors, are discussed, and some peculiarities of the transfer or advection of a fieldf ₀(, ) by are examined. Finally, complementary points are disclosed on the structure of the electronic computer «Temp» which performs automatically the mathematical operations of the basic formula, and a brief report is given of the present state of its construction.

     

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.     

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.     

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.

     

Relations entre le creusement et le comblement des perturbations et leur déplacement

Résumé On commence par définir le creusement et le comblement d'une fonctionp(, t) du tempst et des points (, ) d'une surface régulière fermée en se donnant, sur cette surface, un vecteur vitesse d'advection ou de transfert tangent à . Le creusement (ou le comblement) est la variation dep sur les particules fictives se déplaçant constamment et partout à la vitesse , A chaque vecteur et pour un mêmep(, ,t) correspond naturellement une fonction creusementC (, ,t) admissible a priori; mais une condition analytique très générale (l'intégrale du creusement sur toute la surface fermée du champ est nulle à chaque instant), à laquelle satisfont les fonctions de perturbation sur les surfaces géopotentielles, permet de restreindre beaucoup la généralité des vecteurs d'advection admissibles a priori et conduit à des vecteurs de la forme: , oùT est un scalaire régulier, () une fonction régulière de la latitude , le vecteur unitaire des verticales ascendantes etR/2 une constante. Ces vecteurs sont donc une généralisation naturelle des vitesses géostrophiques attachées à tout scalaire régulier. Dans le cas oùp(, ,t) est la perturbation de la pression sur la surface du géoïde, le vecteur d'advection par rapport auquel on doit définir le creusement est précisément une vitesse géostrophique: on a alors ()=sin etT un certain champ bien défini de température moyenne.On déduit ensuite une formule générale de géométrie et de cinématique différentielles reliant la vitesse de déplacement d'un centre ou d'un col d'un champp(, ,t) à son champ de creusementC (, ,t) et au vecteur d'advection correspondant. Cette formule peut être transformée et prend la forme d'une relation générale entre le creusement (ou le comblement) d'un centre ou d'un col et la vitesse de son déplacement, sans que le vecteur d'advection intervienne explicitement. On analyse alors les conséquences de ces formules dans les cas suivants: 1^o) perturbations circulaires dans le voisinage du centre; 2^o) perturbations ayant, dans le voisinage du centre, un axe de symétrie normal ou tangent à la vitesse du centre; 3^o) évolution normale des cyclones tropicaux.Finalement, on examine les relations qui existent entre le creusement ou le comblement d'un champ, le vecteur d'advection et la configuration des iso-lignes du champ dans le voisinage d'un centre.Ces considérations permettent d'expliquer plusieurs propriétés bien connues du comportement des perturbations dans différentes régions.

---

Summary The deepening and filling (development) of a functionp(, ,t) of the timet and the points (, ) of a regular closed surface is first of all defined, in respect to a given advection or transfer velocity field tangent to , as the variation ofp on any fictitious particle moving constantly and everywhere with the velocity . For a givenp(, ,t) and to any there corresponds a well defined development fieldC (, ,t). All theseC fields are a priori admissible, but a very general analytical condition of the perturbation fields in synoptic meteorology (the integral of the development fieldC (, ,t) on any geopotential surface vanishes at any moment), leads to an important restriction to advection vectors of the form: , whereT is any regular scalar, () any regular function of latitude, the unit vector of the ascending verticals andR/2 a constant. These vectors are a natural generalisation of the geostrophic velocities attached to any regular scalar. Whenp(, ,t) is the pressure perturbation at sea level, its development must be defined in respect to a geostrophic advection vector belonging to the above defined class of vectors with ()=sin andT a well defined mean temperature field.A general formula of the differential geometry and kinematics ofp(, ,t) is then derived, giving the velocity of any centre and col of ap(, ,t) as a function of the advection vector and the corresponding development fieldC (, ,t). This formula can be transformed and takes the form of a general relation between the deepening (and filling) of a centre (or a col) of ap(, ,t) and its displament velocity, the advection vector appearing no more explicitly. A detailed analysis of the consequences of these formulae is then given for the following cases: 1^o) circular perturbations in the vicinity of a centre; 2^o) perturbations having, in the vicinity of a centre, an axis of symmetry normal or tangent to the velocity of the centre; 3^o) normal evolution of the tropical cyclones.Finally, the relations between the developmentC (, ,t) of a fieldp(, ,t), the advection velocity vector and the configuration of the iso-lines in the vicinity of a centre are analysed.These theoretical results give a rational explanation of several well known properties of the behaviour of the perturbations in different geographical regions.

---

---

Communication à la 2ème Assemblée de la «Società Italiana di Geofisica e Meteorologia» (Gênes, 23–25 Avril 1954).     

Radiance and color of the sky at twilight: Perturbations caused by stratospheric haze

Layers of stratospheric aerosol with optical thicknesses as small as 10^–4 cause noticeable perturbations in the monochromatic logarithmic brightness gradient,G, and the color ratio,C, of the twilight sky. Modeling of the twilight's radiant properties shows that definite single-valued relationships exist between maxima inC or minima inG and optical thickness, , physical thickness h, and mean altitude, , of stratospheric layers. It is therefore possible to determine , h and and monitor their variations by performing either single wavelength measurements ofG or two-wavelength spectrophotometric measurements ofC. The presence of haze in the lower troposphere and the occurrence of multiple scattering both have relatively minor influences on the recovery of the stratospheric dust properties, provided that 10< <30 km.Formal mathematical inversions of the single-scattering twilight equations are possible in principle, but difficult in practice because of non-linearities. Inversions incorporating an iterative linearization process with constrained smoothing, successfully recovered the features of the haze layer, but tended to oversmooth the vertical profile and underestimate the mean altitude of the haze layer.     

Anisotropy of Cs137 gamma-ray absorption in rocks

Summary The anisotropy of the Cs¹³⁷ gamma-ray absorption coefficient has been investigated in layered models and several rock samples. It was found that in contrast to other petrophysical parameters the gamma-ray absorption coefficient is characterized by a relatively small coefficient of anisotropy (3–5%). Therefore, the possible effect of the anisotropy can be neglected in the measurements of for rocks as it is usually below current observational errors.     

Propagation of love waves in a non-homogeneous stratum of finite depth sandwitched between two semi-infinite isotropic media

Summary The aim of this paper is to study a problem in which the intermediate layer is non-homogeneous, the rigidity varying exponentially with depth i.e. ₂=₂ v ₀ ² e ^2pz , the density being constant, velocity varies also exponentially with depth according to the law =v ₀ e ^pz . The variability ofKH with the change of phase velocity is shown graphically.     

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     

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      

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

On some problems of seismic crustal phase

ＯｎｓｏｍｅｐｒｏｂｌｅｍｓｏｆｓｅｉｓｍｉｃｃｒｕｓｔａｌｐｈａｓｅＨｕａｎ－ＣｈｅｎｇＧＥ（葛焕称）（ＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒｅａｕｏｆＪｉａｎｇｓｕＰｒｏｖｉｎｃｅ,Ｎａｎｊｉｎｇ２１００１４,Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｉｎｔｈｉｓ...     

The sunspot areas and the relative sunspot numbers

Summary Within each sunspot cycle the yearly means (A) of the daily sunspot areas increase faster than the corresponding sunspot numbers (R) from the minimum to the maximum of solar activity and then decrease also faster than these numbers till the next minimum. Relation (A)=16.7 (R), frequently used so far, is approximately valid only for the years in the vicinity of the sunspot maximum. Instead of that, author gives the relations: for the years preceding the sunspot maximum, for the years following the sunspot maximum, wherea andb are constants,T _R is the time of rise of the corresponding sunspot cycle expressed in years, andk takes the valuek=0 for the year of maximum solar activity andk=1, 2, 3, ... for the first, second, third ... year preceding or following that of maximum solar activity. The monthly means of the daily sunspot areas show a similar variation, but in this case the ratioq=AR varies with a greater amplitude both within each sunspot cycle and from cycle to cycle. The values ofq corresponding to all months of a given year in the sunspot cycle are contained between two limits depending on the time of rise.

---

Résumé Les valeurs moyennes (A) des aires diurnes des taches solaires à chaque année depuis 1878 augmentent plus rapidement du minimum vers le maximum de l'activité solaire que les nombres de Wolf correspondants (R). Elles diminuent aussi plus rapidement que les nombres de Wolf du maximum vers le minimum de l'activité solaire. La relation adoptée (A)=16.7 (R) ne s'applique pas avec une approximation satisfaisante que seulement pour les années voisines celle du maximum de l'activité solaire. L'auteur propose les relations: pour les années qui précédent le maximum, pour les années qui suivent le maximum, oùa, b sont des constantes,T _R le temps d'ascension du cycle correspondant exprimé en années et la parametrek prend la valeurk=0 à l'année du maximum de l'activité solaire etk=1, 2, 3 ... pour les années qui précédent et qui suivent celle du maximum. Les valeurs moyennes des aires diurnes des taches à chaque mois, suivent la même marche mais dans ce cas le rapportq=AR present des larges variations. Il oscille pourtant extre deux limites qui dependent du temps d'ascension.

     

Propagation of Rayleigh waves in an axially symmetric heterogeneous layer lying between two homogeneous halfspaces

Summary The propagation of Rayleigh type waves in an axially symmetric inhomogeneous layer lying between two halfspaces is studied. The halfspaces are supposed to be identical in their elastic properties. The variation of the parameters in the layer is assumed to be of the form where is a constant andz is the distance measured from one interface into the layer. With this assumption, the vector wave equation for the layer is separable. The solution is obtained in terms ofWhittaker's functions and the frequency equation of Rayleigh type waves is derived.     

Time pattern of large earthquake sequences in Europe in 1976–1979

Summary The time development of three large sequences, Friuli 1976, Thessaloniki 1978 and Monte Negro-Albania 1979, was investigated by analysing the characteristics of the time distribution of magnitudes of the events, M(t), of average magnitudes, (t), and of the number of events, n(t), as well as the strain release curves and the magnitude-frequency relations, N(M). The existence of time gaps and of an anomalous decay of the numbers of shocks before major events, as well as the oscillation of ( const b^–1) may be considered as precursory phenomena. The clustering of events has been tested by the Poisson and negative binomial laws.Communication presented at the XVIIth General Assembly of the European Seismological Commission in Budapest, 21–29 August, 1980 under the title Properties of Seismic Regime and Recent Earthquake Sequences in Europe.     

A new dissipation model based on memory mechanism

Summary The model of dissipation based on memory introduced by Caputo is generalized and checked with experimental dissipation curves of various materials.List of symbols unidimensional stress - unidimensional strain - Q ^–1 specific dissipation function - c(t) creep compliance - m(t) relaxation modulus - c ₀ instantaneous compliance - m equilibrium modulus - (t) creep function - relaxation function - () spectral distribution of retardation times - spectral distribution of relaxation times - c ^*() complex compliance - m ^*() complex modulus - tang loss-tangent     

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.     

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.     

Zur Interpretation von seismischen Refraktionsmessungen

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

---

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