On the onset of thermal convection in slowly rotating fluid shells

Abstract

The problem of the removal of the degeneracy of the patterns of convective motion in a spherically symmetric fluid shell by the effects of rotation is considered. It is shown that the axisymmetric solution is preferred in sufficiently thick shells where the minimum Rayleigh number corresponds to degree l = 1 of the spherical harmonics. In all cases with l > 1 the solution described by sectional spherical harmonics Y^l _l (θ,φ) is preferred.     

BAYESIAN ESTIMATION IN SEISMIC INVERSION. PART I: PRINCIPLES1

This paper gives a review of Bayesian parameter estimation. The Bayesian approach is fundamental and applicable to all kinds of inverse problems. Its basic formulation is probabilistic. Information from data is combined with a priori information on model parameters. The result is called the a posteriori probability density function and it is the solution to the inverse problem. In practice an estimate of the parameters is obtained by taking its maximum. Well-known estimation procedures like least-squares inversion or l₁ norm inversion result, depending on the type of noise and a priori information given. Due to the a priori information the maximum will be unique and the estimation procedures will be stable except (in theory) for the most pathological problems which are very unlikely to occur in practice. The approach of Tarantola and Valette can be derived within classical probability theory. The Bayesian approach allows a full resolution and uncertainty analysis which is discussed in Part II of the paper.     

Seismic behaviour of confined masonry walls

The results of tests of plain and confined masonry walls with h/l ratio equal to 1·5, made at 1:5 scale, have been used to develop a rational method for modelling the seismic behaviour of confined masonry walls. A trilinear model of lateral resistance–displacement envelope curve has been proposed, where the resistance is calculated as a combination of the shear resistance of the plain masonry wall panel and dowel effect of the tie-columns’ reinforcement. Lateral stiffness, however, is modelled as a function of the initial effective stiffness and damage, occurring to the panel at characteristic limit states. Good correlation between the predicted and experimental envelopes has been obtained in the particular case studied. The method has been also verified for the case of prototype confined masonry walls with h/l ratio equal to 1·0. Good correlation between the predicted and experimental values of lateral resistance indicates the general validity of the proposed method. © 1997 John Wiley & Sons, Ltd.     

Dynamic stiffness matrices for two circular foundations

A systematic procedure is presented for generating dynamic stiffness matrices for two independent circular foundations on an elastic half-space medium. With the technique reported in References 1–3, the analytic solution of three-dimensional (3D) wave equations satisfying the prescribed traction due to the vibration of one circular foundation can be found. Since there are two analytic solutions for two prescribed tractions due to the vibrations of two circular foundations, the principle of superposition must be used to obtain the total solution. The interaction stresses (prescribed tractions) are assumed to be piecewise linear in the r-directions of both cylindrical co-ordinates for the two circular foundations. Then, the variational principle and the reciprocal theorem are employed to generate the dynamic stiffness matrices for the two foundations. In the process of employing the variational principle, a co-ordinate transformation matrix between two cylindrical co-ordinate systems is introduced. Some numerical results of dynamic stiffness matrices for the interaction of two identical rigid circular foundations are presented in order to show the effectiveness and efficiency of the present method, and some elaborations for its future extensions are also discussed.     

Eddy correlations for water flow in a single fracture with abruptly changing aperture

Fluid flow in single fractures with non‐uniform apertures is an important research subject in many disciplines. The abruptly changing aperture is a special case of such non‐uniformity. This paper simulates water flow in a single fracture with abruptly changing aperture (SF‐ACA) using the Lattice Boltzmann Method (LBM) and the Finite Volume Method (FVM). The flow occurs with the Reynolds number (Re) ranging from 5 to 900 and a ratio of aperture change (E) of 3 (E = D/d, where D and d are the larger and smaller apertures, respectively). For Re values between 5 and 100, both LBM and FVM can successfully simulate the eddy development in the expansion regime of an SF‐ACA. Flow with high Re values (up to 900) is simulated by FVM, which appears to be numerically more stable than LBM for high‐Re flow problems studied here. The flow symmetry in the expansion regime breaks at the Re value between 400 and 500. Our simulation result shows a linear relationship between l₁/d and Re at low Re (5–100) or higher Re (110–900) values, where defined as the length from the location of abrupt expansion to the right edge of the first eddy along the flow direction. If considering the simulation results for the entire simulated range of Re (5–900), the l₁/d–Re relationship is better described by a non‐linear logarithmical function. The l₁/d approaches an asymptotic constant at large Re. Copyright © 2011 John Wiley & Sons, Ltd.     

