Joint conditional simulations and the spectral approach for flow modeling

The use of data to condition single random fields has a well-established history. However, the joint use of data from several cross-correlated random fields is not as well developed. For example, the use of both transmissivity and head data in a steady state 2-d stochastic flow problem is essentially an inverse problem that is very important for both flow and transport predictions. This problem is addressed here by using a combination of numerical simulation and analytical methods and its application illustrated. The type of information conveyed by the different data categories is explored. The results presented are especially interesting in that head and transmissivity each give different information: Head values would appear to constrain the geometry of the paths while transmissivity data yields information about travel times. The linearized model is expanded to an iterative procedure and the true conditional distribution at several locations is compared with the iterative solution.The problem mentioned above is one with a special transfer function specified by the flow equation. In the second part of the paper a Fast Fourier Transform method for generation and conditioning of two or more random fields is introduced. This procedure is simple to implement, fast and very flexible.     

Experimental investigation of turbulent processes in a complex terrain

In the present paper the first results of the international KOPEX-86 experiment are presented. The experiment took place at the Kopisty Atmospheric Observatory of the Institute of Physics of the Atmosphere in Prague as part of a special project of the Commission of the Academy of Sciences in Planetary Geophysics (KAPG) in June and July 1986. Using 4 ultrasonic anemometers at 4 levels up to 80 m, Doppler-SODAR and wind-, temperature- and radiation balance gradients from 2 to 80 m, a complex investigation of the atmospheric boundary layer in an industrial area was made. The authors present the first results of

Turbulent “Dynamo” in the Passive Scalar Problem

A turbulent magnetic dynamo can be considered as the evolution of a vector field in a turbulent fluid flow. The problem of evolution of scalar fields (e.g., number density of small particles) in a turbulent fluid flow is similar to the turbulent magnetic dynamo. The dynamo instability results in generation of magnetic field. The most important effect which can cause a generation of mean magnetic field in a turbulent fluid flow is the -effect: = – (1/3) u · ( × u), where u is the turbulent velocity field with the correlation time . A similar instability in the passive scalar problem results in formation of large-scale inhomogeneous structures in a spatial distribution of particles due to the -effect: = u_p ( · u_p), where u _p is the random velocity field of the particles which they acquire in a turbulent fluid velocity field. The effect is caused by inertia of particles which results in divergent velocity field of the particles. This results in additional turbulent nondiffusive flux of particles. The mean-field dynamics of inertial particles are studied by considering the stability of the equilibrium solution of the derived evolution equation for the mean number density of the particles in the limit of large Péclet numbers. The resulting equation is reduced to an eigenvalue problem for a Schrödinger equation with a variable mass, and a modified Rayleigh-Ritz variational method is used to estimate the lowest eigenvalue (corresponding to the growth rate of the instability). This estimate is in good agreement with obtained numerical solution of the Schrödinger equation. Similar effects arise during turbulent transport of gaseous admixtures (or light noninertial particles) in a low-Mach-number compressible fluid flow. The discussed effects are important in planetary and atmospheric physics (cloud formation, pollutant dynamics, preferential concentration of particles in protoplanetary disks and also planetesimals in them).     

Deconvolution of squared velocity waveform as applied to the study of a noncoherent short-period radiator in the earthquake source

