Models of sub-grid variability in numerical simulations of solute transport in heterogeneous porous formations: three-dimensional flow and effect of pore-scale dispersion

Modeling transport of contaminants in the earths subsurface relies on numerical solutions over grids with blocks larger than Darcys scale. The hydraulic conductivity is homogenized over the grid blocks and the plumes spreading is reduced as a consequence of the wiped-out variability. To compensate for this loss Rubin et al. (1999) proposed to augment mixing by block-effective dispersion coefficients, and Rubin et al. (2003) showed, by means of two dimensional simulations, how this concept can be applied in practice. In this paper, we present new solutions of the block-effective dispersion tensor for an axisymmetric exponential covariance model. In addition, we discuss the influence of pore-scale dispersion in both two- and three-dimensional applications.     

Comparison of electric and elastic parameters of indian precambrian rocks

Summary Examination of the relation between electric resistivity, compressional wave velocity and density of Indian Precambrian rocks is presented.

---

¶rt;m auumu ¶rt; muu nmuu, m n¶rt; uu u nmm ¶rt;uu n¶rt; u ¶rt;uu.

     

Least squares inversion of one-dimensional magnetotelluric data: An assessment of procedures employed by Brown University

The pirmary objective of this paper is to review inversion algorithms employed to interpret magnetotelluric data by the geophysical/electromagnetics group at Brown University. Examples are presented from a variety of tectonic regimes — the Rio Grande rift, the Colorado Plateau, and Iceland. Following a brief tutorial in which the basic relations for the inverse prrolem are developed, we discuss the detectability and resolvability of various elements of the models. The correlation between certain model parameters is stressed, a phenomena well-known to interpreters. In this paper we emphasize the importance of identifying the precise range of values for which two or more parameters are coupled — these ranges are defined by what wer term saturation limits. Generalized inverse theory is employed to arrive at optimal models which are evaluated in terms of resolution and information matrices. The degree to which models fit a given dat base can be assessed using several quantitative approaches. We have found that it is impotant not only to determine the best fitting model, but also to identify those marginal solutions which also fall within the range of uncertainty of the data. The results from our direct invverse method agree favorably with those from Monte Carlo methods.     

Random maze models of flow through porous media

Summary This paper discusses a class of stochastic models of flow through porous media in which the randomness is attached to the structure of the medium rather than to the flow path. These models are obtained by generalizing an earlier model available in the literature where a regular crystal was taken in which bonds (flow channels) were dammed in a random fashion, yielding a random maze. The hydraulic properties of general models of this type are calculated; in particular, it is shown that they exhibit the phenomenon of dispersion whereby the factor of dispersion turns out to be a linear function of the percolation velocity.     

Electrical conductivity of the earth's crust and upper mantle

This review summarizes recent results of electrical resistivity studies of the earth's crust and upper mantle. Where available, the data are discussed in the context of further regional geophysical information. Electrical resistivity is very sensitive to a wide range of petrological and physical parameters, e.g., to carbon, fluids, volatiles and enhanced temperatures, making electrical resistivity methods a powerful tool in crust and upper mantle investigations. Yet, the general increase in resistivity data of the crust and mantle has not ended the battle of explanations for anomalous crustal conductivities.     

A Free Boundary Problem Related to the Location of Volcanic Gas Sources

A mathematical model of general gas emitting systems is derived, and a sample of relevant mathematical results is offered. The present paper indicates that shallow subsurface gas sources in typical volcanic areas can be located if appropriate physico-chemical measurements are made on the Earths surface and put to use.     

Far-field Gravity and Tilt Signals by Large Earthquakes: Real or Instrumental Effects?

