首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 93 毫秒
1.
A distribution-free estimator of the slope of a regression line is introduced. This estimator is designated Sm and is given by the median of the set of n(n – 1)/2 slope estimators, which may be calculated by inserting pairs of points (X i, Yi)and (X j, Yj)into the slope formula S i = (Y i – Yj)/(X i – Xj),1 i < j n Once S m is determined, outliers may be detected by calculating the residuals given by Ri = Yi – SmXi where 1 i n, and chosing the median Rm. Outliers are defined as points for which |Ri – Rm| > k (median {|R i – Rm|}). If no outliers are found, the Y-intercept is given by Rm. Confidence limits on Rm and Sm can be found from the sets of Ri and Si, respectively. The distribution-free estimators are compared with the least-squares estimators now in use by utilizing published data. Differences between the least-squares and distribution-free estimates are discussed, as are the drawbacks of the distribution-free techniques.  相似文献   

2.
A study of the lithogeochemistry of metavolcanics in the Ben Nevis area of Ontario, Canada has shown that factor analysis methods can distinguish lithogeochemical trends related to different geological processes, most notably, the principal compositional variation related to the volcanic stratigraphy and zones of carbonate alteration associated with the presence of sulphides and gold. Auto- and cross-correlation functions have been estimated for the two-dimensional distribution of various elements in the area. These functions allow computation of spatial factors in which patterns of multivariate relationships are dependent upon the spatial auto- and cross-correlation of the components. Because of the anisotropy of primary compositions of the volcanics, some spatial factor patterns are difficult to interpret. Isotropically distributed variables such as CO 2 are delineated clearly in spatial factor maps. For anisotropically distributed variables (SiO 2 ), as the neighborhood becomes smaller, the spacial factor maps becomes better. Interpretation of spatial factors requires computation of the corresponding amplitude vectors from the eigenvalue solution. This vector reflects relative amplitudes by which the variables follow the spatial factors. Instability of some eigenvalue solutions requires that caution be used in interpreting the resulting factor patterns. A measure of the predictive power of the spatial factors can be determined from autocorrelation coefficients and squared multiple correlation coefficients that indicate which variables are significant in any given factor. The spatial factor approach utilizes spatial relationships of variables in conjunction with systematic variation of variables representing geological processes. This approach can yield potential exploration targets based on the spatial continuity of alteration haloes that reflect mineralization.List of symbols c i Scalar factor that minimizes the discrepancy between andU i - D Radius of circular neighborhood used for estimating auto- and cross-correlation coefficients - d Distance for which transition matrixU is estimated - d ij Distance between observed valuesi andj - E Expected value - E i Row vector of residuals in the standardized model - F(d ij) Quadratic function of distanced ij F(d ij)=a+bd ij+cd ij 2 - L Diagonal matrix of the eigenvalues ofU - i Eigenvalue of the matrixU;ith diagonal element ofL - N Number of observations - p Number of variables - Q Total predictive power ofU - R Correlation matrix of the variables - R 0j Variance-covariance signal matrix of the standardized variables at origin;j is the index related tod andD (e.g.,j=1 ford=500 m,D=1000 m) - R 1j Matrix of auto- and cross-correlation coefficients evaluated at a given distance within the neighborhood - R m 2 Multiple correlation coefficient squared for themth variable - S i Column vectori of the signal values - s k 2 Residual variance for variablek - T i Amplitude vector corresponding toV i;ith row ofT=V –1 - T Total variation in the system - U Nonsymmetric transition matrix formed by post-multiplyingR 01 –1 byR ij - U i Componenti of the matrixU, corresponding to theith eigenvectorV i;U i= iViTi - U* i ComponentU i multiplied byc i - U ij Sum of componentsU i+U j - V i Eigenvector of the matrixU;ith column ofV withUV=VL - w Weighting factor; equal to the ratio of two eigenvalues - X i Random variable at pointi - x i Value of random variable at pointi - y i Residual ofx i - Z i Row vectori for the standardized variables - z i Standardized value of variable  相似文献   

3.
The Chayes-Kruskal procedure for testing correlations between proportions uses a linear approximation to the actual closure transformation to provide a null value,p ij , against which an observed closed correlation coefficient,r ij , can be tested. It has been suggested that a significant difference betweenp ij andr ij would indicate a nonzero covariance relationship between theith andjth open variables. In this paper, the linear approximation to the closure transformation is described in terms of a matrix equation. Examination of the solution set of this equation shows that estimation of, or even the identification of, significant nonzero open correlations is essentially impossible even if the number of variables and the sample size are large. The method of solving the matrix equation is described in the appendix.  相似文献   

