Chemical weathering indices are useful tools in characterizing weathering profiles and determining the extent of weathering. However, the predictive performance of the conventional indices is critically dependent on the composition of the unweathered parent rock. To overcome this limitation, the present paper introduces an alternative statistical empirical index of chemical weathering that is extracted by the principal component analysis (PCA) of a large dataset derived from unweathered igneous rocks and their weathering profiles. The PCA analysis yields two principal components (PC1 and PC2), which capture 39.23% and 35.17% of total variability, respectively. The extent of weathering is reflected by variation along PC1, primarily due to the loss of Na2O and CaO during weathering. In contrast, PC2 is the direction along which the projections of unweathered felsic, intermediate and mafic igneous rocks appear to be best discriminated; therefore, PC1 and PC2 represent independent latent variables that correspond to the extent of weathering and the chemistry of the unweathered parent rock. Subsequently, PC1 and PC2 were then mapped onto a ternary diagram (MFW diagram). The M and F vertices characterize mafic and felsic rock source, respectively, while the W vertex identifies the degree of weathering of these sources, independent of the chemistry of the unweathered parent rock.
The W index has a number of significant properties that are not found in conventional weathering indices. First, the W index is sensitive to chemical changes that occur during weathering because it is based on eight major oxides, whereas most conventional indices are defined by between two and four oxides. Second, the W index provides robust results even for highly weathered sesquioxide-rich samples. Third, the W index is applicable to a wide range of felsic, intermediate and mafic igneous rock types. Finally, the MFW diagram is expected to facilitate provenance analysis of sedimentary rocks by identifying their weathering trends and thereby enabling a backward estimate of the composition of the unweathered source rock. 相似文献
The Non-linear lterative Partial Least Squares(NIPALS)algorithm is used in principal componentanalysis to decompose a data matrix into score vectors and eigenvectors(loading vectors)plus a residualmatrix.N1PALS starts with some guessed starting vector.The principal components obtained by NIPALSdepends on the starting vector;the first principal component could not always be computed.Wold hassuggested a starting vector for NIPALS,but we have found that even if this starting vector is used,thefirst principal component cannot be obtained in all cases.The reason why such a situation occurs isexplained by the power method.A simple modification of the original NIPALS procedure to avoid gettingsmaller eigenvalues is presented. 相似文献