samedi 28 mai 2016

Neoproterozoic subduction-related metavolcanic and metasedimentary rocks from the Rey Bouba Greenstone Belt of north-central Cameroon in the Central African Fold Belt: New insights into a continental arc geodynamic setting

Geological sketch map of Cameroon



This paper presents new geochemical and geochronological data for the low to medium grade Rey Bouba Greenstone Belt (RBGB), located in northern Cameroon at the northern margin of the Central African Fold Belt (CAFB), and discusses the maximum age of volcanic activity, the maximum depositional age of metasedimentary rocks and geotectonic implications. Geochemical data on volcanic rocks highlight the predominance of transitional to calc-alkaline magma compositions (Zr/Y = 3.62–26.38) with mostly andesite to basaltic andesite with an unusually high Mg concentration (>5 wt.%, high-Mg andesite), but also basalt and trachy-andesite protoliths. Moreover, chondrite-normalized REE patterns and primitive mantle-normalized spidergrams show enrichment of LREE relative to HREE with flat to depleted elements, and moderate to slight negative Nb–Ta, Ti and Eu anomalies respectively; which are consistent with a continental arc setting related to a subduction zone. U–Pb zircon LA-ICP-MS geochronology on felsic metavolcanic or metatuff fix the maximum age for the volcanic activity in the RBGB at ca 670 Ma. Detrital zircon grains indicate that Neoproterozoic zircon (Ediacaran to Cryogenian) are the main source for the detritus that fed the basin, combined with minor Paleoproterozoic and Mesoproterozoic inputs. The maximum depositional age, corresponding to the youngest graphical age peak controlled by multiple grain ages is consistently constrained between 645 and 630 Ma, whereas the age of low grade metamorphism weakly recorded by overgrowths on detrital zircon in the RBGB basin is around 600 Ma. These results provide new insights into the geodynamic processes during the Neoproterozoic along the northern margin of the CAFB of northern Cameroon, suggesting that the RBGB, where high-Mg andesite magmatism has taken place is consistent with a continental arc-related basin

M. Houketchang Bouyoa,b,∗, Y. Zhaoa,∗∗, J. Penaye b,∗∗, S.H. Zhanga, U.O. Njel b

( note: please, you can find all details on the orginal document. see in the end of the text)

1. Introduction 

The assembly of the Gondwana Supercontinent during the Late Neoproterozoic–Cambrianinvolved closure ofthe intervening Neoproterozoic ocean basins and subduction of a substantial volume of oceanic lithosphere along a number of convergent margins (Collins et al., 2007; Santosh et al., 2009, 2012; Boger, 2011; Ngako and Njonfang, 2011). Major accretionary processes contributing to continental growth in Africa during the Neoproterozoic have also been  identified (Condie, 2003; Penaye et al., 2006; Pouclet et al., 2006; Tchameni et al., 2006; Isseini et al., 2012). The CAFB is a collage of Paleoproterozoic microcontinents and Neoproterozoic plutonic and volcanic arcs attached to the Archaean Congo craton during a Pan-African continental collision at ca 600 Ma (Toteu et al., 2004, 2006; Van Schmus et al., 2008; Bouyo Houketchang et al., 2009, 2013; Nkoumbou et al., 2013). Well-preserved Neoproterozoic magmatic arcs bounded by narrow low to medium grade volcanosedimentary schist belts have been described in northern Cameroon and in southwestern Chad Republic (Fig. 1). These volcanosedimentary sequences include the Poli, Bibemi-Zalbi and Rey Bouba greenstone belts in northern Cameroon, and Zalbi and Goueygoudoum greenstone belts in southwest Chad. In previous works, these sequences were generally interpreted as pre-tectonic back-arc basins intruded by or associated withthe calc-alkaline TTG

∗ Corresponding author at: Centre for Geological and Mining Research, PO Box 333, Garoua, Cameroon. Tel.: +237 696 215 566. ∗∗ Corresponding authors. E-mail address: mbouyo2100@yahoo.fr (M.H. Bouyo).

suite of the Sinassi and Mayo Kebbi Batholiths (Toteu et al., 1984, 1987, 2006; Pinna et al., 1994; Pouclet et al., 2006). In the Cameroonian part, where most of the geological, geochronological and isotopic data have been well-documented in the Poli, but also in the Bibemi-Zalbi belts, the Rey Bouba belt has been poorly studied and very little is known about its geochemistry, age and provenance of detrital and volcanic material. The ages and geochemistry of volcanic and sedimentary rocks from the Rey Bouba basin are therefore crucial to better understand the tectonic setting and geodynamic evolution of the CAFB along its northern margin. We present and discuss here new geochemical and laser inductively coupled plasma mass spectrometry (LA-ICP-MS) U–Pb data for zircon from metavolcanic and metasedimentary rocks to show that the Rey Bouba Greenstone Belt, was likely deposited in a continental arc setting related to a subduction zone.

2. Geological setting

 The NE–SW trending Rey Bouba Greenstone Belt is one of the volcanosedimentary belts of the CAFB in northern Cameroon (Fig. 1 and Table 1). Located in the Western Cameroonian Domain, it defines a narrow greenschist belt extending about 80 × 16 km, closely related to the low to high grade Poli Belt and continues beyond the Chadian border by the greenschist Goueygoudoum Belt and the Bibemi-Zalbi Belt. Furthermore, it is located along the Tcholliré-Banyo shear zone (TBF) interpreted as a major terrain boundary separating the younger Neoproterozoic to Mesoproterozoic Western Cameroonian Domain on the west side from the older reworked Paleoproterozoic Adamawa-Yadé Domain on the east side (Penaye et al., 1989; Pinna et al., 1994; Toteu et al., 2001, 2004; Van Schmus et al., 2008; Bouyo Houketchang et al., 2009).

Fig. 2. Geological map of the Sinassi region showing the volcanosedimentary Rey Bouba Greenstone Belt and the associated granitoids. A schematic stratigraphic column for the RBGB is inserted.

The RBGB mainly consists of greenschist-facies mafic to felsic volcanic, volcanosedimentary and sedimentary rocks associated with a set of pre-, syn- and post-tectonic granitoids and dykes. It has been defined as a back-arc basin related with an oceanic plate subduction below the southeastern continental margin of the Adamaoua-Yadé Domain (Pouclet et al., 2006; Ngako and Njonfang, 2011). However, only very few geochronological data obtained by Pb–Pb minimum ages on single zircon are available for the RBGB indicating ages of 557 ± 17 Ma for the post-tectonic Vaimba granite and 750 ± 20 Ma for the Gatougel dacitic tuff (Pinna et al., 1994).
The Poli Belt defines a pre- to syn-collisional basin developed upon, or in the vicinity of young magmatic arcs. The filling of the basin occurred in a docking-arc/back-arc context (Toteu et al., 2006). It consists of low, medium- to high-grade Neoproterozoic schists and gneisses of volcanic, volcano-sedimentary and sedimentary origin. Metavolcanics are tholeiitic basalt and calcalkaline rhyolite emplaced in an extensional crustal environment (Njel, 1986; Toteu, 1990). The depositional age is constrained between 700–665 Ma; detrital sources comprise ca. 920, 830, 780 and 736 Ma magmatic rocks (Toteu et al., 1987, 2006).
The Bibemi-Zalbi Belt extends continuously within both the Cameroonian and Chadian territories, and is locally called Bibemi and Zalbi Greenstone Belt, respectively. It is well studied in SW Chad, where it is dated around 700 ± 10 Ma on metabasalt (Isseini, 2011) and 777 ± 5 Ma on epiclastite (Doumnang, 2006). Following the assumption of Pinna et al. (1994) of an arc and back-arc basin systemlinked to a subduction beneath theAdamawa-Yadé Domain, Pouclet et al. (2006) interpreted the region in terms of fore-arc basin – volcanic arc – back-arc basin that were accreted eastward to the Adamawa-Yadé Domain, along the Tcholliré-Banyo shear zone. This model implies the subduction and closure of a western oceanic basin that was located between the Western Cameroonian Domain that could have belonged to the Central Saharan Ghost Craton (Black and Liégeois, 1993) and the remobilized Paleoproterozoic Adamawa-Yadé Domain. All these terranes are now included in the southern part of the Saharan Metacraton (Abdelsalam et al., 2002).
At the scale of the CAFB in Cameroon, the tectonic history is complex and summarized by Ngako et al. (2008) and Ngako and Njonfang (2011) in three main tectonic events related to PanAfrican collision and post-collision evolution: (i) crustal thickening (ca 630–620 Ma, and even 600 Ma, Bouyo Houketchang et al., 2009, 2013); (ii) left lateral wrench movements (613–585 Ma); and (iii) right lateral wrench movement (585–540 Ma), successively; the latter being related at the global scale to the final amalgamation of the Gondwana Supercontinent (Alkmim et al., 2001; Meert, 2003; Collins and Pisarevsky, 2005) during Latest Neoproterozoic-Earliest Cambrian.

3. General sample descriptions and procedures for geochronological and geochemical analyses 