Changes in flood regime by use of the modified curve number method

Abstract

This paper describes the use of a simple two stage rainfall-runoff model in which a curve number (CN) principle is used to calculate the soil water content and, subsequently, the rainfall contribution to direct runoff and groundwater flow. The maximum soil water retention, S, is used to express various characteristics of a catchment (infiltration rate, soil cover and land use, as in the CN method) relevant to flood formation. Using historical flood events, the model is calibrated, and the statistical distribution parameters of peak flows determined. With the same historical input data scenarios (rainfall), sets of flood hydrographs are simulated for various values of the parameter S, and corresponding distribution parameters of peak flows are determined. This procedure is used to demonstrate possible changes in flood regime to be expected due to changes of the catchment soil properties and its vegetation cover. A case study is presented for the River Hron catchment, area 582 km², in the mountainous region of central Slovakia.     

THE APPLICATION OF LINEAR FILTER THEORY TO THE DIRECT INTERPRETATION OF GEOELECTRICAL RESISTIVITY SOUNDING MEASUREMENTS*

Koefoed has given practical procedures of obtaining the layer parameters directly from the apparent resistivity sounding measurements by using the raised kernel function H(λ) as the intermediate step. However, it is felt that the first step of his method—namely the derivation of the H curve from the apparent resistivity curve—is relatively lengthy. In this paper a method is proposed of determining the resistivity transform T(λ), a function directly related to H(λ), from the resistivity field curve. It is shown that the apparent resistivity and the resistivity transform functions are linearily related to each other such that the principle of linear electric filter theory could be applied to obtain the latter from the former. Separate sets of filter coefficients have been worked out for the Schlumberger and the Wenner form of field procedures. The practical process of deriving the T curve simply amounts to running a weighted average of the sampled apparent resistivity field data with the pre-determined coefficients. The whole process could be graphically performed within an quarter of an hour with an accuracy of about 2%.     

Modeling Transient Streaming Potentials in Falling‐Head Permeameter Tests

We present transient streaming potential data collected during falling‐head permeameter tests performed on samples of two sands with different physical and chemical properties. The objective of the work is to estimate hydraulic conductivity (K) and the electrokinetic coupling coefficient (C_l) of the sand samples. A semi‐empirical model based on the falling‐head permeameter flow model and electrokinetic coupling is used to analyze the streaming potential data and to estimate K and C_l. The values of K estimated from head data are used to validate the streaming potential method. Estimates of K from streaming potential data closely match those obtained from the associated head data, with less than 10% deviation. The electrokinetic coupling coefficient was estimated from streaming potential vs. (1) time and (2) head data for both sands. The results indicate that, within limits of experimental error, the values of C_l estimated by the two methods are essentially the same. The results of this work demonstrate that a temporal record of the streaming potential response in falling‐head permeameter tests can be used to estimate both K and C_l. They further indicate the potential for using transient streaming potential data as a proxy for hydraulic head in hydrogeology applications.     

On the propagation of topographically trapped waves

Abstract

In the context of ageostrophic theory in a homogeneous ocean, a nondimensional number is determined which corresponds to the Ursell number for long gravity waves. It is defined as Q = NL ²/h ³, where N is the amplitude of the wave travelling along the long length-scale direction, L is its length and h (which for gravity waves is the water depth) is given by h=(l ⁴ f ²/g)^1/3. where l is the short length-scale, f the Coriolis parameter and g the acceleration due to gravity. The physical meaning of Q is as follows: if Q? O(1) the free evolution of the wave is linear and weakly dispersive, if Q = O(1) nonlinear and dispersive effects balance out and finally if Q ?O(1) the evolution is nonlinear and non-dispersive. Expressions for the time scales for the development of dispersive and nonlinear effects are also determined. These results apply to topographically trapped waves, namely barotropic continental shelf and double Kelvin waves travelling along a rectilinear topographic variation.     

Azimut di una geodetica ellissoidica passante per due punti molto lontani fra loro