4.
The magnetic fabrics of 150 samples from 25 sites from the English Lake District were determined using the magnetic susceptibility anisotropy technique. The susceptibility ellipsoid for 22 of the sites is oblate and for the remaining 3 prolate. The maximum intermediate susceptibility planes of the ellipsoids always coincide with the recorded cleavages which have a Caledonian strike and are sub-vertical. The orientations of the principal susceptibility axes are found to be in excellent agreement with those of principal strains determined by researchers using petrofabric methods on samples from near the author's sites. A correlation of the magnetic susceptibility anisotropy with the recorded strains shows that the magnetic anisotropy parameters can be related to those of strain by the empirical relationship: (for i = 1, 2, 3; j = 1, 2, 3 and i ≠ j), where χi and χj are orthogonal principal axes and li and lj are the corresponding orthogonal principal strain axes. The exponent a for the correlated sites from the andesitic tuffs of the Borrowdale volcanic group in the English Lake District is 0.142 ± 0.001. This index corresponds closely to that (0.145 ± 0.005) found for the Cambrian slate belt of North Wales (Rathore, 1979).  相似文献   

5.
Pyrite and arsenopyrite are the predominant ore minerals in the Korean Au-Ag deposits of this study. The X Ni py , X Co py , X Ni apy , and X Co apy values range between 100 and 3,000 ppm, 200 and 6,000 ppm, 200 and 8,200 ppm, and 100 and 10,200 ppm, respectively. Most X Ni py /X Co py values fall in the field lower than values varying 0.16–1.30. Arsenopyrite also tends to prefer cobalt rather than nickel showing X Ni apy /X Co apy values between 0.20 and 1.40. The concentrations of minor elements in ores and gangue minerals vary 1–55 ppm Au and 1–1,120 ppm Ag for the former and 4–57 ppm Ni and 2–45 ppm Co for the latter. The Au/Ag ratio in ore has a good correlation to the Ni/Co ratio of arsenopyrite to gangue. The (Ni/Co)py-(Ni/Co)gangue and (Ni/Co)apy-(Ni/Co)gangue diagrams revealed that the values from the Korean Au-Ag deposits plot in the field lower than 900 °C which is the lowermost temperature determined by previous partitioning experiments.  相似文献   

6.
Strain analysis of the Baraitha conglomerate is attempted by direct measurements on extracted pebbles and by micrometric analysis. The overall deformation is of flattening type, with thek value lower by more than half in the matrix than in the pebbles. The viscosity contrast between pebbles and matrix (μ i/μm) is in the ratio of 2:1 and the bulk deformation appears to be strongly controlled by Ci (concentration of pebbles expressed as percentage). The total shortening (≃35%) in the Baraitha conglomerate is comparable with the shortening accomplished in the folding of the overlying Bijawar Group volcanosedimentary sequence. The bulk strain axesX t, Yt andZ t, as determined from the analysis of the deformed conglomerate, are unsymmetrically oriented with reference to folds formed by oblique flexural-slip with neitherX t norY tcoincident with the fold hinges. The lack of transection of folds by cleavage again suggests flattening deformation. The extension in theY tdirection is greater in the matrix than in the pebbles.  相似文献   