A wide set of dynamics phenomena (i.e., Geodynamics, Post Glacial Rebound, seismicity and volcanic activity) can produce time gravity changes, which spectrum varies from short to long (more than 1 year) periods. The amplitude of the gravity variations is generally in the order of consequently their detection requires instruments with high sensitivity and stability: then, high quality experimental data. Spring and superconducting gravimeters are intensively used with this target and they are frequently jointed with tiltmeters recording stations in order to measure the elasto-gravitational perturbation of the Earth. The far-field effects produced by large earthquakes on records collected by spring gravimeters and tiltmeters are investigated here. Gravity and tilt records were analyzed on time windows spanning the occurrence of large worldwide earthquakes; the gravity records have been collected on two stations approximately 600 km distant. The background noise level at the stations was characterized, in each season, in order to detect a possible seasonal dependence and the presence of spectral components which could hide or mask other geophysical signals, such as, for instance, the highest mode of the Seismic Free Oscillation (SFO) of the Earth. Some spectral components (6.5; 8; 9; 14, 20, 51) have been detected in gravity and tilt records on the occasion of large earthquakes and the effect of the SFO has been hypothesized. A quite different spectral content of the EW and NS tiltmeter components has been detected and interpreted as a consequence of the radiation pattern of the disturbances due to the earthquakes. Through the analysis of the instrumental sensitivity, instrumental effects have been detected for gravity meters at very low frequency.F.S.E. (Fondo Sociale Europeo -European Community -)     

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

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

Large-scale structure of magnetospheric plasma

Recent investigations of magnetospheric plasma structure are summarized under the broad categories of empirical models, transport across boundaries, formation, and dynamics of the plasma sheet. This report reviews work in these areas during the period 1991 to 1993. Fully three-dimensional empirical models and simulations have become important contributors to our understanding of the magnetospheric system. Some new structural concepts have appeared in the literature: the entry boundary and geopause, the plasma sheet region 1 vortices, the low-energy layer, the adiabaticity boundary or wall region, and a region in the tail to which we refer as the injection port. Traditional structural concepts have also been the subject of recent study, notably the plasmapause, the magnetopause, and the plasma sheet. Significant progress has been made in understanding the nature of plasma sheet formation and dynamics, but the acceleration of electrons to high energy remains somewhat mysterious.     

Does the interplanetary magnetic field and the solar activity contribute to the tropospheric circulation?

A significant tendency is shown for both Etesians and (+, –) sector boundaries of the interplanetary magnetic field (IMF) to occur on the same solar rotation days, during the main period of the Etesians effect (July–August). In addition, the solar activity seems to control the Etesians distribution within the IMF sector structure. In the epoch of maximum there is a significant tendency of Etesians to occur during toward IMF days. In contrast, in the epoch of minimum Etesians occur mainly during away IMF days. Finally, in the epoch of intermediate the Etesians are uniformly distributed in away and toward IMF days. Since these conclusions are statistically significant at high confidence levels, it is fair to assume that IMF and solar activity seem to contribute, to some extent, to the Etesians occurrences, as well as to their distribution within the solar rotation and the IMF sector structure; that is, some solar contribution to the tropospheric circulation is implied.     

Review of Ray Theory Applications in Modelling and Imaging of Seismic Data

Throughout the last twenty years, 3D seismic ray modelling has developed from a research tool to a more operational tool that has gained growing interest in the petroleum industry. Various areas of application have been established and new ones are under development. Many of these applications require a modelling system with flexible, robust and efficient modelling algorithms in the core. The present paper reviews the basic elements of such a system, based on the open model concept and the wavefront construction technique. In the latter, ervený's dynamic ray tracing is an intrinsic part. The modelling system can be used for generating ray attributes and synthetic seismograms for realistic 3D surveys with tens of thousands of shots and receivers. Moreover, some other types of application areas are illustrated: Production of Green's functions for prestack depth migration and hybrid modelling (combined ray and finite-difference modelling), attribute mapping and illumination analysis, both for survey planning and interpretation. Finally, the concepts of isochron rays and velocity rays related to seismic isochrons have been introduced recently, with very interesting future applications.     

Numerical modeling of the high-latitude F-layer anomalies