At the local scale, the geological framework of northern Cameroon is dominated by a NE–SW extensive magmatic arc province which includes a heterogeneous and complex plutonic (TTG) and volcanosedimentary sequences (Fig. 2) that underwent a polyphase deformation characterized by sub-vertical foliation, sub-horizontal stretching lineation and folds associated with greenschist to amphibolite facies metamorphism. In this study, we focused on greenschist facies volcanic and sedimentary components of the RBGB from which forty five samples were collected (Fig. 2). From metric to multi-metric thickness, mafic metavolcanic rocks are massive, in elongated bands, flagstones, or blocks, which commonly alternate with felsic metavolcanic, metasandstone, metasiltstone, quartzite, schist and slaty shale locally highly distorted. Conglomeratic layers consisting of angular to rounded clasts of mafic volcanic, granitoid, quartzite and mineral  fragments such as quartz within a finer-grained chlorite-rich matrix are observed in Baba Sara area (Fig. 3a–d). Thirty samples were selected for thin sections, 18 for geochemistry and four for geochronology. In thin section, most of the rocks examined display poorly preserved primary magmatic and sedimentary textures superimposed by metamorphic recrystallisation fabrics (Fig. 3e–h). Major elements except FeO were analyzed on fused glass discs by X-ray fluorescence spectrometry and FeO contents by classical wet chemical analysis at the Analytical Laboratory of the Beijing Research Institute of Uranium Geology. Trace element concentrations were determined using inductively coupled plasma-mass spectrometry (ICP-MS) at the Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou, China. Detailed analytical methods for ICP-MS were described by Chen et al. (2010). For most of the trace elements, analytical precision and accuracy are better than 10%. Major element analyses were recalculated to 100 wt.% anhydrous basis for intercomparisons, where Fe2+ is assumed to be 80% total Fe, with the prior conversion of reported Fe concentration (FeO or Fe2O3) as Fe2O3t from appropriate atomic or molecular weights (Vocke, 1999; Verma and Armstrong-Altrin, 2013). Chondrite and primitive mantle reservoir compositions are those of Sun and McDonough (1989). Mg-number indicating the level of evolution of volcanic rock was defined as 100MgO/(MgO + FeO) in mol per cent. Being highly resistant to chemical and physical influences, zircon is a particularly usefulmineralfor petro-chronological investigations (Corfu et al., 2003; Sircombe, 2004; Vermeesch, 2004; Andersen, 2005). In order to constrain the timing of volcanic activity and depositional age of the RBGB, one felsic metavolcanic tuff, VA-R 127, and three representative metasedimentary samples including, VA-R 193, VA-R 204 and VA-R 222, were collected for detailed study and dating. Zircon grains from each sample were separated from 2 to 5 kg crushed rock samples by conventional magnetic and methylene iodide liquid separation. The separated zircon grains were handpicked and mounted in epoxy resin. The epoxy mounts were polished to expose the mounted minerals, carbon coated, photographed in both transmitted and reflected light, and imaged via Cathodoluminescence (CL) using the JSM- 6510 scanning electron microscope (SEM) at Beijing Geoanalysis Limited Company. U–Pb dating of zircon was conducted by LA-MCICP-MS (Gehrels et al., 2008; Johnston et al., 2009; Jackson et al., 2004) at the Isotopic Laboratory of Tianjin Centre of Geological Survey-China. Laser sampling was performed using a UP193-FX, and a Neptune MC-ICP-MS instrument was used to acquire ionsignal intensities. Heliumwas applied as a carrier gas for the ablated material. Trace element compositions of zircon were calibrated against Nist610 using Si as an internal standard. GJ-1 standard gem zircon from Sydney (Australia), with 207Pb/206Pb age based on 8 TIMS determinations of 609 Ma (Jackson et al., 2004), used as an external standard for U–Pb dating was mounted with samples cleaned in 1 N nitric acid immediately prior to analysis to remove surface Pb contamination, and then measured twice every eight analyses. The counting time during the analyses is 60 s (including the first 20 s for background). Fractionation and instrumental mass bias are corrected by direct calibration against a zircon standard analyzed under carefully matched conditions, using He as the ablation gas to increase the reproducibility of the Pb/U fractionation. Time-resolved data acquisition is employed to evaluate zircon homogeneity and to allow selective integration of signals to minimize common Pb contributions and Pb loss, and thus to maximize concordance. LA-ICP-MS measured values for GJ-1show that zircon is relatively low in Th, with mean U and Th contents of 230 and 15 ppm, respectively. The Calibration errors (2) on the 206Pb/238U, 207Pb/235U and 207Pb/206Pb ratio were 1.9%, 3.0%

Fig. 3. Field outcrop photographs and thin section microphotographs of various rock types of the RBGB. (a) Mafic metavolcanic rock at Mayo Rian-Mayo Vaïmba confluence; (b) Fine-grained metasandstone in Mayo Vaïmba river; (c) Metasandstone at Baba Sara; (d) Metaconglomerate around Baba Sara; (e) Felsic metavolcanic rock (VA-R 235) with abundant very small crystals (suggesting recrystallisation has occurred) of quartz ± plagioclase in the lighter bands, but also few epidote ± sericite ± calcite in darker bands; (f) Felsic metatuff (VA-R 127) with quartz + feldspar embedded in argillaceous matrix; (g) Medium metasandstone (VA-R 191) showing mature sediments with well sorted grains; (h) Metasiltstone (VA-R 166) showing a very fine-grained matrix and concentration of small muscovite flake + quartz intergrown after clasts of feldspar.
Fig. 4. CaO/Al2O3–MgO–SiO2 diagram showing geochemically unaltered nature of most of the metavolcanic rock samples of the RBGB. 

and 2.4%, respectively and propagated through the error analysis. During the experiment, sites for dating were selected on the basis of CL and photomicrograph images, in order to obtain as representative a population as possible. Laser energy of 10–11 J/cm2, a frequency of 8 Hz and a spot size of 35 m were used for the instrument. Off-line selection, integration of background, analyte signals, time-drift correction and quantitative calibration for trace element analyses were performed by ICPMSDataCal (Liu et al., 2008). Concordia diagrams, histograms overlain by probability density distributions and weighted mean calculations were made using Isoplot (Ver3.23) (Ludwig, 2003) to show the distribution of the zircon populations over the full age range of all grains (see text, Tables 3–6 and Figs. 6–8 for details on individual samples).

 4. Geochemical results 

Major and trace element data of the investigated rocks from the RBGB are presented in Table 2 (complete data set may be accessed from the online data repository (see Appendix A)). Postdepositional processes as well as low-temperature metamorphism (metasomatism) are known to affect the mobility of certain elements (Cs, Rb, Li, B, Ba, Sr, Si,K, etc.). Index of alteration also includes the degree of hydration (i.e. Loss on Ignition which it is generally less than 1%). Some of the analyzed samples in the present study have Loss on Ignition (LOI) values higher than expected (Table 2). The data were, therefore, plotted in the CaO/Al2O3–MgO–SiO2 diagram, proposed by Schweitzer and Kröner (1985) to test alteration effects in metavolcanic rocks. As seen, majority of the samples plot in the field of “unaltered rocks” except few samples with highest LOI (VA-R 223 and VA-R 127), which also show lowest (0.35 wt.%) and highest CaO (6.39 wt.%) contents respectively (Fig. 4; Table 2). Both mafic to intermediate and felsic metavolcanic rocks of the RBGB are characterized by magmatic compositions ranging from andesite/basaltic andesite to trachy-andesite. Only one sample clearly shows basaltic composition in the bivariate Zr/Ti versus Nb/Y diagram (Fig. 5). In addition, they show (Zr/Y = 3.62–26.38) ratios generally comparable to those of transitional (Zr/Y = 4.5–7.0) and calc-alkaline volcanic compositions. On the basis of petrographic characteristics, major elements, and chondrite and primitive mantle normalized rare earth and trace element patterns, they have been subdivided into three major groups (Fig. 6; Table 2): mafic to intermediate metavolcanic, felsic metavolcanic and metasedimentary rocks.