7.
We present new experimental data on diffusion of divalent cations in almandine-spessartine diffusion couples in graphite capsules in the P-T range of 14–35 kb, 1100–1200° C. The tracer diffusion coefficients of the major divalent cations, viz. Fe, Mg and Mn, retrieved from the multicomponent diffusion profiles, have been combined with earlier data from our laboratory at 29–43 kb, 1300–1480° C (Loomis et al. 1985) to derive expressions of the P-T dependence of the diffusion coefficients at fO2 approximately corresponding to that defined by equilibrium in the system graphite-O2. We review the conditions, discussed earlier by Cooper, under which the flux of a component in a multicomponent system becomes proportional to its concentration gradient (Fickian diffusion), as if the entire solvent matrix behaves as a single component, and also suggest a method of incorporating the thermodynamic effect on diffusion in the same spirit. Regardless of the magnitude or sign of the off-diagonal terms of the D matrix, it is always possible to define an effective binary diffusion coefficient (EBDC) of a component in a semi-infinite multicomponent diffusion-couple experiment such that it has the property of the Fickian diffusion coefficient, provided that there is no inflection on the diffusion profiles. It is shown that the success of Elphick et al. in fitting the experimental diffusion profiles of all components over a limited concentration range by a single diffusion coefficient is due to fortuitous similarity of the EBDCs of the components (Fe, Mg, Mn and Ca) in their diffusion couple experiments. In common metapelitic garnets showing compositional zoning, the EBDCs of the divalent cations do not differ from each other by more than a factor of 2.5. However, the EBDC of a component changes from core to rim by a factor of 3 to 12, depending on the composition. We suggest a method of volume averaging of the EBDC which should prove useful in approximate calculations of diffusion flux during relaxation of compositional zoning. The EBDC of Mn is found to reduce essentially to D MnMn, the main diagonal term of the D matrix, and consequently can be calculated quite easily. Evaluation of EBDC of Fe, Mg and Mn in garnets from a prograde Barrovian sequence did not reveal any significant dependence on the extent of relaxation of garnet. The diffusion data have been applied to calculate the cooling rate of natural biotite-garnet diffusion couple from eastern Finland and diffusional modification of growth zoning in garnet in early Proterozoic Wopmay orogen, Canada. The results are in good agreement with geochronological and other independent constraints.Symbols and abbreviations a Radius of a spherical garnet crystal - BSE Back-scattered electron imaging - C Column vector of (n-1) independent components - D Diffusion coefficient matrix - D ij An element of the diffusion matrix on the i th row and j th column - D * i Tracer diffusion coefficient of component i - D(i) Effective interdiffusion coefficient (EIC) of various components in a multicomponent solution rich in the component i - D(i-j) Interdiffusion coefficient of components i and j in a binary solution - D i (EB) Effective binary diffusion coefficient of component i in a multicomponent solution - D i (EB:Ideal) D i (EB) under condition of ideal thermodynamic mixing of the diffusing species - D i (EB:thermo) Thermodynamic component of D i(EB) - D O A Interdiffusion coefficient at peak temperature T 0 in the phase A - D 0 Pre-exponential factor in an Arrhenius relation - EBDC Effective binary diffusion coefficient between a solute and a multicomponent solvent matrix - FEC Fixed edge composition model - EIC Effective interdiffusion coefficient - f i Fugacity of component i - HM Hematite-magnetite oxygen fugacity buffer - kb Kilobars - P Pressure - Q Activation energy (enthalpy) of diffusion - Extent of relaxation defined as the difference between core and rim compositions normalized to the same difference in the initial zoning profile - R Gas constant - s Cooling rate - T 0, T Ch Peak temperature and characteristic temperature, respectively - t Time - VEC Variable edge composition model - V + Activation volume - W ij Simple mixture interaction parameter between i and j - W i(EB) Effective simple mixture interaction parameter of a component i in a multicomponent solution - ij Margules interaction parameter between i and j - X i Mole fraction of component i - i Activity coefficient of component i - A dimensionless variable =D t/a 2 - ij Kronecker delta (i=j, ij =1; ij, ij =0) - Zi Charge on the ion i  相似文献   

8.
Simulation of geological surfaces using fractals   总被引：2，自引：0，他引：2  
Methods suggested in the past for simulated ore concentration or pollution concentration over an area of interest, subject to the condition that the simulated surface is passing through specifying points, are based on the assumption of normality. A new method is introduced here which is a generalization of the subdivision method used in fractals. This method is based on the construction of a fractal plane-to-line functionf(x, y, R, e, u), where(x, y) is in[a, b]×[c, d], R is the autocorrelation function,e is the resolution limit, andu is a random real function on [–1, 1]. The simulation using fractals escapes from any distribution assumptions of the data. The given network of points is connected to form quadrilaterals; each one of the quadrilaterals is split based on ways which are extensions of the well-known subdivision method. The quadrilaterals continue to split and grow until resolution obtained in bothx andy directions is smaller than a prespecified resolution. If thex coordinate of theith quadrilateral is in[a i ,b i ] and they coordinate is in[c i ,d i ], the growth of this quadrilateral is a function of(b i –a i ) and(d i –c i ); the quadrilateral could grow toward the positive or negativez axis with equal probability forming four new quadrilaterals having a common vertex.This paper was presented at Emerging Concepts, MGUS-87 Conference, Redwood City, California, 13–15 April 1987.  相似文献   

