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

Chang-Jo F. Chung 《Mathematical Geology》1993,25(7):851-865

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.  相似文献   

硅质碎屑砂岩中胶结物来源的新认识     

杨香华 《地质科技情报》1993,12(4):34-39

分析了浅埋藏砂岩的成岩环境,并从岩石学、矿物学、同位素及质量平衡估算方面论证硅质碎屑砂岩中的胶结物主要来源于砂岩浅埋藏期间的地下水,进而分析浅埋藏带内地下水的流动机制及早期胶结作用对晚期成岩作用的影响。  相似文献   

Estimation for partially observed Markov processes     

M. E. Thompson  T. N. Kaseke 《Stochastic Hydrology and Hydraulics》1995,9(1):33-47

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.  相似文献   

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia     

J. Zahradník  M. Antonini  G. Grünthal  J. Janský  D. Procházková  E. Schmedes  A. Špičak  J. Zedník 《Pure and Applied Geophysics》1990,133(1):53-71

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.  相似文献   

Mobility of heavy metals from tailings to stream waters in local zinc smelting activities contaminated site     

Taoze LIU Yuangen YANG Bangyu LIU Xiangyang BI Feili LI 《中国地球化学学报》2006,25(B08):41-42

It has been found that stream waters were severely contaminated with wastes from a long-time smelting factory in Hezhang, Guizhou, China. The main sources of contamination are the smelting wastes stored in the open air and abandoned in the vicinity of stream. A method of lead isotope was adopted in order to identify relations between tailings and water contamination. Representative samples of tailings and stream sediments were collected. Mineralogical characterizations were conducted using XRD and TEM/SEM, while acid digestion was carried out for determining metal contents. BCR sequential leaching tests were performed in order to assess metal mobility. The tremendous ‘actual＇ and ‘potential＇ mobility of heavy metals indicates that the smelting waste and stream sediments present a considerable threat to the environment. Besides the chemical remobilization of heavy metals from the sediments and the reworking of riverbed sediments act as a secondary source of pollution. Also groundwater and stream water were sampled in specific locations and were measured.  相似文献   

On the Importance of Geolithological Features for the Estimate of the Site Response: The Case of Catania Metropolitan Area (Italy)     

G. Lombardo  H. Langer  S. Gresta  R. Rigano  C. Monaco  G. De Guidi 《Natural Hazards》2006,38(3):339-354

The evaluation of seismic site response in the urban area of Catania was tackled by selecting test areas having peculiar lithological and structural features, potentially favourable to large local amplifications of ground motion. The two selected areas are located in the historical downtown and in the northern part of Catania where the presence of a fault is evident. Site response was evaluated using spectral ratio technique taking the horizontal- to-vertical component ratio of ambient noise. Inferences from microtremor measurements are compared with results from synthetic accelerograms and response spectra computed at all drillings available for this area. Such method is particularly suitable in urban areas where the nature of the outcropping geological units is masked by city growth and anthropic intervention on the surface geology. The microtremor H/V spectral ratios evaluated at soft sites located within the downtown profile tend to be smaller than that usually reported in the literature for such soils. A tendency for amplifications to peaks near 2 Hz is observed only in some sites located on recent alluvial deposits. Evidences for amplifications of site effects (frequency range 4–8 Hz) were observed in the sampling sites located on the fault, with a rapid decrease of spectral amplitude just a few tenths of metres away from the discontinuity. Numerical simulations evidenced the importance of geolithological features at depth levels even greater than 20–30 m. Besides this, the results strongly confirm the importance of the subsurface geological conditions, in the estimate of seismic hazard at urban scale.  相似文献   

New Generation of Probabilistic Seismic Hazard Assessment for the Area Cologne/Aachen Considering the Uncertainties of the Input Data     

G. Grünthal  R. Wahlström 《Natural Hazards》2006,38(1-2):159-176

