SUBSURFACE GRAVITY MEASUREMENTS AND DENSITY MAPPING IN THE MAILARAM COPPER MINES,ANDHRA PRADESH,INDIA*

Gravity measurements were made in the Mailaram copper mines, Andhra Pradesh. The observations were distributed between the two shafts situated about 220 m apart and in the three levels up to a maximum depth of 100 m. Assuming a normal free-air gradient, average densities for the three layers were determined as 2.631, 2.604, and 2.823 g cm^-3, respectively. Upon incorporating the weighted mean density values from measurements on samples, the free-air gradients were found to be 0.315 mGal m^-1 for the second layer (i.e. between the first and second levels) and 0.2978 mGal m^-1 for the third layer (i.e. between the second and third levels). The density variation map obtained from the gravity data, the deduced anomalies, the weighted mean density values from measurements on rock samples, and the varying free-air coefficients all suggest correspondence with the concentration of ore lodes.     

A numerical-physical planetary boundary layer model

A numerical-physical model for the planetary boundary layer has been formulated for the purpose of predicting the winds, temperatures and humidities in the lowest 1600 m of the atmosphere. An application of the model to the synoptic situation of 30 August, 1972, demonstrates its ability to produce useful forecasts for a period of 24 h. Results are illustrated in terms of horizontal maps and time-height sections of winds and temperatures. The model is divided in the vertical direction into three layers that are governed, respectively, by different physical formulations. At the lowest level, which is the surface of the earth, forecasts of temperature and humidity are computed from empirical relations. In the first layer, the surface layer, application is made of the similarity theories of Monin-Obukhov, Monin-Kazanski and Businger’s form of the universal functions. The second layer, the Ekman layer, is 1550 m deep and is governed by diagnostic momentum and time-dependent thermodynamic and humidity equations. External input to the model are large-scale pressure gradients and middle-level cloudiness. Cressman’s objective analysis procedure is applied to conventional surface and upper air data over a horizontal region of about 2500 km by 2500 km, centered about Lake Ontario. With a grid distance of 127 km and a time interval of 30 min, the computer time required on Control Data Cyber 76 for a 24 h forecast for the case study is less than two minutes.     

Wavelet correlation filter for wide-angle seismic data

A new filtering technique for single‐fold wide‐angle reflection/refraction seismic data is presented. The technique is based on the wavelet decomposition of a set of adjacent traces followed by coherence analysis. The filtering procedure consists of three steps. In the first, a wavelet decomposition of traces into different detail levels is performed. In the second, the coherence attributes for each level are evaluated by calculating cross‐correlation functions of detail portions contained in a space–time moving window. Finally, the filtered traces are obtained as a weighted reconstruction of the trace details. Each weight is obtained from the coherence‐attributes distribution estimated in a proper interval. A sequence of tests is then conducted in order to select possible optimum or unsuitable wavelet bases. The efficiency of the filter proposed was assessed by calculating some properly designed parameters in order to compare it with other standard de‐noising techniques. The proposed method produced a clear signal enhancement in high‐density wide‐angle seismic data, thus proving that it is a useful processing tool for a reliable correlation of seismic phases.     

地形对长偏移距瞬变电磁测深的影响研究(英文)

用基于张量格林函数的体积分方程法对三维异常体进行瞬变电磁响应的正演模拟,首先在频率域内计算电磁场分量的频率域响应,然后利用快速数字滤波技术将计算结果转换到时间域。设计和计算了水平电偶极子源激发下层状水平地层模型背景下的常见地形如山谷、山峰地形的模型,并考察分别把源和接收器放于这些地形中的瞬变电磁场响应,详细分析了这些地形对长偏移距瞬变电磁测深(LOTEM)的影响。结果表明,山谷和山峰地形对LOTEM的结果均有不同程度的影响。当电偶极子源放在山谷谷底时,地形对观测异常场的畸变非常严重;当接收器放在山谷中时,接收器处地形的影响强烈但该影响在空间和时间上只是局部的。总体来讲,不论山峰地形位于何处,其对LOTEM的影响相对较小。当地形处于发射源与接收器之间时,地形对LOTEM的影响非常小,表明在进行LOTEM勘探时,选择发射源的放置比接收器的位置更加重要,野外勘探是尽量把发射源选择在开阔的平坦位置。     