Riassunto L'Autore presenta un procedimento per la ricerca dell'azimut di una geodetica ellissoidica passante per due punti molto lontani fra loro. Istituisce dapprima una corrispondenza fra la geodetica ed un arco di cerchio massimo, corrispondenza definita con l'imporre uguali coordinate agli estremi dei due archi e la latitudine sferica normale uguale a quella ellissoidica. Sviluppa poi la longitudine ellissoidicaw in funzione di quella sferical, e poichè l fra gli estremi dell'arco di cerchio deve essere uguale a w, ne deduce una equazione per ricavare l'azimut.

---

Summary The Author outlines a procedure to research the azimuth of an ellipsoidal geodesic, which passes through two points far distant from another. A corrispondence is first instituted between the geodesic and a bow of maximum circle; this correspondence may be defined by imposing the same coordinates to the extreme parts of the two bows and the spherical normal latitude equal to the ellipsoidal one. The ellipsoidal longitudew is next developed as a function of the spherical onel and, as l between the extreme parts of the circle bow must equal w, an equation is hereby deduced to find out the azimuth.

     

Power law or power flaw?

Since its introduction by Svensson in 1959, the power law curve y = ax^b (where x and y are horizontal and vertical direction, respectively) has been widely used in morphological analysis of glacial trough cross-profiles. The numerical constants a and b are obtained by a linear regression analysis of the logarithmic form of the power law curve (ln y = ln a + b ln x). The value b then gives a measure for the form of the cross-section. However, over the years this form of the power law has endured a lot of criticism. This criticism is well founded, since this particular form of the power law is not suitable for curve fitting in morphological analyses. In this paper a general power law is proposed, of the form y − y₀ = a|x − x₀|^b (where x₀, y₀ are the coordinates of the origin of the cross-profile). A unique and unbiased solution for this equation is obtained with a general least-squares method, thereby minimizing the error between the cross-profile data and the curve, and not between the logarithmic transform of the data and its regression line. This provides a robust way to characterize trough cross-section forms. © 1998 John Wiley & Sons, Ltd.     

A class of analytic solutions of the magnetic force-free field equations

Abstract

Formation of electric current sheets in the corona is thought to play an important role in solar flares, prominences and coronal heating. It is therefore of great interest to identify magnetic field geometries whose evolution leads to variations in B over small length-scales. This paper considers a uniform field B ₀[zcirc], line-tied to rigid plates z = ±l, which are then subject to in-plane displacements modeling the effect of photospheric motion. The force-free field equations are formulated in terms of field-line displacements, and when the imposed plate motion is a linear function of position, these reduce to a 4 × 4 system of nonlinear, second-order ordinary differential equations. Simple analytic solutions are derived for the cases of plate rotation and shear, which both tend to form singularities in certain parameter limits. In the case of plate shear there are two solution branches—a simple example of non-uniqueness.     

Canopy transpiration obtained from leaf transpiration,sap flow and FAO‐56 dual crop coefficient method

With a maize seed planting area of about 67 000 hm², Zhangye city supplies the seeds for more than 40% of the maize planting area in China. Irrigation water is often overused to ensure the quality of the maize seeds, leading to serious water shortage problems in recent years. An accurate and convenient estimate of canopy transpiration is of particular importance to ease the problem. In this paper, leaf transpiration and sap flow in a maize field were measured in 2012 using a portable photosynthesis system and a heat balance sap flow system. Based on a large amount of meteorological data and relevant maize plant‐growing parameters, canopy transpiration was up‐scaled from both leaf transpiration (T_l) and sap flow (T_f), and also calculated by the FAO‐56 dual crop coefficient method (T). Comparing these three types of transpiration, T_f was proved to be more reliable than T_l. Taking T_f as a benchmark, the basal crop coefficient (K_cb, the key parameter of FAO‐56 dual crop coefficient method) was further adjusted and verified for the maize plants in this region. In addition, the errors when using up‐scaling methods and FAO‐56 dual crop coefficient method are summarized. Copyright © 2014 John Wiley & Sons, Ltd.     

Maximum Entropy Density Determination for the Earth and Terrestrial Planets