9.
Sodic pyroxene (jadeite content X jd=0.1–0.3) occurs locally as small inclusions within, albite porphyroblasts and in the matrix of hematite-bearing quartz schists in the Sanbagawa (Sambagawa) metamorphic belt, central Shikoku, Japan. The sodic, pyroxene-bearing samples are characteristically free from chlorite and their typical mineral assemblage is sodic pyroxene+subcalcic (or sodic) amphibole+phengitic mica+albite+quartz+hematite+titanite±epidote. Spessartine-rich garnet occurs in Mn-rich samples. Sodic pyroxene in epidote-bearing samples tends to be poorer in acmite content (average X Acm=0.26–0.50) than that in the epidote-free samples (X Acm=0.45–0.47). X Jd shows no systematic relationship to metamorphic grade, and is different among the three sampling regions [Saruta-gawa, Asemi-gawa and Bessi (Besshi)]. The average X Jd of the Saruta-gawa samples (0.21–0.29) is higher than that of the Asemi-gawa (0.13–0.17) and Bessi (0.14–0.23). The P-T conditions of the Asemi-gawa and Bessi regions are estimated at 5.5–6.5 kbar, >360°C in the chlorite zone, 7–8.5 kbar, 440±15°C in the garnet zone and 8–9.5 kbar, 520±25°C in the albite-biotite zone. Metamorphic pressure of the Saruta-gawa region is systematically 1–1.5 kbar higher than that of the Asemi-gawa and Bessi regions, and materials of the Saruta-gawa region have been subducted to a level 3–5 km deeper than materials that underwent metamorphism at equivalent temperatures and are now exposed in the Asemi-gawa and Bessi regions. Pressure slightly increases toward the north (structurally high levels) through the Sanbagawa belt of central shikoku. Two types of zonal structure were observed in relatively coarse-grained sodic pyroxenes in the matrix. One type is characterized by increasing X Jd from core to rim, the other type by decreasing X Jd from core to rim. Both types of zoned pyroxenes show an increase in X Fe 2+[=Fe2+/(Fe2++Mg)] from core to rim. The first type of zoning was observed in a sample from the chlorite zone of lowest grade, whereas the latter occurs in the garnet and albite-biotite zones of higher grade. The contrast in zonal structure implies that dP/dT during prograde metamorphism decreased with increasing metamorphic grade and may have been negative in some samples from the higher-grade zones. The estimated dP/dT of the prograde stage of the chlorite zone is 3.2 kbar/100°C, and that of the garnet and albite-biotite zones is -1.8 to 0.9 kbar/100°C. The variation of dP/dT at shallow and deep levels of a subduction system probably reflects the difference of heating duration and/or change in thermal gradient of the subduction zone by continuous cooling of the surrounding mantle.  相似文献   

10.
A new set of reversal experiments for coexisting ortho- and clinopyroxenes in the system CMAS at conditions between 1,000–1,570° C and 30 to 50 kb is presented and combined with literature data. Pyroxene behaviour, particularly that of clinopyroxene, is very complicated and different styles of Al incorporation into the pyroxene structure for low and high concentrations of Al are indicated, strongly influencing the exchange of the enstatite component between ortho- and clinopyroxene. Thermodynamic modelling of this exchange is problematic because of the large number of unknown coefficients compared to the number of experiments. Thermometry based on such models becomes very dependent on accuracy of experimental data and analyses of small quantities of elements. Despite this complexity very simple empirical thermometric equations are capable of reproducing experimental conditions in the systems CMS and CMAS over a wide range of P, T conditions. We derived the equation which gives a mean error of estimate of 25° C when applied to CMS and CMAS data.Abbreviations Used in the Text cpx clinopyroxene - di diopside, CaMgSi2O6 - en enstatite, Mg2Si2O6 - opx orthopyroxene - px Pyroxene - py pyrope - a i j activity of component i in phase j - activity coefficient - G P,T (A) molar Gibbs free energy difference of reaction (A) at P, T - X i j mole fraction of component i in phase j  相似文献   