We consider an inverse problem of determination of short-period (high-frequency) radiator in an extended earthquake source. This radiator is assumed to be noncoherent (i.e., random), it can be described by its power flux or brightness (which depends on time and location over the extended source). To decide about this radiator we try to use temporal intensity function (TIF) of a seismic waveform at a given receiver point. It is defined as (time-varying) mean elastic wave energy flux through unit area. We suggest estimating it empirically from the velocity seismogram by its squaring and smoothing. We refer to this function as observed TIF. We believe that one can represent TIF produced by an extended radiator and recorded at some receiver point in the earth as convolution of the two components: (1) ideal intensity function (ITIF) which would be recorded in the ideal nonscattering earth from the same radiator; and (2) intensity function which would be recorded in the real earth from unit point instant radiator (intensity Green's function, IGF). This representation enables us to attempt to estimate an ITIF of a large earthquake by inverse filtering or deconvolution of the observed TIF of this event, using the observed TIF of a small event (actually, fore-or aftershock) as the empirical IGF. Therefore, the effect of scattering is stripped off. Examples of the application of this procedure to real data are given. We also show that if one can determine far-field ITIF for enough rays, one can extract from them the information on space-time structure of the radiator (that is, of brightness function). We apply this theoretical approach to short-periodP-wave records of the 1978 Miyagi-oki earthquake (M=7.6). Spatial and temporal centroids of a short-period radiator are estimated.     

The k -2 Rupture Model Parametric Study: Example of the 1999 Athens Earthquake

We employ a realistic finite-extent k ^-2 stochastic source model with k-dependent rise time, introduced by Bernard et al. (1996) and extended by Gallovi and Brokeová (2004), to synthesize the strong ground motions in near source distances. For a given magnitude, this rupture model involves a number of free parameters (fault dimensions, slip roughness, maximum rise time, slip velocity function). We study their influence on the peak ground acceleration (PGA) maps (parametric study) for an example of the 1999 Athens earthquake of M _w = 5.9. The observed macroseismic intensities serve us as a rough constraint on the simulated PGA's to find a set of suitable k ^-2 models. These models agree with empirical source scaling relations.     

Semigroup solutions to stochastic unsteady groundwater flow subject to random parameters

Two methods for the solution of partial differential equations (PDE) for the general case of random in time physical parameters are presented and their application to the solution of unsteady regional groundwater flow equations are illustrated. The first method is the semigroup approach which directly offers a solution without resorting to closure approximations (hierarchy techniques), perturbation techniques, or Montecarlo simulation techniques. The semigroup approach can also handle the general stochastic problem when randomness also appears as initial conditions, boundary conditions or forcing terms. The second method is an approximation scheme to obtain the semigroup solution in complex cases and permits the solution of equations with more than one random coefficient.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

Terrestrial heat flow in Hungarian Permian strata

Summary Latest measurements of terrestrial heat flow in the Hungarian and Russian parts of the Carpathian basin confirm previously measured high flow values between 2.0–3.3 cal/cm² sec. Recent measurement in the Permian anticline structure of the Mecseck. Mts. in Hungary gave 2.4 cal/cm²sec, whereas in the Russian part of the basin, near to the Hungarian border 2.6 cal/cm² sec was measured in Miocene sediments. For more than 100000 km² surface of the Carpathian basin covered by Hungary and parts of Slovakia and Russian the high heat flow is an established fact.     

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

Scaling vesicle distributions and volcanic eruptions

Models of coalescence-decompressive expansion of the later stages of bubble growth predict that for diverse types of volcanic products the vesicle number densities (n(V)) are of the scaling form where V is the volume of the vesicles and B₃ the 3-dimensional scaling (power law) exponent. We analyze cross sections of 9 pumice samples showing that over the range of bubble sizes from 10 m to 3 cm, they are well fit with B₃0.85. We show that to within experimental error, this exponent is the same as that reported in the literature for basaltic lavas, and other volcanic products. The importance of the scaling of vesicle distributions is highlighted by the observation that they are particularly effective at packing bubbles allowing very high vesicularities to be reached before the critical percolation threshold, a process which—for highly stressed magmas—would trigger explosion. In this way the scaling of the bubble distributions allows them to be key actors in determining the rheological properties and in eruption dynamics.Editorial responsibility: D. Dingwell     

Theoretical background for the inversion of seismic waveforms including elasticity and attenuation

To account for elastic and attenuating effects in the elastic wave equation, the stress-strain relationship can be defined through a general, anisotropic, causal relaxation function ^ijkl (x, ). Then, the wave equation operator is not necessarily symmetric (self-adjoint), but the reciprocity property is still satisfied. The representation theorem contains a term proportional to the history of strain. The dual problem consists of solving the wave equation withfinal time conditions and an anti-causal relaxation function. The problem of interpretation of seismic waveforms can be set as the nonlinear inverse problem of estimating the matter density (x) and all the functions ^ijkl (x, ). This inverse problem can be solved using iterative gradient methods, each iteration consisting of the propagation of the actual source in the current medium, with causal attenuation, the propagation of the residuals—acting as if they were sources—backwards in time, with anti-causal attenuation, and the correlation of the two wavefields thus obtained.     

The four primary geodetic parameters

Summary The four primary geodetic parameters defining the geodetic reference system are discussed from the point of view of their physical meaning and current estimation of their actual accuracy. The geopotential scale factor has been treated as the primary geodetic parameter defining the Earth's dimensions.

---

¶rt;am m nu¶rt;uu naama, n¶rt;u¶rt;u um mumu, mu u uu a u mmu. ama amnmuaa ¶rt;am am nu¶rt;u naama, n¶rt; a u.

     

Radioactivity and geochemistry of granulites and some related rocks from the saxonian granulite complex

Summary The distribution of radioactive(Th, U, K), major and selected trace(Rb, Sr, Ba, Y, Zr, V, Cr, Ni) elements of granulites from the Saxonian Granulite Complex was studied. Similarly to the South Bohemian granulites, the Saxonian granulites can be divided according to the contents of their major and trace elements into two main groups, groupA containing mostly acid and subacid granulites (K ₂ O>2.5%, SiO ₂ >68%), and groupB containing mostly intermediate and basic granulites (K ₂ O<2.5%, SiO ₂ <68%). Statistically significant differences between groupsA andB were found for all major oxides and several trace elements(Rb, V, Cr, Ni). The Saxonian granulites follow the same calc-alkaline trend as the South Bohemian, granulitesA being placed mostly in the rhyolite field and granulitesB mostly in the dacite, andesite and basalt fields of this trend. The investigated granulites are characterized by a considerable scatter ofTh andU contents accompanied by very variableTh/U ratios; theTh andU concentrations of granulitesA are substantially lower than is usual for rocks of corresponding acidity.

---

¶rt;a an¶rt;u a¶rt;uamu(Th, U, K) u ua ¶rt;u(Rb, Sr, Ba, Y, Zr, V, Cr, Ni) m aum n¶rt;a aaum na. naa, m u¶rt;aum n uu aam n aaum u ¶rt;u am aua, u u uu. aum n u uu ma a¶rt;um ¶rt; ¶rt;nn; nnA nua¶rt;ama a au¶rt; u au¶rt;aum (K ₂ O>2,5%, Si O ₂ >68%), nnB ¶rt;u u aum (K ₂ O<2,5%, SiO ₂ <68%). ¶rt; muunnau mm mamumuu m au ¶rt; a u u ¶rt; m ¶rt;u m(Rb, V, Cr, Ni). auaum n¶rt;¶rt;m um- m¶rt; a u -uaum;aumA a¶rt;ma a uum n, uaumB a a ¶rt;aum, a¶rt;um u aam n m m¶rt;a. ¶rt;aum — u unnA — aamum uu ¶rt;au da¶rt;uamu mTh uU.

     

Transformation of phosphorus species in settling seston and during early sediment diagenesis

Sequential P extraction was combined with electron microscop and X-ray spectroscopy to characterise various P species and to study their transformation in settling seston and in recent sediment. During early diagenesis most of the particulate P formed in the water was redissolved. No net transformation into species that would resist dissolution was observed.It was shown that

Effective hydraulic conductivity of nonstationary aquifers

Assuming that the ln hydraulic conductivity in an aquifer is mathematically approximated by a spatial deterministic surface, or trend, plus a stationary random noise, we treat the problem of finding what the effective hydraulic conductivity of that aquifer is. This problem is tackled by spectral methods applied to a type of diffusion equation of groundwater flow, together with suitable coordinate transformations. Analytical (exact) solutions in terms of elementary functions are presented for one- and three-dimensional finite and infinite domains. Stability criteria are obtained for the solutions, in terms of a critical parameter, that turns out to involve the product of correlation scale and trend gradient. For the case of finite and symmetrical domains, additional provisions to insure the stability of numerical calculations of effective hydraulic conductivity are provided. Effective hydraulic conductivity is an important property, with potential applications in the calibrations of groundwater and transport numerical models.     

On the general theory of equivalent projection

Summary The paper deals with the general theory of equivalent projection. By specifying the results derived, special cases may be treated. Examples of methods are described which enable free functions to be found in a general solution. A detailed procedure of solving the problem with a numerical example is given, taking into account the optimization according to the global criterton, established for extreme angular distortion.

---

aamuam a mu auu nu; ¶rt;am m am au. nu m¶rt;, nu amu ¶rt; uu uu. umam nmuauau n aamumu au, ma ¶rt; aua uau , u nuam n¶rt; m¶rt; u u nu.

---

---

Dedicated to Professor Frantiek Fiala on the Centenary of his Birthday     

Propagation of surface waves on the surface of a non-homogeneous aeolotropic cylindrical shell

Summary The object of the present paper is to investigate the propagation of surface waves on a non-homogeneous aeolotropic cylindrical shell surrounded by vacuum. The elastic constantsc _ij (i, j=1,2...) and density of the material of the shell are assumed to be of the form and respectively, where _ij, ₀ are constants andk ₁,k ₂ are any integers.     

Single- and inter-station transfer functions for non-synchronous geomagnetic arrays-II. the results of the field data analysis

Summary The theory of methods of computing single- and inter-station transfer functions in both the spectral and time domains was developed in paper[1]. Both approaches are applied to the variation data recorded at field stations along two non-simultaneous profiles traversing the eastern margin of the Bohemian Massif, where a zone of anomalous induction seems to mark an important geological boundary of formations with different histories of development. The results of both analyses are found to coincide within reasonable bounds of 20–30% in the principal induction characteristics.

---

u m¶rt; ama ¶rt;-u -mau n¶rt;am u nma u am a ua am[1]. am nua m am a n¶rt;¶rt;a nu ¶rt; aaua ¶rt;aaum auau aumua a n mau ¶rt; u nu, nau m au aua,¶rt; aa a aa u¶rt;uu. a, u¶rt;u, mamau a¶rt;a ¶rt; ¶rt; ¶rt;uuau au umuu aumu. mam aau nma u am auam a 20–30% ¶rt; u¶rt;u naam.

     

On the transmissivity of the atmosphere in the water vapour absorption (6–6.5 μ) band

Summary The atmospheric transmissivity in the water vapour absorption band (6–6.5 ) is evaluated using the recent Tiros III data and based onGoody's statistical model.Work completed while on leave from the Department of Meteorology, University of California, Los Angeles.     

Intensity of vlf and elf phenomena in the outer ionosphere

Summary Spatial and temporal variability of natural ELF-VLF phenomena field intensity has been studied using data from Interkosmos satellites with apogees below 2000 km. The results show a distinct dependence of the diurnal variation and latitudinal intensity distribution on the magnetic activity. The L-value at which the maximum of plasmaspheric noise occurs, e.g., correlates better with the Dst-index than with other indicies.

---

mama u a auau anmu n mm - u a ua a uu ¶rt;a uu m an u 2000 . mam naam m auum m ¶rt;a u um an¶rt;u uu m aum amumu. uua L-naama, a m nm au naa uu, um Dst-u¶rt; ¶rt;uu u¶rt;auaum amumu.