Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

POPULATION EXPANSION IN CHINA''''S METROPOLISES AND ITS SCALE-CONTROL

Metropolises are centers of politics, economy, culture and information, and also economic growth poles in comparatively large areas. After 1949, China's urban structure have changed greatly, from original "pyramid" shape into "upside -down pyramid" shape, that is, more megalopolises and metropolises, fewer medium-sized and small cities. The most distinguishing feature of urban population growth was the remarkable increase of population in megalopolises and metropolises. The key reason for urban population expansion is to want to obtain the benefits of scale -agglomeration economy of metropolises. The number of urban functions directly determines the quality and quantity of urban population. Metropolises appeal not only agricultural population but also the population of medium-sized and small cities, which put much more population pressure on metropolises. The ideal goals of scale-control of metropolitan population are discussed in terms of economy, society and ecollogy. In order to control the metropoli     

山东省德州市“中华碘泉”矿泉水的发现与开发

在山东省德州市首次探明国内稀有含碘饮用用天然矿泉水－－“中华碘泉”,研究表明,矿泉水是在新生代断陷盆地的构造背景和弱碱氧化的水文地球化学环境条件下,地下水与含水介南长期水－岩反应中经溶滤形成的,矿泉水水量丰富,水质优良,为含碘水苏打型矿泉水,但氟含量偏讥;为合理开发这一富足资源,可选用电渗析,倾析方法降低氟含量,并将其作为规模化含碘矿泉水生产的技术工艺。     

TOWARD ESTABLISHING THE CONCEPT OF PHYSICAL URBAN AREA IN CHINA

I.IntroduchollWiulegeOgraphystudiestheterritoria1systemsabouttherelationsbforeenhumanbeingSandthenature[I4],urbanmpaphyfocusesspecificallyonthespatialoopltizationandprocessofurbanterritory,anditsinteraCtionwiti1Otherterritorialsystems.Usually,wedoresearch…     

On the periodic solutions of a particular Lame equation

We consider the Hill's equation: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca% WGKbWaaWbaaSqabeaacaaIYaaaaOGaeqOVdGhabaGaamizaiaadsha% daahaaWcbeqaaiaaikdaaaaaaOGaey4kaSYaaSaaaeaacaWGTbGaai% ikaiaad2gacqGHRaWkcaaIXaGaaiykaaqaaiaaikdaaaGaam4qamaa% CaaaleqabaGaaGOmaaaakiaacIcacaWG0bGaaiykaiabe67a4jabg2% da9iaaicdaaaa!4973!\[\frac{{d^2 \xi }}{{dt^2 }} + \frac{{m(m + 1)}}{2}C^2 (t)\xi = 0\]Where C(t) = Cn (t, {frbuilt|1/2}) is the elliptic function of Jacobi and m a given real number. It is a particular case of theame equation. By the change of variable from t to defined by: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaqcaawaaOWaaiqaaq% aabeqaamaalaaajaaybaGaamizaGGaaiab-z6agbqaaiaadsgacaWG% 0baaaiabg2da9OWaaOaaaKaaGfaacaGGOaqcKbaG-laaigdajaaycq% GHsislkmaaleaajeaybaGaaGymaaqaaiaaikdaaaqcaaMaaeiiaiaa% bohacaqGPbGaaeOBaOWaaWbaaKqaGfqabaGaaeOmaaaajaaycqWFMo% GrcqWFPaqkaKqaGfqaaaqcaawaaiab-z6agjab-HcaOiab-bdaWiab% -LcaPiab-1da9iab-bdaWaaakiaawUhaaaaa!51F5!\[\left\{ \begin{array}{l}\frac{{d\Phi }}{{dt}} = \sqrt {(1 - {\textstyle{1 \over 2}}{\rm{ sin}}^{\rm{2}} \Phi )} \\\Phi (0) = 0 \\\end{array} \right.\]it is transformed to the Ince equation: (1 + · cos(2)) y + b · sin(2) · y + (c + d · cos(2)) y = 0 where % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGGipm0dc9vqaqpepu0xbbG8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaqcaawaaiaadggacq% GH9aqpcqGHsislcaWGIbGaeyypa0JcdaWcgaqaaiaaigdaaeaacaaI% ZaGaaiilaiaabccacaWGJbGaeyypa0Jaamizaiabg2da9aaacaqGGa% WaaSaaaKaaGfaacaWGTbGaaiikaiaad2gacqGHRaWkcaaIXaGaaiyk% aaqaaiaaiodaaaaaaa!4777!\[a = - b = {1 \mathord{\left/{\vphantom {1 {3,{\rm{ }}c = d = }}} \right.\kern-\nulldelimiterspace} {3,{\rm{ }}c = d = }}{\rm{ }}\frac{{m(m + 1)}}{3}\]In the neighbourhood of the poles, we give the expression of the solutions.The periodic solutions of the Equation (1) correspond to the periodic solutions of the Equation (3). Magnus and Winkler give us a theory of their existence. By comparing these results to those of our study in the case of the Hill's equation, we can find the development in Fourier series of periodic solutions in function of the variable and deduce the development of solutions of (1) in function of C(t).     

城市生态景观的构建研究

运用景观生态学、旅游规划学、城市规划理论来认识城市生态景观并分析其具有和谐性、整体性、多样性、畅达性、安全性、可持续性等特点，阐述以人为本、尊重地域和历史文化、生态、满足市民休闲与游憩需求等构建原则，认为重点抓好自然廊道和斑块等自然生态景观的建设、重视人工廊道和城市建筑密集区的景观效果、完善城市生态景观的空间结构体系等是构建和谐城市生态景观的重要途径。     

西安南郊丰水年秋季土壤水分研究

 通过对西安南郊丰水年秋季土壤含水量的测定，研究了地下0~6 m之间土壤含水量的变化与土壤干层的恢复问题。研究结果表明，丰水年西安南郊人工林下2~3 m的土壤含水量在20％左右，远远高于正常年份的土壤含水量，而且大于表层和深层的土壤含水量。分析得出，土壤干层发育较弱的地区在降水丰富的年份可以得到一定程度的恢复；西安地区基本适合进行人工造林；通过人工措施，增强降水入渗和土壤含水量，对树木成活和长期的生长应当具有明显的作用。     

河南省城区经济发展研究

文章首先介绍了城区经济的概念、特征和类型，在分析河南省城区经济发展状况的基础上，论述了该省城区经济的地位和作用，提出了河南省城区经济发展构想以及对策措施。     

Large eddy simulation of hydrocyclone—prediction of air-core diameter and shape

Hydrocyclones are widely used in the mining and chemical industries. An attempt has been made in this study, to develop a CFD (computational fluid dynamics) model, which is capable of predicting the flow patterns inside the hydrocyclone, including accurate prediction of flow split as well as the size of the air-core. The flow velocities and air-core diameters are predicted by DRSM (differential Reynolds stress model) and LES (large eddy simulations) models were compared to experimental results. The predicted water splits and air-core diameter with LES and RSM turbulence models along with VOF (volume of fluid) model for the air phase, through the outlets for various inlet pressures were also analyzed. The LES turbulence model led to an improved turbulence field prediction and thereby to more accurate prediction of pressure and velocity fields. This improvement was distinctive for the axial profile of pressure, indicating that air-core development is principally a transport effect rather than a pressure effect.     

Cross-formational hydrocarbon fluid flows in the Tertiary deltaic system of the Beaufort–Mackenzie Basin

Molecular and isotopic compositions of crude oils in the Beaufort–Mackenzie Basin confirm three genetic end-member oil groups and suggest extensive cross-formational hydrocarbon fluid flows in the Tertiary deltaic system. Inter- and intra-fractional variations in the geochemistry of the Tertiary-reservoired oils indicate that the oil source/maturity signatures were substantially masked by biomarkers that were picked up along migration pathways. Thus, many of the previously recognized “immature non-marine oils” are in fact thermally mature, probably derived from unpenetrated deeper marine source rocks. Although the effective source rock volumes have not been evaluated and their exact stratigraphic levels remain unknown, the relative timing of oil generation versus trap formation, rather than poor source quality, may be the cause of under-filled traps in the offshore area.