An alternative method for density variation modeling of the crust based on 3-D gravity inversion

In this paper we are proposing an alternative method for determination of density variations of the crust from constrained inversion of the terrestrial gravity data. The main features of the method can be summarized as follows: (i) Constructing a band-pass filter to remove the long and short wavelength signals from the terrestrial gravity data. (ii) Using an iterative method for stabilization and solution of the inverse problem. The mentioned regularization method is first validated by simulated gravity data and next the methodology is used for development of a new regional density variation model of the crust in three layers based on real gravity data in geographical area of Iran. Application of the band-pass filter to the latter data resulted the residual gravitation variations in the range of − 300 to 50 (mGal) which next based on the iterative method resulted following ranges for residual densities: −120 to 40 (kg/m³) in first layer, −40 to 40 (kg/m³) in second layer, and − 40 to 40 (kg/m³) in third layer.     

A Comprehensive Method for Fractionating Soil Organic Matter Not Protected and Protected from Decomposition by Physical and Chemical Mechanisms

The objective of this work was to describe a method for isolating meaningful and measurable soil organic matter (SOM) pools that differ in the mechanisms by which they are protected from decomposition. The proposed method is appropriate for soil C stabilization and sequestration studies. Unlike previous fractionation schemes, this procedure allows free SOM located between aggregates (unprotected C pool) and SOM occluded within both macroaggregates and microaggregates (C weakly and strongly protected by physical mechanisms, respectively) to be recovered separately, freed from the soil mineral matrix and the mineral‐associated SOM pool (C pool protected by chemical mechanisms) and thus well suited to advanced chemical characterization by ¹³C‐NMR. Briefly, free SOM is isolated by an initial density separation. Stable macroaggregates are broken up into stable microaggregates and intra‐macroaggregate SOM, which is then separated by density. Finally, intra‐microaggregate SOM is isolated from mineral‐associated SOM by a third density separation after ultrasonic disruption. The SOM dissolved during the fractionation procedure is also recovered. Results obtained on soil samples with contrasting textures suggested that clay content induces a decrease of the proportion of free organic C and an increase of mineral‐associated organic C content. Free SOM is characterized by a marked presence of undecayed organic material and biologically labile substances, such as carbohydrates and proteins. In contrast, SOM occluded within aggregates, especially within microaggregates, represents a more decomposed fraction, relatively enriched in unsubstituted‐aliphatic material, most probably lipid biopolymers.     

LOW-PASS FILTERING OF NOISY FIELD SCHLUMBERGER SOUNDING CURVES PART II: APPLICATION*

The basic principles of the application of the linear system theory for smoothing noise-degraded d.c. geoelectrical sounding curves were recently established by Patella. A field Schlumberger sounding is presented to demonstrate first their application and validity. To achieve this purpose, firstly it is pointed out that the required smoothing or low-pass filtering can be considered as an intrinsic property of the transformation of original Schlumberger sounding curves into pole-pole (two-electrode) curves. Then we sketch a numerical algorithm to perform the transformation, opportunely modified from a known procedure for transforming dipole diagrams into Schlumberger ones. Finally we show a field example with the double aim of demonstrating (i) the high quality of the low-pass filtering, and (ii) the reliability of the transformed pole-pole curve as far as quantitative interpretation is concerned.     

Interpretation of self-potential anomalies of some simple geometric bodies

