Mourad Essalhi

Stanislas Sizaret
email: stanislas.sizaret@univ-orleans.fr

Luc Barbanson

Yan Chen

France Lagroix

François Demory

José M Nieto

Reinaldo Sáez

M Ángeles Capitán

A CASE STUDY OF THE INTERNAL STRUCTURES OF GOSSANS AND WEATHERING PROCESSES IN THE IBERIAN PYRITE BELT USING MAGNETIC FABRICS AND PALEOMAGNETIC DATING

Keywords: AMS, Iberian Pyrite Belt, gossans, massive sulphides, terraces, magnetic fabric, paleomagnetic dating

In the Rio Tinto district of the Iberian Pryrite Belt of South Spain, the weathering of massive sulphide bodies form iron caps, i.e. true gossans and their subsequent alteration and re-sedimentation has resulted in iron terraces, i.e. displaced gossans. To study the stucture and evolution of both types of gossans, magnetic investigations have been carried out with two foci: (i) the characterisation and spatial distribution of magnetic fabrics in different mineralised settings, including massive sulphides, gossans and terraces, and (ii) paleomagnetic dating. Hematite has been identified as the suceptibility carrier in all sites and magnetic fabric investigation of four gossans reveals a vertical variation from top to bottom, with: (i) a horizontal foliation refered to as "mature" fabric in the uppermost part of the primary gossans; (ii) highly inclined or vertical foliation interpreted as "immature" fabric between the uppermost and lowermost parts, and (iii) a vertical foliation interpreted to be inherited from Hercynian deformation in the lowermost part of the profiles. In terraces, a horizontal foliation dominates and is interpreted to be a "sedimentary" fabric.

Rock magnetic studies of gossan samples have identified goethite as the magnetic remanence carrier for the low temperature component (LTC), showing either a single direction close to the present Earth field (PEF) direction or random directions. Maghemite, hematite and occasionally magnetite are the remanence carriers for the stable high temperature component (HTC) that is characterized by non PEF directions with both normal and reversed magnetic polarities. No reliable conclusion can be yet be drawn on the timing of terrace magnetization due to the small number of samples. In gossans, the polarity is reversed in the upper part and normal in the lower part. This vertical distribution with a negative reversal test suggests remanence formation during two distinct periods.

Remanence in the upper parts of the gossans is older than in the lower parts, indicating that the alteration proceeded from top to bottom of the profiles. In the upper part, the older age and the horizontal "mature" fabric is interpreted to be a high maturation stage of massive sulphides' alteration. In the lower part, the age is younger and the inherited "imature" vertical Hercynian fabric indicates a weak maturation stage. These two distinct periods may reflect changes of paleoclimate, erosion and/or tectonic motion.

Introduction