11.
We carried out reversed piston-cylinder experiments on the equilibrium paragonite = jadeite + kyanite + H2O at 700°C, 1.5–2.5 GPa, in the presence of H2O-NaCl fluids. Synthetic paragonite and jadeite and natural kyanite were used as starting materials. The experiments were performed on four different nominal starting compositions: X(H2O)=1.0, 0.90, 0.75 and 0.62. Reaction direction and extent were determined from the weight change in H2O in the capsule, as well as by optical and scanning electron microscopy (SEM). At X(H2O)=1.0, the equilibrium lies between 2.25 and 2.30 GPa, in good agreement with the 2.30–2.45 GPa reversal of Holland (Contrib Miner Petrol 68:293–301, 1979). Lowering X(H2O) decreases the pressure of paragonite breakdown to 2.10–2.20 GPa at X(H2O)=0.90 and 1.85–1.90 GPa at X(H2O)=0.75. The experiments at X(H2O) = 0.62 yielded the assemblage albite + corundum at 1.60 GPa, and jadeite + kyanite at 1.70 GPa. This constrains the position of the isothermal paragonite–jadeite–kyanite–albite–corundum–H2O invariant point in the system Na2O–Al2O3–SiO2–H2O to be at 1.6–1.7 GPa and X(H2O)~0.65±0.05. The data indicate that H2O activity, a(H2O), is 0.75–0.86, 0.55–0.58, and <0.42 at X(H2O)=0.90, 0.75, and 0.62, respectively. These values approach X(H2O)2, and agree well with the a(H2O) model of Aranovich and Newton (Contrib Miner Petrol 125:200–212, 1996). Our results demonstrate that the presence or absence of paragonite can be used to place limits on a(H2O) in high-pressure metamorphic environments. For example, nearly pure jadeite and kyanite from a metapelite from the Sesia Lanzo Zone formed during the Eo-Alpine metamorphic event at 1.7–2.0 GPa, 550–650°C. The absence of paragonite requires a fluid with low a(H2O) of 0.3–0.6, which could be due to the presence of saline brines.  相似文献   

12.
We report for the first time the evidence for prograde high-pressure (HP) metamorphism preceding a peak ultrahigh-temperature (UHT) event in the northernmost part of the Madurai Block in southern India. Mg–Al-rich Grt–Ged rocks from Komateri in Karur district contain poikiloblastic garnet with numerous multi-phase inclusions. Although most of the inclusion assemblages are composed of gedrite, quartz, and secondary biotite, rare staurolite + sapphirine and spinel + quartz are also present. The XMg (=Mg/[Fe+Mg]) of staurolite (0.45–0.49) is almost consistent with that reported previously from Namakkal district in the Palghat–Cauvery Shear Zone system (XMg = 0.51–0.52), north of the Madurai Block. The HP event was followed by peak UHT metamorphism at T = 880–1040 °C and P = 9.8–12.5 kbar as indicated by thermobarometric computations in the Grt–Ged rock and associated mafic granulite. Symplectic intergrowth of spinel (XMg = 0.50–0.59, ZnO < 1.7 wt.%) and quartz, a diagnostic indicator of UHT metamorphism, probably formed by decompression at UHT conditions. The rocks subsequently underwent retrograde metamorphism at T = 720–760 °C and P = 4.2–5.1 kbar. The P–T conditions and clockwise exhumation trajectory of the Komateri rocks, comparable to similar features recorded from the Palghat–Cauvery Shear Zone system, suggest that the Madurai Block and the Palghat–Cauvery Shear Zone system underwent similar HP and UHT metamorphic history probably related to the continent–continent collision during the final stage of amalgamation of Gondwana supercontinent.  相似文献   

13.
Optimal discrimination among several groups can be achieved by simultaneous diagonalization of pooled within-group, W, and among-group, A, sums of squares and cross-product matrices formed by utilizing axial-ratio sample statistics of quartz grains belonging to different sieve grades. This method maximizes the ratio of among-group cross products to within-group cross product quadratic forms (V'AV/V'WV)and simultaneously yields discriminant scores whose correlation coefficients are zero for group means as well as for within each group. This procedure enables a simple Euclidean distance measure for partitioning the discriminant space for assignment. Although W–1 and Amatrices are symmetric, the W–1 Amatrix needed for multigroup discrimination is asymmetric and hence the eigenstructure of W–1 Ais obtained by simultaneous diagonalization of Wand Amatrices. The first four sample statistics (mean, standard deviation, skewness, kurtosis) of normalized axial-ratios are required for discrimination, although the mean and standard deviation are the most important discriminators.  相似文献   