Summary An entirely new procedure for interpreting selfpotential anomalies of spheres, rods and dipping sheets is presented. The anomaly of a sphere is divided into two parts — the anomaly of odd symmetry and the anomaly of even symmetry — from which the depth can be obtained by fitting them with the master curves. The self-potential anomalies of a finite rod are transformed to the anomalies of a veritcal sheet, for which standard curves are presented. The case of a sheet was divided into three parts; (a) finite line of poles; (b) infinite double line of poles and (c) finite double line of poles. For the first case logarithmic curves were prepared and presented; by their comparison with the field profile, different parameters can be obtained. In the second case, a geometrical construction is provided to obtain the various values. In the third case, the anomalies of finite sheet (finite double line of poles) are transformed into those due to an infinite double line of poles for interpretation.     

全球大洋混合层深度的计算及其时空变化特征分析

本文利用2005-2009年的全球网格化Argo数据,分别采用温度判据和密度判据计算了全球大洋混合层深度(Mixed Layer Depth, MLD),讨论了障碍层(Barrier Layer, BL)和补偿层(Compensated Layer, CL)对混合层深度计算的影响,得到了合成的混合层深度,并研究了其时空变化特征. 研究表明:(1)在赤道西太平洋(10°S -5°N,150°E-150°W),孟加拉湾,热带西大西洋(10°N-20°N,30°W-60°W)是障碍层高发区域. 冬季的北太平洋副热带区域(30°N附近)以及东北大西洋(40°N-60°N,0°-30°W)是补偿层发生的区域. (2) 在各个半球的夏季MLD都比较浅,在各个半球的冬季MLD则普遍比较深. 北太平洋和北大西洋的MLD的分布和变化比较相似,印度洋MLD受季风影响显著,呈现半年周期变化. 太平洋和大西洋的MLD 的经向分布大致呈现出"两端深,中间浅"的拱形特点. (3)混合层深度距平场EOF第一模态时间变化为周期的年信号,北太平洋和北大西洋、南大洋(尤其是南极绕流区)都是MLD变化剧烈的海域,第二模态显示全球大洋混合层深度距平存在着一个半年的变化周期.     

ESTIMATING INFILTRATION PARAMETERS FROM BASIC SOIL PROPERTIES

Infiltration data were collected on two rectangular grids with 25 sampling points each. Both experimental grids were located in tropical rain forest (Guyana), the first in an Arenosol area and the second in a Ferralsol field. Four different infiltration models were evaluated based on their performance in describing the infiltration data. The model parameters were estimated using non-linear optimization techniques. The infiltration behaviour in the Ferralsol was equally well described by the equations of Philip, Green–Ampt, Kostiakov and Horton. For the Arenosol, the equations of Philip, Green–Ampt and Horton were significantly better than the Kostiakov model. Basic soil properties such as textural composition (percentage sand, silt and clay), organic carbon content, dry bulk density, porosity, initial soil water content and root content were also determined for each sampling point of the two grids. The fitted infiltration parameters were then estimated based on other soil properties using multiple regression. Prior to the regression analysis, all predictor variables were transformed to normality. The regression analysis was performed using two information levels. The first information level contained only three texture fractions for the Ferralsol (sand, silt and clay) and four fractions for the Arenosol (coarse, medium and fine sand, and silt and clay). At the first information level the regression models explained up to 60% of the variability of some of the infiltration parameters for the Ferralsol field plot. At the second information level the complete textural analysis was used (nine fractions for the Ferralsol and six for the Arenosol). At the second information level a principal components analysis (PCA) was performed prior to the regression analysis to overcome the problem of multicollinearity among the predictor variables. Regression analysis was then carried out using the orthogonally transformed soil properties as the independent variables. Results for the Ferralsol data show that the parameters of the Green–Ampt and Kostiakov model were estimated relatively accurately (maximum R² = 0.76). For the Arenosol, use of the second information level together with PCA produced regression models with an R² value ranging from 0.38 to 0.68. For the Ferralsol, most of the variance was explained by the root content and organic matter content. In the Arenosol plot, the fractions medium and fine sand explained most of the observed variance.     