The maximum entropy principle of the information theory gives rise to a general regularization strategy for ill-posed inverse problems. The methods based on this principle have become standard in various branches of engineering sciences. Of course, ill-posed problems frequently appear in Earth sciences, too. Nonetheless, the concept of maximum entropy is not very popular here. Therefore, we review the basic approaches employing the principle of maximum entropy in one way or another. We can distinguish at least three different approaches, partly yielding coincident results. One possible area of application is the determination of Earth and planetary models, although the paper cannot treat this in its practical complexity. Most of the discussion is restricted to the determination of the Earth's mass density function from various sources of data. Three sample problems are solved using the principle of maximum entropy: a spherical and an ellipsoidal problem related to the Earth and an ellipsoidal problem related to Mars. This illustrates the numerical procedure, which is non-trivial in many cases. It also shows some results, partly compared to standard solutions. The pros and cons of the approaches are discussed.     

Estimation of the elastic earth parameters using slr data for the low satellites STARLETTE and STELLA

In this paper we present estimated values for the global elastic parameters (k ₂, k ₃) and (h ₂, l ₂) derived from the analysis of Satellite Laser Ranging (SLR) data. We analyse SLR data for two low satellites, STARLETTE and STELLA, collected over a period of two years, from 1 January 2005 to 1 January 2007, from 18 globally distributed ground stations. We carry out a sequential analysis for the two satellites jointly, and study the stability of the estimates as a function of the length of the data set used. The adjusted final values of (k ₂, k ₃) and (h ₂, l ₂) for STARLETTE and STELLA are compared to, and are largely found to support, the estimates we previously published based on data for two high satellites LAGEOS 1 and LAGEOS 2. A major discrepancy between the two solutions was only found for the Shida number l ₂.     

TRAITEMENT AUTOMATIQUE DES SONDAGES ELECTRIQUES AUTOMATIC PROCESSING OF ELECTRICAL SOUNDINGS*

I. Introduction In this section the problem is stated, its physical and mathematical difficulties are indicated, and the way the authors try to overcome them are briefly outlined. Made up of a few measurements of limited accuracy, an electrical sounding does not define a unique solution for the variation of the earth resistivities, even in the case of an isotropic horizontal layering. Interpretation (i.e. the determination of the true resistivities and thicknesses of the ground-layers) requires, therefore, additional information drawn from various more or less reliable geological or other geophysical sources. The introduction of such information into an automatic processing is rather difficult; hence the authors developped a two-stage procedure:

• a) the field measurements are automatically processed, without loss of information, into more easily usable data;
• b) some additional information is then introduced, permitting the determination of several geologically conceivable solutions.