14.
New experimental data on compositions of garnets in two-pyroxene — garnet assemblages in the system CaO –MgO –Al2O3 –SiO2 (CMAS) are presented for conditions between 1,100 and 1,570° C and 30 to 50 kb. Garnets in these assemblages become less calcic with increasing pressure. Garnet-orthopyroxene barometry (Al-solubility-barometry) pertinent to geobarometry for garnet lherzolites has been evaluated with a set of experimental data covering the range 900 to 1,570° C and 15 to 100 kb. Various formulations of this barometer work well to 75 kb. Phase equilibria are not sufficient to positively verify the thermodynamic validity of any of such models. Empirical garnet-orthopyroxene barometry at least in the system CMAS can be formulated to obtain a pressure estimate without previous temperature estimation (P(kb)=34.4-19.175 1n X Al M1 +17.702 1n X Ca M2 ). The potential application of an analogous garnetclinopyroxene equilibrium is limited because the amount of Ca-Tschermaks in natural clinopyroxenes is usually quite small in garnet lherzolites and many eclogites. The Ca-Mg exchange between garnet and clinopyroxene appears however sufficiently sensitive to pressure to allow calibration of a CMAS barometer. The reaction 3CaMgSi2O6+Mg3Al2Si3O12=3Mg2Si2O6+Ca3Al2Si3O12 has a V o of 3.5 cm3. The total pressure dependency of this reaction is however closer to a theoretical V o of about 5 cm3 when excess volume properties of the phases involved are taken into account. We have calibrated such a barometer (mean error of estimate 2.8 kb) for assemblages with pyrope-rich (py>80) garnets and orthopyroxenes. This may provide the basis for a geobarometer for eclogites from kimberlites.Abbreviations Used in the Text CaTs Ca-tschermak's molecule, CaAl2SiO6 - cpx clinopyroxene - di diopside, CaMgSi2O6 - en enstatite, Mg2Si2O6 - gr grossular, Ca3Al2Si3O12 - gt garnet - MgTs Mg-Tschermak's molecule, MgAl2SiO6 - opx orthopyroxene - px pyroxene - py pyrope, Mg3Al2Si3O12 - a i j activity of component i in phase j - activity coefficient - G(I) molar Gibbs free energy difference of reaction (I) at standard state unless specified otherwise - H(I), (H I) molar enthalpy (difference) of phase (reaction) (I) at standard state unless specified otherwise - S (I), (S I) molar entropy (difference) of phase (reaction) (I) at standard state unless specified otherwise - V o, (V I o) molar volume (difference) of phase (reaction) (I) at standard state - X i j mole fraction of component i in phase j  相似文献   