Lateral variability of flow over a sill in a channel of southern Chile

Measurements of velocity and density profiles were used to describe the tidal and mean flow structure across and along a sill in Refugio Channel, a fjord-like inlet in Southern Chile (43.9°S). These are the first oceanographic measurements of any kind effected in Refugio Channel. Current profiles were obtained with a 307.2-kHz acoustic Doppler current profiler during two semidiurnal cycles along a repeated triangular circuit. Two along-channel transects formed the sides of the triangle that crossed the sill and were identified as the western and eastern transects. One cross-channel transect, the base of the triangle, was located on the seaward side of the sill. Density profiles were obtained at the corners of the triangle. The longitudinal mean flow in the western transect showed a two-layer exchange structure over the landward side of the sill. The structure of net seaward flow at the surface and landward flow at depth was disrupted by the sill in such a way that over the seaward side of the sill, only seaward flow was observed throughout the water column. This likely resulted from the blocking of landward net flow by the sill. In the eastern transect, two-layer exchange dominated over most of the transect and was consistent with the observed density profiles. Over the seaward side of the sill, a surface layer, ∼10m deep, flowed landward as a third layer. This feature should have been caused by river input further seaward (to the north) and produced a surface convergence region over the sill. In terms of tidal flows, the greatest tidal current amplitudes were 40cm s⁻¹ over the sill as the flow accelerated through the reduced cross-sectional area of the channel. Near-surface flow convergences were identified over both along-channel transects.     

Fallout radionuclides in the Pacific Ocean: Vertical and horizontal distributions,largely from GEOSECS stations

From GEOSECS stations, largely, the 1974 distributions of Pu and of¹³⁷Cs are described in the Pacific Ocean north of about 20°S latitude. Changes in some of these distributions are described from 1978 cruises by the authors.The Pacific exhibited, everywhere, a shallow subsurface layer of Pu-rich water with its concentration maximum at about 465 m in 1974; over a large portion of the central North Pacific a second layer of Pu-labelled water, less concentrated than the shallow layer, lay just above the bottom. Similar features were not observed in the case of¹³⁷Cs.The inventories of both Pu and¹³⁷Cs in the water column at most 1974 stations are substantially greater than those to be expected from world-wide fallout alone; these inventory excesses appear to be attributable to close-in fallout, but only if the ratio Pu/¹³⁷Cs in this source was much higher than in world-wide fallout. The North Pacific mean ratio of the inventories is 2.2 times that observed in world-wide fallout.Resolubilization of Pu both from sinking particles and from sediments explains peculiarities of its depth distributions.There is little evidence for tracer movement by sliding downward along density surfaces;¹³⁷Cs appears to have moved to depth by downmixing at the edge of the Kuroshio, and then moved horizontally and upward alongσ_t contours. The shallow Pu-rich layer shows no coordination with density, salinity or O₂ isopleths. The deep Pu-rich layer is restricted to a narrow range of O₂ concentrations that confirm its origin in the Aleutian Trench and rapid spread southward and laterally. Near-bottom circulation processes have been much more active than here-to-fore described.     

The equatorial dynamics of a shallow,homogeneous ocean

Abstract

It is shown that the linear equatorial dynamics of a shallow ocean is characterized by two boundary layers of width γ^? L and γL (γ is the Ekman number of the flow, assumed small, and L is a horizontal dimension of the basin). In the γ^? layer stress in the bottom Ekman layer is comparable to that in the surface Ekman layer. In the γ layer vertical friction is important throughout the depth of the ocean. Should the Rossby number ? be so large as to invalidate a linear theory (? > γ^5/3), then inertial effects become important at a distance ?^2/5 L from the equator. The role played in the circulation of the basin by the non-linear equatorial current first studied by Charney (1960) is shown to be similar to that of the γ layer of the linear theory. Though lateral friction is unimportant in a linear model of the flow, shear layers at the equator are found to be a necessary feature of non-linear flow.     