Fig. 5. Zr/Ti versus Nb/Y classification diagram for metavolcanic rocks of the RBGB (Table 2). Metasedimentary rocks are also plotted for intercomparisons. Compositional fields revised by Pearce (1996) after Winchester and Floyd (1977).

 4.1. Mafic to intermediate metavolcanic rocks

 Mafic to intermediate metavolcanic rocks (Fig. 5) are enriched in MgO (2.89 and 6.19 wt.%), moderate to strong fractioned with a variable Mg-number (Mg#) between 38 and 65% (Table 2). They have low to medium SiO2 (50.95–63.00 wt.%) contents and relatively high TiO2% (0.65–1.03), MgO% (2.89–6.19), and CaO% (0.35–8.86) respectively with an average of (0.83, 4.50, and 4.64 wt.%; Table 2). The relatively high concentrations of transition metals V (up to 332 ppm), Cr (249 ppm), Ni (162 ppm), Cu (236 ppm), Sc (28 ppm) and total REE (313 ppm) are noticeable (Table 2). On chondrite and primitive mantle normalized diagrams, they display enrichment in LREE and slight depletion in HREE with no Eu anomaly (Fig. 6a and Table 2); but also erratic distribution in the mobile LILE (Rb, Ba, K) and negative anomalies in Nb-Ta and Ti (Fig. 6b). Inter-element ratios Nb/Ta and (La/Yb)N give values of 0.18–0.58 and 1.92–39.00 (with an average of 10.62), respectively.

4.2. Felsic metavolcanic rocks

 Felsic metavolcanic rocks display the highest SiO2 (69.51–75.37 wt.%) and lowest MgO (0.41–1.46 wt.%) contents (Table 2). They are also characterized by low concentrations oftransition metals: Ti (0.11–0.49 wt.%), Sc (4–13 ppm), V (14–62 ppm), Cr (2–22 ppm), Ni (1–12 ppm) and total REE (33–68 ppm). Felsic metavolcanic rocks exhibit enriched LREE and flat to depleted HREE patterns (Fig. 6c). The trace element patterns of this group are relatively enriched in LILE (Rb, Ba, K) and have negative anomalies in Nb-Ta, P and Ti (Fig. 6d); but variable Eu anomalies (Eu/Eu* = 0.60–1.31; Table 2). Nb/Ta and (La/Yb)N ratios are between (0.13–0.38) and (1.45–28.35 with an average of 12.48) respectively.

4.3. Metasedimentary rocks

Metasedimentary rocks, compared with the previous groups show variable contents of SiO2 (58.61–66.86 wt.%), Mg (2.36–3.71 wt.%); transition metals Ti (0.57–1.00 wt.%), Sc (13–20 ppm), V (88–131 ppm), Cr (65–114 ppm), Ni (34–64 ppm) and total REE (65–157 ppm) (Table 2) between the highest values for mafic metavolcanic rocks and the lowest for felsic metavolcanic rocks and vice versa. However, Ba and Zr have the highest contents, reaching respectively up to 1739 ppm and 212 ppm. The chondrite normalized plots show enrichment in LREE relative to HREE which have a flat pattern and slight fractionated or negative Eu anomaly (Fig. 6e; Table 2). Regarding the primitive mantle normalized diagram, they show enrichment in LILE (Rb, Ba, K) and negative anomalies in Nb-Ta and Ti (Fig. 6f). Their Nb/Ta and (La/Yb)N ratios are respectively (0.33–0.78) and (2.66–8.35 with an average of 5.61).

 Fig. 6. Chondrite and primitive mantle normalized diagrams for (a and b) Mafic to intermediate metavolcanic rocks, (c and d) Felsic metavolcanic rocks and (e and f) Metasedimentary rocks.

5. Geochronological results 

In order to give an adequate, accurate and precise feature of the age distribution ofthe volcanic and sedimentary rocks of RBGB, several hundred grains of zircon were selected randomly from each sample. Among the four samples selected for dating, 318 zircon crystals were individually analyzed in three samples, and only 26 crystals saved for the fourth; while a maximum of 120 representative target spots from each sample were located using CL images. In total, 394 analytical spots were performed on 344 zircon grains.

5.1. Sample VA-R 127: Felsic metatuff 

This sample is interpreted as a felsic metavolcanic rock or metatuff comprising a mixture of fragments and euhedral to subeuhedral fine- to medium-grained quartz and feldspar crystals embedded in a fine-grained matrix. Collected near Vaïmba village (Fig. 2), and generally displaying a grey colour or grey yellowish to brownish weathered surfaces, the rock has very fine-grained illite, sericite and chlorite associated with sub-angular to sub-rounded quartz, and minor feldspar, epidote and oxides. Late stage microfractures or veinlets filled by calcite are also locally observed. Zircon grains from this sample are light-pink to brownish in colour. They are commonly prismatic, but occasionally rounded or fragmentary in shape, possibly related to different populations. The grain size varies from ca 80 to 230 m, the majority being between ca 105 and 115 m. On CL images (Fig. 7a), all the grains display igneous features with oscillatory or sector zoning. In most cases, grains show different colours in grey, with faint, broad zoning, and in fewer cases, they are unzoned. Overgrowths are very rare to absent. A number of 110 points were analyzed in situ on 107 grains. Most of the analyzed zircon grains (92%) show 232Th/238U ratios higher than 0.2 (Table 3, complete data set may be accessed from the online data repository (see Appendix A)), pointing to a magmatic origin (Williams and Claesson, 1987). From the data, two main groups of concordant and discordant zircon (Fig. 8a and Table 3) are observed. The concordant group is very clustered (inset on Fig. 8a) and display three distinct sub-groups or populations on the probability density distribution and histogram diagrams (Fig. 8b), with a high concentration of 206Pb/238U ages at ∼670 Ma (33 analyses; median = 668 ± 8.5 [1] Ma), at ∼700 Ma (27 analyses; median = 697 ± 3.1 [1] Ma) and at ∼760 Ma (5 analyses; median = 758 ± 2.8 [1] Ma). The discordant group shows signifi- cant episodic loss of lead which may be related either to analytical issues or to relatively high U contents of analyzed zircon grains (dark domains, Fig. 7a: 127–3 and Table 3: 127–3: 713 ppm, 22: 371 ppm, 85: 476 ppm, 94: 253 ppm, 98: 451 ppm, etc.). However they yield an upper intercept of 1431 Ma and lower trajectory to ca 0 Ma (Fig. 8a).

5.2. Sample VA-R 193: Schistose sandstone 

RockVA-R193 is a light- tomedium-grey, schistose sandstone or metasiltstone from Baba Sara, showing a very fine-grained groundmass of silica and clayey minerals in association with clasts of quartz, pyrite and iron oxides. The unit is locally affected by recrystallisationfabrics due to intense deformation, and interbedded with felsic metavolcanic rock and quartzite bands, and commonly contains lithic fragments giving it the appearance of conglomeratic volcanosedimentary schist. Detrital grains of zircon range in size mostly from ca 110 to 200 m, with pale pink to yellowish colour. Grains are euhedral to sub-euhedral, mostly displaying oscillatory or sector zoning typical of igneous rocks (Fig. 7b). Moreover, some individual crystals show strong variations in the development of zoned domains, where one large uniform central zone is succeeded by much finer oscillatoryzoned bands. A total of 120 spots on 105 zircon grains were analyzed, with 94% indicating 232Th/238U ratios higher than 0.2 (Table 4, complete data set may be accessed from the online data repository (see Appendix A)). The U–Pb isotope data obtained from the zircon grains show a wide range of ages almost all concordant (115 analyses) between 529 ± 3 Ma and 906 ± 5 Ma (Fig. 8c) with a peak at 645 Ma representing a dominantly Neoproterozoic source (Fig. 8d). The presence of Late Mesoproterozoic and Paleoproterozoic components is also indicated by two zircon ages at 1016 ± 6 Ma (193–32), and 2077 ± 12 Ma (193–119) respectively. Five discordant data (193–10, 12, 48, 52 and 104; Table 4) are not considered.

5.3. Sample VA-R 204: Chl-schistose sandstone 

Sample VA-R 204 is a greenish grey to dark green, fine-grained chlorite-bearing schistose sandstone or volcanoclastic schist, collected near the river Mayo Godi, about 6 km NW from Baba Sara. The rock consists essentially of quartz and chlorite in a cryptocrystalline matrix. It is locally strongly deformed and interbedded with well-foliated slaty schists. Zircon grains are euhedral to sub-euhedral, and commonly vary in size from ca 135 to 180 m and up to 235 m in places. They are mostly pale pink and display morphologies and internal structures characterized by more or less well-preserved magmatic growth zoning (Fig. 7c). These grains for the most part have less blunted or eroded terminations, suggesting a short distance of transportation prior to deposition, while very few show thin overgrowths. Some 44 laser ablation analyses were performed on 26 crystals. 232Th/238U ratios range from 0.03 to 1.67 with the majority (80%) higher than 0.2 (Table 5, complete data set may be accessed from the online data repository (see Appendix A)). Data are plotted in Fig. 8e. Most analyses are concordant defining three main groups with 206Pb/238U ages ranging from 585 ± 4 to 876 ± 6 Ma (36 analyses, insert Fig. 8e) for the first group, one analysis at 1259 ± 10 Ma for the second group, and the third group at ca. 2000 Ma (four analyses). A predominant Neoproterozoic source is shown by a broad probability density distribution from Ediacaran to Tonian with maximum peaks at 663 and 700 Ma, but also at 682, 630 and 600 Ma (Fig. 8f). This sample also confirms the presence of older Mesoproterozoic and Paleoproterozoic sources for the sediments.

5.4. Sample VA-R 222: Ferruginous metasandstone

 The rock specimen is brown to ocher, probably related to its dominantly ferruginous cementingmaterial.Itis afine- tomediumgrained ferruginous metasandstone, collected upstream of the Mayo Beoulpir river. Minerals identified in thin section include principally sub-angular to sub-rounded quartz, but also plagioclase, calcite, pyrite, iron oxides, epidote, illite, chlorite and indistinct grains of brownish altered ferromagnesian mineral. Folded veinlets of quartz and calcite are also observed. Detrital crystals of zircon are prismatic to round, and pink, yellowish, or brown in colour. Crystals are sub-euhedral to euhedral, range in size from ca 75 to 200 m, with the majority between 110 and 150 m, and mostly show a well-developed oscillatory to sector zoning (Fig. 7d). The differences in morphology and colour among zircon may reflect distinct sources, and the small degree of abrasion suggests a proximal provenance. A whole of 120 measurements were performed in situ on 106 zircon grains, with 99% displaying 232Th/238U ratios higher than 0.2 (Table 6, complete data set may be accessed from the online data repository (see Appendix A)). Concordia diagrams (Fig. 8g) show scattered discordant points with relatively high U contents (dark domains, Fig. 7d: 222–13 and 26; Table 6: 222–13: 172 ppm, 26:

Fig. 7. Cathodoluminescence images of zircon from representative metavolcanic and metasedimentary samples of the RBGB basin: (a) felsic metatuff, VA-R 127; (b) schistose sandstone, VA-R 193; (c) chlorite-bearing schistose sandstone, VA-R 204 and (d) ferruginous metasandstone, VA-R 222. The figure also shows circles indicating analyzed spots (not at scale), with numbers representing 206Pb–238U age and spot position listed respectively in Tables 3–6.

227 ppm, 29: 156 ppm, 30: 240 ppm, 65: 327 ppm, 117: 326 ppm, etc.) possibly due either to analytical issues or to common lead or extreme lead loss; and concordant points (56 analyses, diagram insert) all yielding Neoproterozoic 206Pb/238U ages from 574 ± 4 to 893 ± 6 Ma. The combined probability density distribution and histogram diagrams on concordant points show a maximum peak at 660 Ma, but also at 740, 795 and 848 Ma (Fig. 8h).

6. Discussion

6.1. Age of the volcanic activity

U–Pb zircon LA-ICP-MS geochronology on a felsic metavolcanic rock or metatuff (sample VA-R 127; insert Fig. 8a) is used to constrain the timing of magmatic activity within the RBGB. On CL images (Fig. 7a), zircon grains from VA-R 127 show typical oscillatory growth or sector zoning and are characterized by high 232Th/238U ratios (Table 3) indicating a magmatic origin. Concordant zircon grains (65 analyses) reported on Fig. 8 a (insert) and b define three peaks (ca 760, 700 and 670 Ma) which we interpreted as the times of different magmatic sources in the region. They indicate that the magmatic activity in the RBGB is prominent and took place in the Neoproterozoic. The peak at ∼760 Ma reflects an earlier inherited magmatic event which is consistent with the Pb/Pb age of 750 ± 20 Ma on Gatouguel dacitic tuff (Pinna et al., 1994), whereas the youngest peak at ∼670 Ma fix the maximum age for the volcanic activity in the RBGB. From previous studies in the Poli Belt by Toteu (1990) and Toteu et al. (2001, 2006), a wide range of Neoproterozoic detrital zircon (Tonian–Cryogenian) at 700, 780, 830, and 920 Ma have been recorded (Table 1), some of which correspond to well identified magmatic rocks originating from neighbouring areas (e.g., Goldyna metarhyolite, medium-grade schists and underlying volcanic sequence). In the Chadian part of the Bibémi-Zalbi Belt, Isseini (2011) constrains a crystallization U–Pb age on metabasalt at 700 ± 10 Ma, quite different from the previous U–Pb age of 777 ± 5 Ma derived from an epiclastite by Doumnang (2006). The distribution of ages implies at least three major periods or episodes, representing formation times of volcanic rocks at the northern margin of the CAFB: (1) an early volcanic activity from ca 920 to 830 reported only in the Poli Belt,(2) a second from ca 780 to 700 reported in the Poli and Bibémi-Zalbi belts, and (3) a third volcanic event around 670 Ma or younger reported in the RBGB belt. These results show that the volcanic activity is older in the Poli and Bibémi-Zalbi belts and younger in the RBGB.

 6.2. Provenance of the sediments

 LA-ICP-MS U–Pb detrital zircon dating presented here provides new constraints onthe geological evolutionoftheRBGB, dominated by Neoproterozoic zircon with some Paleoproterozoic, Mesoproterozoic, and Early Cambrian inputs. Over 70% of zircon grains have Neoproterozoic (Ediacaran to Tonian) ages ranging from ca 574 to ca 906 Ma (samples VA-R 193, VA-R 204 and VA-R 222), with major

Fig. 8. U–Pb concordia and probability density distribution-histogram diagrams for detrital zircon analyses from representative metavolcanic and metasedimentary samples of the RBGB basin: (a and b) felsic metatuff, VA-R 127; (c and d) schistose sandstone, VA-R 193; (e and f) chlorite-bearing schistose sandstone, VA-R 204 and (g and h) ferruginous metasandstone, VA-R 222.

peaks at 645, ca 660 and 700 Ma (Fig. 8d, f and h) but also at 600, 630, ca 680, 740, ca 800 and ca 850 Ma (Fig. 8f and h); indicating that the Neoproterozoic (Ediacaran to Cryogenian) was the main source for the Rey Bouba basin. These ages correspond well with the U–Pb and Ndmodel ages ofthe juxtaposed NeoproterozoicWestern Cameroonian Domain (granitoids and volcanosedimentary sources at ca 595, 600, 615, 630, 660, 693, 736, 750, 780, 830, and 920 Ma; Table 1; Toteu et al., 2001, 2004, 2006; Bouyo Houketchang et al., 2009; Dawaï et al., 2013) as well as the Pan-African Mayo Kebbi Batholith or magmatic arc of SW Chad (ca 640, 665 and 740 Ma; Table 1; Penaye et al., 2006; Isseini et al., 2012). Looking at the high proportion (more than 90%) of 232Th/238U ratios higher than 0.2, this implies that the predominant sources of Neoproterozoic detritus that fed the RBGB basin were derived from erosion of the adjacent magmatic arc and Western Cameroon Domain or Block. In this regard, our data clearly demonstrate that northern Cameroon volcanosedimentary basins at the northern margin of CAFB are younger than the magmatic arc, and not intruded by the latter as suggested by previous works (Pinna et al., 1994; Penaye et al., 2006; Pouclet et al., 2006). However, there is also minor input of Paleoproterozoic (Orosirian to Rhyacian) zircon at ca 2000 Ma and 2077 ± 12 Ma, recorded on samples VA-R 193 (Fig. 8c) and VA-R 204 (Fig. 8e). Rare Mesoproterozoic (Stenian to Ectasian) zircon at 1016 ± 6 Ma (sample VA-R 193; Fig. 8c) and 1259 ± 10 Ma (sample VA-R 204; Fig. 8e), and very rare (less than 1%) Early Cambrian (Fortunian) zircon at 529 ± 3 Ma (sample VA-R 193; insert Fig. 8c) also occur. Paleoproterozoicdetrital zirconreportedfrommetavolcanosedimentary schist of the RBGB may be correlated with the nearby gneisses of Adamawa-Yadé Domain (e.g. Mayo Makat and Mayo Kout gneisses at 2218 ± 14 Ma, and also Mbé gneiss at 2014 ± 7 Ma, Penaye et al., 1989; Mayo Kout orthogneiss at 2039 ± 26 Ma, Bouyo Houketchang et al., 2009; Table 1). The Mesoproterozoic zircon may be regarded as evidence of eroded rocks that remain unknown or hidden, since the area is not mapped in detail. The Neoproterozoic ages at ca 600 Ma (sample VA-R 204, Fig. 8f) agree with tectonometamorphic ages previously reported from volcanosedimentary basins within the CAFB (zircon U–Pb at 620 ± 10 Ma, Yaoundé: Penaye et al., 1993; composite Sm–Nd isochron at 616 Ma, Bivouba-Yaoundé: Toteu et al., 1994; Pb–Pb zircon evaporation at 611 ± 20 Ma, Yaoundé: Stendal et al., 2006; zircon U–Pb at ca 600 Ma, Lom: Toteu et al., 2006; Sm–Nd garnet-whole rock at 628 ± 68 Ma, Bafia: Numbem Tchakounté et al., 2007; zircon U–Pb and Sm–Nd garnet-whole rock at ca 600 Ma, Banyo and Tcholliré: Bouyo Houketchang et al., 2009, 2013; monazite U–Th–Pb at 622 ± 43 Ma, Boumnyebel-Yaoundé: Yonta-Ngouné et al., 2010; at 613 ± 33 Ma, Yaoundé: Owona et al., 2011). The Early Cambrian ages are likely to be minimum estimates for the formation of RBGB or may represent post-metamorphic cooling.

6.3. Age of metamorphism 

Although very rare to absent on CL images (Fig. 7), younger metamorphic overgrowth rims around older magmatic cores are locally present and characterized by 232Th/238U ratios lower than 0.2 (less than 9% of the entire samples; Tables 3–6). Except for a few very young ages (e.g. 288 Ma or 432 Ma, VA-R 127–3 and 127-24 respectively; Table 3), which may be due to calibration problems or lead loss after deposition, about 67% of this metamorphic population yield Late Neoproterozoic ages, most often discordant. However, the youngest and best example to constrain the age of metamorphism is given by sample VA-R 204-24 with a concordant 206Pb/238U age of 600 ± 4 Ma, and 232Th/238U ratio of 0.03 (Table 5). This result is consistent with the ca 600 Ma age of metamorphism obtained from high pressure granulites (Bouyo Houketchang et al., 2009, 2013) in Banyo and Tcholliré (Poli Belt) 0 10 20 30 40 50 60 500 600 700 800 900 1000 1100 Age (Ma) Number All of the metasedimentary rocks; n=207 645

 Fig. 9. Synthesis diagram displaying the combined probability density distributionhistogram for the geochronological data from the three metasedimentary samples. within the CAFB of north-central Cameroon (Table 1). This is the most direct evidence for the timing of continental collision. 

6.4. Maximum depositional age of the RBGB basin

The maximum depositional age for Precambrian stratigraphic units lacking preserved biostratigraphic age control such as fossils is commonly constrained by the youngest U–Pb ages of zircon grains in populations of detrital zircon (Jones et al., 2009). As mentioned in previous sections, Neoproterozoic ages around 600 Ma and Early Cambrian ages are best interpreted as tectonometamorphic ages and post metamorphic cooling respectively, and therefore, they were not considered in the detrital zircon populations. The maximum depositional age for each sample is defined here as the 206Pb/208U age of the youngest dated 99.5% to 100.5% concordant zircon at 621 ± 5 Ma in the case of sample VA-R 204- 23 (Table 5; Fig. 8e); 631 ± 4 Ma in the case of sample VA-R 222-16 (Table 6; Fig. 8g), and 626 ± 4 Ma in the case of sample VA-R 193-13 (Table 4; Fig. 8c). Our results show somewhat different youngest ages from one sample to another, and so individually are not yet considered to represent actual depositional ages. Dickinson and Gehrels (2009) tested the research strategy of using youngest U–Pb ages of detrital zircon to constrain the maximum depositional ages of strata containing zircon grains by comparing a total of 5365 concordant or nearly concordant U–Pb ages of detrital zircon in 58 samples of Mesozoic sandstone from the Colorado Plateau and adjacent areas with depositional ages known independently from biostratigraphy. Their analysis confirms that using youngest detrital zircon ages to constrain maximum depositional age is a valid procedure, but indicates that results vary somewhat depending upon how the youngest grain age is specified. Moreover, they add that in general, the youngest-age measurements based on multiple grain ages are more consistently compatible with depositional ages. In our approach to constrain the maximum depositional age, we consider the geochronological data from the three metasedimentary samples as a composite sample of detrital zircon, in order to do statistical comparisons of similarity between zircon populations within the RBGB. From a total of 284 spots on 237 grains analyzed, 207 are considered significant in terms of concordance and have been used to construct a probability density distribution-histogram diagram (Fig. 9; Tables 4–6). From Fig. 9, a single peak controlled by multiple ages is defined at 645 Ma,representing the main provenance source of the detrital zircon, as well as the youngest probability density distribution peak. In light of Figs. 8f and 9, the maximum depositional age can be consistently constrained between 645–630 Ma, whereas 600 Ma represents the age of low grade metamorphism recorded by the detrital zircon overgrowths from the RBGB basin.

 6.5. Petrogenesis of volcanic and sedimentary rocks 

Knowledge of the tectonic setting of an ancient basin is very important for understanding the geodynamic evolution and associated mineral resources, as well as the paleogeographic environment. Although there are many discrimination diagrams for ultrabasic and basic tectono-magmatic, and sedimentary (Bhatia, 1983; Roser and Korsch, 1986; Armstrong-Altrin and Verma, 2005) settings in the literature with success rate sometimes controversial for a specific geodynamic context, very few or almost none exist for intermediate composition rocks. In this work, we constrain the RBGB basin tectonic setting based on new sets of multidimensional tectonic discrimination diagrams for intermediate magmas from five tectonic settings (island arc-IA, continental arc-CA, continental rift-CRand ocean island-OI combined together, and collision-Col) of Verma andVerma (2013), but also for sedimentary rocks from three tectonic settings (arc, continental rift, and collision) of Verma and Armstrong-Altrin (2013); all based on worldwide examples with high success rates of about 63–100%. In two different sets of five discrimination diagrams for intermediate rocks designed to discriminate among five tectonic settings based on major-elements (Fig. 10A, accessible from the online data repository (see Appendix A)) and immobile trace-elements (Fig. 10B, accessible from the online data repository (see Appendix A)), about 80% of our samples (mafic to intermediate metavolcanic rocks) fall in the arc tectonic setting field consistent with the dominantly transitional to calc-alkaline andesite to basaltic andesite geochemical signatures of the studied rocks (Fig. 5). Similarly, in the new multi-dimensional major-element based diagrams for tectonic discrimination of siliciclastic sediments from three main tectonic settings, the RBGB metasedimentary rocks mostly plot in the arc tectonic setting field for high-silica (Fig. 11a, accessible from the online data repository (see Appendix A)) and low-silica (Fig. 11b, accessible from the online data repository (see Appendix A)) contents. Moreover, the chondrite and primitive mantle normalized diagrams show respectively enrichment in LREE relative to HREE and moderate to slight negative Nb-Ta, Ti and Eu anomalies for all analyzed samples (Fig. 6); compatible with a subduction zone. Thus an arc setting related to convergence and subduction can be clearly inferred for the northern margin of the CAFB during the Neoproterozoic (645–630 Ma), and therefore indicates the depositional environment of the RBGB basin prior to collision at ca 600 Ma.

6.6. Geodynamic interpretation 

As mentioned above, mafic to intermediate metavolcanic rocks from the RBGB show geochemical features similar to those of wellstudied high-Mg andesite (HMA) from Mitsuiyama (Okoppe area of North Hokkaido in Japan, with MgO contents ranging from 4.9 to 7.9 wt.%, Ayabe et al., 2012). These include: andesitic composition (Fig. 5) with high contents of MgO (2.89–6.19% or Mg-number between 38–65% with an average of 52; Table 2), high concentrations of transition metals V (up to 332 ppm), Cr (249 ppm), Ni (162 ppm), Cu (236 ppm) and Sc (28 ppm), relatively low (La/Yb)N ratios between 1.92–39.00 with an average of 10.62, depleted rare earth elements (Fig. 6), and also high field strength elements (Nb, Ta, Ti), which are indicators of arc volcanic rocks developed in continental setting (Figs. 5, 10 and 11). This is supported by the detrital sources that fed the RBGB which are derived predominantly from erosion of the adjacent Neoproterozoic Mayo-Kebbi and Sinassi juvenile magmatic arcs of the Western Cameroonian Domain and from the Paleoproterozoic Adamawa-Yade Domain. Therefore, the RBGB may correspond to an extensional basin developed upon or behind a Neoproterozoic magmatic arc. Fig. 12 (accessible from the online data repository (see Appendix A) illustrates a subductionrelated geodynamic setting model for the RBGB (Fig. 12A) in the network of Neoproterozoic magmatic arcs (Fig. 12B) of northcentral Cameroon in the CAFB. Although the andesitic rocks from the RBGB are Mg-rich, their geotectonic setting isquitedifferentfromthat oftheprimitivehighMg andesite as defined by Wood and Turner (2009) or by Ayabe et al. (2012) in the Okoppe volcanic field of North Hokkaido, Japan. Finally, it seems evident that there are different types of HMA in subduction-related settings; some are primitive and result from interaction of a melt derived from subducted oceanic basaltic crust and the overlying mantle wedge peridotite (Shiraki et al., 1978; Meijer, 1980; Crawford et al., 1981; Tatsumi and Ishizaka, 1982; Xu et al., 2000) and others aremore evolved and resultfrominteraction of melt derived from subducted oceanic crust and the continental magmatic arc.

 7. Conclusions 

The main objective of our work was to acquire new geochemical and geochronological data from the low to medium grade volcanosedimentary Rey Bouba Greenstone Belt (RBGB) of northern Cameroon at the northern margin of the Central African Fold Belt (CAFB), in order to contribute for a better knowledge and understanding of tectonic settings and geodynamic evolution during its Precambrian history. The overall results of our work show that: - Petrogenesis of mafic to intermediate metavolcanic rocks from the RBGB is consistent with the dominantly transitional to calcalkaline andesite to basaltic andesite, indicating geochemical features that are similar to typical high-Mg andesite including: andesitic composition with high contents of MgO (2.89–6.19% or Mg-number between 38 and 65%), high concentrations of transition metals V (up to 332 ppm), Cr (249 ppm), Ni (162 ppm), Cu (236 ppm) and Sc (28 ppm), relatively low (La/Yb)N ratios between 1.92 and 39.00 with an average of 10.62, depleted rare earth elements, and also in high field strength elements (Nb, Ta, Ti), which are suitable signatures of arc volcanic rocks; - U–Pb zircon LA-ICP-MS geochronology on felsic metavolcanic or metatuff fix the maximum age for the volcanic activity in the RBGB at ∼670 Ma, thus younger than in the neighbouring Poli and Bibémi-Zalbi belts; - the dating of detrital zircon indicate that the main source for the detritus that fed the RBGB is Neoproterozoic (Ediacaran to Cryogenian) with minor Paleoproterozoic and Mesoproterozoic inputs from the neighbouring Mayo Kebbi-Sinassi magmatic arc and Adamawa-Yade Domain; - the age of low grade metamorphism weakly recorded by younger metamorphic overgrowth rims around older magmatic cores, characterized by 232Th/238U ratios lower than 0.2 is 600 Ma; - the maximum depositional age of the RBGB, corresponding to the youngest graphical age peak controlled by multiple grain ages is consistently constrained between 645 and 630 Ma. Our data therefore provide new insights into the geodynamic processes during the late Neoproterozoic, suggesting thattheRBGB, where high-Mg andesite magmatism has taken place is consistent with a continental arc setting related to a subduction zone.

Source:
M. Houketchang Bouyo, Y. Zhao, J. Penaye, S.H. Zhang, U.O. Njel, 2015.Neoproterozoic subduction-related metavolcanic and metasedimentary rocks from the Rey Bouba Greenstone Belt of north-central Cameroon in the Central African Fold Belt: New insights into a continental arc geodynamic setting
Precambrian Research 261:40-51.



Barrières géologiques : mieux comprendre leur comportement

Du nanomètre au kilomètre… Un saut d’échelle qu’il faut maîtriser pour prédire le comportement des barrières naturelles ou ouvragées. L’objectif des études dites "multi-échelles" est, à terme, de mieux comprendre les propriétés de transfert de ces barrières afin de mieux maîtriser leurs usages, et ce, pour une meilleure protection environnementale.
Étudier les propriétés d’une barrière géologique naturelle, ou  fabriquée par l’homme (ouvragée), permet d’optimiser ses conditions d’utilisation, et ce dans un but de protection de l’environnement. Un défi, vu le caractère "multi-échelle" des objets étudiés et de leurs propriétés ! C’est tout l’enjeu de travaux récents de spécialistes en environnement du BRGM. Objectif : mieux comprendre à petite échelle pour mieux prédire à plus grande échelle et sur le long terme.

Applications : stockage de déchets et séquestration dans le sous-sol

Ce défi n’est pas des moindres. Il peut s’agir, par exemple, d’étudier les propriétés de transfert d’une barrière qui s’étend sur  plusieurs kilomètres carrés et une centaine de mètres d’épaisseur. Et ce pour prédire son  évolution sur des périodes allant jusqu’à cent mille ans !
L’intérêt  de  tels  travaux est  majeur.  Les applications de stockage de déchets (toxiques ou  radioactifs) ou  de séquestration (du CO2) en milieu géologique reposent, en  grande partie, sur  les  propriétés des  barrières qui empêchent ou retardent la migration des produits confinés. Dans  le cadre de concepts d’utilisation énergétique du sous-sol, pour le stockage de chaleur, de pression ou de gaz, ces mêmes propriétés doivent assurer, sur  toute la durée d’exploitation, la conservation quantitative et qualitative du vecteur énergétique entre les périodes de production et de consommation.
Observation de la microstructure d'un matériau argileux (illite en marron) ayant subi une perturbation par précipitation in situ de célestine (sulfate de strontium en bleu). © BRGM

Prédire le comportement des barrières

Clé du succès de ces applications : notre capacité à prédire le comportement à moyen et long termes des propriétés des  barrières, en fonction des  perturbations physiques et  chimiques qu’elles subiront durant  leur période d’exploitation.
C’est un vrai  challenge scientifique. En effet,  ce type de  prédiction nécessite de  modéliser l’évolution  des systèmes de barrières de la façon la plus fiable possible, ce qui demande d’appréhender la complexité multi-échelle des  processus physico-chimiques mis  en  jeu.  Et ce, du nanomètre, la taille  représentative des  composants des barrières, minéraux argileux ou phases cimentaires, au kilomètre, la taille représentative des formations géologiques étudiées.  Mais  également de la nanoseconde, l’échelle de  temps sondée pour quantifier certains processus fondamentaux comme la diffusion moléculaire,  jusqu’au million d’années, l’échelle utile pour le stockage des déchets radioactifs.
Pour cela, le BRGM, dans le cadre de ses activités en partenariat avec l’Andra, a développé de nouvelles approches, qui  intègrent expérimentation et  modélisation, dans le but de décrire les propriétés des  matériaux poreux étudiés sur un continuum d’échelles spatiales. Leur originalité repose en partie sur  l’amélioration des méthodes permettant  d’observer directement la structure et  la minéralogie des  matériaux telles qu’elles se présentent in situ, avec des résolutions d’ordre nano à micrométrique. Des observations rendues possibles grâce  à l’utilisation intensive des  outils de  dernière génération en  microscopie électronique, et qui  nous permettent de mieux contraindre et fiabiliser nos modèles à grandes échelles. En suivant cette démarche, le BRGM espère gagner en confiance dans ses prédictions à long terme, en réduisant autant que possible le caractère empirique de la détermination des paramètres d’entrée des modèles.
Source: www.brgm.fr

lundi 23 mai 2016

Geomagnetic detection of the sectorial solar magnetic field and the historical peculiarity of minimum 23–24


Risultati immagini per Geomagnetic detection of the sectorial solar magnetic field[1] Analysis is made of the geomagnetic-activity aa index covering solar cycle 11 to the beginning of 24, 1868–2011. Autocorrelation shows 27.0-d recurrent geomagnetic activity that is well-known to be prominent during solar-cycle minima; some minima also exhibit a smaller amount of 13.5-d recurrence. Previous work has shown that the recent solar minimum 23–24 exhibited 9.0 and 6.7-d recurrence in geomagnetic and heliospheric data, but those recurrence intervals were not prominently present during the preceding minima 21–22 and 22–23. Using annual-averages and solar-cycle averages of autocorrelations of the historical aa data, we put these observations into a long-term perspective: none of the 12 minima preceding 23–24 exhibited prominent 9.0 and 6.7-d geomagnetic activity recurrence. We show that the detection of these recurrence intervals can be traced to an unusual combination of sectorial spherical-harmonic structure in the solar magnetic field and anomalously low sunspot number. We speculate that 9.0 and 6.7-d recurrence is related to transient large-scale, low-latitude organization of the solar dynamo, such as seen in some numerical simulations. Citation: Love, J. J., E. Joshua Rigler, and S. E. Gibson (2012), Geomagnetic detection of the sectorial solar magnetic field and the historical peculiarity of minimum 23–24, Geophys. Res. Lett., 39, L04102, doi:10.1029/2011GL050702. 1. Introduction [2] Recurrent geomagnetic activity changes over the course of each solar cycle. This can be understood in terms of the different phases of the Sun’s dynamo cycle [e.g., Solanki et al., 2006] and the oscillatory exchange of energy between the Sun’s toroidal and poloidal magnetic field ingredients. At solar maximum, the Sun’s magnetic field is primarily a toroidal quadrupole. Buoyancy brings toroidal field up to the solar surface, and with its emergence through the photosphere, sunspot groups are formed. During the declining phase of the solar cycle, energy shifts from the toroidal field to a poloidal dipolar field. This leads to a diminishment in sunspot number, an increase in broad regions of open coronal magnetic field lines and outflowing solar wind, corresponding to each end of the strengthening dipole, and the development of an organized heliospheric current sheet, corresponding to the dipolar magnetic equator. If the dipole is tilted with respect to the heliographic axis [e.g., Suess, 2008, Figure 7.1], then with 27.0-d synodic Carrington rotation of the Sun, the solar wind forms a heliospheric current sheet whose intersection with the heliographic equator is an Archimedean spiral [e.g., Smith, 2008]. This current sheet can potentially cross the Earth twice per solar rotation. Solar wind from the north and south heliospheric sectors drives low-level geomagnetic activity and small magnetic storms having 13.5-d lagged recurrence [e.g., Mursula and Zieger, 1998]. More often, however, the Sun is not quite so symmetrical. Either one dipolar end emits more solar wind than the other, or high-speed streams of solar wind are emitted from semi-isolated coronal holes that can persist for several months. As a result, the dominant interval for geomagnetic-activity recurrence is 27.0 d [e.g., Tsurutani et al., 2006]. At sunspot minimum, the Sun’s poloidal dipole field reaches its greatest strength and is roughly aligned with the Sun’s rotational axis. During this time, recurrent geomagnetic activity is mild but detectable. For review of the toroidal-poloidal decomposition, see Chandrasekhar [1961, Appendix III]. [3] In contrast to this idealized description, the declining phase of cycle 23 and the depth of the ensuing minimum 23–24 were unusual [e.g., Russell et al., 2010]. There were fewer sunspots in 2008 than in any year since 1913, minimum 14–15. Solar-wind data collected by the ACE satellite in 2005 record semi-persistent 9.0 and 6.7-d recurrence intervals [Temmer et al., 2007], corresponding to the third and fourth harmonics of the synodic solar-rotational period. These same harmonics have been identified in geomagnetic activity and thermospheric density [Lei et al., 2008; Thayer et al., 2008], in auroral electrons [Emery et al., 2009], and in relativistic radiation belt electrons [Gibson et al., 2009], each for years near minimum 23–24. These observations have been interpreted in terms of low-latitude coronal holes [de Toma, 2012] and multipolar ingredients in the solar magnetic field [Abramenko et al., 2010; DeRosa et al., 2010] that gave structural complexity to the heliospheric current sheet [McComas et al., 2006; Hathaway and Suess, 2008]. [4] Curiously, none of the above studies show prominent 9.0 and 6.7-d recurrence in data covering minima 21–22 and 22–23, leading to the perception that minimum 23–24 was unusual. This motivates our study. To establish whether or not minimum 23–24 was truly unusual in terms of recurrence, analysis of data covering many earlier solar cycles is required. Sunspots are not directly useful here, since they do not record recurrent phenomena. Fortunately, measures of geomagnetic activity, obtained from ground-based geomagnetic observatories [e.g., Love, 2008], are continuous in time from the middle of the 19th century to the present. Using autocorrelation methods, we analyze 13 solar-cycle minima of recurrent geomagnetic activity, 1868–2011, cycle 11 through to the beginning of 24. On the basis of comparisons between the historical geomagnetic results and other modern .


1 Geomagnetism Program, U.S. Geological Survey, Denver, Colorado, USA. 2 High Altitude Observatory, NCAR, Boulder, Colorado, USA.

data recording related solar-terrestrial phenomena for recent solar cycles, we can draw inferences about the long-term behavior of the solar magnetic field and better understand just how unusual minimum 23–24 actually was. 2. Data [5] We analyze 5 different time series. (1) The three-hour geomagnetic-activity aa index, for years 1868–2011, from cycle 11 to the beginning of 24, is derived from a pair of ground-based observatories [Mayaud, 1980]. This index forms the basis of most of the analysis presented here. (2) Sunspot group numbers G [Hoyt and Schatten, 1998], 1976–2011, are a qualitative measure of variable solar activity. (3) Solar-wind velocity V and (4) the radial BX component of the interplanetary magnetic field, for 2006 and 2008, were measured by the ACE [Stone et al., 1999] satellite, 1.5 million km from the Earth and toward the Sun on the Sun-Earth line. (5) Wilcox Observatory potential-field models of the coronal magnetic field, 1976–2011, are fitted to magnetogram data with radial boundary conditions at the photosphere (1.0R⊙) and in the corona (2.5R⊙) [Altschuler and Newkirk, 1969; Wang and Sheeley, 1992]. These models consist of spherical-harmonic coefficients, but they can also be shown as synoptic maps. 3. Geomagnetic Autocorrelation 1991–2011 [6] To measure geomagnetic-activity recurrence, we calculate Pearson autocorrelations r(l) of the aa index as a function of integer-day lag l. The computer algorithm [Press et al., 1992, Chapter 14.5] is applied to 100.0-d over-lapping time-series segments of daily averages of aa, thus identifying recurrence that is persistent over a 100.0-d duration of time, or slightly more than three Carrington rotations. In Figure 1, we show annual averages of the autocorrelations for 1991–2011. During solar-cycle rise and maximum, such as 1998–2001, many magnetic storms result from the sporadic occurrence of coronal-mass ejections, and so there are few obvious features in the autocorrelation curves. But during solar-cycle decline and minimum, such as 1993–1996, cycles 22–23, 27.0-d recurrence is seen as distinctive peaks. Smaller peaks represent 13.5-drecurrence, for example, during 1995. There is a hint of a 6.7-d interval in 1992, but it does not persist for lags much greater than 27.0 d. [7] Of more interest, here, are autocorrelations for the declining phase of cycle 23 and minimum 23–24. The year 2008 (blue) shows 9.0-d geomagnetic-activity recurrence intervals, corresponding to the third harmonic of synodic solar rotational, and for 2006 (blue), there is a 6.7-d recurrence interval, the fourth harmonic of solar rotation; these can be compared with the power spectra of Thayer et al. [2008, Figure 2]( see on the original document) . It is fair to say that minimum 23–24 was different from 22–23. This assessment is consistent with that of Emery et al. [2009, Figure 5], who identified differences in the harmonic content of solar wind data between these two minima. Annual-average autocorrelation plots for years 1868– 1993 are in the auxiliary material; a panoramic inspection of aa autocorrelation across many solar cycles is made in Section 6.1 4. Examples of 9.0-d and 6.7-d Recurrence [8] In Figure 2( see on the original document)  we plot 6 Carrington rotations of solarterrestrial data for 2008 showing 9.0-d recurrence. The Wilcox synoptic maps of the radial coronal magnetic field at 2.5R⊙ show a warped “heliomagnetic equator”, corresponding to the heliospheric current sheet that divides solar magnetic hemispheres of opposite polarities. For days 165–290, with Figure 1. Annual averages of Pearson autocorrelation r(l) of the geomagnetic-activity aa index as a function of integer-day lag l, 1991–2011. Results for 2006 and 2008 are shown in blue and should be compared with Figures 2 and 3. The amplitude scale is given in the upper right-hand corner, and the horizontal gray line for each autocorrelation shows its zero-level baseline.
1 Auxiliary materials are available in the HTML.

each solar rotation, outward flowing solar wind on either side of the kink in the magnetic equator is measured by ACE as peaks in V separated by 9.0 d, followed by 18.0-d gaps of slower solar wind. The changing sign of ACE interplanetary magnetic field BX indicates passage from one sector to another. After day 290, a slightly more regular, every-9.0-d pattern emerges. As for geomagnetic activity measured by aa, a corresponding 9.0 and 18.0-d recurrence results from solar wind-magnetospheric coupling. [9] Similar data are shown in Figure 3. Here, the current sheet is especially scalloped, having a tidy 90°-sectorial structure. During this time, and especially for days 238–319, with each solar rotation, 6.7-d recurrence is seen in solar wind velocity V and geomagnetic activity aa; 13.5-d recurrence is seen in interplanetary magnetic field BX. These observations, and those for Figure 2, are not of a heliosphere that is just generically non-axisymmetric, with a current sheet having random warps here and there. Instead, they show that the heliosphere can be organized in its asymmetry. Extending the inference made by Mursula and Zieger [1998], who focussed on the tilted dipole, the heliosphere near solar-cycle minimum is shaped by the Sun’s low-degree sphericalharmonic poloidal field, and this can be detected in geomagnetic-activity recurrence intervals. 5. Sectorial Solar Magnetic Field [10] We refer these observations to solar-cycle variation of the solar magnetic field. For each degree l and order m, we denote the Wilcox coronal magnetic field, radial at 2.5R⊙, as Blm 2.5. Spherical integration gives an energy spectrum,

 Elm ¼ 1 4p I 4p B2:5 lm ⋅B2:5 lm sinqdqdf;                          (1)

where q is colatitude and f is longitude. From this we estimate the relative proportion of energy per degree and order,

Plm ¼ Elm ∑jk Ejk :          (2)

The spherical-harmonic ingredients contributing the greatest amount of nonaxisymmetry in the solar magnetic field are sectorial, for which m = l. For example, a sectorial quadrupolar field has four equatorial patches of open field (2 of each sign), and this can drive 6.7-d geomagnetic-activity recurrence. In Figure 4 ( see on the original document)  we show time series of the proportion of dipolar P11 and quadrupolar P22 sectorial energy, 1976–2011; we also show, as a superposition, annual-average sunspot number G. Figure 4a is related to plots of dipole tilt seen in many papers, and it is equivalent to Hoeksema [2009], his Figure 2c divided by his Figure 2a. ( see on the original document). Otherwise, the sectorial quadrupole energy shown in Figure 4b( see on the original document)  is different from that shown in other work (contrast with Abramenko et al. [2010, Figure 5] and DeRosa et al. [2010, Figure 4]). [11] In Figure 4,( see on the original document)  it is important to note that the correlation between both P11 and P22 with G is high from cycle minimum 20–21 (1976) until maximum 23 (2000). Afterwards and into minimum 23–24 (2008), there is a departure, with sectorial energy remaining elevated during the declining phase of cycle 23, while sunspot numbers diminish and reach minimum. This combination is a phenomenological basis for recent 9.0 and 6.7-d aa recurrence: (1) non-axisymmetric, and especially sectorial, ingredients in the solar magnetic field give a non-axisymmetric heliosphere that drives geomagnetic recurrence, and (2) low sunspot numbers correspond to relatively few coronal-mass ejections and relatively little sporadic geomagnetic activity that would otherwise obscure measures of periodic recurrence. 6. Secular Change 1868–2011 [12] To put the preceding observations into a long-term context, in Figure 5 we show solar-cycle averages of geomagnetic aa autocorrelation r(l). Each average is taken over a duration extending from one sunspot maximum to the next, where, at each maximum, the axial dipole has a strength of approximately zero. Thus each average encloses a period of geomagnetic recurrence during solar-cycle decline and minimum. In the same figure, we also show the long-term average autocorrelation taken across all cycles (orange). Again, distinctive peaks correspond to 27.0-d recurrence, but we see, now, that the amplitude of recurrence has slowly changed from one cycle to another [Sargent, 1985]. In contradiction to Rangarajan [1991], 13.5-d recurrence is not always present. It is, however, seen for averages 15–16 and 16–17; anomalous autocorrelation, defined as the difference between an individual average and the long-term average, is shown for 16–17 (red). We have inspected the annualaverage autocorrelations for each year since 1868 (auxiliary material); 13.5-d recurrence is seen for some isolated years, such a 1895, 1922, and 1942, and, prominently, for the consecutive years of 1929 and 1930 [Newton, 1931]. Apparently, some solar polarity transitions are accomplished in a way that includes a tilting of the poloidal dipole, while other transitions result more from the diminishment of the axial dipole and its reappearance with the opposite polarity. With respect to historically quiet years, 1901 had an annualaverage sunspot number of 2.5 and extremely low geomagnetic activity levels; aa autocorrelation shows a very faint 6.7-d recurrence for 1901, not nearly of the amplitude for 2006. Perhaps the 1901 solar magnetic field did not have very prominent sectorial ingredients. 7. Discussion and Dynamo Context [13] With respect to the recent minimum 23–24 (2000– 2011), Figure 5 and the material in the auxiliary material clearly show that it was unusual – even, “peculiar”. None of the 12 preceding minima for cycles 11–23 show prominent 9.0 and 6.7-d aa recurrence. It is the nonaxisymmetric heliosphere that drives such recurrence, and since the heliosphere is controlled by the Sun, it is reasonable to conclude that the recent minimum 23–24 was distinguished by the solar dynamo obtaining a state of unusual asymmetry. Since solar convection is highly supercritical [e.g., Miesch and Toomre, 2009], a wide and continuous range of turbulent lengthscales, all shorter than the radius of the convection zone, is normal. But recent numerical simulations show that small-scale stellar convection cells can have larger-scale, low-latitude organization [Dikpati and Gilman, 2005; Brown et al., 2008] that can be described in terms of low-degree, sectorial spherical harmonics. We speculate that the transient development of such large-scale nonaxisymmetric organization within the Sun leads to nonaxisymmetric structure in the solar poloidal field and in the heliosphere. This might be a physical explanation for sectorial structure in the solar magnetic field that can lead to 9.0 and 6.7-d geomagnetic-activity recurrence, such as seen during minimum 23–24. [14] Acknowledgments. We thank: (1) The British Geological Survey and Geoscience Australia for observatory data, (2) the Wilcox Solar Observatory for coronal magnetic field models, (3) the ACE Science Center for solar wind and interplanetary magnetic field data, and NASA’s OmniWeb team for making the data available, (4) NOAA’s National Geophysical Data Center for archiving sunspot group number. We thank C. A. Finn, J. L. Gannon, M. S. Miesch, and K. Mursula for reviewing a draft manuscript, and V. Courtillot, J. T. Hoeksema, and G. de Toma for help and conversations. This work was supported by the US Geological Survey and the National Center for Atmospheric Research, which is supported by the National Science Foundation. [15] The Editor thanks Edward Smith and an anonymous reviewer for their assistance in evaluating this paper.

Jeffrey J. Love,1 E. Joshua Rigler,1 and Sarah E. Gibson2
Source: www.geology.usgs.gov

mercredi 11 mai 2016

Earthquake north of Kumamoto: M6.2 - 7km SW of Ueki, Japan

Tectonic Summary

The April 14, 2016 M 6.2 earthquake north of Kumamoto, on the island of Kyushu in southwest Japan, occurred as the result of strike-slip faulting at shallow depth. Focal mechanisms for the earthquake indicate slip occurred on either a left-lateral fault striking to the northwest, or on a right-lateral fault striking northeast. While the earthquake occurred several hundred kilometers northwest of the Ryukyu Trench, where the Philippine Sea plate begins its northwestward subduction beneath Japan and the Eurasia plate, the shallow depth and faulting mechanism of this earthquake indicate it occurred on a crustal fault within the upper Eurasia plate. At the location of this event, the Philippine Sea plate converges with Eurasia towards the northwest at a velocity of 58 mm/yr.
Moderate-to-large, shallow earthquakes in Kyushu are infrequent – most seismicity in the region is related to the subduction of the Philippine Sea plate at depth. Thirteen M 5+ earthquakes have occurred at shallow depths (< 50 km) within 100 km of the April 2016 event over the past century. In January and April of 1975, two shallow events with magnitudes of M 5.8 and M 6.1 - 40 km and 65 km to the northwest of the April 2016 earthquake, respectively – caused injuries, but no known fatalities. A shallow M 6.6 earthquake in March 2005, just off the north coast of Kyushu and 110 km north of the April 2016 event, caused over 1000 injuries and at least one fatality.
Mapped faults in the region generally trend east-west or northeast-southwest, in agreement with the right-lateral plane of preliminary focal mechanisms, and the trend of early aftershocks. In the first three hours after the M 6.2 event (12:26:36 UTC), 7 aftershocks have been located, the largest of which is a M 6.0 event at 15:03:47 UTC.

Seismotectonics of the Philippine Sea and Vicinity

The Philippine Sea plate is bordered by the larger Pacific and Eurasia plates and the smaller Sunda plate. The Philippine Sea plate is unusual in that its borders are nearly all zones of plate convergence. The Pacific plate is subducted into the mantle, south of Japan, beneath the Izu-Bonin and Mariana island arcs, which extend more than 3,000 km along the eastern margin of the Philippine Sea plate. This subduction zone is characterized by rapid plate convergence and high-level seismicity extending to depths of over 600 km. In spite of this extensive zone of plate convergence, the plate interface has been associated with few great (M>8.0) ‘megathrust’ earthquakes. This low seismic energy release is thought to result from weak coupling along the plate interface (Scholz and Campos, 1995). These convergent plate margins are also associated with unusual zones of back-arc extension (along with resulting seismic activity) that decouple the volcanic island arcs from the remainder of the Philippine Sea Plate (Karig et al., 1978; Klaus et al., 1992).
South of the Mariana arc, the Pacific plate is subducted beneath the Yap Islands along the Yap trench. The long zone of Pacific plate subduction at the eastern margin of the Philippine Sea Plate is responsible for the generation of the deep Izu-Bonin, Mariana, and Yap trenches as well as parallel chains of islands and volcanoes, typical of circum-pacific island arcs. Similarly, the northwestern margin of the Philippine Sea plate is subducting beneath the Eurasia plate along a convergent zone, extending from southern Honshu to the northeastern coast of Taiwan, manifested by the Ryukyu Islands and the Nansei-Shoto (Ryukyu) trench. The Ryukyu Subduction Zone is associated with a similar zone of back-arc extension, the Okinawa Trough. At Taiwan, the plate boundary is characterized by a zone of arc-continent collision, whereby the northern end of the Luzon island arc is colliding with the buoyant crust of the Eurasia continental margin offshore China.
Along its western margin, the Philippine Sea plate is associated with a zone of oblique convergence with the Sunda Plate. This highly active convergent plate boundary extends along both sides the Philippine Islands, from Luzon in the north to the Celebes Islands in the south. The tectonic setting of the Philippines is unusual in several respects: it is characterized by opposite-facing subduction systems on its east and west sides; the archipelago is cut by a major transform fault, the Philippine Fault; and the arc complex itself is marked by active volcanism, faulting, and high seismic activity. Subduction of the Philippine Sea Plate occurs at the eastern margin of the archipelago along the Philippine Trench and its northern extension, the East Luzon Trough. The East Luzon Trough is thought to be an unusual example of a subduction zone in the process of formation, as the Philippine Trench system gradually extends northward (Hamburger et al., 1983). On the west side of Luzon, the Sunda Plate subducts eastward along a series of trenches, including the Manila Trench in the north, the smaller less well-developed Negros Trench in the central Philippines, and the Sulu and Cotabato trenches in the south (Cardwell et al., 1980). At its northern and southern terminations, subduction at the Manila Trench is interrupted by arc-continent collision, between the northern Philippine arc and the Eurasian continental margin at Taiwan and between the Sulu-Borneo Block and Luzon at the island of Mindoro. The Philippine fault, which extends over 1,200 km within the Philippine arc, is seismically active. The fault has been associated with major historical earthquakes, including the destructive M7.6 Luzon earthquake of 1990 (Yoshida and Abe, 1992). A number of other active intra-arc fault systems are associated with high seismic activity, including the Cotabato Fault and the Verde Passage-Sibuyan Sea Fault (Galgana et al., 2007).
Relative plate motion vectors near the Philippines (about 80 mm/yr) is oblique to the plate boundary along the two plate margins of central Luzon, where it is partitioned into orthogonal plate convergence along the trenches and nearly pure translational motion along the Philippine Fault (Barrier et al., 1991). Profiles B and C reveal evidence of opposing inclined seismic zones at intermediate depths (roughly 70-300 km) and complex tectonics at the surface along the Philippine Fault.
Several relevant tectonic elements, plate boundaries and active volcanoes, provide a context for the seismicity presented on the main map. The plate boundaries are most accurate along the axis of the trenches and more diffuse or speculative in the South China Sea and Lesser Sunda Islands. The active volcanic arcs (Siebert and Simkin, 2002) follow the Izu, Volcano, Mariana, and Ryukyu island chains and the main Philippine islands parallel to the Manila, Negros, Cotabato, and Philippine trenches.
Seismic activity along the boundaries of the Philippine Sea Plate (Allen et al., 2009) has produced 7 great (M>8.0) earthquakes and 250 large (M>7) events. Among the most destructive events were the 1923 Kanto, the 1948 Fukui and the 1995 Kobe (Japan) earthquakes (99,000, 5,100, and 6,400 casualties, respectively), the 1935 and the 1999 Chi-Chi (Taiwan) earthquakes (3,300 and 2,500 casualties, respectively), and the 1976 M7.6 Moro Gulf and 1990 M7.6 Luzon (Philippines) earthquakes (7,100 and 2,400 casualties, respectively). There have also been a number of tsunami-generating events in the region, including the Moro Gulf earthquake, whose tsunami resulted in more than 5000 deaths.
Source: http://geology.usgs.gov/