(U-Th)/He Dating of Supergene Iron (Oxyhydr-)Oxides of the Nefza-Sejnane District (Tunisia): New Insights into Mineralization and Mammalian Biostratigraphy

: The mining district of Nefza-Sejnane (Tunisia) encloses numerous ores and raw material deposits, all formed in relation with successive Fe-rich ﬂuids of meteoric and/or hydrothermal origins. Here, for the ﬁrst time in Tunisia, (U-Th)/He ages were obtained on supergene goethite from various localities/deposits of the district highlight direct dating of signiﬁcant weathering episodes during late Tortonian and late Pleistocene. These weathering events are most likely associated with favorable conditions that combine (i) wet climate displaying sufﬁcient meteoric water/ﬂuid; and (ii) regional exhumation, due to large-scale vertical lithospheric movements enhancing the percolation of ﬂuids. Matched with previous works, these results reﬁne the stratigraphic frame for the polymetallic mineralization and clay deposits in the district, conﬁrming the inﬂuence of meteoric ﬂuids circulation during the late Cenozoic. As a consequence of the new (U-Th)/He data, we moreover propose a taxonomic and stratigraphic revision of the well-known mammalian fauna from the Fe-rich Douahria locality, suggesting an early Tortonian age for the fossils, i.e., prior to the ﬁrst episode of meteoric event in the area. and discussing their biostratigraphic outputs, in order to further discuss GHe ages and paleontological dating.