A design procedure for buildings equipped with energy dissipation devices using nonclassical damping and iso-performance curves

This article describes a design procedure for elastic buildings equipped with linear and nonlinear energy dissipating devices. The objective is to achieve a design that responds to a target building performance following a simple and robust step-by-step algorithm. The proposed procedure identifies first the modal significance of key design performance indicators and controls the modal properties by solving a singular two-parameter eigenvalue problem. For that purpose, a new modal significance metric is proposed, and a target frequency shift and damping ratio for the complete structure are obtained from the so-called iso-performance design curves. The design algorithm employs linear-equivalent stiffness and damping properties, which are then transformed into parameters characterizing inelastic force-deformation constitutive models corresponding to physical devices. The design algorithm leads to an optimal damper distribution corresponding to the minimum global amount of supplemental equivalent damping needed to achieve a maximum modal perturbation. The design procedure is first demonstrated using a five-story building example and then a real and complex 22-story free-plan building with two towers of rhomboid-shape plan with a very singular dynamic behavior.     

Hypolimnetic density currents traced by sulphur hexafluoride (SF6)

The artificial tracer sulphur hexafluoride (SF₆) has been used to study the density-driven deep water exchange between two sill-separated basins of Lake Lucerne, Gersauersee and Urnersee. The sources of the density gradients between the two basins are (1) salinity differences between the major inlets due to the different geology of their drainage areas, and (2) temperature differences due to spatial variation of wind forcing. Wind speeds are generally larger in Urnersee, especially in spring during the so-called Föhn events, when winds blow from the south. In contrast, Gersauersee is protected form these winds. In spring 1989, a total of 630 g of SF₆ was released at 80 to 120 m depth in the small Treib Basin located between Urnersee and Gersauersee. During about 100 days the distribution of SF₆ in the lake was determined by gaschromatography. Two models are used to quantify the exchange flow, (1) a one-box mass balance model for SF₆ in the deep part of Treib Basin, and (2) a one-dimensional diffusion/advection model describing the temporal and vertical temperature variation in Urnersee. According to the first model, the flow into the deep hypolimnion of Urnersee, decreases from 21·10⁶ m³·d^?1 at the end of March to about 8·10⁶ m³·d^?1 in late April. The second model yields similar flow rates. The decrease of the flow rate during spring, confirmed by both approaches, is consistent (1) with the decreasing strength of the density gradient above the sill during spring and early summer, and (2) with hydrographic information collected in Lake Lucerne during other years.     

Fission track dating of kimberlitic zircons

The only reliable method for dating kimberlites at present is the lengthy and specialized hydrothermal procedure that extracts²⁰⁶Pb and²³⁸U from low-uranium zircons. This paper describes a second successful method by fission track dating of large single-crystal zircons, 1.0–1.5 cm in dimension. The use of large crystals overcomes the limitations imposed in conventional fission track analysis which utilizes crushed fragments. Low track densities, optical track dispersion, and the random orientation of polished surfaces in the etch and irradiation cycle are effectively overcome.Fission track ages of zircons from five African kimberlites are reported, from the Kimberley Pool (90.3 ± 6.5 m.y.), Orapa (87.4 ± 5.7 and 92.4 ± 6.1 m.y.), Nzega (51.1 ± 3.8 m.y.), Koffiefontein (90.0 ± 8.2 m.y.), and Val do Queve (133.4 ± 11.5 m.y.). In addition we report the first radiometric ages (707.9 ± 59.6 and 705.5 ± 61.0 m.y.) of crustal zircons from kimberlites in northwest Liberia. The fission track ages agree well with earlier age estimates. Most of the zircons examined in this study are zoned with respect to uranium but linear correlations are established (by regression analysis) between zones of variable uranium content, and within zones of constant uranium content (by analysis of variance). Concordance between the fission track method and the U/Pb technique is established and we concluded that track fading from thermal annealing has not taken place. Kimberlitic zircons dated in this study, therefore, record the time of eruption.     