The results of a detailed numerical study of the behaviour of the convective polar ionosphere are presented. The developed theoretical model produces three-dimensional distribution of electron density, electron and ion temperature. The effects of auroral particle precipitation on the density and temperature structures are studied for winter and low solar activity conditions. The high-latitude ionospheric features, such as the tongue of ionization, the main trough the polar ionization peak, the auroral ionization peak, the high-latitude ionization hole, the tongue of electron temperature, the high latitude minima of electron, temperature, and the ion temperature hot spot are obtained from calculations.Numerically obtained results are used for determining the HF propagation paths in the polar ionosphere. The effects of ionospheric irregularities on high frequency ionospheric radio waves are investigated by using a three-dimensional ray-tracing computer program. Ray-path trajectories are presented for different values of the elevation angle of transmission. From our study, it was found that large-scale irregularity structures of the high-latitude ionosphere, in the presence of the earth's magnetic field, significantly affect high frequency radio wave propagation.     

Solution of stochastic groundwater flow by infinite series,and convergence of the one-dimensional expansion

This paper investigates analytical solutions of stochastic Darcy flow in randomly heterogeneous porous media. We focus on infinite series solutions of the steady-state equations in the case of continuous porous media whose saturated log-conductivity (lnK) is a gaussian random field. The standard deviation of lnK is denoted . The solution method is based on a Taylor series expansion in terms of parameter , around the value =0, of the hydraulic head (H) and gradient (J). The head solution H is expressed, for any spatial dimension, as an infinite hierarchy of Green's function integrals, and the hydraulic gradient J is given by a linear first-order recursion involving a stochastic integral operator. The convergence of the -expansion solution is not guaranteed a priori. In one dimension, however, we prove convergence by solving explicitly the hierarchical sequence of equations to all orders. An infinite-order stochastic solution is obtained in the form of a -power series that converges for any finite value of . It is pointed out that other expansion methods based on K rather than lnK yield divergent series. The infinite-order solution depends on the integration method and the boundary conditions imposed on individual order equations. The most flexible and general method is that based on Laplacian Green's functions and boundary integrals. Imposing zero head conditions for all orders greater than one yields meaningful far-field gradient conditions. The whole approach can serve as a basis for treatment of higher-dimensional problems.     

Results of concentration and mobility measurements of positively and negatively charged particles taken by a rocket-borne parachuted aspiration (Gerdien) probe in the height region from 72 to 39 km

Summary With a self-neutralizing, rocket-borne parachute aspiration probe concentration and mobility measurements of positively and negatively charged carriers were made between 72 and 40 km over Sardinia. Between 72 km and 60 km a heavy and a light group of positively charged particles were found. The heavy group was more abundant than the hight group in the greater heights. For the negatively charged particles only one mobility group was found which roughly corresponded to the mobility of the light group of the positively charged carriers. The probe was designed so that no electrons could enter the analyzer section of the aspiration system. By this the electron concentration was given by the difference in concentration between the positively and negatively charged particles.This paper was read byH. Dolezalek, in an abbreviated form supplied by the authors.     

A Local Least-Squares Modification of Stokes’ Formula

The combination of Stokes formula and an Earth Gravity Model (EGM) for geoid determination has become a standard procedure. However, the way of modifying Stokes formula vary from author to author, and numerous methods of modification exist. Most methods are deterministic, with the primary goal of reducing the truncation bias committed by limiting the area of Stokes integration around the computation point, but there are also some stochastic methods with the explicit goal to reduce the global mean square error of the geoid height estimator stemming from the truncation bias as well as the random errors of the EGM and the gravity data. The latter estimators are thus, at least from a theoretical point of view, optimal in a global mean sense, but in a local sense they may be far from optimality.Here we take advantage of the error variance-covariance matrices of the EGM and the terrestrial gravity data to derive the modification parameters of Stokes kernel in a local least-squares sense. The solution is given for the unbiased type of modification of Stokes formula of Sjöberg (1991).     

