A la vista <strong>de</strong> la naturaleza <strong>de</strong> los ceros existentes, en nuestro <strong>estudio</strong> aplicaremos elreemplazamiento multiplicativo en función <strong><strong>de</strong>l</strong> umbral <strong>de</strong> cada parte (Martín-Fernán<strong>de</strong>z,2003). La cantidad <strong>de</strong> ceros presente en las diferentes partes y los patronesobservados nos sugieren <strong>de</strong>sarrollar nuestro análisis estadístico contemplandodiferentes subcomposiciones. Naturalmente, realizaremos el <strong>estudio</strong> incluyendo las 12partes <strong><strong>de</strong>l</strong> conjunto <strong>de</strong> datos. A continuación, repetiremos el análisis consi<strong>de</strong>randotodas las partes menos las tres tierras raras, es <strong>de</strong>cir, analizando la sub<strong>composición</strong> <strong>de</strong>9 partes que no incluya las tierras raras. Finalmente, el elevado número <strong>de</strong> cerospresentes nos aconseja repetir el <strong>estudio</strong> únicamente consi<strong>de</strong>rando las partes quecontienen una cantidad nula o <strong>de</strong>spreciable <strong>de</strong> ceros. En este caso analizaremos unasub<strong>composición</strong> <strong>de</strong> 6 partes.3.2. Principios iniciales <strong>de</strong> la interpretación geo<strong>química</strong> <strong>de</strong> los mo<strong><strong>de</strong>l</strong>os<strong>geoestadístico</strong>s.Nuestra etapa <strong>de</strong> análisis <strong>de</strong>scriptivo se inicia representando el biplot en los treshiperespacios correspondientes según la sub<strong>composición</strong> que analicemos. Tenemoshiperespacios <strong>de</strong> tres tipos: 12-dimensional, 9-dimensional y 6-dimensional. En losbiplots <strong>de</strong> las Figuras 3, 5 y 6, cada elemento químico o parte se i<strong>de</strong>ntifica con un rayo<strong><strong>de</strong>l</strong> hiperespacio <strong>de</strong> los datos clr-transformados. Distancias pequeñas entre los vértices<strong>de</strong> los rayos sugieren un gran grado <strong>de</strong> semejanza entre la variación <strong>de</strong> las partescorrespondientes. Las longitu<strong>de</strong>s <strong>de</strong> los rayos muestran la magnitud <strong>de</strong> la variabilidad<strong><strong>de</strong>l</strong> parte correspondiente. Los valores asociados a las distancias entre los vértices <strong><strong>de</strong>l</strong>os rayos ⎯Var(ln(X i /X j ))⎯ y asociados las longitu<strong>de</strong>s <strong>de</strong> los rayos ⎯Var(clr(X i ))⎯ semuestran en la Tabla 4.Tabla 4. Matriz <strong>de</strong> variabilidad relativa Var(ln(X i /X j )). La última columna muestra la variabilidad totalasociada a cada componente X iXjVar. totalXi Mn Fe Cu Zn Ga Rb Ba La Ce Nd Sr Pb XiMn 0 0,606 0,923 0,934 0,807 0,786 0,529 0,807 0,849 0,804 0,725 1,022 0,366Fe 0,606 0 0,416 0,651 0,292 0,325 0,173 0,259 0,303 0,274 0,209 0,986 0,187Cu 0,923 0,415 0 0,744 0,694 0,663 0,581 0,576 0,607 0,613 0,603 1,082 0,313Zn 0,934 0,651 0,744 0 0,925 0,928 0,852 1,041 1,088 0,984 0,736 1,282 0,424Ga 0,807 0,292 0,694 0,925 0 0,305 0,230 0,295 0,356 0,348 0,284 1,389 0,247Rb 0,786 0,325 0,663 0,928 0,305 0 0,208 0,339 0,365 0,436 0,139 1,013 0,229Ba 0,529 0,173 0,581 0,852 0,230 0,208 0 0,204 0,227 0,252 0,177 0,960 0,183La 0,807 0,259 0,576 1,041 0,295 0,339 0,204 0 0,079 0,126 0,345 1,139 0,217Ce 0,849 0,303 0,607 1,088 0,356 0,365 0,227 0,079 0 0,166 0,371 1,100 0,230Nd 0,804 0,274 0,613 0,984 0,348 0,436 0,252 0,126 0,166 0 0,392 1,120 0,230Sr 0,725 0,209 0,603 0,737 0,284 0,139 0,177 0,346 0,372 0,392 0 1,110 0,212Pb 1,022 0,986 1,081 1,282 1,389 1,0129 0,960 1,139 1,100 1,120 1,110 0 0,50820
3.3. Biplot en el clr-hiperespacio <strong>de</strong> 12-dimensiones3.3.1. Factores mineralógicosEn la Figura 3 se muestra el biplot composicional que se obtiene utilizando las 12partes <strong><strong>de</strong>l</strong> conjunto <strong>de</strong> datos. Este gráfico biplot está confeccionado en el espaciobidimensional <strong>de</strong> los dos primeros ejes asociados a la <strong>de</strong>s<strong>composición</strong> dvs. Estos dosprimeros ejes recogen un 52% <strong>de</strong> la variabilidad total, sugiriendo que la representaciónbiplot no tiene una calidad alta. En consecuencia, todas las características especialesasociadas al comportamiento <strong>de</strong> las partes y todas las relaciones entre ellas que seansugeridas por este biplot <strong>de</strong>berán ser tenidas en consi<strong>de</strong>ración <strong>de</strong> una manera muycuidadosa y, en ningún caso, <strong>de</strong> una manera concluyente.clr(Pb)clr(Mn)clr(Ce)clr(La)clr(Nd)clr(Ba)clr(Rb)clr(Cu)clr(Fe)clr(Sr)clr(Ga)clr(Zn)Figura 3. Biplot <strong>de</strong> los datos clr-transformados con 12 variables (Proporción explicada <strong><strong>de</strong>l</strong>a varianza: 52%).Tal como se aprecia en esta Figura 3 y en la Tabla 4 los vértices asociados a las partesLa y Ce están muy cercanos ⎯Var(ln(La/Ce))=0.079⎯ y, por lo tanto, sugiere que sonpartes aproximadamente proporcionales. En la región <strong>de</strong> investigación las tierras rarasno son contaminantes <strong><strong>de</strong>l</strong> medioambiente, no tienen origen tecnogénico, su presenciase pue<strong>de</strong> ciertamente <strong>de</strong>terminar cómo una causa mineralógica, es <strong>de</strong>cir laspropieda<strong>de</strong>s geo<strong>química</strong>s <strong>de</strong> los suelos condicionan su <strong>composición</strong> mineralógica ynada más. Los minerales más conocidos, que contienen La y Ce, son la bastnaesita[(Ce, La)(CO 3 )F] y la monacita [(Ca, Ce, La, Th…)PO 4 ]. En la presentación <strong><strong>de</strong>l</strong> este<strong>estudio</strong> se ha comentado la presencia <strong>de</strong> monacita en las rocas arenosas parentales<strong>de</strong> la región <strong>de</strong> investigación. El rayo <strong>de</strong> la tercera tierra rara, Nd, también aparece(Figura 3) muy cerca <strong>de</strong> los rayos <strong>de</strong> La y Ce. Las propieda<strong>de</strong>s geo<strong>química</strong>s <strong>de</strong> esteelemento son un poco distintas <strong>de</strong> las <strong>de</strong> Ce y La, en particular, su contenido en losminerales es esencialmente menor. A<strong>de</strong>más el comportamiento en procesosgeoquímicos <strong>de</strong> Ce y La, <strong>de</strong> un lado, y el Nd, <strong><strong>de</strong>l</strong> otro, es también bastante distinto. Lavariabilidad relativa entre el neodimio y la <strong>de</strong> los lantano y cerio (Tabla 4; 0.126 y21