Cross-correlation estimates of telluric and magnetotelluric parameters

Procedures using cross-correlation functions to analyze telluric and magnetotelluric field data can be designed which, in certain applications, are more efficient than conventional techniques using fast Fourier transforms. One such application, involving the processing ofband-limited data, is presented here. The linear coupling relations between fields in a limited frequency band are estimated from transient time series by minimizing, in a least squares sense, the residuals between observed and predicted values of the frequency coefficients. The resulting normal equations contain integral averages over the continuous auto- and cross-energy spectra which are efficiently evaluated as Fourier transforms of windowed auto- and cross-correlation functions in the time domain. The method is outlined in general terms, then illustrated with a specific example involving the analysis of 30–80 second pulsation data. The procedure involves three stages:

Stage I: Data sections approximately 2000 seconds long are digitally sampled at 1 second intervals, filtered at 50 seconds and decimated to 500 point series with 4 second sampling intervals.

Stage II: The correlation functions are formed for 16 lags (of 4 secs. each) on either side of zero and multiplied by a Gaussian window.

Stage III: These modified correlation functions are Fourier transformed at the single period, 50 seconds; the band-averaged energy spectra which result are used to solve the desired field relations for either the telluric or the magnetotelluric coupling coefficients.

Several built-in advantages are demonstrated. Since the window is concentrated around the origin, the correlation functions only need to be calculated out to relatively small lags. Decimation of the data during filtering further improves efficiency by dramatically reducing the number of points summed over. Visual display of intermediate results throughout the analysis not only provides added insight, but improves reliability by pinpointing problems early. For the examples considered here, linear coupling coefficients are shown to be stable within about 5 percent over several data sets. The method gives results within 1 percent of those determined by fast Fourier transform techniques while using half the computer time.     

Reflected and refracted waves from a linear transition layer — A periodic source

Summary The problem of a transiton layer lying between two homogeneous liquid media has been discussed in a previous paper. There, some of the terms in the integrands in the expressions for the displacement potentials due to a periodic point source, lying in the upper layer, have been evaluated along the Sommerfeld loops. Both first order and second order discontinuities in density and bulk modulus at the boundaries of the transition layer have been considered. In the present paper, the integrals for the displacement potentials have been evaluated by the method of stationary phase and the results obtained are compared with those obtained previously. Also the first order discontinuities in density and velocity at the two ends of the transition layer are removed and the work is extended to an explosion type pulse. It is found that for a pulse source, the decrease with distance is the same as for a periodic point source. Also the frequency dependence for the refraction arrivals, when frequency is extremely high, is ^1/3 as compared to ^–1 in the case when there are discontinuities at the boundaries of the transition layer.     

Ray tracing in elliptical anisotropic media using the linear traveltime interpolation (LTI) method applied to traveltime seismic tomography

The linear traveltime interpolation (LTI) method is a suitable ray‐tracing technique for modelling first‐arrival times in isotropic media. LTI is extended to elliptical anisotropic media and applied to a tomographic inversion procedure. A theoretical formulation is first derived and then LTI implementation is discussed in terms of source–receiver arrays and cell size. The method is then combined with the tomographic inversion procedure adopted. The matching of the ray tracing with inversion in elliptical anisotropic media posed a double non‐linear problem. Thus two assumptions were made: the velocity in each cell is uniform and the main directions of anisotropy are known. To take into account the geometrical characteristics of the area under investigation (depth and velocity of the weathering, and thickness of the inner media), cells of varying size were considered. No hypothesis was made on anisotropy weakness.
The algorithm was first tested on synthetic models and then applied to a field survey. On comparing the results of the synthetic models and the field survey with those obtained with a linear raypath approximation, it was found that there were fewer data misfits.