The final interpretation remains with the geophysicist who has to adjust the results of the processing to all the specific conditions of his actual problem. II. Principles of the procedure In this section the fundamental idea of the procedure is given as well as an outline of its successive stages. Since the early thirties, geophysicists have been working on direct methods of interpreting E.S. related to a tabular ground (sequence of parallel, homogeneous, isotropic layers of thicknesses h_i and resistivities ρ_i). They generally started by calculating the Stefanesco (or a similar) kernel function, from the integral equation of the apparent resistivity: where r is the distance between the current source and the observation point, S₀ the Stefanesco function, ρ(z) the resistivity as a function of the depth z, J₁ the Bessel function of order 1 and λ the integration variable. Thicknesses and resistivities had then to be deduced from S₀ step by step. Unfortunately, it is difficult to perform automatically this type of procedure due to the rapid accumulation of the errors which originate in the experimental data that may lead to physically impossible results (e.g. negative thicknesses or resistivities) (II. 1). The authors start from a different integral representation of the apparent resistivity: where K₁ is the modified Bessel function of order I. Using dimensionless variables t = r/2h₀ and y(t)=ζ (r)/ρ₁ and subdividing the earth into layers of equal thicknesses h₀ (highest common factor of the thicknesses h_i), ø becomes an even periodic function (period 2π) and the integral takes the form: The advantage of this representation is due to the fact that its kernel ø (function of the resistivities of the layers), if positive or null, always yields a sequence of positive resistivities for all values of θ and thus a solution which is surely convenient physically, if not geologically (II.3). Besides, it can be proved that ø(θ) is the Fourier transform of the sequence of the electric images of the current source in the successive interfaces (II.4). Thus, the main steps of the procedure are: a) determination of a non-negative periodic, even function ø(θ) which satisfies in the best way the integral equation of apparent resistivity for the points where measurements were made; b) a Fourier transform gives the electric images from which, c) the resistivities are obtained. This sequence of resistivities is called the “comprehensive solution”; it includes all the information contained in the original E.S. diagram, even if its too great detail has no practical significance. Simplification of the comprehensive solution leads to geologically conceivable distributions (h, ρ) called “particular solutions”. The smoothing is carried out through the Dar-Zarrouk curve (Maillet 1947) which shows the variations of parameters (transverse resistance R_i= h_i.ρ_i–as function of the longitudinal conductance C_i=h_i/ρ_i) well suited to reflect the laws of electrical prospecting (principles of equivalence and suppression). Comprehensive and particular solutions help the geophysicist in making the final interpretation (II.5). III. Computing methods In this section the mathematical operations involved in processing the data are outlined. The function ø(θ) is given by an integral equation; but taking into account the small number and the limited accuracy of the measurements, the determination of ø(θ) is performed by minimising the mean square of the weighted relative differences between the measured and the calculated apparent resistivities: minimum with inequalities as constraints: where t_l are the values of t for the sequence of measured resistivities and p_l are the weights chosen according to their estimated accuracy. When the integral in the above expression is conveniently replaced by a finite sum, the problem of minimization becomes one known as quadratic programming. Moreover, the geophysicist may, if it is considered to be necessary, impose that the automatic solution keep close to a given distribution (h, ρ) (resulting for instance from a preliminary interpretation). If φ(θ) is the ø-function corresponding to the fixed distribution, the quantity to minimize takes the form: where: The images are then calculated by Fourier transformation (III.2) and the resistivities are derived from the images through an algorithm almost identical to a procedure used in seismic prospecting (determination of the transmission coefficients) (III.3). As for the presentation of the results, resorting to the Dar-Zarrouk curve permits: a) to get a diagram somewhat similar to the E.S. curve (bilogarithmic scales coordinates: cumulative R and C) that is an already “smoothed” diagram where deeper layers show up less than superficial ones and b) to simplify the comprehensive solution. In fact, in arithmetic scales (R versus C) the Dar-Zarrouk curve consists of a many-sided polygonal contour which múst be replaced by an “equivalent” contour having a smaller number of sides. Though manually possible, this operation is automatically performed and additional constraints (e.g. geological information concerning thicknesses and resistivities) can be introduced at this stage. At present, the constraint used is the number of layers (III.4). Each solution (comprehensive and particular) is checked against the original data by calculating the E.S. diagrams corresponding to the distributions (thickness, resistivity) proposed. If the discrepancies are too large, the process is resumed (III.5). IV. Examples Several examples illustrate the procedure (IV). The first ones concern calculated E.S. diagrams, i.e. curves devoid of experimental errors and corresponding to a known distribution of resistivities and thicknesses (IV. 1). Example I shows how an E.S. curve is sampled. Several distributions (thickness, resistivity) were found: one is similar to, others differ from, the original one, although all E.S. diagrams are alike and characteristic parameters (transverse resistance of resistive layers and longitudinal conductance of conductive layers) are well determined. Additional informations must be introduced by the interpreter to remove the indeterminacy (IV.1.1). Examples 2 and 3 illustrate the principles of equivalence and suppression and give an idea of the sensitivity of the process, which seems accurate enough to make a correct distinction between calculated E.S. whose difference is less than what might be considered as significant in field curves (IV. 1.2 and IV. 1.3). The following example (number 4) concerns a multy-layer case which cannot be correctly approximated by a much smaller number of layers. It indicates that the result of the processing reflects correctly the trend of the changes in resistivity with depth but that, without additional information, several equally satisfactory solutions can be obtained (IV. 1.4). A second series of examples illustrates how the process behaves in presence of different kinds of errors on the original data (IV.2). A few anomalous points inserted into a series of accurate values of resistivities cause no problem, since the automatic processing practically replaces the wrong values (example 5) by what they should be had the E.S. diagram not been wilfully disturbed (IV.2.1). However, the procedure becomes less able to make a correct distinction, as the number of erroneous points increases. Weights must then be introduced, in order to determine the tolerance acceptable at each point as a function of its supposed accuracy. Example 6 shows how the weighting system used works (IV.2.2). The foregoing examples concern E.S. which include anomalous points that might have been caused by erroneous measurements. Geological effects (dipping layers for instance) while continuing to give smooth curves might introduce anomalous curvatures in an E.S. Example 7 indicates that in such a case the automatic processing gives distributions (thicknesses, resistivities) whose E.S. diagrams differ from the original curve only where curvatures exceed the limit corresponding to a horizontal stratification (IV.2.3). Numerous field diagrams have been processed (IV. 3). A first case (example 8) illustrates the various stages of the operation, chiefly the sampling of the E.S. (choice of the left cross, the weights and the resistivity of the substratum) and the selection of a solution, adapted from the automatic results (IV.3.1). The following examples (Nrs 9 and 10) show that electrical prospecting for deep seated layers can be usefully guided by the automatic processing of the E.S., even when difficult field conditions give original curves of low accuracy. A bore-hole proved the automatic solution proposed for E.S. no 10, slightly modified by the interpreter, to be correct.     