A new earthquake catalogue for central, northern and northwestern Europe with unified M_w magnitudes, in part derived from chi-square maximum likelihood regressions, forms the basis for seismic hazard calculations for the Lower Rhine Embayment. Uncertainties in the various input parameters are introduced, a detailed seismic zonation is performed and a recently developed technique for maximum expected magnitude estimation is adopted and quantified. Applying the logic tree algorithm, resulting hazard values with error estimates are obtained as fractile curves (median, 16% and 84% fractiles and mean) plotted for pga (peak ground acceleration; median values for Cologne 0.7 and 1.2 m/s² for probabilities of exceedence of 10% and 2%, respectively, in 50 years), 0.4 s (0.8 and 1.5 m/s²) and 1.0 s (0.3 and 0.5 m/s²) pseudoacclerations, and intensity (I₀ = 6.5 and 7.2). For the ground motion parameters, rock foundation is assumed. For the area near Cologne and Aachen, maps show the median and 84% fractile hazard for 2% probability of exceedence in 50 years based on pga (maximum median value about 1.5 m/s²), and 0.4 s (>2 m/s²) and 1.0 s (about 0.8 m/s²) pseudoaccelerations, all for rock. The pga 84% fractile map also has a maximum value above 2 m/s² and shows similarities with the median map for 0.4 s. In all maps, the maximum values fall within the area 6.2–6.3° E and 50.8–50.9° N, i.e., east of Aachen.  相似文献   

Spontaneous Structural Reform in Australian Local Government: the case of the Gilgandra Co-operative Model in New South Wales     

Brian Dollery  Warrick Moppett  Lin Crase 《The Australian geographer》2006,37(3):395-409

Dire financial constraints and the threats of forced structural reform have had the unforeseen effect of forcing municipal councils across Australia to reconsider their operational activities and organisational arrangements. With considerable ingenuity, numerous municipalities have proposed and sometimes adopted new structural formations that embody various forms of co-operative service provision. This remarkable development has unfortunately been largely ignored in the scholarly literature on Australian local government. In a modest effort aimed at remedying this neglect, the present paper seeks to outline the small, rural New South Wales Gilgandra Shire Council's (2004) ‘Co-operative/Local Government Service Company’ model, place it in the broader context of alternative models of local governance suitable for Australian conditions, and evaluate its characteristics.  相似文献   

Re-evaluation of focal depths and source mechanisms of selected earthquakes in the Afar depression   总被引：1，自引：0，他引：1  

L. Hagos  H. Shomali  R. Roberts 《Geophysical Journal International》2006,167(1):297-308

We present a stepwise inversion procedure to assess the focal depth and model earthquake source complexity of seven moderate-sized earthquakes  (6.2 > M _w > 5.1)  that occurred in the Afar depression and the surrounding region. The Afar depression is a region of highly extended and intruded lithosphere, and zones of incipient seafloor spreading. A time-domain inversion of full moment tensor was performed to model direct P and SH waves of teleseismic data. Waveform inversion of the selected events estimated focal depths in the range of 17–22 km, deeper than previously published results. This suggests that the brittle–ductile transition zone beneath parts of the Afar depression extends more than 22 km. The effect of near-source velocity structure on the moment tensor elements was also investigated and was found to respond little to the models considered. Synthetic tests indicate that the size of the estimated, non-physical, non-isotropic source component is rather sensitive to incorrect depth estimation. The dominant double couple part of the moment tensor solutions for most of the events indicates that their occurrence is mainly due to shearing. Parameters associated with source directivity (rupture velocity and azimuth) were also investigated. Re-evaluation of the analysed events shows predominantly normal faulting consistent with the relative plate motions in the region.  相似文献   

Evidence of reworking and resuspension of carbonates during last glacial maximum and early deglacial period along the southwest coast of India     

Kamlesh Verma  M Sudhakar 《Journal of Earth System Science》2006,115(6):695-702

A gravity core collected from the upper slope of southwest of Quilon at a water depth of 776 m (Lat: 8°12′263″N, Long: 76°28′281″E) was analysed for texture (carbonate free), calcium carbonate and organic carbon. Variation in silicic fraction seems to be controlled by silt, i.e., enrichment from 15 ka BP to 10 ka BP and then constant in Holocene. Below 15 ka BP, the silicic fraction gets depleted compared to the Holocene section with a minimum around 21 ka BP. Clay content remains nearly constant except in the Holocene where it shows an enrichment. Carbonate content of less than 63 micron when computed by subtracting coarse fraction content from the total carbonate suggests that the total carbonates are mainly concentrated in the finer fraction. All these carbonate phases show an inverse relationship with silicic fraction except in Holocene. Below 15 ka BP, CaCO₃ dominates in sediments comprising more than 65%, such an increase is also seen in the coarse fraction. Coarse fraction from these sections contains abundant nodular type aggregates encrusting small forams. This period is marked by a high sedimentation rate comparable to Holocene. These parameters suggest that the productivity and precipitation have increased in the Holocene due to the intensification of the southwest monsoon. During the last glacial maximum and early deglacial period the high sedimentation rate indicates redeposition of the carbonates from the existing carbonate lithofacies situated between Quilon and Cape Comorin probably due to the slope instability.  相似文献