versión impresa ISSN 0016-7169
Geofís. Intl v.48 n.4 México oct./dic. 2009
Paleomagnetic and rockmagnetic survey of Brunhes lava flows from Tancitaro volcano, Mexico
R. Maciel Peña1*, A. Goguitchaichvili1, V. H. Garduño Monroy2, V. C. Ruiz Martinez3, B. Aguilar Reyes1, J. Morales1 L. AlvaValdivia4, C. Caballero Miranda4 and J. UrrutiaFucugauchi4
1 Laboratorio Interinstitucional de Magnetismo Natural, Instituto de Geofísica, Sede Michoacán, Universidad Nacional Autónoma de México, Campus Morelia, Mexico. *Corresponding author: email@example.com.
2 Departamento de Geología y Mineralogía, Instituto de Investigaciones Metalúrgicas, Universidad Michoacana, San Nicolás de Hidalgo, Morelia, Michoacán, Mexico.
3 Departamento de Física de la Tierra, Astronomía y Astrofísica I, Universidad Complutense de Madrid, 28040, Madrid, Spain.
4 Laboratorio de Paleomagnétismo y Geofísica Nuclear, Instituto de Geofísica, Universidad Nacional Autónoma de México, Ciudad Universitaria, Del Coyoacán 04510 México City, Mexico.
Received: March 31, 2009
Accepted: May 13, 2009
Este estudio presenta los resultados paleomagnéticos del Volcán Tancitaro, (perteneciente al Campo Volcánico Michoacán Guanajuato), en el oeste de México, como una contribución a la base de datos de promediados en el tiempo del campo global. Se realizó el estudio paleomagnético y de magnetismo de roca en once flujos de lava. Se colectaron 120 núcleos orientados en el volcán Tancitaro y sus áreas aledañas. Todos los sitios fueron previamente fechados por 40Ar 39Ar (Ownby et al., 2007,) y van desde 793 ka hasta el presente. Se realizaron experimentos de magnetismo de roca, curvas de susceptibilidad e histéresis magnética y en la mayoría de los casos la remanencia magnética fue por la presencia de titanomagnetita pobre en Titanio con estructura magnética de pseudodominio simple. Ocho de los once flujos de lava dieron polaridad normal, mientras que los tres restantes produjeron paleodirecciones inconsistentes y no se tomaron en cuenta para el análisis y se debió probablemente por efectos de rayos. La paleodirección principal obtenida de los ocho flujos es Inc=39.5°, Dec= 356.4°, k=29, αc95=9.1° lo cual corresponde a la posición del polo con Plat=84.4°, Plong=219.9°, K=33 y αc95=8.5° Esta dirección es prácticamente consistente con la paleodirección esperada del Plio Cuaternario, derivada del polo de referencia para el Craton Norteamericano. La variación paleosecular coincide con otros estudios de la misma latitud y con modelos estadísticos recientes. La inclinación media cae dentro de la incertidumbre del Eje Dipolar Geomagnético más 5% de contribuciones cuadrupolares.
Palabras clave: Paleomagnetismo, variación paleosecular, promedio de tiempo de campo, epoca Brunhes, Cinturón Volcánico TransMexicano.
This study presents paleomagnetic results from Tancitaro volcanism in the Michoacan Guanajuato Volcanic Field, western Mexico, as a contribution to the time averaged field global database. Detailed paleomagnetic and rockmagnetic studies were carried out on eleven independent lava flows; 120 oriented, standard paleomagnetic cores were collected from Tancitaro volcano and surrounding areas. All sites were dated by means of 40Ar39Ar (Ownby et al., 2007) as originating from 793 ka to present. Rockmagnetic experiments included continuous susceptibility and hysteresis measurements. Remanence is carried mostly by Tipoor titanomagnetite of pseudosingledomain magnetic structure. Eight out of eleven flows yield normal magnetic polarities while three sites yielded inconsistent paleodirections most probably due to lightning. Mean paleodirection from eight flows is Inc=39.5°, Dec=356.4°, k=29, αc95=9.1° which corresponds to a pole position with Plat=84.4°, Plong=219.9°, K=33 and A95=8.5°, practically undistinguishable from expected PlioQuaternary paleodirections, for the North American Craton. Paleosecular variation is compatible with other studies at the same latitude bands and with recent statistical models. The mean inclination falls within the uncertainties of the Geomagnetic Axial Dipole plus 5% quadrupolar contributions.
Key words: Paleosecular variation, time averaged field, Brunhes chron, TransMexican Volcanic Belt.
Fluctuations of the geomagnetic field with time are essential for understanding the geodynamo, and conditions in the Earth's liquid core and at the coremantle boundary. Paleosecular variation may indicate modulation of geodynamo action in the core by the convective state of the lower mantle. Thus, determinations of these variations are decisive for understanding the processes in the core that give rise to the geomagnetic field and how and why the Earth's magnetic field reverses polarity.
The fundamental factors in the variability of the Earth Magnetic Field are the Time Averaged Field (TAF) and paleosecular variation (PSV). The TAF initiative has begun to update the database of geomagnetic observations over the last five million years (Mejia et al., 2005; Tauxe et al., 2004, Lawrence et al., 2006; Lawrence et al., 2009). A simple method to estimate the PSV consists in observing of the angular standard deviation (ASD) of virtual geomagnetic poles (VGPs) for a given locality. Several combinations of dipole and nondipole components predict the ASD characteristic of PSV with latitude (McFadden et al., 1988, 1991; Lawrence et al., 2006; Johnson et al. 2008)).
Johnson et al. (2008) reported a detailed synthesis of a new generation of paleomagnetic data compilations. The latitudinal dependence of VGP (virtual geomagnetic poles) scatter for these data appears much less important. The data at low latitudes seems to be more scattered than those at high latitude. This issue depends critically on a set of data of latitude 20°N. Published directional data from Trans Mexican Volcanic Belt (TMVB) meeting some very basic selection criteria (HerreroBervera et al., 1986; RuizMartinez et al., 2000; Osete et al., 2000, Morales et al., 2001, AlvaValdivia et al., 2001; Mejia et al., 2005; Petronille et al., 2005, RodriguezCeja et al., 2006) show an unusually high degree of scatter, which might be due to genuine geomagnetic field behavior.
The present study reports time averaged field global database and paleosecular investigations at low latitudes, based on a detailed rockmagnetic and paleomagnetic study of eleven lava flows associated to Tancitaro Volcano in the Michoacan Guanajuato Volcanic Field (MGVF). All sites were recently dated by means of Ar40Ar39 sy stematics (Ownby et al., 2007). The available ages range from 793 Ka to present.
Available age and sampling
The TransMexican Volcanic Belt contains stratovolcanoes, cinder cone fields and silicic caldera complexes. Volcanism in this region dates from late Miocene related to the subduction of Cocos and Rivera plates at the Middle America trench. The western TMVB is associated with subduction of ~ 9 Ma Rivera plate, whereas the 1218 Ma Cocos plate subducts under North America to the east (Klitgord and Mammerickx, 1982).
The MichoacanGuanajuato volcanic field is located in the westcentral part of TMVB (Fig. 1), east of Chapala Lake. The geographic boundaries of the MGVF are 18°45' N to 20° 15' N and 100°25' W to 100°45' W, with an area of 40,000 km2, It contains over 2,000 small monogenetic volcanoes including cinder cones (90%), maars, tuff rings, lava domes and lava flows with hidden vents all are predominant calcalkaline but some alkaline and transitional rocks are also found Silica content varies from 47% to 70% for olivine basalt and basaltandesite rocks (Hasenaka, 1994; Hasenaka et al., 1994).
Volcán Tancítaro (VT) is a large andesitic, composite volcano in the MichoacánGuanajuato Volcanic Field in westcentral Mexico. With a height of 3840 m (Ownby et al., 2007), VT is the dominant feature in the MGVF. Twentysix new Ar40Ar39 ages indicate that Tancítaro became active >793+22 ka and that the most recent effusive activity occurred at 237± 4 ka. Two catastrophic sectorcollapse events are identified and dated; the first one on the west flank between 694 and 571 ka, and the second on the east flank between 261 and 238 ka (Ownby et al., 2007).
Our sampling strategy was largely conditioned by Ownby et al. (2007) who gave 26 new ArAr incremental heating ages for Tancitaro and adjacent sorounding areas. We sampled only sites with available radiometric dating information (Table 1 and Fig. 2 trough 4), of easy access and yielding fresh, apparently unanltered outcrops. In total, 120 oriented samples belonging to 11 individual lava flows were collected. The samples were distributed throughout each flow both horizontally and vertically. All lava flows sampled were horizontal (dip less than 4°). In general, samples were obtained at the bottom of flows with the hope of collecting samples with the finest grained material. Cores were sampled with a gasolinepowered portable drill, and oriented in most cases with both magnetic and sun compasses.
Summary of magnetic experiments
In order to obtain the directions of characteristic remanent magnetization and to identify the principal magnetic carriers, following experiments were carried out: a) Thermal and alternating field demagnetizations, b) measurements of continuous thermomagnetic curves (low field susceptibility versus temperature), and c) hysteresis experiments.
The remanent magnetizations of seven to nine samples from each lava flow (Table 1) were measured with a JR6 (AGICO LtD) spinner magnetometer (nominal sensitivity ~ 109 Am2) at the Laboratorio Interinstitucional de Magnetismo Natural (LIMNA) in Morelia (Mexico). Both alternating field (AF) demagnetization (using a molspin AFdemagnetizer) and stepwise thermal demagnetization up to 560 °C (using a ASC TD48 furnace) were carried out. During thermal demagnetization, the lowfield susceptibility at room temperature was measured after each step with a Bartington susceptibility meter.
Twocomponent magnetizations were systematically recognized for most of studied units (Fig. 4). The secondary components are sometimes much stronger than primary ones (samples 08T012A, 08T012B and 08T093A). The characteristic magnetizations components are isolated after applying 40 mT peak alternating field. It should be noted that AF treatments proved to be more efficient than thermal demagnetization. This is illustrated at Fig. 4. Samples 08T012B and 08T012A belong to the same core. While thermal treatment is unable to isolate primary remanence, the alternating fields could reveal the primary, characteristic magnetization at last steps of demagnetization procedure. We believe that the origin of this strong secondary overprint is due to the lightning effect. This is in agreement with the fact that the NRM (natural remanent magnetization) directions show huge dispersion on the equalarea projection. We note that no ChRM directions were obtained from sites TAN1, TAN4 and TAN5.
For remaining sites, relatively small, secondary components, probably of viscous origin were detected and were easily removed applying 10 mT (Fig. 4, sample 08T107A). The greater part of remanent magnetization, in most cases was removed at temperatures between 520 and 560°C, which indicate to lowTi titanomagnetites as responsible for magnetization. The median destructive fields (MDF) range mostly from 20 to 40 mT, suggesting pseudosingle domain grains as remanent magnetization carriers (Dunlop and Özdemir, 1997).
A characteristic magnetization direction was determined by the least squares method (Kirschvink, 1980), 4 to 10 points being taken in the principal component analysis for this determination. The obtained directions were averaged by unit and the statistical parameters calculated assuming a Fisherian distribution.
Susceptibility vs. Temperature
Lowfield susceptibility measurements (kT curves) under air were carried out using AgicoKapabrdige susceptibility meter equipped with furnace in Saint Maur (France) IPGP laboratory. One sample from each site were heated up to about 630°C at a heating rate 20°C/ min and then cooled at the same rate. Curie temperature was determined by the Prévot et al's (1983) method. Alternatively, lowtemperature (from about 185°C to room temperature) susceptibility was recorded using the same apparatus.
LowT susceptibility experiments (Fig. 5) show a rather monotonic increase from about 185°C to room temperature with no indication of Verwey transition. Some titaniumpoor titanomagnetite may be responsible for remanent magnetization. As showed by Özdemir et al. (1993), the Verwey transition may be largely suppressed for the titanomagnetites with variable titanium content. Alternatively, similar behavior may also belong to nonstoichiometric (partially oxidized) magnetite. Corresponding highT susceptibility experiments (Fig. 6) also indicate the presence of Tipoor titanomagnetites. However, the cooling and heating curves are not perfectly reversible, probably because of low initial value of magnetic susceptibility. This may also be due to some moderate mineralogical alteration at high temperatures. Both Tirich and Tipoor titanomagnetites seem to coexist in few lava flows (samples 08T107 and 08T004A). These curve yields apparently two different thermomagnetic phases during heating. The lower Curie point ranges between 300400°C, while the highest one is about 570°C.
Magnetic hysteresis measurements were performed at room temperature on a specimen from all sampled sites at IPGP (Saint Maur) laboratory apparatus in fields up to 0.8 Tesla. The histeresis parameters were calculated after correction for the paramagnetic contribution. The coercivity of remanence (Hcr) was determined by applying progressively increasing backfield after saturation. Typical hysteresis plot are reported in Fig. 7a. The representative curve is simple, symmetrical and reflects very restricted ranges of the coercivities (Tauxe et al., 1996). Judging from the ratios of hysteresis parameters (Fig. 7b), it seems that all samples fall in the pseudosingle domain (PSD) grain size region (Day et al., 1977). This may also indicate a mixture of multidomain (MD) and a significant amount of single domain (SD) grains (Dunlop and Özdemir, 1997; Dunlop, 2002).
Results and discussion
Beside strong lightning effect, the average unit directions are rather precisely determined for 8 independent lava flows out of 11 collected (Table 1, Fig. 8a). All α95 are less than 10°. All flows yielded normal polarity magnetization as may be expected for the Bruhnes age rocks. We consider the paleodirections determined in this study to be of primary origin. Thermomagnetic curves show that the remanence is carried in most cases by Tipoor titanomagnetite, resulting of oxyexsolution of original titanomagnetite during the initial flow cooling, which indicates that the primary magnetization is a TRM (thermoremanent magnetization). Moreover, unblocking temperature spectra and relatively high coercivities point to pseudosingle domain magnetic structure grains as responsible for remanent magnetization. The mean paleodirection obtained from eight flows is Inc=39.5°, Dec=356.4°, k=29, α95=9. 1° which corresponds to the pole position Plat=84.4°, Plong=219.9°, K=33 and A95=8.5°. These directions are practically undistinguishable (Fig. 8a and b) from both the spin axis and the expected PlioQuaternary paleodirections, as derived from reference poles for the North American craton (Besse and Courtillot, 2002). This may indicates that no major regional tectonic rotation occurred in the area since about 1 My. The mean inclination overlaps within the uncertainties to those derived from the GAD (Geomagnetic Axial Dipole) and GAD plus a 5% quadrupolar contributions.
The formula was used for estimating paleosecular variation in this study where, ST is the total angular dispersion [Cox, 1969], N the number of sites used in the calculation, δi the angular distance of the ith virtual geomagnetic pole (VGP) from the axial dipole, SW the within site dispersion (following McEllhinny and McFadden, 1997) and, n the average number of sample per site. All new VGPs obtained on this study yield lesser colatitudes (maximum value is 28.3°) than generally adopted 45° cutoff angle (Johnson et al., 2008). We obtained SF=16.0 with SU=25.5 and SL=8.1 (upper and lower limits respectively) which reasonably agree to the model G of McFadden et al. (1988, 1991) fit to the McElhinny and McFadden (1997) and Johnson et al. (2008) databases for the last 5 Myr (Fig. 8b and 9).
The combination of our data with previously published results from Central and Western Mexico (Conte, 2004) yields that the amplitude of the secular variation is consistent with values obtained from other worldwide scattered sites. Thus, the hypothesis of the Pacific dipole window (Doell and Cox, 1971) may be rejected. This supports the findings of McElhinny et al. (1996) and RuizMartínez, (2004).Withintheuncertainties,the. Within the uncertainties, the results obtained in this study agree with the PSV values for Hawaii (+20° latitude) and Reunion (20° latitude). However, the amplitude of dispersion found is much higher in Mexico with respect to other places at the same latitude bands in agreement with findings of Lawrence et al., 2006. More high quality studies and reevaluation of old sites are needed in order to estimate whether this is a genuine characteristic of geomagnetic field.
This study was supported by UNAMDGAPA grant n° 102007 CONACYT grant n° 54957. We thank Bernard Henry and Maxim le Goff for help during rockmagnetic measurements in Saint Maur. VCRM is grateful to the financial support given by MEC (Spain) "José Castillejo Program" ref. JC200700314.
AlvaValdivia L. M, A. Goguitchaichvili and J. UrrutiaFucugauchi, 2001. Further Constraints for the PlioPleistocene Geomagnetic Field Strength: new results from Los Tuxtlas volcanic field (Mexico). Earth, Planets and Space 53, 9, 873881. [ Links ]
Besse, J. and V. Courtillot, 2002. Apparent and true polar wander and the geometry of the geomagnetic field over the last 200 Myr. J. Geophys. Res., 107, B11, 1029/2000JB000050. [ Links ]
Conte, G., 2004. Estudio integrado de rocas volcánicas PlioCuaternarias de México y Italia, PhD Thesis, National University of México, 154pp. [ Links ]
Cox, A., 1969. Confidence limits for the precision parameter k. Geophys. J. R. astr. Soc. 18, 545549. [ Links ]
Day, R., M. Fuller and V. A. Schmidt, 1977. Hysteresis properties of titanomagnetites: grain size and compositional dependence. Phys. Earth Planet. Int., 13, 260267. [ Links ]
DelgadoGranados, H., 1993. Late Cenozoic tectonics offshore western Mexico and its relation to the structure and volcanic activity in the western TransMexican volcanic belt. Geofísica Internacional 32, 543559. [ Links ]
Doell, R. and A. Cox, 1971. Pacific geomagnetic secular variation. Science, 71, 248254. [ Links ]
Dunlop, D. and Ö. Özdemir, 1997. Rock magnetism: fundamentals and frontiers, Cambridge University Press, 573 pp. [ Links ]
Dunlop, D., 2002. Theory and application of the Day plot (Mrs/Ms versus Hcr/Hc) 1. Theoretical curves and tests using titanomagnetite data. J. Geophys. Res., 107, 10.1029/2001JB000486. [ Links ]
Johnson, C. L., C. G. Constable, L. Tauxe, R. W. Barendregt, L. L. Brown, R. S. Coe, P. Layer, V. Mejia, N. D. Opdyke, B. S. Singer, H. Staudigel, and D. Stone,2008. Recent investigations of the 05 ma geomagnetic field recorded in lava flows. Geochemistry, Geophysics, Geosystems, 9:Q04032, doi:10.1029/2007GC001696. [ Links ]
Hasenaka, T., 1994. Size, distribution, and magma output rate for shield volcanoes of the MichoacánGuanajuato volcanic field, Central Mexico. Journal of Volcanology and Geothermal Research, 63, 1331. [ Links ]
Hasenaka, T., M. Ban and H. DelgadoGranados, 1994. Contrasting volcanism in the MichoacánGuanajuato Volcanic Field, central Mexico: Shield volcanoes vs. cinder cones. Geofísica Internacional, 33, 1, 125138. [ Links ]
HerreroBervera E., J. UrrutiaFucugauchi, A. Martin del Pozzo, H. Bohnel and J. Guerrero, 1986. Normal amplitude Brunhes paleosecular variation at lowlatitudes: A paleomagnetic record from the TransMexican Volcani Belt. Geophysical Research Letters, 13, 14421445. [ Links ]
Kirschvink, J. L. 1980. The leastsquare line and plane and analysis of palaeomagnetic data. Geophys. J. R. Astron. Soc., 62, 699718. [ Links ]
Klitgord, K. D. and J. Mammerickx, (1982). Northern East Pacific RiseMagnetic anomaly and bathymetric framework. J. Geophys. Res., 87, 67256783. [ Links ]
Lawrence, K. P., C. G. Constable, C. L. Johnson, 2006. Paleosecular variation and the average geomagnetic field at ±20° latitude. Geochem. Geophys. Geosyst., doi:10.1029/2004GC000871. [ Links ]
Lawrence, K. P., L. Tauxe, H. Staudigel, C. G. Constable, A. Koppers, W. McIntosh and C. L. Johnson, 2009. Paleomagnetic field properties near the southern hemisphere tangent cylinder. Geochem. Geophys. Geosys., 10, Q01005, doi: 10.1029/2008GC002072. [ Links ]
McEllhinny, M. W., P. L. McFadden and R. Merrill, 1996, The myth of the Pacific dipole window. Earth Planet Sci. Lett., 143,1322. [ Links ]
McElhinny, M. W. and P. L. McFadden, 1997. Palaeosecular variation over the past 5 Myr based on a new generalized database. Geophys. J. Int., 131, 240252. [ Links ]
McFadden, P., T. Merrill and W. McElhinny, 1988, Dipole/Quadrupole Family Modeling of Paleosecular Variation. Journal of Geophysical Research, 93, 11,58311,588. [ Links ]
McFadden, P. L., R. Merrill, M. W. McEllhinny and S. Lee, 1991. Reversals of the Earth's magnetic field and temporal variations of the dynamo families. J. Geophys. Res., 96, 39233933. [ Links ]
Mejia, V., H. Bohnel, M. A. OrtegaRivera, J. Lee and J. ArandaGomez, 2005. Paleosecular variation and timeaveraged field recorded in Late PlioceneHolocene lava flows from Mexico. Geochem. Geophys. Geosyst., 6, doi : 10.1029/2004GC000871. [ Links ]
Morales, J., A. Goguitchaichvili and J. UrrutiaFucugauchi, 2001. A rockmagnetic and paleointensity study of some Mexican volcanic lava flows during the Latest Pleistocene to the Holocene. Earth Planets Space, 53, 839902. [ Links ]
Osete, M. L., V. C. RuizMartinez, C. Caballero, C. Galindo, J. UrrutiaFucugauchi and D. Tarling, 2000, Southward migration of continental volcanic activity in the Sierra de Las Cruces, Mexico: Palaeomagnetic and radiometric evidence. Tectonophysics, 318, 201215. [ Links ]
Ozdemir, O., D. Dunlop and B. M. Moskowitz, 1993. The effect of oxidation on the Verwey transition in magnetite. Geophys. Res. Lett., 20, 16711674. [ Links ]
Ownby, S. E., H. DelgadoGranados, R. A. Lange and C. M. Hall, 2007. Volcan Tancitaro, Michoacan, Mexico, 40ar/39ar constraints on its history of sector collapse. J. Volcanology and Geothermal Res., 161, 114. [ Links ]
Petronille, M., A. Goguitchaichvili, B. Henry, L. AlvaValdivia, J. RosasElguera, M. Rodríguez Ceja and M. CalvoRathert, 2005. Paleomagnetism of ArAr dated lava flows from the CeborucoSan Pedro volcanic field (western Mexico): Evidence for the MatuyamaBrunhes transition precursor and a fully reversed geomagnetic event in the Brunhes chron. J. of Geophys. Res., 110, b08101, doi:10.1029/2004jb003321. [ Links ]
Prévot, M., E. A. Mankinen, S. Grommé and A. Lecaille, 1983. High paleointensities of the geomagnetic field from thermomagnetic study on rift valley pillow basalts from the midAtlantic ridge. J. Geophys. Res., 88,23162326. [ Links ]
RodriguezCeja, M., A. Goguitchaichvili, M. CalvoRathert, J. MoralesContreras, L. AlvaValdivia, J. R. Elguera, J. U. Fucugauchi and H. D. Granados, (2006). Paleomagnetism of Tequila Volcanic Field: Western Mexico. Earth Planets Space, 58,10, 13491358. [ Links ]
RuizMartinez, V. C., M. L. Osete, R. Vegas, J. I. Nunez Aguilar, J. UrrutiaFucugauchi and D. H. Tarling, 2000. Palaeomagnetism of Late Miocene to Quaternary volcanics from the eastern segment of the TransMexican Volcanic Belt. Tectonophysics, 318, 217233. [ Links ]
RuizMartínez, V. C., 2004. Estudio paleomagnético del cinturón volcánico transmexicano: implicaciones tectónicas, PhD thesis, Universidad Complutense de Madrid, ISBN 8466917616, 288 pp., Madrid. [ Links ]
Tauxe, L., T. A. T. Mullender and T. Pick, 1996. Potbellies, waspwaists and superparamagnetism in magnetic hysteresis. J. Geophys. Res., 95, 1233712350. [ Links ]
Tauxe, L., C. Luskin, P. Selkin, P. Gans and A. Calvert, 2004. Paleomagnetic results from the Snake River Plain: Contribution to the timeaveraged field global database. Geochem. Geophys. Geosyst., 5, Q08H13, doi:10.1029/2003GC000661. [ Links ]