The equatorial F-layer: progress and puzzles

This work reviews some aspects of the ionospheric F-layer in the vicinity of the geomagnetic equator. Starting with a historical introduction, brief summaries are given of the physics that makes the equatorial ionosphere so interesting, concentrating on the large-scale structure rather than the smaller-scale instability phenomena. Several individual topics are then discussed, including eclipse effects, the asymmetries of the equatorial trough, variations with longitude, the semiannual variation, the effects of the global thermospheric circulation, and finally the equatorial neutral thermosphere, including superrotation and possible topographic influences.     

A note on the general concept of wave breaking for Rossby and gravity waves

A recently proposed general definition of wave breaking is further discussed, in order to deal with some points on which misunderstanding appears to have arisen. As with surface and internal gravity waves, the classification of Rossby waves into breaking and not breaking is a generic classification based on dynamical considerations, and not a statement about any unique signature or automatically recognizable shape. Nor is it a statement about passive tracers uncorrelated with potential vorticity on isentropic surfaces. A strong motivation for the definition is that proofs of the nonacceleration theorem of wave, mean-flow interaction theory rely, explicitly or implicitly, on a hypothesis that the waves do not break in the sense envisaged.The general definition refers to the qualitative behaviour of a certain set of material contours, namely those, and only those, which would undulate reversibly, with small slopes, under the influence of the waves' restoring mechanism, in those circumstances for which linearized, nondissipative wave theory is a self-consistent approximation to nonlinear reality. The waves' restoring mechanism depends upon the basic-state vertical potential density gradient in the case of gravity waves, and upon the basic-state isentropic gradient of potential vorticity in the case of Rossby waves. In the usual linearized theory of planetary scale Rossby waves on a zonal shear flow, the relevant material contours lie along latitude circles when undisturbed.     

Tilt Observations in the Normal Mode Frequency Band at the Geodynamic Observatory Cueva de los Verdes,Lanzarote

We have conducted observations with the aid of a seismo-tiltmeter station, which is based on the Ostrovsky pendulum and installed at the Geodynamic Observatory Cueva de los Verdes at Lanzarote Island since 1995. In this station the signal is separated into two frequency bands – tidal tilts (from 0 to 5 mHz) and ground oscillations in the frequency range of free Earths normal modes (from 0.2 to 5 mHz). The later band, called accelerometer channel, has additional amplification. We analyzed the background records in the frequency range of Earths free oscillations from August 2000 to September 2001, as well as, Earths normal modes after strong earthquakes. We found several distinctive persistent peaks in the spectra of background oscillations. Both amplitudes of distinguished peaks and noises have seasonal variations. We found that spectra of background oscillations are different in the frequency interval between 1.4 and 2.5 mHz for North- South and East-West components.     

Generators in atmospheric electricity

Summary The definition of a generator in atmospheric electricity is considered, and various phenomena are discussed as to whether they can be described as generators.     

Geophysical Investigations Of Seismically Induced Surface Effects: Case Study Of A Landslide In The Suusamyr Valley, Kyrgyzstan

In summer 1998, a geophysical survey including seismic profilesand electrical tomography has been carried out in the Suusamyr valley, Kyrgyzstan. The scope wasto investigate surface effects induced by the Ms = 7.3 Suusamyr earthquake, the 19th of August, 1992. Inthis paper, special attention is paid to the case study of a debris slide triggered by the earthquake.Seismic data are analysed by P-wave refraction technique and by surface wave inversion. Electrical tomographicprofiles are processed by 2D-inversion.Using geotechnical and geological information, P-velocity modelsand resistivity sections are interpreted in terms of geological materials, in order to build a geological3D model. On the basis of the latter, we carried out static finite element computations as well as staticand pseudo-static calculations with Janbus method. Newmark displacement was computed, considering or notthe influence of the shallow soft deposits. The results are compared to the real displacementobserved in the field and conclusions are drawn about the mechanism of the landslide.