Iron cap is an iron oxi-hydroxides-rich superficial material formed by oxidation, weathering or decomposition of former rocks. In metallogeny, a gossan is an iron cap formed by intensive weathering and oxidation of sulphides that is located on the upper part of sulphide mineralization [START_REF] Kosakevitch | Chapeaux de fer : problème de définition et de nomenclature pratique[END_REF][START_REF] Wilhelm | Chapeaux de fer[END_REF]1979). A classic gossan is composed mainly of neoformed iron oxi-hydroxides with residual unaltered minerals such as quartz. In addition to manganese and iron oxi-hydroxides, various precious metals and minerals, such as gold and silver and various sulfate and silicate minerals can occur in gossans. Commonly, gossans appear as a red stain to the country rocks and soil due to the abundance of oxidized iron and they may form elevated plateaus due to the abundance of erosion-resistant quartz and iron oxides. Gossans are considered also to be an important marker for the exploration of ore deposits [START_REF] Kosakevitch | Chapeaux de fer : problème de définition et de nomenclature pratique[END_REF][START_REF] Wilhelm | Chapeaux de fer[END_REF]Assiri and Moussa 2008). Previous studies on gossans are focused mostly on the geochemistry, environmental impacts and life in extreme environments (e.g. [START_REF] Boyle | Geochemistry and genesis of the Murray Brook precious metal gossan deposit, Bathurst Mining Camp[END_REF][START_REF] Scott | The geochemistry, mineralogy and maturity of gossans derived from volcanogenic Zn-Pb-Cu deposits of the eastern Lachlan Fold Belt, NSW, Australia[END_REF][START_REF] Al | Effects of acid-sulfate weathering and cyanide-containing gold tailings on the transport and fate of mercury and other metals in Gossan Creek: Murray Brook mine, New Brunswick, Canada[END_REF][START_REF] Egal | Iron isotopes in acid mine waters and iron-rich solids from the Tinto-Odiel Basin (Iberian Pyrite Belt, Southwest Spain)[END_REF]. Studies of gossan formation and evolution are less common. Moreover, paleomagnetic studies to date gossan formation, such as the Murray Brook and Heath Steele Au-Ag gossan deposits are rare [START_REF] Symons | Pliocene-Pleistocene genesis for the Murray Brook and Heath Steele Au-Ag gossan ore deposits, New Brunswick, from paleomagnetism[END_REF], and to our knowledge, no investigation on magnetic fabrics of a gossan has ever been performed.

However, these fabrics and their associated sulphides may help to understand the processes involved in gossan formation. Anisotropy of Magnetic Susceptibility (AMS) is the principal technique used in this work. AMS has proven to be an efficient technique to quantify the preferred orientation of magnetic fabrics with a high sensitivity and therefore is considered to be a powerful tool to study the texture of rocks (e.g. [START_REF] Rochette | Rock magnetism and the interpretation of the anisotropy of magnetic susceptibility[END_REF][START_REF] Tarling | The Magnetic Anisotropy of Rocks[END_REF][START_REF] Sizaret | Magnetic fabrics and fluid flow directions in hydrothermal systems. A case study in the Chaillac Ba-F-Fe deposits (France)[END_REF]. AMS has been applied in fields such as structural geology and rock petrofabrics (e.g. [START_REF] Hrouda | Magnetic anisotropy of rocks and its application in geology and geophysics[END_REF][START_REF] Tarling | The Magnetic Anisotropy of Rocks[END_REF], granitic pluton emplacement (e.g. [START_REF] Jover | Magnetic mineralogy of some granites from the French Massif Central: origin of their low-field susceptibility[END_REF][START_REF] Bouillin | Granite emplacement in an extensional setting: an AMS study of the magmatic structures of Monte Capanne (Elba, Italy)[END_REF][START_REF] Talbot | Pull-apart emplacement of the Margeride granitic complex (French Massif Central). Implications for the late evolution of the Variscan orogen[END_REF], magmatic flow direction (e.g. [START_REF] Nomade | Paleomagnetic study of French Guyana Early Jurassic dolerites: hypothesis of a multistage magmatic event[END_REF], rock deformation (e.g. [START_REF] Borradaile | The influence of deformation mechanisms on magnetic fabrics in weakly deformed rocks[END_REF][START_REF] Evans | Strain partitioning of deformation mechanisms in limestones: examining the relationship of strain and anisotropy of magnetic susceptibility (AMS)[END_REF], and recently to track paleoflow circulation (e.g. [START_REF] Sizaret | Anisotropie de susceptibilité magnétique (ASM) et chimie des traces : une nouvelle méthodologie pour démêler processus hydrothermaux et supergènes. Application au gisement à Ba-Fe-F de Chaillac (Indre, France)[END_REF]Sizaret et al. 2006a and b;[START_REF] Essalhi | Track of fluid paleocirculation in dolomite host rock at regional scale by the Anisotropy of Magnetic Susceptibility (AMS): An example from Aptian carbonates of La Florida, Northern Spain[END_REF][START_REF] Sizaret | Estimating the local paleo-fluid flow velocity: New textural method and application to metasomatism[END_REF]. Although the AMS method has already been applied to weathered soils (e.g. [START_REF] Borgne | Susceptibilité magnétique anormale du sol superficiel[END_REF][START_REF] Özdemir | A preliminary study of soil samples from west-central Minnesota[END_REF] and lateritic profiles (Cf. [START_REF] Mathé | Volumetric changes in weathered profiles: iso-element mass balance method questioned by magnetic fabric[END_REF][START_REF] Théveniaut | Paleomagnetism applied to lateritic profiles to assess saprolite and duricrust formation processes: the example of Mont Baduel profile (French Guiana)[END_REF], no AMS study has been devoted so far to gossans.

Actively explored since the Chalcolithic era, the Iberian Pyrite Belt (IPB) is one of the world's oldest and most important sulphide mining districts with original reserves estimated at > 1700 Mt [START_REF] Sáez | Geological constraints on massive sulphide genesis in the Iberian Pyrite Belt[END_REF]) and a mining history of > 4500 years [START_REF] Leblanc | 4.500-year-old mining pollution in southwestern Spain: long-term implications for modern mining pollution[END_REF][START_REF] Sáez | The extractive metallurgy of copper from Cabezo Juré, Huelva, Spain: Chemical and mineralogical study of slags dated to the Third Millenium B.C. Can[END_REF]. These deposits are located in an arc-shaped belt of 230 km in length and 40 km in width, extending from Seville in Spain to south of Lisbon in Portugal.

Several structural, petrological and geochemical studies have been published on the IPB (e.g. [START_REF] Leistel | The volcanic-hosted massive sulphide deposits of the Iberian Pyrite Belt : Review and preface to the Thematic Issue[END_REF][START_REF] Sáez | The Iberian type of volcanosedimentary massive sulphide deposits[END_REF]. The sulphide deposits were formed during the Devono-Carboniferous period and folded by the Variscan orogeny. Percolation of meteoritic water within the massive sulphides and associated stockworks induced the formation of gossans above the water table.

The process of gossan formation starts with the oxidation of pyrite and other sulphides to form Fe oxi-hydroxides plus sulfuric acid. The initial product of weathering is typically an amorphous ferric hydroxide expressed as Fe(OH) 3 that gradually crystallizes to goethite and/or hematite. Other important associated minerals with the Fe oxi-hydroxides are jarosite, gypsum, clay minerals and silica (e.g. [START_REF] Assiri | The geochemistry of gossans associated with Sarcheshmeh porphyry copper deposit, Rafsanjan, Kerman, Iran: Implications for exploration and the environment[END_REF].

The aim of the present study is to describe the magnetic fabric and to date terrace and gossan formation, in order to better understand the sulphide weathering processes in the Spanish part of the Iberian Pyrite Belt.

Geological setting of the Iberian Pyrite Belt

The IPB with the Baja-Acebuches ophiolitic complex and Pulo do Lobo antiform constitute the South Portuguese Zone (Fig. 1a). The IPB stratigraphic succession is composed of three main units that rest on an unknown substratum (Figs. 1a and1b). The lowest unit is the Late Devonian Phyllite-Quartzite (PQ) formation, which is about 1000 to 5000 m thick and composed mainly of alternating sandstone and shale. This formation is characterized by typical stable epicontinental platform features [START_REF] Moreno | Evidence for catastrophism at the Famennian-Dinantian boundary in the Iberian Pyrite Belt[END_REF]. The boundary between the PQ formation and the intermediate unit is defined by quartz-arenite layers [START_REF] Moreno | Evidence for catastrophism at the Famennian-Dinantian boundary in the Iberian Pyrite Belt[END_REF]. The middle unit is the Late Famennian to Late Visean Volcano-Sedimentary Complex (VSC; [START_REF] Van Den Boogaard | Conodonts of the Upper Devonian and Lower Carboniferous age from Southern Portugal[END_REF][START_REF] Oliveira | South Portuguese Zone: introduction. Stratigraphy and syn-sedimentary tectonism in the South Portuguese Zone[END_REF]) that is composed of three acidic volcanic cycles separated by two basic volcanic cycles [START_REF] Lécolle | Successions lithologiques et stratigraphiques dans la Province de Huelva Espagne; position des minéralisations magnésifères et pyriteuses[END_REF][START_REF] Routhier | La ceinture sud-iberique à amas sulfurés dans sa partie espagnole médiane[END_REF] with mudstone and some chemical sediments. The VSC varies in thickness between 100 and 1300 m. The IPB massive sulphide deposits are hosted in the VSC. The youngest unit in the IPB is the Culm Group [START_REF] Moreno | Postvolcanic Paleozoic of the Iberian Pyrite Belt: an example of basin morphologic control on sediment distribution in a turbidite basin[END_REF], which ranges from Late Visean to Middle-Late Pennsylvanian [START_REF] Oliveira | South Portuguese Zone: introduction. Stratigraphy and syn-sedimentary tectonism in the South Portuguese Zone[END_REF]). The Culm group is a detrital unit of synorogenic turbidites up to 3000 m thick. The whole succession of the IPB suffered slight regional metamorphism during the Hercynian Orogeny [START_REF] Munhá | Blue amphiboles, metamorphic regime and plate tectonic modelling in the Iberian Pyrite Belt[END_REF][START_REF] Munhá | Metamorphic evolution of the south Portuguese/Pulo do Lobo zone[END_REF][START_REF] Abad | The phyllosilicates in diagenetic-metamorphic rocks of the South Portuguese Zone, Southwestern Portugal[END_REF]. Deformation in the IPB has been extensively discussed and a thin-skinned tectonic model is widely accepted that describes the IPB as a southward-verging thrust and fold belt [START_REF] Schermerhorn | Folded overthrusts at Aljustrel[END_REF][START_REF] Schermerhorn | An outline stratigraphy of the Iberian Pyrite Belt[END_REF][START_REF] Ribeiro | Structure of the South Portuguese Zone[END_REF][START_REF] Silva | Structural outline of the south Portuguese zone[END_REF][START_REF] Soriano | Tectonica del cabalgamientos en la Faja Piritica Iberica (Zona Sur Portuguesa): la lamina de cabalgamientode Sanlucar de Guadiana y el antiform de Puebla de Guzman[END_REF][START_REF] Quesada | A reappraisal of the structure of the spanish segment of the Iberian Pyrite Belt[END_REF][START_REF] Onézime | Structural evolution of the southernmost segment of the West European Variscides: the South Portuguese Zone (SW Iberia)[END_REF].

Massive sulphide metallogeny in the IPB Three sulphide deposit styles coexist in the IPB: massive ore bodies, stockworks and disseminations. Devono-Carboniferous mineralization took place during the Variscan orogeny at about 350 Ma (Re-Os isotopic method; [START_REF] Mathur | Age and sources of the ore at Tharsis and Rio Tinto, Iberian Pyrite Belt, from Re-Os isotopes[END_REF][START_REF] Nieto | Estudio isotópico con el sistema Re-Os de las mineralizaciones de sulfuros de la Faja Pirítica Ibérica[END_REF]. Palynology has provided a similar age in several sectors of the IPB, such as Tharsis, Aznalcóllar, Sotiel-Coronada or Neves-Corvo [START_REF] González | Ore genesis age of the Tharsis Mining District (Iberian Pyrite Belt): a palynological approach[END_REF][START_REF] Sáez | Synchronous deposition of massive sulphide deposits in the Iberian Pyrite Belt: New data from Las Herrerías and La Torerera ore-bodies[END_REF], which suggest that a single mineralization event was responsible for the genesis of all massive sulphide deposits in the IPB.

The massive sulphides of the IPB are associated with volcanic (Río Tinto), volcanoclastic (La Zarza) or purely sedimentary rocks (Tharsis). These massive sulphides constitute an intermediate genetic type between the Clastic-dominated sediment hosted and the Volcanogenic Massive Sulphides (VMS) types. For this reason, metallogenists have given it a specific name: the Iberian Type [START_REF] Sáez | Geological constraints on massive sulphide genesis in the Iberian Pyrite Belt[END_REF][START_REF] Sáez | The Iberian type of volcanosedimentary massive sulphide deposits[END_REF].

The massive sulphides in the IPB are composed mainly of pyrite (all deposits) sometimes with chalcopyrite (e.g. Aznalcóllar, El Carmen), sphalerite (e.g. Tharsis, San Platón), galena (e.g. Los Silos, Neves-Corvo), as well as small quantities of pyrrhotite (e.g. Tharsis, Sotiel) (Marcoux et al. 1996). Accessory minerals are described by Marcoux and Leistel (1996).

The gossans in the IPB

In the IPB, two types of sulphide-related iron-oxided formations can be distinguished: (i) primary or "in place" gossans (true gossans) developed immediately over the sulphide bodies (Filón Sur in Tharsis, La Joya) and (ii) secondary or "remobilized" gossans (false gossans)

formed by subsequent supergene weathering, transportation and sedimentation (La Zarza, Alto de la Mesa). To avoid confusion, the latter type of gossan is called a "terrace" (e.g. [START_REF] Capitán | Caracterización textural y mineralógica del gossan del Filón Sur (Tharsis, Huelva)[END_REF][START_REF] Nieto | Beudantite: a natural sink for As and Pb in sulphide oxidation processes[END_REF]. [START_REF] Wilhelm | Chapeaux de fer[END_REF] and [START_REF] Viallefond | The massive sulphide deposits of the South Iberian Pyrite Province: geological setting and exploration criteria[END_REF] show that true gossans generally contain high amounts of Cu, Pb, Ag, As, Sb, Bi, Au and Sn, erratic values in Mo and Co, and low values in Mn, Ni, and B, whereas terraces are rather poor in metals. Several gossans in the IPB have been exploited for gold, such as thos at Río Tinto and in the Filón Sur of Tharsis. In the latter gold is enriched up to about ≈ 6 times with an ore content varying between 1 and 3 g/t [START_REF] Capitán | Caracterización textural y mineralógica del gossan del Filón Sur (Tharsis, Huelva)[END_REF][START_REF] Nieto | Beudantite: a natural sink for As and Pb in sulphide oxidation processes[END_REF]. A large number of IPB gossans have been well studied through mineralogical and textural analysis. According to [START_REF] Capitán | Caracterización textural y mineralógica del gossan del Filón Sur (Tharsis, Huelva)[END_REF], the Filón Sur gossan was formed in two main stages: first by oxidation of primary sulphide mineralization above the water table, and second by subsequent weathering. These authors described this gossan as dominated, from bottom to top, by jarosite, goethite and hematite. Hematite is considered to represent the final stage of oxidation.

Method and results

Sampling

Cores were drilled with a portable gasoline drill and oriented using both magnetic and sun compasses when possible to correct local magnetic field anomalies [START_REF] Tarling | The Magnetic Anisotropy of Rocks[END_REF]. A total of 214 standard oriented cores were collected from 27 sites with at least 4 samples per site in the RioTinto district of the IPB. Sampling was done of massive sulphides and true gossans at Corta Atalaya, San Miguel, Poderosa and Angostura and from terraces at Alto de la Mesa (Fig. 1c). The gossans and underlying massive pyrite were sampled along vertical profiles to compare the magnetic signals of the primary sulphides and their true gossans. The San Miguel gossan is well developped with massive iron oxide about 15 m thick. Locally the oxides form stalactites in the upper part of the gossan. At Poderosa sampling occurs in a pinched structure very regular from the top of the gossan to the bottom within the massive sulphide zone. In the Angostura pit, samples have been taken in a silicified pinched structure at the top of the gossan and in a massive body in the lower part. At Corta Atalya, the samples were taken from the stockwork veins and from close to the oxidation front. Finally, the terrace of Alto de la Mesa in Rio Tinto city was sampled in the fine oxihydroxide ribbons where the detrital fraction is low.

Mineralogy of the gossan and massive sulphide

The massive sulphides are mainly composed of pyrite with inclusions of pyrrhotite and chalcopyrite (Fig. 2a and2b). In the gossans the dominant minerals are iron oxy-hydroxides, mainly goethite and hematite, with ilmenite and quartz (Fig. 2c-2h). Spherolites of goethite crystallized from colloidal solution forms mamillated textures (Figs 2c and2d). These ribbons sometimes cover pyrite boxworks or ilmenite. The open spaces between the spherulites are filled by quartz, goethite and hematite fragments. Ilmenite is commonly altered to rutile (Fig. 2g).

Magnetic fabrics study (AMS)

Anisotropy of Magnetic Susceptibility

Magnetic susceptibility (K) is the capacity of a material to be magnetized in a given magnetic field i.e. K=M/H where M and H are the induced and applied magnetic fields, respectively.

The susceptibility is described as a symmetric second rank tensor, which can be depicted by an ellipsoid defined by three perpendicularly oriented axes K 1 , K 2 and K 3 (with K 1 ≥ K 2 ≥ K 3 ).

K 1 and K 3 are considered to be the magnetic lineation and the pole of the magnetic foliation, respectively. A number of parameters can be calculated to characterize the AMS ellipsoid, but the most widely used [START_REF] Tarling | The Magnetic Anisotropy of Rocks[END_REF] are: the mean magnetic susceptibility (Km), the anisotropy degree (P'), and the shape parameter (T) with T>0 for an oblate ellipsoid and T<0 for a prolate shape [START_REF] Jelínek | Statistical processing of anisotropy of magnetic susceptibility measured on groups of specimens[END_REF]. The definitions of Km, P' and T are given in Table 1.

Magnetic susceptibility carrier

Magnetic fabrics are usually dominated by ferromagnetic (s.l) phases that have higher magnetic susceptibility and are the main susceptibility carrier(s). All studied samples show a positive and high magnetic susceptibility with 95% of having magnetic susceptibilities exceeding 10 -4 [SI] (Fig. 3) and a clear dominance of ferromagnetic (s.l.) minerals with the probable participation of some paramagnetic ones [START_REF] Rochette | Rock magnetism and the interpretation of the anisotropy of magnetic susceptibility[END_REF][START_REF] Tarling | The Magnetic Anisotropy of Rocks[END_REF].

X-ray diffraction was performed on 12 representative samples from different parts of the studied profiles (gossans, terraces and massive sulphides). The gossans are mainly composed of goethite and hematite, while the massive sulphides contain quartz, chlorite and pyrite. It is worth noting that the diffraction spectra of the primary gossans and terraces are similar.

Curie temperature measurements were performed on 11 representative samples. The different paths for heating and cooling, especially in massive sulphides (Fig. 4b), indicate a mineralogical transformation during heating. In the gossans, four susceptibility drops record Curie temperatures at about 100°C, 300°C, 560°C and 670°C (Fig. 4a). Heating experiments up to 350°C show a cooling curve below the heating path and suggests the destabilization of maghemite. The successive drops could be interpreted as goethite, maghemite destabilization, magnetite (or titanomagnetite) and hematite Curie temperatures [START_REF] O'reilly | Rock and mineral magnetism[END_REF][START_REF] Dekkers | Grain-size dependence of the magnetic behavior of pyrrhotite during its low temperature transition at 34 K[END_REF]Dunlop and Özdemeir 1997). Some neoformed magnetite and hematite in the gossans appears after heating with the significant increase of the magnetic susceptibility.

Therefore, at least some of these results, reflect the transformation of preexisting minerals. In the massive sulphides, a mineralogical transformation is well evidenced by the appearance of magnetite and hematite on the cooling curve (Fig. 4b). During cooling neoformed pyrrhotite has been well defined and probably resulted from the thermal transformation of pyrite. This strong increase of the magnetic susceptibility during cooling at about 300°C is thus due to the abundant formation of hexagonal and monoclinic pyrrhotites. This observation is confirmed by X-ray diffraction measurements carried out on the same samples before and after heating.

Isothermal remanent magnetization (IRM) has been measured for 14 samples of gossan, massive sulphide and terrace samples from different localities. The samples are magnetized in an increasing direct magnetic field and the resulting IRM saturation is characteristic of the different ferromagnetic minerals. The samples from the gossans do not show any magnetic saturation at an applied magnetic field up to 1 Tesla (Fig. 4c), revealing high coercivity magnetic minerals, such as goethite and/or hematite (Dunlop and Özdemeir 1997). Samples of massive sulphides show a different behavior i.e. a rapid increase of the induced magnetic moment at weaker applied fields, but without reaching total saturation by 1 T (Fig. 4d). This IRM behavior corresponds to an iron sulphide such as pyrrhotite (O'Reilly 1984).

Hysteresis analysis has been performed on 11 representative samples to test the ferrimagnetic components. The loop from ferromagnetic (s.l.) minerals within the gossan (Fig. 4e) is open and shows clearly the coexistence of ferromagnetic and antiferromagnetic minerals (probably with paramagnetic ones also). This behavior indicates ferromagnetic minerals with a high coercive field (Hc) of 248.8 mT in the gossan. In the massive sulphides, the hysteresis loop shows either a paramagnetic or antiferromagnetic behavior (Fig. 4f). Hysteresis thus suggests the presence of goethite and/or hematite in the gossans, and rather paramagnetic or perhaps hexagonal pyrrhotite within the massive sulphides [START_REF] Dunlop | Rock magnetism : Fundamentals and Frontiers[END_REF].

From a magnetic susceptibility point of view, gossans contain goethite and hematite with possibly a trace of magnetite. Though its susceptibility is much weaker than that of magnetite, hematite is probably the principal magnetic fabric carrier because of its great abundance and its high magnetic susceptibility compared to goethite [START_REF] Rochette | Rock magnetism and the interpretation of the anisotropy of magnetic susceptibility[END_REF]). In the massive sulphides, the principal magnetic susceptibility carrier is pyrrhotite in spite of its low abundance.

Selected powdered samples were weighed and subjected to two low temperature magnetism experiments. First, at room temperature (300°K) saturation isothermal remanent magnetization (RT-SIRM) was induced in a 2.5 T field and then its temperature dependence monitored in zero field on cooling from 300°K to 10°K at a 5°K/min sweep rate. Secondly, at 10°K, the sample was subjected once again to a 2.5 T magnetic field acquiring a lowtemperature (LT) SIRM. The temperature dependence of the LT-SIRM was monitored in zero field on warming from 10° to 300°K at a 5°K/min sweep rate. The latter experiment is analogous to a zero-field cooled (ZFC) experiment after the necessary pre-treatment. The monitoring during cooling will be referenced as RT-SIRM hereafter. The example ZFC curve from a sample of massive sulphide shows a 77% decrease of its initial LT-SIRM by 35°K, then remains constant up to 115°K, and decreases once again by 20% between 115°K and 125°K (Fig. 5a). The initial remanence loss (Tp) is attributed to pyrrhotite as previously observed by [START_REF] Dekkers | Grain-size dependence of the magnetic behavior of pyrrhotite during its low temperature transition at 34 K[END_REF] and [START_REF] Rochette | Magnetic transition at 30-34 Kelvin in pyrrhotite: insight into a widespread occurrence of this mineral in rocks[END_REF]. The second remanence loss coincides with the Verwey transition (T V ) of magnetite (Özdemir et al. 1993) where Tv is an electronordering transition occurring in a mixed-valent system that results in an ordering of formal valence states in the low-temperature phase. T V is also observed over the same temperature range on cooling of the RT-SIRM and is equally sharp. Lastly, we observe that only 2% of the ZFC remanence survives thermal demagnetization to 300°K, while 38 % of the RT-SIRM remains following low-temperature demagnetization to 10°K. The sample from gossan (Fig. 5b) does not show any indication of a Verwey transition in either ZFC or RT-SIRM experiments, suggesting the absence of magnetite. It shows a sharp decrease up to 35°K of their ZFC remanence. The amplitude of this decrease is about 30 % of the initial ZFC and much less than the 77 % decrease observed for the massive sulphide sample. Beyond 35°K, the ZFC remanence decreases continuously to a low of 34%. The RT-SIRM increases continuously on cooling such that the remanence at 10°K is 85 % greater than the initial remanence. This low-temperature behavior of the RT-SIRM is reminiscent of goethite (e.g. [START_REF] Guyodo | Rock magnetic, chemical and bacterial community analysis of a modern soil from Nebraska[END_REF].

Magnetic fabrics

Principal susceptibility axes and other AMS parameters were calculated with the SUSAR program of AGICO. Mean tensors were calculated using [START_REF] Jelínek | Statistical processing of anisotropy of magnetic susceptibility measured on groups of specimens[END_REF] statistics; the mean orientations of the principal axes of the ellipsoid of magnetic susceptibility for each site are calculated with "bivariate statistics" (Henry and Le Goff 1995) and the stereograms were plotted using the software Anisoft 4.2 (Chadima and Jelínek 2008).

The AMS data and their associated principal parameters are given in Table 1. In the P'-T plot (anisotropy degree-shape parameter; Fig. 6), the majority of sites display a positive shape parameter, indicating a dominance of magnetic foliation. The degree of anisotropy (P') is mostly lower than 1.02 in the gossans and between 1.02 and 1.05 in the massive sulphides.

These low P' values of ≤ 1.05 indicate the absence of deformation in the studied profiles [START_REF] Tarling | The Magnetic Anisotropy of Rocks[END_REF].

The results of magnetic fabrics along vertical profiles sampled in the primary gossans and terraces are presented in Figures 7 and8, respectively. The profile measured from San Miguel (Figs. 7a-7d; 7h-7l and 7s) is the most representative of a primary gossan because it was rooted in massive sulphides and well developed at its top. Along the studied profiles, the general characteristics of magnetic fabrics are a follows. The poles of the magnetic foliation (K 3 ) give a well-defined average direction. The foliation orientation varies along the profile, being vertical or highly inclined in the massive sulphides (e.g. Figs. 7v-7x) and in the lower parts of the primary gossans (e.g. Figs. 7l, 7o, 7q and 7t), but horizontal in the upper parts of the primary gossans (e.g. Figs. 7a-7e and 7g). Some exceptions to this general trend were observed where the lower part of the gossan had a horizontal foliation (e.g. Fig. 7k). These exceptions are probably due to the heterogeneity within the gossans. In terraces, the horizontal magnetic foliation dominates the magnetic fabric (e.g. Figs. 8a-8d).

Paleomagnetic study

Thermal demagnetization was applied using a laboratory-built furnace (ISTO). The magnetic remanence was measured using JR5 spinner magnetometer (ISTO) to remove unwanted modern viscous remanent magnetization components from the magnetic remanence of the specimens. The magnetic remanence directions were isolated by principal component analysis [START_REF] Kirschvink | The least squares line and the analysis of paleoamgnetic data[END_REF]) and the mean directions were computed by spherical statistics [START_REF] Fisher | Dispersion of a sphere[END_REF]) using a paleomagnetic software packages of by [START_REF] Cogné | PaleoMac: A Macintosh™ application for treating paleomagnetic data and making plate reconstructions[END_REF] and R. Enkin (unpublished).

Progressive magnetic remanence cleaning was done in 16 steps from 20 to 700°C. Intensity drops in magnetic remanence intensity-temperature diagrams show the presence of goethite at about 120°C, maghemite (destabilization) at 320°C and a high quantity of hematite at 680°C

with rare magnetite at 585°C (Figs. 9a, 9c, 9e and 9g). Figures 9b,9d, 9f and 9h show representative samples in orthogonal projection on demagnetization (Zijderveld 1967). Two components are usually isolated: a lower temperature component (LTC) from 0 to150°C, and a higher temperature one (HTC) from about 250 to 680°C that normally decays toward the origin of an orthogonal diagram. The LTC is probably carried by goethite, while the HTC is mainly carried by hematite and pyrrhotite with, sometimes, trace magnetite. The LTCs often exhibit either a normal polarity close to the PEF direction (Figs. 9b, 9d and 9f) or a random direction (Fig. 9h) that probably records the recent geomagnetic field or drilling heating.

Nevertheless, the HTCs are stable and differ from the PEF. Both normal (Figs. 9b and9d) and reversed (Figs. 9f and9h) polarities for the remanence is indicative of a characteristic remanent magnetization.

As discussed above, the high temperature component is principally carried by hematite and pyrrhotite. The similar magnetic directions recorded by two different minerals indicate that their magnetizations are probably acquired at the same time, at least at a specimen level.

Figure 10 presents stereoprojections of magnetic directions obtained from the four primary gossans (Fig. 10a) and the terrace (Fig. 10b). The magnetic directions of each polarity are well grouped and fit a Fisherian distribution [START_REF] Fisher | Dispersion of a sphere[END_REF]. The reversal test of McFadden and McElinny (1990) was applied to each site, and the test is negative for all four gossans, but positive for the terrace (Table 2). In the gossans the reverse polarity is found in the upper part, as is well displayed in the San Miguel gossan where the magnetic fabric foliation is horizontal, and also in the Corta Atalaya oxidised stockwork in the same area with its horizontal foliation. Combining the reversed remanence directions of all gossan samples gives a direction with a declination of 188.2° and an inclination of -55.8°. The normal polarity samples are mainly associated with the lower part of the gossans where the magnetic foliation is partially overprinted, and its declination is 358° and the inclination 51.8°. In the terrace the normal and reversed polarities are mixed and the statistical calculation gives a global mean direction with a declination of 4.1° and an inclination of 62.3°.

Discussion

Magnetic fabrics in gossan from top to bottom AMS study of several gossan and terrace profiles identifies four different types of magnetic fabrics:

(i) In massive sulphides, the magnetic fabric is characterized by an E-W to NW-SE vertical foliation (e.g. Fig. 7v). This foliation is also observed in the country rock shales where two sites also give an E-W magnetic foliation (Fig. 7f) that probably corresponds to a tectonic structure inherited from the Hercynian orogenic phase (e.g. [START_REF] Quesada | A reappraisal of the structure of the spanish segment of the Iberian Pyrite Belt[END_REF]). Although massive sulphide samples have higher P' values of ≤ 1.05 than the gossans samples of ≤ 1.02 (Fig. 6 and Table 1), these values from the intensely deformed Paleozoic ore deposits are still much lower than those cited for highly deformed rock of >1.2 [START_REF] Tarling | The Magnetic Anisotropy of Rocks[END_REF]. This suggests either that the massive sulphides have not completely preserved the fabric acquired during deformation, or that they underwent some mineralogical transformations that modified their primary Hercynian deformation fabrics, although the E-W to NW-SE vertical magnetic foliation is still termed a "Hercynian" fabric.

(ii) In the uppermost parts of the primary gossans, the magnetic fabric is characterized by a well-defined horizontal magnetic foliation (e.g. Figs. 7a-7e and 7g) with a weak degree of anisotropy (Fig. 6 and Table 1). The horizontality of the magnetic fabric can be explained by compaction of the uppermost zones of the gossan due to leaching of labile elements by meteoric water. This transition between a vertical foliation with relatively stronger P' in sulphide protore to a horizontal foliation with weaker P' in weathered zones is also observed in the lateritic profiles [START_REF] Mathé | Volumetric changes in weathered profiles: iso-element mass balance method questioned by magnetic fabric[END_REF][START_REF] Théveniaut | Paleomagnetism applied to lateritic profiles to assess saprolite and duricrust formation processes: the example of Mont Baduel profile (French Guiana)[END_REF]. This typical fabric is called a "mature" fabric (see below).

(iii) Between the root of massive sulphides and the horizontal fabric at the top of the primary gossans, an intermediate transitional zone can be defined that is a combination of the two fabrics (e.g. Figs. 7h-7t). In this zone, the magnetic fabric is characterized by a vertical foliation similar to the massive sulphide fabric but with variable orientations and a very weak anisotropy degree P', as in the upper parts of the gossans. In this transitional zone, the magnetic response is a mixture of "Hercynian" and "mature" fabric that is termed an "immature" fabric (see below).

(iv) In terraces, the AMS data display a well defined horizontal magnetic foliation (e.g. Figs.

8a-8b) that reveals the internal stratification. Such classical "sedimentary fabric" is extensively described in the AMS literature [START_REF] Tarling | The Magnetic Anisotropy of Rocks[END_REF][START_REF] Charreau | Magnetostratigraphy and rock magnetism of the Neogene Kuitun He section (northwest China): implications for Late Cenozoic uplift of the Tianshan mountains[END_REF].

The variation of the anisotropy degree (Fig. 6 and Table 1) probably reflects variable degrees of sediment compaction.

The vertical zonality identified by the magnetic fabrics can be correlated to the mineralogical zoning defined by [START_REF] Capitán | Caracterización textural y mineralógica del gossan del Filón Sur (Tharsis, Huelva)[END_REF] and [START_REF] Nieto | Beudantite: a natural sink for As and Pb in sulphide oxidation processes[END_REF] in the Filón Sur of Tharsis gossan where: (i) the upper zones, dominated by hematite, would correspond to their so-called "hematitic gossan"; and, (ii) the lower parts, dominated by goethite and jarosite with scarce hematite, would correspond to their so-called "goethitic and jarositic gossans". The distinction found between the two magnetic fabrics in gossans is also observed in two sites sampled in oxidized stockworks that show the two different fabrics in their corresponding topographic positions. When the stockwork is exposed to intense alteration, it shows a horizontal magnetic foliation (Fig. 7e). Conversely, when it is exposed to a weak alteration, an overprinted magnetic foliation is observed (Fig. 7o). [START_REF] Capitán | Caracterización textural y mineralógica del gossan del Filón Sur (Tharsis, Huelva)[END_REF] and [START_REF] Nieto | Beudantite: a natural sink for As and Pb in sulphide oxidation processes[END_REF] showed that hematite is the typical mineral formed in the last stage of alteration. Thus a fairly good correlation can be drawn between the mineralogical and magnetic fabric zonalities that can be interpreted as the grade of maturation of a gossan. (i) The higher degree of maturation is observed in the upper part of the profile. It is characterized by a higher concentration in hematite, weak values of anisotropy degree (P'), and a horizontal magnetic foliation. The higher maturation grade is defined by a complete obliteration of the fabric inherited from the Hercynian orogeny and the resultant "mature fabric". (ii) The lower degree of maturation is observed in the deeper parts of the gossans. Its main features are (1) the dominance of goethite with respect to hematite, (2) relatively lower values of anisotropy degree P' than in the underlying massive sulphides, and (3) a strongly inclined magnetic foliation that is interpreted to be a relic of the Hercynian orogenic deformation. The lower maturation grade in the lower part of a gossan corresponds to an "immature fabric".

The boundary between the "Hercynian" and the "immature" fabrics corresponds to an oxidation front (Fig. 12), where the sulphides are oxidized to iron oxi-hydroxides. Also, according to [START_REF] Capitán | Caracterización textural y mineralógica del gossan del Filón Sur (Tharsis, Huelva)[END_REF], [START_REF] Nieto | Beudantite: a natural sink for As and Pb in sulphide oxidation processes[END_REF] and this study, the transition zone between "immature" and "mature" fabrics in a gossan corresponds to a dehydration and recrystallization front (Fig. 12).

A scenario of sulphide alteration, formation of the primary gossans and their later development of associated terraces can be proposed as follows. Massive sulphides, initially horizontal (Fig. 12a), were folded during the Hercynian orogeny and acquired an E-W vertical fabric (Fig. 12b). Later, alteration took place and the primary gossans were formed immediately on top of the massive sulphides above the water table. Based upon mineralogical and fabric considerations, two grades of alteration can be distinguished: a higher maturation grade in the upper part of the gossans (Fig. 12c) and a lower maturation grade in the lower part above the massive sulphides (Fig. 12d). The exposed parts of primary gossans may then have undergone chemical weathering and erosion. The fragments of gossan were transported and deposited with iron oxi-hydroxides in the form of terraces, showing a horizontal magnetic foliation (Fig. 12d) and other characteristics of sedimentary environments.

Age of the massive sulphides oxidation and chemical remanent magnetization acquisition

As described above, the high temperature components are principally carried by pyrrhotite and hematite with occasional traces of magnetite. These minerals exhibit a common, nonpresent Earth field direction as a stable magnetic remanence with both normal and reverse magnetic polarities. Therefore the isolated directions record a characteristic remanent magnetization for the gossans and terraces. To a first order, if the slight angular difference between two polarities of gossans is considered as negligible, two individual paleomagnetic poles are consequently calculated ( Reversal tests on the gossan's remanence directions data reveal a statistically significant negative response. Table 2 lists the paleomagnetic poles of both polarities (Fig. 11b). The slight angular difference between the two gossan poles suggests two distinct periods of magnetization. Close examination of the polarity distribution along the vertical profiles, suggests a polarity distribution that generally follows the vertical logic of the magnetic fabrics: i.e. the samples carrying the reversed polarity are located in the upper part in the mature part of the gossan and those wih normal polarity are in the lower immature part of the gossan except for some mixed polarities in the middle part of the Angostura profile. This vertical polarity distribution accords with the hypothesis that the magnetization has been acquired during two different periods. The paleomagnetic pole of normal polarity is close to the 5-10 Ma sector of the APWP, however, the cone of 95% confidence of the reversed polarity pole intersects the cones of 15-25 Ma poles of the APWP. Although this comparison remains aproximate due to the small number of studied samples and probably also to the shorter time of the remanence acquisition, the normal polarity remanence seems to present a younger age than that of the reversed remanence. This observation indicates that the alteration of the massive sulphides took place from top to depth and that the process may involved two main stages. that have been recognized and correlated around the World (e.g. [START_REF] Zachos | Trends, rhythms, and aberrations in global climate 65 Ma to present[END_REF][START_REF] De Man | Late Oligocene warming event in the southern North Sea Basin: benthic foraminifera as paleotemperature proxies[END_REF][START_REF] Mosbrugger | Cenozoic continental climatic evolution of Central Europe[END_REF][START_REF] Alegret | Effects of the Oligocene climatic events on the foraminiferal record from Fuente Caldera section (Spain, western Tethys)[END_REF]. In order to trace the paleoenvironmental evolution across the Oligocene epoch, [START_REF] Alegret | Effects of the Oligocene climatic events on the foraminiferal record from Fuente Caldera section (Spain, western Tethys)[END_REF] carried out a quantitative analysis of benthonic and planktonic foraminiferal assemblages from the Oligocene section of the Fuente Caldera (Southern Spain). According to these authors, the occurrence of some planktonic foraminiferal species is directly linked to Late Oligocene warming, which induced a major drop in sea level. This event could be the cause of a water table drop that provoked a massive sulphide oxidation (e.g. [START_REF] Alegret | Effects of the Oligocene climatic events on the foraminiferal record from Fuente Caldera section (Spain, western Tethys)[END_REF].

The second alteration stage in the gossans was probably initiated by a water table drop related to the Messinian Salinity Crisis. The crisis resulted from the drying up of the Mediterranean Sea between 5.96 and 5.33 Ma [START_REF] Krijgsman | Chronology, causes and progression of the Messinian salinity crisis[END_REF], and it is one of the most dramatic Cenozoic climatic events to have occurred on Earth [START_REF] Hsü | Late Miocene desiccation of the Mediterranean[END_REF]. Westward Eocene and Oligocene movement of the Alboran block into the gap between Iberia and northern Africa caused a slab roll-back of Tethys oceanic lithosphere and a closure of the Tethys marine corridor linking the Mediterranean Sea with the Atlantic Ocean [START_REF] Duggen | Magmatic evolution of the Alboran region: The role of subduction in forming the western Mediterranean and causing the Messinian Salinity Crisis[END_REF]). Isolation of the Mediterranean Sea resulted in its desiccation in latest Miocene times, causing the Messinian Salinity Crisis [START_REF] Duggen | Deep roots of the Messinian salinity crisis[END_REF]. The resulting Mediterranean Sea level fall would have caused a water table drop and renewed the alteration of the massive sulphides. Intense erosion has also been inferred to have occurred at this time [START_REF] Haq | Chronology of fluctuating sea levels since the Triassic[END_REF][START_REF] Kosakevitch | Contrôles climatique et geéomorphologique de la concentration de l'or dans les chapeaux de fer de Río Tinto (Province de Huelva, Espagne)[END_REF]).

The two-stage alteration identified in the vertical variation of the paleomagnetic signal can be correlated to the AMS observations of this study. The "mature fabric" corresponds to the older age estimated from the reversed magnetic directions and the age of the "immature" lower part of profiles is rather younger. This change in vertical fabric is accompanied by a weak increase in magnetic anisotropy that may be explained by the lower maturation grade of the alteration. Paleomagnetic and AMS data are coherent with sulphide oxidation and gossan formation in two stages. Besides these geophysical studies based on the rock magnetism, mineralogical investigations of a neighboring gossan (Filón Sur, Tharsis district) reached similar conclusions [START_REF] Capitán | Caracterización textural y mineralógica del gossan del Filón Sur (Tharsis, Huelva)[END_REF], Nieto et al. 2003). According to these authors, the Tharsis gossan was formed in two main stages; the first stage of which was associated with the oxidation of primary sulphide mineralization above the water table, while the second stage was related to the evolution of the gossan after complete oxidation of the superficial part of the sulphide deposit.

Gossan formation

Sulphide oxidation and evolution of gossans are known to progress from the top (oxidized zone) to the base (protore) with a well distinguished vertical zonation (e.g. Park and MacDiarmind 1964;[START_REF] Jébrak | Géologie des ressources minérales[END_REF]. According to our previous discussion on age the process of gossan formation may be described as follows. The oxidation of pyrite and other sulphides starts to form iron oxi-hydroxides and sulfuric acid. The initial product of weathering is typically an amorphous ferric hydroxide expressed as Fe(OH) 3 that gradually crystallizes to goethite and hematite. Other important minerals associated with the Fe oxides are jarosite, gypsum, clay minerals and silica (e.g. Atapour and Aftabi 2007). The oxidation of the massive sulphides starts at the top and progresses downward to the deeper zones as the oxygenated meteoric waters percolate downward. This process occurred in two main periods, during Late Oligocene to Early Miocene and during Late Miocene to Early Pliocene. AMS fabrics suggest that the transformation in the upper part of a gossan forms an iron cap with a "mature" magnetic fabric similar to a laterite, whereas in their lower part the immature AMS fabric is younger with a mixed magnetic fabric influenced by Hercynian deformation (Fig. 12c). An alteration phase then evolved in two stages; the first associated with the vertical magnetic foliation and the second with a horizontal foliation (Fig. 12).

The terraces are invariably characterized by a "sedimentary" fabric with mixed magnetic polarities. They may have formed by erosion/alteration, transportation and deposition of primary gossan in a continuous process throughout the two stages of massive sulphide oxidation (Fig. 12).

Conclusions

This study of gossans and their associated massive sulphides inthe Spanish section of the Iberian Pyrite Belt (IPB) using Anisotropy of Magnetic Susceptibility (AMS) characterizes the magnetic fabrics in different mineralized settings and distinguishes different maturation grades within a primary gossan. Vertical zoning has been evidenced by three types of magnetic fabrics with intrinsic geological or/and mineralogical characteristics. These magnetic fabrics are: (i) a "Hercynian fabric" found in the massive sulphides and inherited from the Hercynian orogeny (Fig. 12b); (ii) a "mature fabric" observed in the uppermost part of primary gossans, produced by weathering involving meteoric water, and having a higher maturation grade (Fig. 12c); and (iii) an "immature fabric" found between the above two fabrics in the lower part of primary gossans (Fig. 12d). Terraces formed by alteration/erosion, transportation and deposition of gossans are characterized by a "sedimentary" fabric (Fig. 12d). The results of this magnetic investigation agree with those of earlier mineralogical studies [START_REF] Capitán | Caracterización textural y mineralógica del gossan del Filón Sur (Tharsis, Huelva)[END_REF][START_REF] Nieto | Beudantite: a natural sink for As and Pb in sulphide oxidation processes[END_REF].

As the values for the anisotropy parameter P' are relatively greater in the high maturation zone than in the immature zone, the anisotropy parameter could possibly be an estimate of the grade of maturation in different parts of a primary gossan. Further work is needed to test the reliability of P' as a useful parameter for mapping the grade of maturation and its potential as a mining exploration tool in gossans.

Dating gossans in the Iberian Pyrite Belt reveals a two stages formation. The paleomagnetic signal of gossans is vertically distributed in two zones: (i) a reversed remanence polarity in the upper part and a polarity one in the lower part. Combining this vertical distribution of magnetic polarities with a negative reversal test, it is reasonable to hypothesize that the magnetizations of each polarity were acquired during two distinct periods with an older alteration age at the top and a younger age at the bottom. This observation indicates that the alteration processed from the top of the gossans to depth as for climatic lateritization. These two paleomagnetic dating can be easily correlated to the AMS results. The older age correlates to the horizontal "mature" fabric in the upper part of the gossan and presents a high maturation stage for alteration of the massive sulphides. The younger age correlates to the vertical "immature" fabric in the lower part and indicates a more recent and weaker maturation stage. These two distinct periods may reflect the changes of paleoclimate and/or erosion inducing a water table drop.

Although this comprehensive magnetic study on gossans has proved its efficiency, more detailed work is needed to improve our knowledge of the alteration process in massive sulphides and the formation of gossans. Further studies should focus on the timing relationship between gossans and terraces, and on the possible spatial variation of alteration age with paleoclimatic or local tectonic changes. Moreover, some questions remain that are more difficult to answer, such as, what was the duration of alteration and magnetization acquisition? [START_REF] Henry | Application de l'extension bivariate de la statistique Fisher aux données d'anisotropie de susceptibilité magnétique: intégration des incertitudes de mesure sur l'orientation des directions principales[END_REF]. , and correspond to K 1 , K 2 and K 3 , respectively. Magnetic foliation is defined by the K 1 -K 2 plane, and the magnetic lineation corresponds to the K 1 axis. Ellipses around poles correspond to the interval of confidence at the 95% level. Lithology of sites is given in Table 1.

( ) 3
3 2 1 K K K K m + + =
, D: declination, I: inclination, K 1 , K 2 , and K 3 :

the axes of the maximum, intermediate and minimum of the ellipsoid of magnetic susceptibility, respectively.

( ) Jelinek, 1978). Values with star (*) are not retained for our interpretation, and double star (**) are calculations requiring more samples to be reliable. relative vertical profiles, from the upper part (a to f)and lower part (g to t) of the gossans and the massive sulfides (u to y) profiles. Mean orientations of the principal axes were calculated with bivariate statistics [START_REF] Henry | Application de l'extension bivariate de la statistique Fisher aux données d'anisotropie de susceptibilité magnétique: intégration des incertitudes de mesure sur l'orientation des directions principales[END_REF]. , and correspond to K 1 , K 2 and K 3 , respectively. Magnetic foliation is defined by the K 1 -K 2 plane, and the magnetic lineation corresponds to the K 1 axis. Ellipses around poles correspond to the interval of confidence at the 95% level. Lithology of sites is given in Table 1. , T = [2(ln K 2 -ln K 3 ) (ln K 1 -ln K 3 )] -1 (Jelinek, 1978). precision parameter; α95: 95% confidence circle around the mean direction (pole); P Lat. and P Lon .: latitude and longitude of pole, respectively, calculated at 37.72°N and 6.64°E.

2 1 3 1 2 ln( 2 exp ' ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = ∑ = i moy i K K P , T = [2(lnK 2 -lnK 3 ) (lnK 1 -lnK 3 )] -1 (

Figures

Values with star (*) are not retained for our interpretation, and double star (**) are calculations requiring more samples to be reliable.

Figure captions :

 captions Figure captions:

Figure 2 .

 2 Figure 2. Textures of massive sulphides (a-b) and associated gossans (c-h): bo: boxwork,

Figure 3 .

 3 Figure 3. Interval of bulk magnetic susceptibility of all measured specimens.

Figure 4 .

 4 Figure 4. Magnetic mineral analyses of representative samples from a gossan (left) and a

Figure 5 .

 5 Figure 5. RT-SIRM cooling branch (filled circles) and ZFC (open circles) curves for (a) a

Figure 6 .

 6 Figure 6. Graphic presentation of the shape parameter T versus the corrected anisotropy

Figure 7 .

 7 Figure 7. Stereographic projections in lower hemisphere of AMS responses from the San

Figure 8 .

 8 Figure 8. Magnetic fabric measured in the terrace located at Alto de la Mesa in Rio Tinto

Figure 9 .

 9 Figure 9. Relative remanent magnetization intensity vs. temperature (left) and orthogonal

Figure 10 .

 10 Figure 10. Equal-area projection of specimen remanence directions in stratigraphic

Figure 11 .

 11 Figure 11. Paleomagnetic poles calculated from: all gossan samples (star), normal-polarity

Figure 12 .

 12 Figure 12. Sketch of the evolution of the Iberian pyrite belt gossans from their deposition

Figure 1

 1 Figure 1. a. Simplified geological map of the major structural units and tectonostratigraphic

Figure 2 .

 2 Figure 2. Textures of massive sulfides (a-b) and associated gossans (c-h): bo: boxwork, cpy:

Figure 3 .

 3 Figure 3. Interval of bulk magnetic susceptibility of all measured specimens.

Figure 4 .

 4 Figure 4. Magnetic mineral analyses of representative samples from a gossan (left) and a

Figure 5 .

 5 Figure 5. RT-SIRM cooling branch (filled circles) and ZFC (open circles) curves for (a) a

Figure 6 .

 6 Figure 6. Graphic presentation of the shape parameter T versus the corrected anisotropy

Figure 7 .

 7 Figure 7. Stereographic projections in lower hemisphere of AMS responses from the San

Figure 8 .

 8 Figure 8. Magnetic fabric measured in the terrace located at Alto de la Mesa in Rio Tinto city;

Figure 9 .

 9 Figure 9. Relative remanent magnetization intensity vs. temperature (left) and orthogonal

Figure 10 .

 10 Figure 10. Equal-area projection of specimen remanence directions in stratigraphic

Figure 11 .

 11 Figure 11. Paleomagnetic poles calculated from: all gossan samples (star), normal-polarity

Figure 12 .

 12 Figure 12. Sketch of the evolution of the Iberian pyrite belt gossans from their deposition until

  declination, I: inclination, K 1 , K 2 , and K 3 : the axes of the maximum, intermediate and minimum of the ellipsoid of magnetic susceptibility,

Table 1

 1

	): Lat. = 86.6 °, Long. = 157.7°, α 95 = 2.2° for

comparison is almost impossible as the paleomagnetic poles from this microplate are too dispersed (see the details in

[START_REF] Schott | Palaeomagnetism of Permo-Triassic red beds in the western Pyrenees: evidence for strong clockwise rotations of the Palaeozoic unit[END_REF][START_REF] Osete | The Late Carboniferous to Late Triassic segment of the apparent polar wander path of Iberia[END_REF]

. This leads us to compare our newly obtained results with Stable Europe paleopoles as the kinematic history between Iberia and Stable Europe has been well established (e.g.

[START_REF] Osete | The Late Carboniferous to Late Triassic segment of the apparent polar wander path of Iberia[END_REF][START_REF] Rosenbaum | Relative motions of Africa, Iberia and Europe during Alpine orogeny[END_REF][START_REF] Symons | Dating of the Reocin MVT Deposit, Spain, by Paleomagnetism[END_REF]

). The comparison with Stable Europe

[START_REF] Besse | Apparent and true polar wander and the geometry of the geomagnetic field over the last 200 Myr[END_REF]

can be done directly (Fig.

11a

) without any kinematic pre-reconstruction because a Miocene age has been attributed to the gossans by paleontology (Phillips 1881;

[START_REF] Moreno | Edad mínima del gossan de Las Cruces: Implicaciones sobre la edad del inicio de los ecosistemas extremos en la Faja Pirítica Ibérica[END_REF]

, and moreover, the Iberian microplate has been solidly accreted to Stable Europe since 45 Ma

[START_REF] Symons | Dating of the Reocin MVT Deposit, Spain, by Paleomagnetism[END_REF]

. As only 10 samples were sampled from the terrace, no reliable conclusion could be made though its paleomagnetic pole seems closer to the present Earth's field. It is noteworthy that the positive reversal test may indicate that the terrace deposits may have been formed in a continuous process. The pole of the gossans falls near the 5-15 Ma sector of the apparent polar wander path (APWP) for stable Europe. This first estimate is consistent with previous age estimates (Phillips 1881;

[START_REF] Moreno | Edad mínima del gossan de Las Cruces: Implicaciones sobre la edad del inicio de los ecosistemas extremos en la Faja Pirítica Ibérica[END_REF]

  Two stage alteration of the massive sulphidesCompared to the APWP of Stable Europe[START_REF] Besse | Apparent and true polar wander and the geometry of the geomagnetic field over the last 200 Myr[END_REF], the ages of the two paleomagnetic poles obtained from this study may correspond to a Late Oligocene to Early Miocene age for the older zones in the upper part of the gossans and a Late Miocene to Early Pliocene age for the younger zones in the lower part of the gossans. This two-stage alteration is probably due to water table changes in response to climatic variation in erosion rates and/or in uplift tectonics.The first oxidation step occurred in the late Oligocene. This period is marked by large

and abrupt climate changes characterized by strong short-term fluctuations of warmer (e.g. the Late Oligocene Warming Event) and colder intervals (e.g. the Oligocene Glacial Maximum)

Table captions Table 1 .

 captions1 AMS data and associated parameters from this study. PG: primary gossan, MS: massive sulphides, OS: oxidized stockworks, SR: surrounding rocks, Te: terrace, n:

	number of measured samples, P': corrected anisotropy degree, T: shape parameter (T>0
	for oblate shape, and T<0 for prolate ones, Jelínek 1978), K m : bulk magnetic
	susceptibility,

Table 2 .

 2 Mean paleomagnetic directions and corresponding poles for studied gossans and

terrace of the Río Tinto District; n: number of samples; Dec.: declination; Inc.: inclination; k: precision parameter; α95: 95% confidence circle around the mean direction (pole); P Lat. and P Lon .: latitude and longitude of pole, respectively, calculated at 37.72°N and 6.64°E.

Table captions

 captions Table1. AMS data and associated parameters from this study. PG: primary gossan, MS: massive sulfides, OS: oxidized stockworks, SR: surrounding rocks, Te: terrace, n: number of measured samples, P': corrected anisotropy degree, T: shape parameter (T>0 for oblate shape, and T<0 for prolate ones, Jelínek 1978), K m : bulk magnetic susceptibility,

	MS	14	7x 6	925	1.019 -0.695 288.7 36.7 121.8 52.5	23.5	6.3
	Depth Fig. n Alto de la Messa (Gossan) 37°44'55'', 6°39'13'' Kv P' T	K1		K2	K3
	Te		(m) 0.3	8a 11	x10 -6 [SI] 1200	D (°) 1.002 0.894 309.5 11.5 218.4 I (°) D (°)	I (°) 5.5 103.4 77.2 D (°) I (°)
	Te	San Miguel (Gossan), 37°45'42'', 6°44'58'' 1.3 8b 15 963 1.005 0.559 297.8	4.8 207.1	9.1	55.6 79.7
	PG Te	3 2	7a 7 8c 10	1820 1200	1,019 1.003 0.792 0.475 268.3 15.3 21.8 14.1 287.3 17.4 148.6 67.3 3 16.9 138.5 66.9
	PG	3.5 7b 6	699	1.005	0.515	92.6	27.2 194.6 22.2 318.3 53.7
	PG	4	7c 8	1290	1.007	0.398	12.3	2.4 102.6	7.3 263.8 82.3
	PG	6	7d 8	1210	1.006 -0.102 107.8 24.7 213.7 30.8 346.4 48.7
	PG	7	7h 5	596	1.007	0.314 218.1 18.6 126.1	6.1	18.7 70.3
	PG	8	7i	7	1660	1.011	0.432 136.3	2.2 227.5 29.7	42.4 60.2
	PG	9	7j	7	549	1.047 -0.383 310.1 66.4 187.2 13.3	92.5 19.1
	PG	12	7k 7	1100	1.006 -0.331 111.4 36.7 290.1 53.3	21	0.6
	PG	14	7l	8	457	1.003	0.815 163.7	3.3 254	3.9	33.6 84.9
	MS	16	7t	9	362	1.048	0.575 111.2	5.6	16.9 37.8 208.3 51.7
	Corta Atalaya (Gossan) 37°42'12'', 6°35'51''
	PG	10	7m 7	5960	1.001 -0.526 188.2 19.7	96.7	4	355.8 69.9
	PG	10.2 7n 8	374	1.001	0.398 161.5	2.3	71.2	7	269.5 82.6
	MS	10.7 7o 9	986	1.039	0.958	93.7 37.3 316.4	44	202.4 22.8
	OS	11.4 7u 3	124	1.018	0.653	17.2 57.1 146.8	22.4 246.8 22.8
	OS		7e 4	355	1.003	0.635 311.8	5.9 220.7	10.4	70.8 78
	SR			7f	7	872	1.002 -0.128 213.3 75.6	73.5	11.3 341.6	9.2
	Poderosa (Gossan) 37°44'55'', 6°39'13''
	PG	3	7g 7	7050	1.006 -0.161	81	3.1	173	33.2 346.6 56.6
	PG	3.5	7p 10	465	1.005	0.189 306.5 12.5 107.8 76.8 215.6	4.1
	MS	4	7v 12	13	1.038 -0.561	6.3	72.2 258.6	5.6 166.9 16.8
	Angostura (Gossan) 37°46'01'', 6°43'11''
	PG	3	7q 11	439	1.003 -0.069	22.9 34.2 223.3 54.1 119.6	9.7
	PG	5	7r 11	530	1.008	0.881 298.1 35.3 176.4 36.6	56.5 33.8
	PG	9	7s 6	407	1.004	0.237 203.4 74.8	20.8 15.2 110.9	0.7
	PG	13	7w 5	433	1.005 -0.016	6.9 54.6 116	13.1 214.5 32.2

Table 2 .

 2 Mean paleomagnetic directions and corresponding poles for studied gossans and terrace of the Río Tinto District; n: number of samples; Dec.: declination; Inc.: inclination; k:

	Zone	n	Dec.	Inc.	k	α95 PLat.	PLon. α95
			(°)	(°)		(°)	(°N)	(°E)	(°)
	San Miguel (Gossan)	22*	2.4	49.2	34.7	5.2	83.1	167.5	4.9
	N	12	359.5	44.7	28.6	8.3	79.6	190.1	7.0
	R	10	186.5 -54.5	71.7	5.7	84.0	115.0	6.3
	Corta Atalaya (Gossan)	9*	7.7	53.4	122.1	4.7	82.6	122.9	5.8
	N	4	0.4	50.2	377.6	4.7	83.3	182.7	5.8
	R	5	194.2 -55.6	162.1	6.0	78.5	199.0	7.1
	Poderosa (Gossan)	6	9.3	57.5	415.4	3.3	82.5	90.7	4.3
	Angostura (Gossan)	18*	352.7 55.9	47.7	5.1	84.3	275.9	5.9
	N	16	352.1 55.2	44.8	5.6	83.5	269.0	6.4
	R	2** 180.1 -62.2	146.8 20.8	83.8	6.9 32.6
	Alto de la Mesa
	(Terrace)	9	4.1	62.3	102.1	5.1	82.7	30.8	6.5
	N	5*	7.0	63.6	113.8	7.2	80.5	36.1	9.7
	R	4*	180.8 -60.6	80.4 10.3	85.3	17.2 13.5
	All Gossan	55*	1.0	53.2	43.7	2.9	86.8	169.5	3.1
	Normal polarity	38	358.0 51.8	38.8	3.8	85.5	208.1	8.9
	Reversed polarity	17	188.2 -55.8	84.7	3.9	83.2	101.0	4.5

Acknowledgements

We present our thanks to David Symons, Fernando Tornos, Patrick Wiliams and anonymous reviewer for their constructive suggestions to improve the manuscript.