A NEW PARAMETER FOR THE INTERPRETATION OF INDUCED POLARIZATION FIELD PROSPECTING (time-domain) *

In a previous paper it has been shown that we can relate the transient IP electric field E_p , existing in a rock after a step wave of polarizing current, with the steady-state current density J_ss during the current step wave as follows: E_p =ρ' J_ss . This relation may be interpreted as a generalized Ohm's law, valid in linear cases, in which ρ’(fictitious resistivity) is defined as the product of the true resistivity ρ with the chargeability m. Supposing E _p=— grad U_p and applying the divergence condition div J_ss = o, one can, for a layered earth, obtain a general expression for the depolarization potential U_p as a solution of Laplace's equation ?²U_p= o. Since the mathematical procedure for the solution of this last equation is identical to that used in resistivity problems, we propose now the introduction of an apparent fictitious resistivity ρ'_a (defined as the product of the apparent resistivity ρ_a with the apparent chargeability m_a) as a new parameter for the interpretations of IP soundings carried out over layered structures with a common electrode array. The most general expression of ρ'_a as a function of the electrode distance turns out to be mathematically identical to the general expression of ρ'_a. Therefore it is possible to interpret a ρ'_a field curve using the same standard graphs for resistivity prospecting with the usual method of complete curve matching. In this manner a great deal of work is saved since there is no need to construct proper m_a graphs for the interpretation of IP soundings, as it has been done up to now. Finally some field examples are reported.     

Two notes on the modern theory of inverse problems in geophysics

On the basis of comparison of the approaches to the solution of inverse problems in information theory and geophysics, it is shown that results, obtained in information theory, are suitable to supplement the theory of geophysical inverse problems. The conditions of the existence and uniqueness of the solutions of inverse problems in their practical discrete statement are specified. The terms of ɛ-entropy H _ɛ and informational capacity C _ɛ, characterizing “volumes” of unknown and observed data, are introduced. It is shown, that the instability of the solution of the inverse problem decreases with increase in H _ɛ, (increase in the “complexity” of studied section), if the relation H _ɛ ≤ C _ɛ is maintained.     

TIME DOMAIN ELECTROMAGNETIC RESPONSE OF A SHIELDED CONDUCTOR*

Velekin and Bulgakov (1967) in an interesting model experiment while studying the transient electromagnetic response of a conductive sphere placed below a thin conductive sheet found that at the earlier stages of the transience, the composite system response corresponded to the response due to the overlying sheet alone and at the later stages, it corresponded to that of the sphere alone. To examine whether such a separation of responses due to individual components can be analytically studied and applied to other source configurations, we have analyzed an idealized model consisting of two spherical shells. We find that in corroboration with the above results, the general nature of the curve consists of two humps representing the responses dominated by the outer shell and the inner shell respectively. In addition, however, we find that the two humps gradually disappear to yield a smooth decay curve for increasing values of the ratio σd₁b/σ₂d₂a (where σ₁, σ₂ are the conductivities, d₁, d₂ are the thicknesses of the outer and inner shells respectively, and b and a are their respective distances from the centre) and the effect of inner shell on the composite system response is considerably reduced.