15.
Temperature and H2O activity can be determined with high precision using metamorphic mineral assemblages that define both a dehydration equilibrium and a temperature-sensitive cation-exchange equilibrium. Such determinations are obtained by applying the Gibbs method and then integrating two resulting differential equations, as illustrated here for the assemblage garnet-chlorite-quartz. The first equation, a geothermometer that monitors temperature based upon Fe–Mg exchange between garnet and chlorite, was calibrated using rocks at Pecos Baldy, New Mexico: 0=0.05 P(bars)–19.02 T(K)+4607 ln K D+24,156 with errors of ±8°C based upon analytical precision. The second equation monitors differences in the activity of water between specimens (1) and (2): 0=(0.1 X Mg–chl, 1 – 2.05)(P 2 – P 1) +[–33.02+5.96 ln(X Fe–chl, 1/X alm, 1)][T 2–T 1 –2.67 RT 1ln[a(H2O)2/a(H2O)1] +5.96 T 1ln(X Fe–chl, 2 X alm, 1/X Fe–chl, 1 X alm, 2).For samples equilibrated at the same pressure and temperature, microprobe analytical errors of 1% limit precision to ±0.01 a(H2O). For samples equilibrated at the same pressure but variable temperature, uncertainty of ±8°C limits precision to ±0.06 a(H2O). Extreme presure sensitivity requires that the H2O-barometer be applied only to rocks where pressure gradients are absent or well-constrained. The geothermometer gives temperatures in agreement with two other garnet-chlorite geothermometers (Dickenson and Hewitt 1986; Ghent et al. 1987) and with garnet-biotite geothermometry (ferry and Spear 1978) over the temperature range 350–520°C. Application of the relative H2O barometer shows variations in the activity of water approaching 0.30 in several study areas. Either pelitic schists commonly equilibrate with a fluid that is not pure H2O, or some pelitic rocks undergo metamorphism in the absence of a free fluid phase.  相似文献   

16.
Summary At the northeastern flank of Gebel Yelleq, northern Sinai, pure limestones of Upper Cretaceous age were subjected to a thermal overprint, caused by a c. 80m thick Tertiary olivine dolerite sill. Metasomatic supply of Si, Al, Fe, Mg and Ti was greater to the c. 7m wide upper than to the c. 25m wide lower thermal aureole. The greater width of the lower aureole is possibly due to a longer duration of the thermal overprint at this contact. Mineral assemblages in both aureoles are (from the contact outward):(i) clinopyroxene + garnet ± wollastonite + calcite(ii) garnet ± wollastonite + calcite;(iii) wollastonite + calcite.In places, late stage xenoblasts of apophyllite and witherite overgrow these assemblages. Garnets are grandites to melanites with Grs56–86Adr14–42Sch0–2Sps0–0.2Prp0 in the lower, and Grs29–94Adr5–64Sch0–12Sps0–0.2Prp0–1.7 in the upper aureole. Close to the upper contact, clinopyroxene is virtually pure diopside with X Mg = Mg/(Mg + Fe2+) = 0.97–1.0, whereas clinopyroxenes farther away from the upper contact and in the lower aureole have X Mg-values of 0.49 and 0.53, respectively.The minimum temperatures reached during contact metamorphism in the upper and lower aureole are defined by the lower stability limit of wollastonite. The temperatures are inferred with a calculated T-X(CO2) projection in the system CMASCH and are estimated at c. 290 °C and 380 °C for X(CO2) values of 0.05 and 0.25, respectively. A pressure of roughly 100 bar is estimated for the lower dolerite-limestone contact. As indicated by one-dimensional thermal modelling, a maximum temperature of 695 °C was attained at this contact, assuming a magma temperature of 1150 °C. Further modelling results indicate (i) wollastonite, which occurs first 13 m away from the lower contact, formed at a maximum temperature of c. 575 °C, (ii) there, wollastonite formation lasted for approximately 170 years and, (iii) at the outer rim of the lower aureole, the maximum temperature reached was 480 °C, and temperatures sufficient for wollastonite formation lasted for about 140 years.  相似文献   

17.
Minerals of the triphylite-lithiophilite, Li(Fe, Mn)PO4, and the triplite-zwieselite-magniotriplite series, (Mn, Fe, Mg)2PO4F, occur in the late stage period of pegmatite evolution. Unfortunately, neither are the genetic relationships between these phosphates fully understood nor are thermodynamic data known. Consequently, phosphate associations and assemblages from 8 granitic pegmatites — Clementine II, Rubicon II and III, and Tsaobismund (Namibia); Hagendorf-Süd and Rabenstein (Germany); Valmy (France); Viitaniemi (Finland) — have been tested for compositional zoning and intercrystalline partitioning of main elements by electron microprobe techniques. Although the selected pegmatites display varying degrees of fractionation, and the intergrowth textures indicate different genetic relationships between the phosphates, the plots of mole fractions X Fe=Fe/(Fe+Mn+Mg+Ca), X Mn=Mn/(Fe+Mn+Mg+Ca), and X Mg=Mg/(Fe+Mn+Mg+Ca) can be fitted relatively well with smooth curves in Roozeboom diagrams. Their deviations from symmetrical distribution curves are mainly dependent upon X Mg or X Ca, and upon non-ideal solutions. Surprisingly small differences between the partition coefficients were detected for intergrowths of different origin. However, the partitioning of shared components among coexisting phases is clearly dependent upon the conditions of formation. Compositional zoning is observed only when both Fe–Mn phosphates are intergrown mutually or with other Fe–Mn–Mg mineral solid-solutios. Thus, the zoning does not seem to be due to continuous crystallization, but to later diffusion processes. The triplite structure has preference for Mn, Mg, and Ca, while Fe prefers minerals of the triphylite series. A quantification of main element fractionation between minerals of the triphylite and the triplite series is possible in the cases where diffusion can be excluded. For the Fe/(Fe+Mn) ratios of core compositions an equation with a high correlation coefficient (R=0.988) was determined: Fe/(Fe+Mn)Tr=[Fe/(Fe+Mn)Li]/{2.737-(1.737)[Fe/(Fe+Mn)Li]} (Tr=triplite series, Li=triphylite series). Consequently, the Fe/(Fe+Mn) ratio of the triplite series can now also be used in the interpretation of pegmatite evolution, just like that of the triphylite series which has been successfully applied in the past.  相似文献   

18.
A novel experimental cell was developed for in situ measurements of transport phenomena in porous media using Fourier-Transform Infrared (FTIR) Spectroscopy. The technique was employed at ambient pressure in the temperatures range of 11–44 °C to study the H2O → D2O exchange between water-saturated weathered feldspars (bulk porosity of 5–19 vol% for feldspar) from granitic saprolites and a surrounding aqueous liquid. Such measurements are an important step for understanding internal weathering reactions of feldspars in soils and aquifers. Effective diffusion coefficients Deff for water in water-saturated porous feldspars were determined assuming one-dimensional diffusion in a quasi-homogeneous medium. The values of Deff vary from 7.2 × 10−10 to 1.9 × 10−11 m2/s and are 1–2 orders of magnitude lower than the diffusion coefficients (D) of protons and molecular H2O in liquid water. The activation energy for the H2O → D2O exchange process in porous feldspars ranges from 7.8 to 18.8 kJ/mol.The results imply that the effective diffusivity of water is mainly controlled by physical properties of the feldspars like porosity, pore connectivity, pore geometry and distribution. Perthitic feldspars with homogeneous pore distribution in the albitic lamellas have diffusional tortuosity factors X = D/Deff between 3 and 10 while alkali feldspars with inhomogeneously distributed and disconnected pores have much higher X values up to 129. Diffusion anisotropy has been verified for a vein perthite with diffusion perpendicular to the lamellas being faster by 0.3–0.5 log units than within the lamellas. It has to be emphasized that the study is based only on few selected feldspars, including perthitic feldspar, and additional work on samples with different weathering stages is needed to test the importance of the different parameters controlling diffusive transport in the pore system.  相似文献   

19.
Summary This paper describes corundum formation in a metasomatic reaction zonation around an ultramafic body within a metapelitic sequence. The investigated body is about 100 m in diameter and is located in the Saualpe of the Austroalpine nappe complex in the Eastern Alps. The body is surrounded by a 10 m wide reaction zone (here called zone 3) containing the paragenesis garnet–staurolite–biotite–margarite–chlorite–corundum. Beyond a further metasomatic transition zone (here called zone 2), there are undisturbed metapelitic host rocks (zone 1) that have the metamorphic peak paragenesis garnet–biotite–plagioclase–staurolite–muscovite–quartz. It is shown that reaction zonation formed around 7.2kbar and 615°C during regional metamorphism, just above the serpentine breakdown reactions in the system MgO–SiO2–H2O. Detailed analysis of the whole rock compositions shows that the reaction zonation formed by infiltration metasomatism that caused significant mass loss in the two alteration zones. These zones are particularly depleted in SiO2, Na2O and possibly K2O. An XNa2O–XSiO2 thermodynamic pseudosection is presented that includes the parageneses of both the unaltered metapelitic host rock and the corundum-bearing parageneses. This suggests that the metasomatic process can be explained by the transfer of SiO2 and Na2O alone. We interpret that the process is driven by water liberated from the previously serpentinised ultramafic body during prograde dehydration during regional Eo-Alpine metamorphism. This fluid flowed outwards from the ultramafic body depleting the surrounding pelites in silica and causing margarite formation from plagioclase and muscovite. This interpretation of the driving mechanism is consistent with our knowledge of the low water activities of the Saualpe during the Eo-Alpine orogenic cycle.  相似文献   

20.
The variance-based cross-variogram between two spatial processes, Z1 (·) and Z2 (·), is var (Z1 ( u ) – Z2 ( v )), expressed generally as a bivariate function of spatial locations uandv. It characterizes the cross-spatial dependence between Z1 (·) and Z2 (·) and can be used to obtain optimal multivariable predictors (cokriging). It has also been called the pseudo cross-variogram; here we compare its properties to that of the traditional (covariance-based) cross-variogram, cov (Z1 ( u ) – Z1 ( v ), Z2 ( u ) – Z2 ( v )). One concern with the variance-based cross-variogram has been that Z1 (·) and Z2 (·) might be measured in different units (apples and oranges). In this note, we show that the cokriging predictor based on variance-based cross-variograms can handle any units used for Z1 (·) and Z2 (·); recommendations are given for an appropriate choice of units. We review the differences between the variance-based cross-variogram and the covariance-based cross-variogram and conclude that the former is more appropriate for cokriging. In practice, one often assumes that variograms and cross-variograms are functions of uandv only through the difference u – v. This restricts the types of models that might be fitted to measures of cross-spatial dependence.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号