Introduction
The mining district of Nefza-Sejnane contains numerous ores and raw material deposits, including the Zn-Pb mineralization of Sidi Driss, the vertebrate-bearing Douahria basin that hosts a SEDEX deposit [1], the halloysite-bearing iron deposit of Tamra Formation [2][3][4][5], the Fe-LREE-U enriched breccia of Oued Belif [5], and kaolinitic clays of the late Oligocene Numidian Flysch Formation [6,7]; (Figures 1 and 2). From the Miocene onwards, these deposits experienced successive circulations of Fe-rich fluids of mixed meteoric and hydrothermal origins. The hydrothermal activity resulted from a relatively high geothermal gradient and Miocene magmatic rocks of the area [8]. Using geochemical geothermal gradient and Miocene magmatic rocks of the area [8]. Using geochemical a guments, Decrée et al. (2014) demonstrated that emplacement of magmatic rocks in th Nefza district has enhanced hydrothermal fluid circulations, leading to the deposition o polymetallic mineralization [9]. In addition, significant periods of weathering are con firmed by (i) pedogenesis on sediments within small basins such as in Tamra basin an the deposition of the subaerial Tamra Formation [2]; and (ii) late supergene fluids lead ing to the formation of botryoidal or euhedral goethite [4,5]. The nature and role of m teoric circulations in metal enrichment of the area has already been stressed by Decrée al. [1,2] and discussed by Dekoninck et al. [4] using δD and δ 18 O on halloysite-kaolini and Fe (oxyhydr-)oxides, and by Decrée et al. [10] using iron isotopes and chemic composition of Fe-(oxi)hydroxides. However, ages of these meteoric influences in th area are still lacking, jeopardizing the regional integration of the ore formation. Fu thermore, the Douahria Fe-deposit contains the famous "Pikermian" vertebrate faun including a rhinocerotid, a giraffid, and an anthracotheriid originally described by Ro man and Solignac [11] in a ferruginous crust. Here we provide, for the first time in Tun sia, new (U-Th)/He ages on goethite (GHe ages) associated with weathering processes i the vertebrate-bearing Nefza-Sejnane polymetallic mining district. Matched with prev ous data, these new investigations allow us to discuss (i) the influence of meteoric fluid in the formation of ore deposits in this important mining district of Tunisia, and (ii) the ro of geodynamic and climatic parameters in weathering processes at a regional scale. More ver, we propose a global revision of the mammalian assemblage, aiming at providing upda ed taxonomic assignments for the concerned remains and discussing their biostratigraph outputs, in order to further discuss GHe ages and paleontological dating.  [4,[12][13][14].  [4,[12][13][14]. Minerals 2021, 11, x FOR PEER REVIEW 3 of 19 Figure 2. Simplified geologic map of the Nefza-Sejnane district, Tunisia (modified after reference [5], showing the location and (U-Th)/He corrected ages of the supergene goethite. G-CS = Ghardimaou Cap Serrat sinistral fault zone.

Geological Setting
The Nefza-Sejnane mining district is located in the Tellian "Nappe Zone" of Northern Tunisia (Figures 1 and 2). The district is characterized by the occurrence of late Cenozoic mineral deposits, overlying a substratum comprising folded marls Albian to late Priabonian in age, and Chattian-Burdigalian Numidian Flysch Formation [15], which are present as thrust sheets in the Tellian "Nappe Zone" [16]. The regional felsic subvolcanic rocks were emplaced between 12.9 ± 0.5 Ma and 8.2 ± 0.4 Ma (Serravallian to Tortonian; e.g., [8]), while basaltic flows of the Zouara basin ( Figure 2) were dated between 8.4 ± 0.4 Ma and 6.4 ± 0.2 Ma [17]. The Ghardimaou-Cap Serrat (G-CS) sinistral fault zone, a major regional structure guiding Alpine tectonics at the tectonic plate scale, and the regional N80°E lineaments controlled the magmatism and subsequent mineralization (Figures 1 and 2; [14]). The first notable mineralization event in the area was associated with a major brecciation episode that led to the formation of the ring-shaped Oued Belif LREE-U hematite breccia [5]. The origin of this breccia is most likely related to the mixing of magmatic-hydrothermal fluids and basinal (evaporite) brines. A K-Fe alteration event, slightly predating the main mineralization, is dated at 9.2 ± 0.3 Ma using K-Ar method on K-feldspars [5]. The main Fe-LREE-U mineralization is coeval with a major episode of brecciation, most likely associated with the emplacement of the Oued Belif rhyolite intrusion. The remaining voids are filled with goethite, which, in some areas, is mixed with kaolinite and Ba-Mn oxide. The two SEDEX deposits (Sidi Driss and Douahria; Figure 2) are formed through thermally-driven fluid circulations related to magmatism and reactivation of shear zones during late Miocene extension in small basins (Decrée et al., 2016 [14]). The Pb-Zn ore of Sidi Driss is correlated with the stratified Fe-ore of Douahria, which yielded a large mammalian assemblage, late Miocene in age (originally spelled Douaria; e.g., [18][19][20]). The Tamra basin, unconformably overlying the Sidi Driss basin, is the host for iron mineralization in the form of a 50 m thick succession (Tamra Formation) of iron hydroxide-stained sediments together with kaolinitehalloysite lenses [1,3,4]. Previous datings on hollandite [Ba(Mn 4+ ,Mn 2+ )8O16] and coro-

Geological Setting
The Nefza-Sejnane mining district is located in the Tellian "Nappe Zone" of Northern Tunisia (Figures 1 and 2). The district is characterized by the occurrence of late Cenozoic mineral deposits, overlying a substratum comprising folded marls Albian to late Priabonian in age, and Chattian-Burdigalian Numidian Flysch Formation [15], which are present as thrust sheets in the Tellian "Nappe Zone" [16]. The regional felsic sub-volcanic rocks were emplaced between 12.9 ± 0.5 Ma and 8.2 ± 0.4 Ma (Serravallian to Tortonian; e.g., [8]), while basaltic flows of the Zouara basin ( Figure 2) were dated between 8.4 ± 0.4 Ma and 6.4 ± 0.2 Ma [17]. The Ghardimaou-Cap Serrat (G-CS) sinistral fault zone, a major regional structure guiding Alpine tectonics at the tectonic plate scale, and the regional N80 • E lineaments controlled the magmatism and subsequent mineralization (Figures 1 and 2; [14]). The first notable mineralization event in the area was associated with a major brecciation episode that led to the formation of the ring-shaped Oued Belif LREE-U hematite breccia [5]. The origin of this breccia is most likely related to the mixing of magmatic-hydrothermal fluids and basinal (evaporite) brines. A K-Fe alteration event, slightly predating the main mineralization, is dated at 9.2 ± 0.3 Ma using K-Ar method on K-feldspars [5]. The main Fe-LREE-U mineralization is coeval with a major episode of brecciation, most likely associated with the emplacement of the Oued Belif rhyolite intrusion. The remaining voids are filled with goethite, which, in some areas, is mixed with kaolinite and Ba-Mn oxide. The two SEDEX deposits (Sidi Driss and Douahria; Figure 2) are formed through thermally-driven fluid circulations related to magmatism and reactivation of shear zones during late Miocene extension in small basins (Decrée et al., 2016 [14]). The Pb-Zn ore of Sidi Driss is correlated with the stratified Fe-ore of Douahria, which yielded a large mammalian assemblage, late Miocene in age (originally spelled Douaria; e.g., [18][19][20]). The Tamra basin, unconformably overlying the Sidi Driss basin, is the host for iron mineralization in the form of a 50 m thick succession (Tamra Formation) of iron hydroxide-stained sediments together with kaolinite-halloysite lenses [1,3,4]. Previous datings on hollandite [Ba(Mn 4+ ,Mn 2+ ) 8 [2]. The dated K-Mn oxides relate to mineralizing events closely connected with hydrothermal circulation and leaching of underlying mineralization of the Sidi Driss Pb-Zn SEDEX type deposit [2]. At Tamra, the δD and δ 18 O of halloysite-kaolinite are not in equilibrium with the potential Miocene-Pliocene meteoric fluids, confirming that they could be related to fluid-rock interactions with the underlying ores, marls, and/or skarns [4]. Late mineralization of galena, pyrite, and siderite in sub-vertical fractures cut the whole sequence of the district and are associated with a new hydrothermal event, still active in the area with numerous hot springs from 35 • C up to 70 • C and regional thermal gradients up to 100 • C/km [8]. Recent circulation of hydrothermal fluids is moreover confirmed by K-Ar dating of adularia at c. 0 Ma in the Oued Belif breccia [5].

Materials and Methods
Samples of iron (oxyhydr-)oxides were collected in the Tamra Fe-mine, the Boukhchiba Fe-mine, the Oued Belif Fe-LREE-U breccia, and the Ouchtata area-Zouara basin ( Figure 2; Table 1) during field campaigns in 2013 and 2014. Sample 14MnTAM was collected in the upper part of the Tamra mine, within the manganiferous zone (see field pictures and log in reference [2]. It contains pure, well-crystallized goethite with botryoidal texture ( Figure 3A,B) that formed during the pedogenesis affecting the Tamra Formation [4]. The 14TAM07 sample, also collected in the Tamra mine (see reference [2]), shows a mixture of goethite and hematite, with potential influences of both weathering and hydrothermal origins, as attested by the occurrence of hematite and Mn(Pb)-oxides ( Figure 3C,D). The 13NEFZA sample, also collected in the Tamra mine, corresponds to Fe-indurations of the intermediate part of the sequence defined by Decrée et al. [2]. This sample contains hematite and displays hydrothermal influences according to δD and δ 18 O data from Dekoninck et al. [4] (their samples CP04-3 and CP04-35 are very similar to our sample 13NEFZA). The 14BK11 sample comes from a mineralized vein in the Boukhchiba Fe-mine ( Figure 2), clearly cross-cutting the Fe-LREE-U-rich breccia described in Decrée et al. [5] and field pictures herein). The sample contains very pure, well-crystallized, and homogeneous goethite ( Figure 3E,F), which precipitated in equilibrium with meteoric waters [4]. The 14OU03 and the 13OULI samples were collected in the Ouchtata area ( Figure 2). They consist of several generations of goethite pseudomorphosing grains of pyrite of the Triassic sediments. "Ouchtata" sample in reference [4], similar to 13OULI sample, was associated with meteoric fluids by the latter authors. Sample 14OB02 was collected in a mineralized vein of the Oued Belif locality. It contains goethite, mixed with euhedral quartz and Kfeldspars ( Figure 3G,H). The goethite of this sample was considered as supergene (formed below the water table) by Decrée et al. [5]. The isotopic data of this sample are consistently away from the meteoric equilibrium water line [4], due to mixing with quartz and Kfeldspars. Ferruginous crust from the mammal-bearing Douahria site was also sampled in order to perform (U-Th)/He dating. The latter, however, contains large amounts of mixed hydrothermal Fe-Mn phases, and could thus not be analyzed with the (U-Th)/He method. * nccSTP = nano cubic-centimeters-at-standard-temperature-and-pressure; ** error at 1 sigma; *** eU = effective uranium. The mineralogy of samples was characterized using an X-ray Panalytical X'P (Malvern Panalytical Ltd., Palaiseau, Essone, France) and a diffractometer PW3710 (Philips, Amsterdam, Holland, the Netherlands) (CuKα radiation), and 7500 (Jeol Ltd., Akashima, Tokyo, Japan) scanning electron microscope (SEM) cou an energy dispersive electron spectrometer (EDS). The mineralogy of samples was characterized using an X-ray Panalytical X'Pert Pro (Malvern Panalytical Ltd., Palaiseau, Essone, France) and a diffractometer Philips PW3710 (Philips, Amsterdam, Holland, the Netherlands) (CuKα radiation), and a JEOL 7500 (Jeol Ltd., Akashima, Tokyo, Japan) scanning electron microscope (SEM) coupled to an energy dispersive electron spectrometer (EDS).
The (U-Th)/He dating method relies on the 4 He production, ejection, and accumulation inside the crystal structure during alpha decay of the radioactive U-Th chains [21]. The (U-Th)/He age is determined by measuring the total He content of the samples and in a second step, the total U and Th content. (U-Th)/He ages were obtained in the GEOPS laboratory (University Paris-Saclay, Orsay, France) for three aliquots per sample, except for 14BK11 where only two aliquots were analyzed. For samples composed of several iron (oxyhydr-)oxide generations, the purest zones were selected to avoid erroneous and meaningless ages. First, fragments manually extracted from the hand-sized samples were cleaned using an ultrasonic bath and rinsed with milli-Q water and ethanol. Aliquots were handpicked under a binocular microscope, weighed, and encapsulated into niobium envelopes (purity 99.9%). Fragments with a length exceeding 0.5 mm and a weight around 0.5 mg for 14OU03 sample and ranging between 6 and 18 mg for other samples were selected. Each encapsulated sample was degassed using a diode Ytterbium laser under vacuum for thirty minutes at temperatures below 1000 • C. This procedure was repeated until all 4 He was degassed. The 4 He gas was mixed in the purification line with a known amount of 3 He, and purified. Helium isotopes ( 3 He and 4 He) were measured with a Pfeiffer Prisma Quadrupole mass spectrometer to ensure the purity of the analytical gas (see reference [22] for more details). After degassing, samples were extruded from the envelopes directly into 5 mL Teflon-capped vials for complete dissolution. Firstly, 100 µL of 5 M HNO 3 containing a known amount of 230 Th and 235 U was introduced in each vial followed by 400 µL of 30% HCl, and a few drops of 38% HF. Vials were closed and heated up to 100 • C overnight. Solutions were evaporated at 100 • C to concentrate the sample. If samples were not completely dissolved, the procedure of adding HCl and HF was repeated. Secondly, 5 mL of 5 M HNO 3 was added and heated for one hour at 100 • C. The solution was diluted with 1 M HNO 3 to set the iron content at 100 ppm. The U and Th measurements were undertaken by using a quadrupole ICP-QMS series II CCT Thermo-Electron at LSCE (Laboratoire des Sciences du Climat et de l'Environnement-Gif sur Yvette, France). The analytical error for He and U-Th measurements is 5% at 1σ.
The fossil specimens illustrated here, collected by local miners in 1934 at the Douahria Fe deposit, are stored in the Paleontological Collection of the University Lyon 1-Villeurbanne (FSL). Table 1 shows the (U-Th)/He data obtained on goethite collected in the Nefza-Sejnane district. Most of the prepared samples are pure goethite. Two samples (14TAM07 and 13NEFZA) are made of a mixture of goethite with other phases such as hematite and a Mn-Pb oxide (see Figure 3D). Another sample (14OU02) is composed of a mixture of goethite, quartz, and K-feldspar ( Figure 3H). All samples present relatively homogeneous, effective uranium (eU = U + 0.24Th) contents ranging from 0.5 to 7 ppm. In addition, most samples have quite similar Th/U ratios around 0.01-0.2. Only one sample is characterized by a significantly higher Th/U ratio of 1.7 ± 0.1 (14TAM07, Table 1). The variability in the Th/U ratio, which may be used to trace differences in the iron oxide populations or the presence of inclusions [23,24], confirms the relatively high content of hematite in sample 14TAM07. The raw (U-Th)/He ages range from 0.5 ± 0.03 to 8.1 ± 0.4 Ma and are mostly homogeneous per sample. No correction of alpha ejection has been considered here, because the mean alpha stopping distance of the ejected He in iron oxides and hydroxides ranges from 14 to 16 µm [25]. These values are small compared to the sample's size (~0.5 mm). On the other hand, as He can be lost by diffusion due to the polycrystalline nature of the samples (see reference [26] for a review), He losses are corrected by applying a 10% correction on the He age and a 10% error (Table 1) (2014) applied a 5% to 30% correction for goethite (U-Th)/He ages based on diffusion data measured directly on the dated crystals [23,27]. In this study, we used a 10% correction on the raw ages and applied an additional error of 10%, allowing to estimate He losses by diffusion, which represents the sum of the analytical error and the error on the He diffusion loss estimation, similarly to reference [22] (Table 1). This correction will not significantly change the age interpretation but will be used in the following discussion since it better reflects crystallization ages. The reproducibility of He ages obtained for Tamra and Boukchiba goethite is better than the 10% error estimated to account for He diffusion processes, which gives good confidence in the present data. However, a 10% error is preferred to take into account a potential He diffusion. With this He correction, the oldest episode, at about 8.5 Ma, is recorded in Tamra and Boukchiba. The latest episode, younger than ca. 1 Ma, is recorded in Ouchtata, Oued Belif, and Tamra (Table 1). In summary, most of the present (U-Th)/He ages are reproducible (Table 1), highlighting two main episodes of supergene goethite formation, with mean ages at 8.6 ± 0.9 Ma and 0.8 ± 0.2 Ma (Table 1 and Figure 4). 10%, allowing to estimate He losses by diffusion, which represents the sum of the analytical error and the error on the He diffusion loss estimation, similarly to reference [22] ( Table 1). This correction will not significantly change the age interpretation but will be used in the following discussion since it better reflects crystallization ages. The reproducibility of He ages obtained for Tamra and Boukchiba goethite is better than the 10% error estimated to account for He diffusion processes, which gives good confidence in the present data. However, a 10% error is preferred to take into account a potential He diffusion. With this He correction, the oldest episode, at about 8.5 Ma, is recorded in Tamra and Boukchiba. The latest episode, younger than ca. 1 Ma, is recorded in Ouchtata, Oued Belif, and Tamra (Table 1). In summary, most of the present (U-Th)/He ages are reproducible (Table 1), highlighting two main episodes of supergene goethite formation, with mean ages at 8.6 ± 0.9 Ma and 0.8 ± 0.2 Ma (Table 1 and Figure 4).

(U-Th)/He Age Significance Using Petrographic Examination of Goethite
Most of the studied samples show pure and well-crystallized supergene goethite, allowing us to obtain, for the first time in Tunisia, direct robust dating of two episodes of major meteoric fluid circulation at 8.6 ± 0.9 Ma and 0.8 ± 0.2 Ma. The ages between these two main episodes reflect the mixture of goethite with other mineral phases. The 14TAM07 sample clearly shows a mixture of goethite and hematite, with several phases of crystallization and potential contamination by Mn-Fe phases ( Figure 3D). This mixture most likely explains the large variation observed in (U-Th)/He ages, from 3.3 ± 0.3

(U-Th)/He Age Significance Using Petrographic Examination of Goethite
Most of the studied samples show pure and well-crystallized supergene goethite, allowing us to obtain, for the first time in Tunisia, direct robust dating of two episodes of major meteoric fluid circulation at 8.6 ± 0.9 Ma and 0.8 ± 0.2 Ma. The ages between these two main episodes reflect the mixture of goethite with other mineral phases. The 14TAM07 sample clearly shows a mixture of goethite and hematite, with several phases of crystallization and potential contamination by Mn-Fe phases ( Figure 3D). This mixture Minerals 2021, 11, 260 9 of 19 most likely explains the large variation observed in (U-Th)/He ages, from 3.3 ± 0.3 Ma to 6.3 ± 0.6 Ma, in between the two main episodes (~8.6 and~0.8 Ma). A similar interpretation is proposed for the 13NEFZA sample, with (U-Th)/He ages ranging from 1.0 ± 0.1 Ma to 2.6 ± 0.3 Ma. The sample 14OB02, which is made of goethite, quartz, and K-feldspar, provides (U-Th)/He ages (0.6 ± 0.1 to 1.1 ± 0.1 Ma) that overlap with the other samples. This suggests that (U-Th)/He method is not affected by potential mixing of goethite with quartz and K-feldspars, the latter typically having a relatively poor content in U, Th, Sm, He compared to goethite.

Age of the Meteoric Fe-Fluids Circulation in the Nefza-Sejnane District
The new (U-Th)/He data highlight two episodes of Fe-rich meteoric fluids circulation in the Nefza-Sejnane district, at 8.6 ± 0.9 Ma and 0.8 ± 0.2 Ma. 5.2.1. Late Tortonian Weathering Event (8.6 ± 0.9 Ma) The oldest (U-Th)/He ages (mean age of 8.6 ± 0.9 Ma) are obtained on goethite samples from the Tamra and Boukchiba localities. The Tamra Formation, deposited during a regional extensional tectonic regime [13], consists of successive sedimentary sequences [2]. Each sedimentary sequence results from (1) sedimentation of iron-rich inherited material, and (2) in situ pedogenesis and iron concentration without evidence of a precursor hydrothermal event. This clearly confirms that Fe-(oxi)hydroxide formation resulted (at least partly) from the involvement of meteoric fluids in the area. In this context, extensional regime has provided the geological drains (fractures, faults, . . . ), allowing the fluids to percolate within sediments, as it has been previously suggested for numerous weathering profiles in various areas (e.g., [28,29]). The 14MnTAM sample, showing GHe ages between 8.9 ± 0.9 Ma and 8.4 ± 0.8 Ma (Table 1), is clearly related to weathering, as suggested by (i) its botryoidal texture ( Figure 3A) and (ii) the δD and δ 18 O data highlighting its formation in equilibrium with meteoric waters [4]. This relatively old age may also be associated with pedogenesis recognized during the deposition of the sequence [2]. This episode of weathering is likewise recognized at Boukchiba, where the GHe ages of 14BK11 range from 8.7 ± 0.9 Ma to 8.2 ± 0.8 Ma, which are consistent with those of the sample 14MnTAM. Sediments at Boukchiba also experienced weathering, as attested by the δD and δ 18 O data obtained onto the same sample of goethite. These data suggest that goethite formed in equilibrium with meteoric waters during late Tortonian period, but was not affected by synsedimentary pedogenesis of the subaerial Tamra Formation [5]. The sample 14BK11 ( Figure 3E,F), which is made of pure goethite, is located in a vein cross-cutting the Oued Belif breccia and associated to Fe-LREE-U mineralization. Then, the late Tortonian weathering period is recorded in the synsedimentary pedogenesis of the Tamra Formation (sample 14MnTAM) and in the surrounding brecciated rocks, suggesting at least a local event of weathering.

Late Pleistocene Weathering Event
The 14OU03, 13OULI, and 14OB02 samples, located in various localities in the district (see Figure 2 and Table 1), present consistent GHe ages, ranging from 1.0 ± 0.1 Ma to 0.6 ± 0.1 Ma. These ages highlight meteoric Fe-rich fluid circulation during a late episode around 0.8 ± 0.2 Ma in the Nefza-Sejnane district.

Stratigraphic Frame of the Nefza-Sejnane District
These new (U-Th)/He ages allow refining the stratigraphic frame of the Nefza-Sejnane district, especially from Tortonian to Pleistocene times ( Figure 5).  The entire nappe pile that forms a large part of the regional substratum is cut by shallow-level, Serravallian to lower Tortonian intrusive rocks. The felsic rocks of the whole district are dated between 12.9 ± 0.2 Ma and 8.2 ± 0.4 Ma (see a detailed review in reference [14]). In the Zouara basin (Figure 2), the transitional basalt sills exposed at Mokta el-Hadid have been dated at 8.4 ± 0.4 Ma [30], and at 6.9 ± 0.3 and 6.4 ± 0.2 Ma [17]. Some of these ages are discussed/questioned by Decrée et al., 2016 [14]. Anyway, the widespread middle-late Miocene bimodal magmatism is largely recognized at the scale of the whole district.
The age of the Tamra Formation was not precisely known but was supposedly Miocene-Pliocene or younger according to the 1:50,000 geological map of Nefza [16,17]. Our new dating refines the age of the Tamra succession. The deposition of the Tamra sediments is globally coeval with the oldest GHe age associated with the synsedimentary pedogenesis, at 8.6 ± 0.9 Ma. This suggests that the Tamra Formation is late Tortonian in age. At the time, meteoric waters may have played a role in the formation of halloysite and kaolinite and/or destabilization of primary clays in the Tamra ore. However, these meteoric fluids would have been mixed at some stage with deep hot saline fluids related to a thermally-driven circulation [4]. It confirms that the Tamra Formation experienced polyphased and mixed hydrothermal/meteoric fluid circulations, as already stressed by Decrée et al., 2008b [2].
The Tamra iron-stained sediments overlay the Sidi Driss sediments ( Figure 5). The latter can be correlated with the Douahria sediments, based on (i) regional lithological similarities; (ii) the presence of volcanoclastic intercalations in the Douahria and Sidi Driss sequences, including cinerites [1,9]; and (iii) their similar gastropod fauna [31,32]. The Sidi Driss and Douahria sediments were attributed to the late Miocene (Messinian?) by reference [1]. Our data suggest a minimum Tortonian age for these sediments, close to or at the end of the magmatic activity in the area ( Figure 5). This is in good agreement with the isotopic results [9], demonstrating that the metals in the Nefza-Sejnane deposits were supplied by the regional magmatic through hydrothermal fluids. Interestingly, Decrée et al. [2] observed that the Sidi Driss and Douahria sulfide ore deposits are hosted within carbonate lenses initially composed of Fe-Mn-enriched dedolomite. The dedolomitization is a well-known process that usually takes place when dolomite interacts with meteoric fluids (e.g., [33]). Concurrently with an Mg loss, a Fe-Mn gain led to formation of Fe-Mn carbonates or pure calcite (in the absence of Fe-Mn gain). At Sidi Driss and Douahria, dedolomitization is observed in the first stage of the ore sequence, before the early replacement by barite and celestite [1]. Consequently, dedolomitization could be related to the late Tortonian episode of weathering recognized at Tamra ( Figure 5).
In the Nefza-Sejnane district, the late Tortonian episode of weathering was followed by circulation(s) of hydrothermal fluids, leading to the precipitation of nontronite and talc observed in several localities of the district and dated at 6.2 ± 0.4 Ma [34]. Neoformation of nontronite and talc is observed mainly in (i) the basaltic flows of Mokhta El Hadid (Zouara basin); and (ii) the Triassic material and Miocene rhyodacite of Oued Belif (Figures 2 and 5). These hydrothermal circulations are most likely involved in the mineralization in Pb-Zn sulfides, barite, celestite, and fluorite at Sidi Driss and Douahria. Mn oxides constitute another milestone for the district. Decrée et al. [2] stressed that hydrothermal enrichment was superimposed on the weathering conditions found in the Tamra sequence, leading to the formation of first, hollandite, romanechite, and Sr-cryptomelane, dated at 4.7 ± 0.1 Ma by 40 Ar- 39 Ar dating, and then coronadite, chalcophanite, and amorphous Mn-oxides at 3.4 ± 0.1 Ma. Hydrothermal circulations were confirmed by Dekoninck et al. [4], studying the δD and δ 18 O isotopes of kaolinite-halloysite not in equilibrium with meteoric waters. Our new data show that a late Tortonian weathering phase predates the Messinian-Pliocene hydrothermal episodes leading to the formation of Mn oxides, Pb-Zn sulfides, barite, celestite, and fluorite.
A late Pleistocene weathering episode has been established by our new (U-Th)/He analyses. This recent meteoric fluid circulation is most probably responsible for the late occurrence of anglesite and iron oxide at Sidi Driss and Douahria [1]. The meteoric fluids associated with this recent weathering event affected the whole district. They are most likely responsible for the numerous weathering profiles observed in the area, including in the superficial part of the Albian-Priabonian marls and in the Numidian Flysch Formation, which hosts kaolinitic industrial clays [3,6,7].

Significance of the Weathering Episodes at a Regional Scale
The suitable conditions for the development of significant weathering profiles generally result from a combination of slow large-scale tectonic movement and (sub)tropical climate. At a regional scale, the new ages obtained contribute to unravel the role of (i) the geodynamic/tectonic setting, since exposure to weathering depends on exhumation (uplift) creating a hydrodynamic gradient for the percolation of meteoric fluids, and (ii) the climate, considering that meteoric water/fluid is needed. Our new data highlight two episodes of weathering, during the late Tortonian and the late Pleistocene, in the Nefza-Sejnane district. Interestingly, weathering events were already known during these periods in Tunisia. Based on C-O stable isotopes of non-sulfide Zn-Pb Tunisian deposits, Garnit et al. [35] suggested that Tunisian supergene deposits resulted from weathering during the middle to late Miocene interval and the Pliocene-Quaternary period. Studying the karstic system in Aïn Khamouda (Central Tunisia), Bruyère et al. [36] concluded that karstification took place in Tortonian or post-Tortonian period, as for other comparable karsts located in the Algerian Constantinois. Our new (U-Th)/He data on goethite refine the previously suggested episodes of weathering in Tunisia.
Climate has been regarded as an important parameter of weathering periods because sufficient water supply is needed for weathering processes (e.g., [37]). Beauvais and Chardon [38] and Chardon et al. [39] in West Africa propose that long-term eustatic sea level fall and climate play a significant role in weathering profiles. The late Pleistocene is marked by wet climate in North Africa ("African Humid Period"; e.g., [40]). Tortonian climate was presumably humid to subhumid around the Mediterranean Sea, while the hinterland was dryer ( [41][42][43]), as is the case today. During the late Tortonian, North Africa experienced a temporary turnover from dry to wet climate, initiating the "Zeit Wet Phase" of the Messinian (e.g., [44]). Conversely, our GHe ages do not match Cenozoic warmest climates in the studied area. Comparing these climatic conditions with the weathering periods obtained by our new data supports that warm climate is a second-order control, whereby wet climatic conditions are necessary to ensure the existence of meteoric fluids. The relatively high precipitation during the late Tortonian and late Pleistocene provided water to the weathering system, when differential uplift movements generated the weathering gradient for the downward migration of meteoric fluids.
Large-scale, vertical, positive movements were able to produce protore (including potential hypogene sulfides) uncapping and weathering, by creating a gradient for the percolation of meteoric fluids. In the last decade, the temporal link between uplift-related exhumation and weathering periods at a regional scale has been documented by several studies, in various geological contexts. Using 40 Ar- 39 Ar method on cryptomelane in Mn ores of Democratic Republic of Congo, De Putter et al. [28], Fontaine et al. [45] and De Putter and Ruffet [46] highlighted successive episodes of weathering since Late Cretaceous and concluded that tectonics accounts for most of the recognized supergene ore formation episodes, which are controlled by vertical lithospheric movements that are ultimately responsible for alternating stages of landscape stability and erosion. In the Ardenne (Western Europe), a recent synthesis showed that the uplift of the area was most probably the main factor triggering the development of weathering mantles, together with the increasing seasonality and high precipitation [47], illustrating the role of geodynamics and climate in the weathering.
In the same way, our results may be interpreted in relation with the regional Mediterranean geodynamics. In Morocco, weathering episodes are most likely associated with some of the four main Atlasic geodynamic events, as attested in recent geodynamic and ore deposits studies [48][49][50]. In Tunisia, both periods of weathering highlighted by our new (U-Th)/He data correspond to two distinct Alpine tectonic pulses that produced the exhumation of sulfides ores [35]. Alpine compressional movements may have caused marked slow and large-scale (at least regional) uplifts, preventing deposition of coeval sediments, except in the local extensional basin, such as in Tamra and Zouara basins. Then, the studied area must be regarded as an elevated zone during late Tortonian and late Pleistocene times, favoring the formation of weathering profiles. In the Iberian Pyrite Belt (Southern Spain), preliminary ages of different parts of the gossans using 40 Ar- 39 Ar techniques show that at least two stages of oxidation took place at 7-8 Ma and about 1-2 Ma [51], quite coeval with the new GHe ages in Tunisia. The late Tortonian meteoric episode corresponds to a turning point in Mediterranean basin dynamics with the disruption of the Mediterranean-Atlantic connection, an important tectonic rearrangement due to rollback of the Gibraltar slab [52,53]. This may suggest that the late Tortonian and late Pleistocene time intervals are periods of weathering at the Mediterranean region scale, partly driven by (at least regional) large-scale uplift-related exhumation associated to Alpine events. Further studies in the Atlasic and Mediterranean domains are, however, required to refine the nature, scale, and role of this uplift-related weathering, assuming the available age record is complete.
The stratified Fe ore of Douahria is stratigraphically correlated with the SEDEX-type Pb-Zn ore of Sidi Driss related to thermally-driven fluid circulations [1]. From biostratigraphic and biogeographic perspectives, the Douahria gastropod assemblage has a late Neogene signature at a Maghreb scale (late Miocene-Pliocene; [31,32]), without further precision. Only a late Miocene age of these fossils is in agreement with our (U-Th)/He dating of the Fe-rich fluids (ferruginous crust) of the area, suggesting that the first episode of ferruginous crust containing the fossils is late Tortonian in age (8.6 ± 0.9 Ma-see above).
A systematic revision of the mammalian fauna from Douahria may provide a better age constraint. The giraffid from Douahria was first referred to as "Helladotherium duvernoyi", a short-necked sivatheriine species originally described by Gaudry and Lartet [56] at Pikermi, a classic late Miocene vertebrate locality from Greece (dated around the Tortonian-Messinian transition, at 7.5-7 Ma; for review, see [57]). The dental and postcranial remains were recently reassigned to the longer-ranged sivatheriine Sivatherium hendeyi [57], documenting the late Miocene-early Pliocene interval in Africa (11.6-4 Ma; [20,58]). The Douahria anthracotheriid was originally assigned by Roman and Solignac [11] to the Asian species "Merycopotamus dissimilis", of late Miocene-earliest Pleistocene range ). In their revision of African anthracotheriids, Holroyd et al. [60] and Lihoreau et al. [61] consider that the referred specimens from Africa belong to Libycosaurus instead. Dental features and moderate dimensions further closely match those of Libycosaurus anisae [62], known to range the latest middle-early late Miocene interval in East and North Africa (~12-10 Ma [60,61]). Aside from Douahria, this species has been recognized in Central Tunisia (Bled Douarah, Beglia Formation, early late Miocene [61]) and Southern Tunisia (El Ziz, late Miocene [63]). The two-horned rhinocerotine rhinocerotid Diceros douariensis was named by Guérin [18] based on the abundant craniomandibular and postcranial material from Douahria. Guérin [64] proposed an age of 9.5 Ma for the latter locality, based on biostratigraphic correlations with the Djebel Krechel el Artsouma (Central Tunisia) and other circum-Mediterranean occurrences of D. douariensis. Geraads [19] later reassigned this species to Ceratotherium and considered that specimens from As Sahabi, Northern Libya (latest Miocene; 6.2-5.33 Ma) were also probably assignable to Ceratotherium douariense, thus significantly extending upward its stratigraphic range (i.e., 9.5-5.33 Ma). From a morpho-anatomical perspective, this species has the closest affinities with Ceratotherium neumayri [54], a conspicuous element of eastern Mediterranean faunas throughout the late Miocene (e.g., [55,65]).
Direct re-examination of the specimens in the Lyon University collection led to the unexpected recognition of a second rhinocerotid taxon at Douahria, as already mentioned by Geraads [19]. Most rhinocerotid remains do belong to C. douariense, as typified by the adult skull FSL 16749 (holotype; Figure 6A,B). Nevertheless, the juvenile skull FSL 16752 ( Figure 6C-F), previously defined as the paratype of "Diceros douariensis Guérin, 1966", has cranial and dental characteristics discarding any referral to a two-horned rhinocerotine, such as the presence of a concave dorsal profile and a rounded nasal notch ( Figure 6C), of laterally-projected orbits ( Figure 6D), of a lingual wall on upper premolars ( Figure 6E), of sagittally-elongated molars with an ectoloph undulated in occlusal view, and the abundance of cement filling the valleys of cheek teeth ( Figure 6F). These features are diagnosing middle and early late Miocene elasmotheriine rhinocerotids from Eurasia and Africa, i.e., the ones retaining non-ever-growing teeth [66][67][68]. The juvenile skull FSL 16752 cannot be assigned to any known elasmotheriine species, notably as documented in Africa (i.e., Kenyatherium bishopi [69]; Ougandatherium napakense [70]; Victoriaceros kenyensis [71]; Samburuceros ishidai [68]; and references therein). It is most likely documenting a new species, the formal description of which is far beyond the scope of the current work. Accordingly, no precise biostratigraphical age can be formally inferred from this occurrence, as usual regarding new taxa. Yet, given its evolutionary stage among elasmotheriine rhinocerotids [66,67], a late middle to early late Miocene age can be suspected for the first appearance of this taxon. Independently, the Douahria assemblage can be refined at~11.6-10 Ma, thanks to the concurrent range of the bothriodontine anthracotheriid Libycosaurus anisae [62] and of the sivatheriine giraffid Sivatherium hendeyi [20] documenting the 12-10 and 11.6-4 Ma interval, respectively (Figure 7). This stratigraphical range as inferred for Douahria is further compatible with its rhinocerotid components. A Tortonian age is thus supported here for the Douahria mammalian assemblage. Accordingly, deposition of the Douahria fossils may have occurred in (early) Tortonian times, just prior to the circulation of the meteoric fluids, here dated late Tortonian (8.6 ± 0.9 Ma) by (U-Th)/He method (Figure 7). No long gap is necessarily to be suspected between both the events.      [19,20,60,61]. AFLMA = African Land Mammal Ages [72].

Conclusions
The new (U-Th)/He ages obtained on goethite in the mineralized and vertebratebearing sediments of the Nefza-Sejnane district (Tunisia) highlight that: • meteoric fluids circulated during late Tortonian (8.6 ± 0.9 Ma) and late Pleistocene (0.8 ± 0.2 Ma) time intervals, confirming the polyphased and mixed (hydrothermal/meteoric) fluid circulations in the district and their role in the genesis of part of the mineralization, • these events allow refining the ages previously suggested for weathering episodes in Tunisia, confirming that (U-Th)/He on rigorously sampled and well-characterized (using petrography and δD-δ 18 O analyses) goethite is a very powerful method to decipher weathering periods, • these episodes are coeval with geodynamic events and wet climates in North Africa, confirming a temporal link between weathering events and both regional tectonics and humid periods, as already stressed by other studies in other areas, • matched with previous stratigraphic data, the chronological frame of the district is refined, allowing the integration of the numerous regional deposits (Pb, Zn, Fe, Mn, LREE, U, clays), • consequently, the deposition of the Douahria fossils in a ferruginous crust, here taxonomically revisited, occurred no longer prior to the late Tortonian circulation of the meteoric Fe-fluids, here dated by (U-Th)/He method.