ON TAKING MEANS OF DENSITY DISTRIBUTIONS IN THE EARTH

Summary. The error is investigated when the mean of two solutions of a non-linear second-order differential equation which gives the variation of density with depth in certain layers of the Earth is used as a solution.
A solution (Model M) of the equation for layer E is given, which has a density of 9.385 g/m³ at the top of layer E.
Identifying, at the top of layer E, the mean of two numerically estimated solutions with a similarly estimated solution, it is shown that the error is a maximum at the bottom of layer E, where it is of the order of 0.04 per cent.     

EARTH MODELS WITH CONTINUOUS DENSITY DISTRIBUTION

Summary. Eight different Earth models have been set up, all with the property that the density p varies continuously from just below the crust to the centre. The distributions of the pressure p , gravity g , incompressibility k and rigidity are also given; and values of a parameter equal to (k/p)dp/dp , indicate the deviations from (chemical) homogeneity in the lower mantle and outer core. The models are designed to provide a numerical background towards testing the view that there is no density jump between the Earth's mantle and core.
A discussion shows that this view is difficult to reconcile with a homogeneous core unless an implausibly low value is assumed for the density just below the crust.     

引黄灌溉与黄河断流

引黄灌溉平均每年为黄河下游地区增产粮食４９．１×１０８ｋｇ;山东省引黄灌区粮食产量年增长速率为６．５％。黄河断流,若以两年一次的洛口—利津河段同时断流统计,则每年至少损失粮食１７．２１×１０８ｋｇ;平均每断流一天,直接经济损失达３９３０×１０４元。因此,研究、缓解黄河断流,保证、发展引黄灌溉,是黄河下游农业持续发展的关键。     

Secular variation of the poloidal magnetic field at the core–mantle boundary

A region of enhanced conductivity at the base of the mantle is modelled by an infinitesimally thin sheet of uniform effective conductance adjacent to the core–mantle boundary. Currents induced in this sheet by the temporally varying magnetic field produced by the geodynamo give rise to a discontinuity in the horizontal components of the poloidal magnetic field on crossing the sheet, while the radial component is continuous across the sheet. Treating the rest of the mantle as an insulator, the horizontal components of the poloidal magnetic field and their secular variation at the top of the core are determined from geomagnetic field, secular variation and secular acceleration models. It is seen that for an assumed effective conductance of the sheet of 10⁸  S, which may be not unrealistic, the changes produced in the horizontal components of the poloidal field at the top of the core are usually ≤10 per cent, but corrections to the secular variation in these components at the top of the core are typically 40 per cent, which is greater than the differences that exist between different secular variation models for the same epoch. Given the assumption that all the conductivity of the mantle is concentrated into a thin shell, the present method is not restricted to a weakly conducting mantle. Results obtained are compared with perturbation solutions.     

Seismic Modelling of Lateral Heterogeneity At the Base of the Mantle

Summary. Results from several recent studies suggest that there are lateral heterogeneities of up to a few per cent in the lowermost 150–200 km of the mantle (Bullen's D " region). Inferred anomaly sizes span the range from less than 50 km to greater than 1000 km.
In this study differences in the velocity structure among regions at the base of the mantle were inferred from an analysis of amplitude ratios of PKPAB and PKPDF for given earthquake-station pairs at distances greater than 155° (Sacks, Snoke & Beach). We distinguish two kinds of regions: A (anomalous) regions in which the mean, median and spread in AB/DF amplitude ratios are significantly higher (> 50 per cent) than for a reference radial earth model and N (normal) regions in which the distribution of the amplitude ratios is as expected.
The AB branch has near-grazing incidence to the core and therefore maximum sensitivity to velocity structure compared to the near-normal incident DF phases. Using an iterative, forward-modelling approach, we have determined general characteristics of the velocity structure for regions at the base of the mantle which can produce amplitude-ratio distributions similar to those for an A region. Agreement between model and data is obtained over the period range from 0.5 s to greater than 10 s using a laterally heterogeneous model for the D " region. the model consists of cells which are 200 km in lateral extent with velocity variations of up to ±1 per cent. This structure is modulated by a region-wide (1000km) perturbation which increases smoothly from zero at the edges of the region to a negative 1 per cent at the centre. Small cells (∼40 km) cannot produce anomalously large amplitude, long-period AB arrivals, and larger cells (∼1000km) cannot match the observed scatter. the ∼200 km scale anomalies could be small-scale convection cells confined to the D " region.     

Upper mantle shear structure of North America

Summary. The waveforms and travel times of S and SS phases in the range 10°–60° have been used to derive upper mantle shear velocity structures for two distinct tectonic provinces in North America. Data from earthquakes on the East Pacific Rise recorded at stations in western North America were used to derive a tectonic upper mantle model. Events on the north-west coast of North America and earthquakes off the coast of Greenland provided the data to investigate the upper mantle under the Canadian shield. All branches from the triplications due to velocity jumps near 400 and 660 km were observed in both areas. Using synthetic seismograms to model these observations placed tight constraints on heterogeneity in the upper mantle and on the details of its structure. SS–S travel-time differences of 30 s along with consistent differences in waveforms between the two data sets require substantial heterogeneity to at least 350 km depth. Velocities in the upper 170 km of the shield are about 10 per cent higher than in the tectonic area. At 250 km depth the shield velocities are still greater by about 4.5 per cent and they gradually merge near 400 km. Below 400 km no evidence for heterogeneity was found. The two models both have first-order discontinuities of 4.5 per cent at 405 km and 7.5 per cent at 695 km. Both models also have lids with lower velocities beneath. In the western model the lid is very thin and of relatively low velocity. In the shield the lid is 170 km thick with very high elocity (4.78 km s^-1); below it the velocity decreases to about 4.65 km s^-1. Aside from these features the models are relatively smooth, the major difference between them being a larger gradient in the tectonic region from 200 to 400 km.     

P-wave velocity anisotropy in crystalline rocks

Summary. Compressional wave velocities and anisotropy coefficients determined at high hydrostatic pressures are compiled from the data published for the main types of crystalline rocks. The crack-free elastic anisotropy of igneous crustal rocks is generally very low, between 1 and 3 per cent on average. The anisotropy of metamorphic rocks is higher (up to 22 per cent), but very variable. The average anisotropy coefficients in schists and amphibolites are about 10 per cent, in gneisses between 3 and 7 per cent, and in granulites less than 3 per cent. The average anisotropy of olivine ultramafites is between 7 and 12 per cent, whereas in pyroxenites and eclogites it is usually less than 4 per cent. A comparison of ranges of average velocities and average anisotropies for the individual rock groups suggests that, whereas in the crust the lateral velocity variations are mainly due to compositional changes, in the olivine of the uppermost mantle the velocity variations due to anisotropic structures could be of the same magnitude as the variations due to inhomogeneities.     

Locating some types of random errors in Digital Terrain Models

The increasing use of Geographical Information System applications has generated a strong interest in the assessment of data quality. As an example of quantitative raster data, we analysed errors in Digital Terrain Models (DTM). Errors might be classified as systematic (strongly dependent on the production methodology) and random. The present work attempts to locate some types of randomly distributed, weakly spatially correlated errors by applying a new methodology based on Principal Components Analysis. The Principal Components approach presented is very different from the typical scheme used in image processing. A prototype implementation has been conducted using MATLAB, and the overall procedure has been numerically tested using a Monte Carlo approach. A DTM of Stockholm, with integer-valued heights varying from 0 to 59 m has been used as a testbed.The model was contaminated by adding randomly located errors, distributed uniformly within 4 m and 4m. The procedure has been applied using both spike shaped (isolated errors) and pyramid-like errors. The preliminary results show that for the former, roughly half of the errors have been located with a Type I error probability of 4.6 per cent on average, checking up to 1 per cent of the dataset. The associated Type II error of the larger errors (of exactly 4m or 4 m) drops from an initial value of 1.21 per cent down to 0.63 per cent. By checking another 1 per cent of the dataset, such error drops to 0.34 per cent implying that about 71 per cent of the 4m errors have been located; Type I error was below 11.27 per cent. The results for pyramid-like errors are slightly worse, with a Type I error of 25.80 per cent on average for the first 1 per cent effort, and a Type II error drop from an initial value of 0.81 per cent down to 0.65 per cent. The procedure can be applied both for error detection during the DTM generation and by end users. It might also be used for other types of quantitative raster data.     

Sensitivity and inversion of borehole flexural dispersions for formation parameters

Propagation characteristics of borehole flexural waves are functions of five borehole, fluid and formation parameters. Uncertainties in any of these model parameters may introduce errors in the estimate of formation shear-wave speed. Other sources of error in the estimate may be caused by deviations from the assumed circular borehole cross-section or heterogeneity in the material properties of formation in the vicinity of a borehole from the usually assumed homogeneous properties. The influence of uncertainties in model parameters on borehole flexural dispersion has been calculated from a general perturbation model based on Hamilton's principle. A sensitivity analysis of the flexural dispersion to small variations in the model parameters shows that the formation shear speed has by far the most dominant influence in a slow formation. In contrast, the flexural dispersion in a fast formation is significantly influenced by three of the five model parameters: the formation shear speed, the borehole fluid compressional speed and the borehole diameter. The frequency dependence of these sensitivity functions indicates that the inversion of flexural dispersion for formation shear speed is optimal in the range 2–4 kHz for a borehole of diameter 25.4 cm. The range of validity of the perturbation model has been estimated by comparing results of concentric annuli of different thicknesses, and shear-wave speeds different from that of the formation with those from exact numerical solutions from a modal search program. Generally, the perturbation results for altered annuli of thicknesses up to 15 cm are accurate to within 1 per cent for a shear-wave speed 10 per cent lower than that of the formation. This difference between the perturbation and exact results increases to approximately 1 to 3 per cent for a shear-wave speed 20 per cent lower than that of the formation.     

冬小麦水分耗散特性与农业节水

本文根据中国科学院禹城综合试验站１９８６年－１９９６年蒸渗仪农田水分模拟试验资料统计；冬小麦全生育期耗水量可达４８２．５ｍｍ，缺水率可达６９．３％；冬小麦全生育不分耗散过程有两个明显的需水峰区和３个关键需水期，为实施节水灌溉提供了实验依据；冬小麦耗水量与环境因子有明显的相关关系，其统计规律可供地下水浅埋区区域灌溉预测参照应用。     

Travel-time anomalies in the mantle under the North Atlantic

Summary. Lateral heterogeneity exists in the Earth's mantle, and may result in seismic velocity anomalies up to several per cent. If convection cells and plumes extend down to the core, then these features may be associated with local inhomogeneities observed in the lower mantle.
Published data for direct and core-reflected P -wave residuals are used to delineate velocity anomalies in the middle—lower mantle under the North Atlantic. Differential ( PcP — P ) residuals indicate travel-time anomalies near the core—mantle transition, and may be due to core topography or lateral variations in velocity. It is assumed that the anomalies occur near the midpoints of the ray paths. The main source of error in the data set may arise from phases which have been identified incorrectly. Hence trend-surfaces are fitted to the residual data to show only the large-scale trends in anomaly values, with wavelengths of the order of 1000 km.
The Azores and Colorado hot spots occur in a region covered by the data. A possible interpretation of the trend maps is that an anomalous zone extends from a relatively fast region at the core boundary at 35° N, 50° W up to these hot spots, at about 30 degrees from the vertical. This may agree with the suggestion of Anderson that plumes are chemical rather than thermal in origin. If inclⁱned plumes do exist, the deviation from the ideal vertical plume or convection cell boundary may imply that lateral shear or other distortion effects exist in the mantle.     

The additivity of partial thermal remanent magnetization in magnetite

Summary. Experiments were done to test the additivity of partial thermal remanent magnetizations (PTRMs) for prepared samples containing magnetite particles whose sizes range from SD (single domain) to MD (multidomain). The experiments compare the sum of two PTRMs with total-TRM, all produced by the same external field of 0.47 oe. The most significant conclusion of this paper is that, to first order, the additivity of PTRMs is obeyed for the magnetites of this study regardless of particle size. However, small, higher order deviations from additivity occur such that ΣPTRM > TRM by an average of about 1 per cent. Though small, these departures from additivity are significant at the 99 per cent confidence level, and they can be understood in terms of Néel's theory for SD particles. The small departures from additivity are intrinsic to the experimental procedure in which some particles acquire remanence twice, in each of the two PTRM steps. In the limit of small inducing fields additivity should be obeyed exactly for the magnetites of this study and for samples of interest in palaeomagnetism. The deviations from additivity should have no effect on palaeointensity determinations by the Thelliers' version of the Thellier palaeointensity method. For palaeointensity determinations by Coe's version of the Thellier method the effects of deviations from additivity would be very small, less than 4 per cent on the average for a worst-case experimental configuration, and these effects can be minimized by producing PTRMs parallel to the original NRM and by using weak laboratory fields.     

The salinity of hypersaline brines: Concepts and misconceptions

Some widespread misconceptions about salinity are pointed out, explained and commented on. Definitions, units, and methods of salinity determination recommended for brines of moderate salinities may become unsuitable for highly saline brines. Various salinity units are examined and commented on, and density at a precisely monitored temperature is chosen as the preferable “substitute parameter” for salinity determination of hypersaline brines. It is shown that ambient air density during pycnometric measurements must be known as accurately as required of the brine's density, and that its value must be added to the direct pycnometric measurement. Variations of atmospheric pressure at the time of pycnometry must be taken into account if greater than 15 mB. For a salinity accuracy of ±0.02 per mille, as required for some physical and chemical studies, the temperature of the sample during salinity measurement must be monitored with an accuracy of at least 0.04°C. A definite curve in the ρ-S plane, corresponding to the conversion of salinity to density at a fixed temperature, andvice-versa, does not exist if the brine is saturated in one of its salts. Non-linearity in the equation of state is shown to affect the surface level drop due to evaporation; the effect is negligible at low salinities, but a 10 per cent correction is already required at the salinity found in the northern Red Sea (41 per mille). Reliable reference points for the conversion of absolute salinity to density of highly saline brines are not known, at present, as accurately as desirable; pending the accurate determination of such reference points, a substantial dilution of hypersaline brines (for the purpose of adapting them to conductivity measurements) would multiply the error range by more than a hundred fold. On the other hand, a minute dilution of hypersaline brine samples to prevent salt crystal formation is possible and would not bring about any major increase of the error range in their density (salinity) determination.     

A three-dimensional radially anisotropic model of shear velocity in the whole mantle

We present a 3-D radially anisotropic S velocity model of the whole mantle (SAW642AN), obtained using a large three component surface and body waveform data set and an iterative inversion for structure and source parameters based on Non-linear Asymptotic Coupling Theory (NACT). The model is parametrized in level 4 spherical splines, which have a spacing of ∼ 8°. The model shows a link between mantle flow and anisotropy in a variety of depth ranges. In the uppermost mantle, we confirm observations of regions with   V_SH > V_SV   starting at ∼80 km under oceanic regions and ∼200 km under stable continental lithosphere, suggesting horizontal flow beneath the lithosphere. We also observe a   V_SV > V_SH   signature at ∼150–300 km depth beneath major ridge systems with amplitude correlated with spreading rate for fast-spreading segments. In the transition zone (400–700 km depth), regions of subducted slab material are associated with   V_SV > V_SH   , while the ridge signal decreases. While the mid-mantle has lower amplitude anisotropy (<1 per cent), we also confirm the observation of radially symmetric   V_SH > V_SV   in the lowermost 300 km, which appears to be a robust conclusion, despite an error in our previous paper which has been corrected here. The 3-D deviations from this signature are associated with the large-scale low-velocity superplumes under the central Pacific and Africa, suggesting that   V_SH > V_SV   is generated in the predominant horizontal flow of a mechanical boundary layer, with a change in signature related to transition to upwelling at the superplumes.     

黄淮海平原农业综合开发的效益和粮食增产潜力

黄淮海平原经过1988-1993年农业综合开发,取得了明显的经济、生态和社会效益。改造中低产田409万ha,开发区粮食年平均增长率达4.0%,其增产量占全区增产总量的42.5%,农业开发的贡献几乎与传统农业的贡献相等。预测1995-2000年,即使是农业综合开发的规模与前六年相同,但是,由于水利工程效益更充分地发挥和农业的科技含量持续增长,开发区粮食完全可能增产达90亿kg.     

南极威尔克斯地SGR钻孔处雪冰转化过程中的密度变化特征

本文着重描述了SGR钻孔处冰盖上的积雪在密实化过程中的特征变化,并对该过程进行了分段的和全面的回归分析.结果表明,冰盖密度随深度增大,但增长幅度随深度减小.作者提出密度变化减小度的概念.计算得出的所研究冰芯钻取点的密度变化减小度为-0.15kg/m~3·m~2,粒雪成冰前的密实速率平均值为4.08kg/m~3·a.本文得到的冰盖密度变化“临界点”与以往报道的有所不同.分析这一现象时,作者强调当积雪还在活动层时冰盖温度的影响,并以此解释密度剖面的异常变化以及离差的回升.特别指出,积雪的密度变化具有气候学意义,它在一定程度上能够反映出积雪形成及变化过程中气候变化的某些信息.本文由密度变化确定的钻孔点雪冰转化深度为50米.     

田间土壤含水率的空间结构及取样数目确定

本文以田间试验资料为基础,研究田间土壤含水率的空间变异性与空间结构。对1m土层不同深度的含水率与1m土层平均含水率作了统计、自相关与半方差分析,表明田间土壤存在着空间变异性。1m土层平均含水率的变差系数Cv小于各不同深度含水率的Cv值。资料还表明含水率在平面位置上相关性不显著。将田间含水率作为空间变异的随机变量,可计算其符合要求的置信水平P1和估值精度的合理取样数目。本文提供的取样数目表亦可供确定合理取样数目时参考。     

Pn anisotropy studies in northern Britain and the eastern and western United States

Summary. Available seismic refraction data from three different continental areas, northern Britain and the eastern and western United States, has been studied for possible P_n , velocity anisotropy using the methods described by Bamford. There are various deficiencies in the time—distance data used in each case but, while the uppermost mantle beneath northern Britain and the eastern United States seems to be isotropic within the limits of measurement error, there is a small but significant anisotropy beneath the western United States.
Both the amount (up to 3 per cent) and the direction (70–80° east of north) of this anisotropy are very similar to the results obtained in the Pacific Ocean off California. We tentatively conclude that this anisotropy is present as a consequence of the subduction of oceanic lithosphere beneath the western United States.     

南水北调东线水质污染和湖泊水生经济的系统分析

本文讨论了我国南水北调东线水质污染的产生及其对调水水质的影响,同时对调水区域内工业排水与治理经济模式和湖泊水生经济系统的最优化分析进行了探论.     

Application of the boundary element method to a terrestrial heat flow problem

b
In this paper, we apply the boundary element method (BEM) to the 2-D steady state heat flow problem of what would be the perturbation to the regional temperature gradient, and hence heat flow density, determined from temperatures measured in a borehole that passes close to, but does not penetrate, a body of anomalous thermal properties. This type of problem with an infinity boundary is particularly well suited to the BEM.
The results have been compared with those obtained from analytical solutions for bodies of simple shape; it is found that for the worst case of a close approach to a boundary of small radius of curvature, a numerical modelling error of less than 1 per cent can still be obtained provided the length of each element is less than the shortest distance between the calculation point and the object.