M Wildman

R Brown

D Chardon

D Rouby

A N Kouamelan

Dall ' Asta

Contrasting Thermal Evolution of the West African Equatorial and Central Atlantic Continental Margins

Keywords: Thermochronology, Apatite fission-track, Transform Margin, Tectonics, West Africa

The landscape of the West African continental margins is the product of tectonic, thermal and surface processes acting in concert during and following the breakup of Gondwana.

Central

Atlantic opening was marked by the emplacement of the Central Atlantic Magmatic Province (CAMP) and continental breakup proceeded through Late Jurassic and Early Cretaceous divergent tectonics while opening of the Equatorial Atlantic was dominated by Early and mid-Cretaceous transform movement. The onshore erosional response to these events is poorly constrained yet is a crucial component of our understanding of topographic evolution and sediment transfer across continental margins. We present new apatite fission-track (AFT) data for 24 samples from Guinea and 11 samples from Ivory Coast, and thermal histories from inverse modelling. Our data and thermal histories show the following: the thermal effect of the CAMP across Guinea and Ivory Coast; rapid cooling along the coast during the early to mid-Cretaceous related to erosion of short-wavelength rift-shoulders; moderate cooling across longer wavelengths reflecting a pattern of erosion across flexural margin upwarps; and low cooling rates from the start of the Cenozoic to present day, consistent with low magnitudes of erosion inferred by onshore geomorphological data. We present our results alongside the published regional AFT dataset and draw inferences on the thermal and tectonic evolution of the onshore margin.

Introduction

Divergent margins are first-order features of global plate tectonics. Constraining the erosional history across the onshore domain of these margins is critical for understanding the response of the Earth's surface to tectonic and climatic processes over million-year timescales [START_REF] Summerfield | Geomorphology and Global Tectonics[END_REF][START_REF] Séranne | South Atlantic continental margins of Africa: a comparison of the tectonic vs climate interplay on the evolution of equatorial west Africa and SW Africa margins[END_REF]. The Guinean-Ivory Coast sector of the West African continental margins (Fig. 1) presents a complex case-study due to the multi-phase history of rifting and interaction of a predominantly Jurassic extensional system in the Central Atlantic along the northern Guinean margin [START_REF] Biari | Opening of the central Atlantic Ocean: implications for geometric rifting and asymmetric initial seafloor spreading after continental breakup[END_REF] and an Early Cretaceous transform dominated system in the Equatorial Atlantic along the Ivory Coast margin (Basile et al., 2005).

Advances in understanding rifting chronology and the development of offshore sedimentary basins along the margins has revealed links between the margins' tectonic and geodynamic evolution and vertical motions (i.e., uplift and subsidence) across both the onshore and offshore domains [START_REF] Attoh | Seismic stratigraphic record of transpression and uplift on the Romanche transform margin, offshore Ghana[END_REF][START_REF] Burke | The African erosion surface: a continental-scale synthesis of geomorphology, tectonics, and environmental change over the past 180 million years[END_REF][START_REF] Grimaud | Cenozoic sediment budget of West Africa and the Niger delta[END_REF][START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF]. Recent work involving mapping of dated lateritic palaeo-landforms (e.g., [START_REF] Beauvais | Modes, tempo, and spatial variability of Cenozoic cratonic denudation: The West African example[END_REF][START_REF] Grimaud | Very long-term incision dynamics of big rivers[END_REF][START_REF] Chardon | Stabilization of large drainage basins over geological time scales: Cenozoic West Africa, hot spot swell growth, and the Niger River[END_REF] and inverting river profiles [START_REF] Lodhia | Observation and Simulation of Solid Sedimentary Flux: Examples From Northwest Africa[END_REF] has revealed very low magnitude denudation and spatially variable uplift from c. 45 Ma and has been interpreted within the context of basin-and-swell mantle driven uplift, long-term sea-level fall, and climate changes during this time. This data has been linked to offshore sediment accumulation data to constrain source-to-sink relationships across the margin during the Cenozoic [START_REF] Grimaud | Cenozoic sediment budget of West Africa and the Niger delta[END_REF][START_REF] Lodhia | Observation and Simulation of Solid Sedimentary Flux: Examples From Northwest Africa[END_REF].

However, the paucity of Mesozoic geological markers means that the tectonic and geomorphic evolution of West Africa from the onset of Central Atlantic rifting in the late Jurassic [START_REF] Labails | An alternative early opening scenario for the Central Atlantic Ocean[END_REF] to the Cenozoic remains poorly understood. Moreover, there remains a lack of thermochronological data that, in the absence of these markers, can provide insights into periods of thermal or tectonic activity, enhanced erosion, and/or burial by constraining rock thermal histories.

The apatite fission-track (AFT) thermochronometer has been used globally to constrain denudation and burial along continental margins (see [START_REF] Wildman | Fission-track thermochronology applied to the evolution of passive continental margins[END_REF] for a review) by providing data on the timing and rate of rock cooling. Apatite fission-track data in West Africa has been limited to onshore studies in Burkina Faso [START_REF] Gunnell | Radiometric ages of laterites and constraints on long-term denudation rates in West Africa[END_REF], Ghana [START_REF] Lisker | Denudation and uplift across the Ghana transform margin as indicated by new apatite fission track data[END_REF][START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF], Benin [START_REF] Wildman | Long-term evolution of the West African transform margin: estimates of denudation from Benin using apatite thermochronology[END_REF], and from offshore wells along a marginal ridge in the Ivory Coast Basin [START_REF] Clift | Apatite fission track analysis of Sites 959 and 960 on the transform continental margin of Ghana, West Africa[END_REF][START_REF] Bigot-Cormier | Denudation of the Côte d'Ivoire-Ghana transform continental margin from apatite fission tracks[END_REF]. We present new AFT data from Guinea and Ivory Coast (Fig. 1b, Fig. 2) and obtain thermal histories for these regions using Bayesian inverse modelling to investigate the thermal response of the upper crust along the Guinean and Ivory Coast margin segments in response to superimposed rifting. We also discuss our data alongside previously published AFT datasets from West Africa to draw more regional conclusions on the tectonic evolution of the region.

2.

West African geological setting 2.1 Summary of rifting and break-up along the West African margin.

The West African continental margin is segmented by transform faults with the extensional Central Atlantic Ocean segment North of the Guinea fracture zone (FZ) [START_REF] Biari | Opening of the central Atlantic Ocean: implications for geometric rifting and asymmetric initial seafloor spreading after continental breakup[END_REF] (Fig. 1a). The Guinea-Liberia margin, Ivory Coast-Ghana margin and Togo-Benin sub-margins are transtensional systems creating pull-apart oblique normal faults, splay faults, and step-over and extensional duplex structures (Nemcok et al. 2016) bounded by the major St. Paul FZ, Romanche FZ and Chain FZ transform faults, respectively [START_REF] Ye | Paleogeographic and structural evolution of northwestern Africa and its Atlantic margins since the early Mesozoic[END_REF][START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF][START_REF] Basile | Phanerozoic geological evolution of the Equatorial Atlantic domain[END_REF][START_REF] Basile | Transform continental margins-Part 1: Concepts and models[END_REF][START_REF] Moulin | A new starting point for the South and Equatorial Atlantic Ocean[END_REF][START_REF] Heine | Kinematics of the South Atlantic rift[END_REF]. Transpressional structures have also been mapped in offshore Ghana and Ivory Coast [START_REF] Tari | Traditional and new play types of the offshore Tano Basin of Cote D'Ivoire and Ghana, West Africa[END_REF]. The structural variations across each of these segments and the onshore and offshore geological record are testament to the prolonged and multi-phase history of rifting [START_REF] Ye | Paleogeographic and structural evolution of northwestern Africa and its Atlantic margins since the early Mesozoic[END_REF][START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF].

In the Central Atlantic, the emplacement of the extensive Central Atlantic Magmatic Province (CAMP) in the early Jurassic is linked to thermal and tectonic processes preceding and during continental rifting [START_REF] Nomade | Chronology of the Central Atlantic Magmatic Province: implications for the Central Atlantic rifting processes and the Triassic-Jurassic biotic crisis[END_REF][START_REF] Greenroyd | Demerara plateau-The structure and evolution of a transform passive margin[END_REF][START_REF] Buiter | A review of Wilson Cycle plate margins: A role for mantle plumes in continental break-up along sutures?[END_REF]. As observed along other margins (e.g., [START_REF] Clemson | Structural segmentation and the influence of basement structure on the Namibian passive margin[END_REF][START_REF] Gibson | Pre-existing basement structure and its influence on continental rifting and fracture zone development along Australia's southern rifted margin[END_REF][START_REF] Peace | The role of pre-existing structures during rifting, continental breakup and transform system development, offshore West Greenland[END_REF], reactivation of pre-existing structures likely focused the location of rifting [START_REF] Attoh | The role of Pan-African structures in intraplate seismicity near the termination of the Romanche fracture zone, West Africa[END_REF][START_REF] Mercier De Lépinay | Transform continental margins-Part 2: A worldwide review[END_REF]. Rifting was predominantly symmetrical within the Central Atlantic, involving normal faulting and the formation of rift basins [START_REF] Withjack | Development of the passive margin of eastern North America: Mesozoic rifting, igneous activity, and breakup[END_REF][START_REF] Biari | Opening of the central Atlantic Ocean: implications for geometric rifting and asymmetric initial seafloor spreading after continental breakup[END_REF] until seafloor spreading commenced around 180-200 Ma [START_REF] Labails | An alternative early opening scenario for the Central Atlantic Ocean[END_REF][START_REF] Mercier De Lépinay | Transform continental margins-Part 2: A worldwide review[END_REF]. While the CAMP coverage was extensive across West Africa and Northeast America [START_REF] Marzoli | Synchrony of the Central Atlantic magmatic province and the Triassic-Jurassic boundary climatic and biotic crisis[END_REF], intercontinental deformation at this time was minor [START_REF] Biari | Opening of the central Atlantic Ocean: implications for geometric rifting and asymmetric initial seafloor spreading after continental breakup[END_REF].

Along the Equatorial Atlantic are alternating zones of divergent and transform tectonics [START_REF] Basile | Transform continental margins-Part 1: Concepts and models[END_REF]. Oblique transform rifting via an eastward propagating system of strike-slip faults initiated in the Early Cretaceous (c. 140 Ma) forming fault bounded pull-apart basins [START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF]. The Guinean Plateau, which sits at the transition between Central and Equatorial Atlantic systems, exhibits evidence of crustal thinning related to the initial Jurassic extension and deformation linked to later Cretaceous transform movement (Mercier de [START_REF] Mercier De Lépinay | Transform continental margins-Part 2: A worldwide review[END_REF]. This suggests the Guinea-Liberia margin initially began rifting during Central Atlantic extension followed by the main transform rift phase during the Early Cretaceous [START_REF] Bennett | Regional 2D seismic interpretation and exploration potential of offshore deepwater Sierra Leone and Liberia, West Africa[END_REF][START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF].

Movement along E-W to ENE-WSW intra-continental transform faults continued through the Valanginian to Aptian (140-113 Ma) [START_REF] Heine | Kinematics of the South Atlantic rift[END_REF]. The end of continental rifting of the Equatorial Atlantic is marked by a Late Albian Breakup Unconformity (107 -100 Ma) mapped and correlated across diachronous oceanic crust from the Sierra Leone segment to the Togo-Benin segment of the margin [START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF]. At the time of break-up, South American plate rotation shifted the tectonic regime of the Ivory Coast-Ghana margin segment from transtension to transpression [START_REF] Attoh | Seismic stratigraphic record of transpression and uplift on the Romanche transform margin, offshore Ghana[END_REF]) and an angular breakup unconformity formed beneath Late Albian-Cenomanian Marine sediments [START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF][START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF].

The transition to seafloor spreading is marked by a so-called 'active-transform' stage from the Turonian to Santonian (c. 94-84 Ma) as oceanic spreading continues along transform faults [START_REF] Basile | Phanerozoic geological evolution of the Equatorial Atlantic domain[END_REF]Nemcok et al., 2016;[START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF][START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF].

The tectonic processes that occurred during the formation of the West African Transform Margin has made this sector the type-example for this form of continental break-up [START_REF] Basile | Transform continental margins-Part 1: Concepts and models[END_REF]. The margin exhibits a classic geometry with an outer corner, where the continental margin extends out towards the oceanic accretion axis, and inner corner, where the margin extends inwards towards the continent, created by the divergent margin segments being intersected and bound by major transform faults (Fig. 1a). There is also the presence of a transform fault parallel marginal ridge along the western side of the Romanche FZ where maximum strike-slip deformation has occurred (Fig. 1a, [START_REF] Basile | Transform continental margins-Part 1: Concepts and models[END_REF] and transform marginal plateaus in the Guinean, Liberian and Côte d'Ivoire-Ghana offshore segments [START_REF] Loncke | Transform marginal plateaus[END_REF].

Regional geology

The Guinea-Ivory Coast sector of West Africa is dominated by the Archean to Palaeoproterozoic granitoids and greenstones of the Leo-Man Shield (Fig. 1b). In Guinea, Archean gneisses in the west are separated from the younger Palaeoproterozoic rocks in the east by a series of N-S and NE-SW trending shear zones [START_REF] Rollinson | Archaean crustal evolution in West Africa: A new synthesis of the Archaean geology in Sierra Leone, Liberia, Guinea and Ivory Coast[END_REF]. The trend of these shear zones reflects the prevailing structural trend in Ivory Coast, with basement structures accommodating transcurrent displacements during the Palaeoproterozoic orogeny [START_REF] Lompo | Paleoproterozoic structural evolution of the Man-Leo Shield (West Africa). Key structures for vertical to transcurrent tectonics[END_REF]. A more complex structural fabric is observed in southeast Guinea where numerous faults crosscut one another and are the product of pre-Cambrian and Pan-African deformation accommodating mostly strike-slip movements [START_REF] Guiraud | Phanerozoic geological evolution of Northern and Central Africa: an overview[END_REF] (Fig. 1a,b).

Across West Africa, Pan-African basin inversion and thrusting is manifested as two distinct events [START_REF] Villeneuve | Paleozoic basins in West Africa and the Mauritanide thrust belt[END_REF]. The first Pan-African event occurred at c. 660 Ma and formed the Mauritanides and Bassarides Belt, with the Mauritanides being reworked during Late Palaeozoic Hercynian deformation [START_REF] Guiraud | Phanerozoic geological evolution of Northern and Central Africa: an overview[END_REF]. The second event at c. 550-500 Ma produced the NNW-SSE trending Rokelides Belt, which extends through southwestern Guinea, beneath the Bove Basin, and runs adjacent to the coastline through Sierra Leone and Liberia (Villeneuve and Cornee, 1994;[START_REF] Guiraud | Phanerozoic geological evolution of Northern and Central Africa: an overview[END_REF][START_REF] Deynoux | Pan-African tectonic evolution and glacial events registered in Neoproterozoic to Cambrian cratonic and foreland basins of West Africa[END_REF] (Fig. 1a).

The Neoproterozoic to Cambrian Madina-Kouta Basin, comprised of alternating continental and shallow marine sandstones, siltstones and shales, overlies the Leo-Man shield in northern Guinea and extends eastward through Mali and Burkina Faso along the southern margin of the larger Taoudeni Basin [START_REF] Villeneuve | Review of the orogenic belts on the western side of the West African craton: the Bassarides, Rokelides and Mauritanides[END_REF]Ennih and Liégeois, 2008) (Fig 1b).

The Bové Basin in northern Guinea (Fig. 1a) sits on top of the pre-Cambrian basement and Pan-African belts and is comprised of Cambrian-Devonian fluvial-deltaic sandstones and conglomerates, Silurian marine shales with interbedded sandstones, and alternating sequences of Devonian marine sandstones and shales [START_REF] Villeneuve | Paleozoic basins in West Africa and the Mauritanide thrust belt[END_REF][START_REF] Deynoux | Pan-African tectonic evolution and glacial events registered in Neoproterozoic to Cambrian cratonic and foreland basins of West Africa[END_REF]. Northwest of Guinea, in Guinea-Bissau, Upper Cretaceous sediments are present and belong to the larger Senegal Basin (Fig. 1a,b), which extends north along the Central Atlantic margin and west into the offshore domain [START_REF] Brownfield | Assessment of the undiscovered oil and gas of the Senegal Province, Mauritania, Senegal, the Gambia, and Guinea-Bissau, Northwest Africa[END_REF]. These rocks are primarily marine shales, siltstones and sandstones, which form the primary hydrocarbon reservoirs and seals in the Senegal Basin [START_REF] Brownfield | Assessment of the undiscovered oil and gas of the Senegal Province, Mauritania, Senegal, the Gambia, and Guinea-Bissau, Northwest Africa[END_REF][START_REF] Davison | Central Atlantic margin basins of North West Africa: geology and hydrocarbon potential (Morocco to Guinea)[END_REF].

Emplacement of the CAMP (Fig. 1) is believed to have occurred during a short-lived (c. < 1 Myr) period of peak magmatism at c. [START_REF] Deckart | Geochemistry and Sr, Nd, Pb isotopic composition of the Central Atlantic Magmatic Province (CAMP) in Guyana and Guinea[END_REF]. Lava flows can be found and traced across sedimentary basins; however, they are thin (<500 m) and rarely preserved [START_REF] Marzoli | The Central Atlantic magmatic province (CAMP): a review[END_REF]. Outcrop in Guinea is predominantly in the form of layered intrusions forming caps on elevated regions. Estimated pre-erosional volumes of the CAMP have been given in the range of 2-3x10 6 km 3 [START_REF] Mchone | Volatile emissions from Central Atlantic Magmatic Province basalts: Mass assumptions and environmental consequences[END_REF][START_REF] Svensen | Gondwana Large Igneous Provinces: plate reconstructions, volcanic basins and sill volumes[END_REF]) (Fig. 1b).

Topography and geomorphology

The topography in Guinea is dominated by the Guinea Rise topographic massif [START_REF] Chardon | Stabilization of large drainage basins over geological time scales: Cenozoic West Africa, hot spot swell growth, and the Niger River[END_REF] that trends NW-SE through Guinea, reaching elevations >1000 m capped by lowrelief plateaus (Fig. 2). Inward from the coast, a series of stepped escarpments mark rapid increases in elevation (up to 700 m) in northern Guinea. Across Sierra Leone and Liberia, a 100-200 km wide low-relief, coastal plain with a gentle elevation increase up to c. 200 m is observed. Along the Guinea sector of the Central Atlantic margin, the coastal plain is much narrower < 100 km, with a region of high relief and elevations up to c. 600 m, following the structural trend of the Rokelide Belt, separating the coastal region from the Guinea Rise. Four hundred kilometres inland, the topography drops to a low-relief, gently inland sloping region with elevations of c. 400 m (Fig. 2). In the southern part of the Guinea Rise, the elevated plateau is fragmented with the eastern crest of the plateau dissected in a series of N-S trending ridges following the cratonic shear zone. Further south, the relief, comprised of peaks and scattered plateaus, is controlled by the NE trending structural grain with maximum elevations reaching 1948m on Mont Loma in NE Sierra Leone. Maximum elevations across Ivory Coast are lower than along the Central Atlantic margin. Elevation gently increases from sea-level at the coast to c. 400-600 m at 500 km inland.

The Guinean Rise is one of three major topographic features (including the Hoggar swell in the north and Jos Plateau to the east), that control the large-scale regional drainage of West Africa [START_REF] Chardon | Stabilization of large drainage basins over geological time scales: Cenozoic West Africa, hot spot swell growth, and the Niger River[END_REF]. These features form part of the continental-scale basin-and-swell topography [START_REF] Burke | The African plate[END_REF], which is proposed to be a product of long-term surface processes in response to the break-up of Gondwana, post-rift tectonics and intraplate mantle driven dynamic uplift [START_REF] Burke | The African plate[END_REF][START_REF] Chardon | Stabilization of large drainage basins over geological time scales: Cenozoic West Africa, hot spot swell growth, and the Niger River[END_REF]. The drainage is divided by [START_REF] Grimaud | Cenozoic sediment budget of West Africa and the Niger delta[END_REF] into 4 main domains: (i) Senegambia (ii) Short Atlantic drainages, (iii) Long Atlantic drainages, and (iv) the Niger catchment (Fig. 2). The existing drainage network is suggested to have been relatively stable since the early Oligocene (29 Ma) or possibly earlier at the end of the Eocene (34 Ma) [START_REF] Chardon | Stabilization of large drainage basins over geological time scales: Cenozoic West Africa, hot spot swell growth, and the Niger River[END_REF].

Reconstruction of the Cenozoic drainage evolution has been achieved through a series of studies using 40 Ar-39 Ar dating of K-rich Mn oxides and regionally correlated remnant lateritic palaeo-landforms across West Africa (e.g., [START_REF] Vasconcelos | 40Ar39Ar analysis of supergene jarosite and alunite: Implications to the paleoweathering history of the western USA and West Africa[END_REF][START_REF] Beauvais | Chemical and physical erosion rhythms of the West African Cenozoic morphogenesis: the 39Ar-40Ar dating of supergene K-Mn oxides[END_REF][START_REF] Beauvais | Modes, tempo, and spatial variability of Cenozoic cratonic denudation: The West African example[END_REF][START_REF] Grimaud | Very long-term incision dynamics of big rivers[END_REF][START_REF] Chardon | Stabilization of large drainage basins over geological time scales: Cenozoic West Africa, hot spot swell growth, and the Niger River[END_REF]. Several weathering sequences were identified and include the bauxitic "African Surface", which developed under humid conditions from the Late Cretaceous to c. 45 Ma, the ferricrete capped "Intermediate" surface, which developed between 29-24 Ma, and three Neogene pediment systems abandoned at c. 11, 6 and 3 Ma, respectively, known as the High, Middle and Low glacis systems [START_REF] Beauvais | Modes, tempo, and spatial variability of Cenozoic cratonic denudation: The West African example[END_REF][START_REF] Chardon | Stabilization of large drainage basins over geological time scales: Cenozoic West Africa, hot spot swell growth, and the Niger River[END_REF][START_REF] Chardon | West African lateritic pediments: landform-regolith evolution processes and mineral exploration pitfalls[END_REF][START_REF] Grimaud | Very long-term incision dynamics of big rivers[END_REF][START_REF] Grimaud | Cenozoic sediment budget of West Africa and the Niger delta[END_REF].

Dissection of these landforms have enabled facilitated estimates of dissection driven denudation to be made with rates being extremely low since the Eocene (c. 2-10 m/Myr, mean values of 5-7 m/Myr) [START_REF] Beauvais | Modes, tempo, and spatial variability of Cenozoic cratonic denudation: The West African example[END_REF][START_REF] Grimaud | Very long-term incision dynamics of big rivers[END_REF][START_REF] Grimaud | Cenozoic sediment budget of West Africa and the Niger delta[END_REF].

2.4

Previous thermochronology studies.

The new AFT data we present here compliments and adds to existing thermochronology data across the West African Equatorial margin (Fig. 1a). Gunnel (2003) presented central AFT ages from Burkina Faso, ranging from 175 ± 10 to 218 ± 7 Ma, and thermal histories implying slow monotonic cooling since the Mesozoic across the West African craton. In Ghana, [START_REF] Lisker | Denudation and uplift across the Ghana transform margin as indicated by new apatite fission track data[END_REF] and [START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF] present AFT data from samples along the coastline and continental interior, respectively. AFT ages from basement samples along the southern Ghanian coastline range from 130 to 415 Ma [START_REF] Lisker | Denudation and uplift across the Ghana transform margin as indicated by new apatite fission track data[END_REF]. Thermal history modelling of this data infers two stages of cooling during the Palaeozoic, following the Pan-African Orogeny, and Cretaceous, related to exhumation following continental rifting. Inland, central AFT ages range from 65 ± 11 to 298 ± 18 Ma [START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF] (2005) are younger than the samples Lower Cretaceous stratigraphic ages, they advocate significant sedimentary burial of these samples followed by erosion during the 'activetransform' stage.

Apatite fission track thermochronology

AFT Methods

Apatite fission-track ages were obtained from 24 samples from Guinea and 11 from the Ivory Coast using the zeta-calibration external detector method [START_REF] Hurford | The zeta age calibration of fission-track dating[END_REF].

Apatite fission-track ages were combined with horizontal confined track length measurements to constrain thermal annealing of tracks through the AFT partial annealing zone (PAZ, 60-110 ± 10°C) [START_REF] Gleadow | Fission track lengths in the apatite annealing zone and the interpretation of mixed ages[END_REF]. The angle of measured track lengths to the crystallographic caxis was measured to correct for anisotropic annealing [START_REF] Ketcham | Improved modeling of fission-track annealing in apatite[END_REF] and Dpar (length of the etch pit formed by tracks intersecting the polished grain surface) was measured as a proxy for the bulk compositional influence on track annealing [START_REF] Carlson | Variability of apatite fission-track annealing kinetics: I. Experimental results[END_REF][START_REF] Donelick | Apatite fission-track analysis[END_REF]. Radial plots and track length distributions for all samples can be found in Supplementary Data, Fig. S1 and Fig. S2, respectively.

Apatite Fission-Track Results

Guinea

Twenty-four samples were collected from Guinea (Table 1; Supplementary Data, Table S1)

and can be divided into two regional groups (Fig. 2). In the northwest, 15 samples (Guinea North) form a coast-perpendicular transect that extends from the coast, through the elevated Guinean Rise to c. 400 km inland. Nine samples (Guinea South) cover the south-eastern region of the Guinean Rise at c. 800-1000 m elevation. Central AFT ages across the entire Guinea dataset range from 99.1 ± 13.0 to 204.2 ± 34.3 Ma. Nine samples fail the P(χ 2 ) age homogeneity test (i.e., P(χ 2 ) < 0.05) with four of these samples showing excess single grain age dispersion (>20%) than that expected for a single population (see radial plots in Supplementary Data, Fig. S1). In most of the Guinea samples, Dpar measurements were collected across the entire sample (i.e., not on individual dated grains) and so the relationship between single grain age and Dpar cannot be assessed. In those samples where Dpar was measured on dated grains (Supplementary Data, Fig. S1) no relationship between Dpar and age is apparent suggesting excess scatter cannot be solely attributed to apatite composition influence annealing. Given the high number of tracks for the samples in the Guinean dataset, it is possible that the chi-sq test is penalising the data for being too precise (e.g., [START_REF] Vermeesch | On the treatment of discordant detrital zircon U-Pb data[END_REF]. Mean track lengths (MTL), corrected for their c-axis orientation (c-MTL) range from 13.57 ± 0.12 to 14.78 ± 0.08 with standard deviations of 0.78 to 1.47.

The Ivory Coast dataset contains 11 samples (Fig. 2) (Table 1; Supplementary Data, Table S1) forms a broadly coast perpendicular transect from the coast to c. 500 km inland with a gradual elevation gain from sea-level to 500 m. Approximately 350 km from the coast, two samples (YOR-1 and YAL-1) fall to the west of the main transect along an outer ridge of the Guinean rise. Central AFT ages in Ivory Coast range from 114.7 ± 13.1 to 188.9 ± 22.5 Ma, c-MTLs range from 13.99 ± 0.14 to 14.39 ± 0.07 μm with unimodal track length distributions and MTL-SD of 0.74-1.03. Five samples fail the P(χ 2 ), however, only one of these samples show excessive dispersion (>20%). As suggested above for the Guinean samples, apatite composition variability may have an influence on creating excess single grain age dispersion but the high track counts in these grains may be a contributing factor causing the P(χ 2 ) to fail.

Our West African data set yields a relatively narrow range in MTLs 13.57 ± 0.12 to 14.78 ± 0.08 over a reasonably large age range 99.1 ± 13.0 to 204.2 ± 34.3 Ma, which does not show any clear correlation (Fig. 3a). The negative correlation between MTL and MTL St. Dev. (Fig. 3b) reflects shorter MTLs having a broader track length distribution.

Inland from the coast, AFT ages in the Ivory Coast and the northern Guinea increase gradually to a distance of c. 200 km (Fig. 3c). In northern Guinea, this distance corresponds to the region of increasing topography from sea-level at the coast to the highest sampled elevations in the Guinean Rise (i.e., samples GN-06 to GN-08) (Fig. 2). The relationship to topography is less clear in Ivory Coast due to the modest elevation range across the entire country (Fig. 2). However, in general this c. 200 km region corresponds to the low-lying (< 200 m) coastal strip.

Beyond 200 km inland, AFT ages gradually decrease again (Fig. 3c). This is best observed in the Guinea North dataset, with three samples from Ivory Coast also plotting on this trend. In northern Guinea, the decreasing trend is coincident with a general decrease in elevation from peak elevation >1000m in the Guinean rise to c. 400 m inland (Fig. 2). Although some of the Ivory Coast data agree with the lateral trend in AFT age, the elevation beyond 200 km inland gradually increases to c. 400 m. Four samples from the Ivory Coast profile also appear to form a negative correlation from c. 350 to 500 km inland that is offset from the main data trend (Fig. 3c). This type of offset may be attributable to fault offsetting with down-to-the-northeast displacement and differential denudation across the fault. However, more samples would be required to better relate the spatial relationship of the AFT data to fault displacement. Samples from the Guinea South dataset, situated between 300 and 400 km inland, are fairly scattered.

Five samples appear to have AFT ages consistent with the main AFT age-distance from coast trend, while two seem more consistent with the offset trend defined by the Ivory Coast samples.

The Guinea North dataset shows the strongest AFT age-elevation correlation (Fig. 3d), which connects increasing and decreasing elevations across the Guinean Rise to the AFT agedistance from the coast relationship. The Guinea South data are more scattered but while they do not show a strong correlation, they appear to be consistent with the Guinea North AFT ageelevation trend. Several of the Ivory Coast samples fall on the trend defined by the Guinea

North data but at 350-400 m there is more scatter in the Ivory Coast AFT ages. The offset samples here are also the samples that form the offset part of the AFT age-distance from the coast trend for Ivory Coast.

4.

Thermal History Modelling

Modelling approach

Thermal histories for all samples were acquired using QTQt, which incorporates a Bayesian transdimensional Markov Chain Monte Carlo (MCMC) approach to data inversion [START_REF] Gallagher | Transdimensional inverse thermal history modeling for quantitative thermochronology[END_REF]. Apatite fission-track data were modelled using the multi-kinetic fissiontrack annealing model of [START_REF] Ketcham | Improved modeling of fission-track annealing in apatite[END_REF] with anisotropic annealing and compositional influences on the annealing rate taken into account using c-axis projected track lengths and sample average Dpar measurements, respectively. We present the expected thermal history (ExTH) and 95% credible intervals in figure 4 and summarise the modelling results below. We provide a more detailed discussion of the modelling strategy, preliminary testing, alternative model scenarios and plots of the observed data vs. model predicted data in Supplementary Data, Data S1, Table S2, Fig. S3-S9. All QTQt input files are provided in Supplementary Data, Dataset S1.

Our modelling strategy involved initially modelling all samples individually and without any additional constraints to generate thermal histories with complexity solely driven by the data of one sample (Supplementary Data, Fig. S3). We then look to combine samples with the aim of reinforcing mutually consistent features observed in individual sample thermal histories and making these parts of the thermal history more robust. We take this step because due to the relatively low precision of AFT thermochronology, the data for a single sample may not robustly constrain detail in the thermal history the model (i.e., timing of maximum heating or inflection points in the thermal history path). Moreover, thermal histories from single samples taken from locations near to one another may result in apparent regional variability in thermal histories between adjacent samples. These variations in thermal histories could be interpreted as having a geological meaning (e.g., migrating patterns of erosion or fault movement causing differential denudation). Before this conclusion can be drawn, it should be assessed whether data from several samples within a region can be explained by a single thermal history that is consistent with all the data and their uncertainties. A multi-sample thermal history may also identify a sample or samples that are inconsistent with data from adjacent samples, and this may indicate a more complex spatial-thermal history relationship.

We take a simplistic approach to grouping samples, grouping first into 100 km bins (starting from the coastline), and assessing whether the joint model suitably reproduces all of the grouped data. If some samples are not reproduced, we group the data into 50 km bins, and then, if necessary to fit the data, resort to the individually modelled sample (see Supplementary Data, Data S1, Fig. S4 for further details). We assign a temperature offset based on an assumed geothermal gradient of 25±5°C/km. However, it should be appreciated that these samples do not form a pseudo-vertical profile or are not taken from a borehole profile where the definition of a thermal offset would have more relevance. The grouped samples are typically at similar elevations, have the same present-day surface temperatures and assumed to have had similar palaeo-temperatures. The specified geothermal gradient and the range on this value is used to build some uncertainty into this assumption.

We do not explicitly define an initial starting constraint for the models and in some cases the output thermal history appears to begin at temperatures somewhere much colder than the base of the PAZ. In these circumstances it is inferred that rapid cooling from temperatures hotter than the PAZ has occurred immediately before the starting point time.

Thermal History Modelling Results

Guinea North

Samples GM01, GM02, GN1, GN3 and GN4 reside on the low elevation coastal margin zone that extends 150-200 km inland from the coastline and are modelled together (Fig. 4a).

The bauxite capped plateaus and pervasive Middle glacis relics in that area suggest that this part of the landscape experienced approximately 400 and 1000 m of denudation between 45 and 6 Ma, followed by negligible denudation to present-day. Assuming a geothermal gradient of 25±5°C/km from 45 Ma, and a surface temperature of 20°C, this equates to cooling from a temperature in the range of 28 to 50°C.

The geomorphological information can be used as additional constraints to force cooling from 45 and 6 Ma To prevent additional time-temperature points being proposed after 45 Ma, the lower time range of the prior was reduced to 45 Ma. The early portion of the ExTH (>70 Ma) shows the main cooling episode observed at c. 110 Ma. Heating of c. 20°C is then inferred between 70 and 45 Ma, before finally cooling to surface temperatures at present-day as constrained by the geomorphology (Fig. 4b). In supplementary information (Supplementary Data, Fig. S5) we explore the necessity for the Early Cenozoic cooling event be forbidding reheating and observe that the data can still be reproduced without reheating. However, as this is a very restrictive constraint to impose, and the inferred reheating is fairly minor we continue our discussion with model presented in Fig. 4a.

While some ambiguity remains in the ExTH model, we can have confidence in the consistent and recurrent parts of the history. Specifically, this includes an initial period of protracted cooling through the Late Jurassic and Early Cretaceous followed by a phase of rapid cooling in the mid-Cretaceous (110 -90 Ma). Some burial in the Late Cretaceous to Early Cenozoic prior to the samples cooling to the surface through the Late Cenozoic is possible, but not unequivocally required by the observed data. Regardless of the complexity of the Cretaceous and Early Cenozoic history, the samples remained at temperatures <50°C since the mid-Late Cretaceous (c. 90 Ma).

The thermal history for samples GN6, GN7 and GN8 shows cooling from 100°C at 190

Ma to 35°C at 100 Ma. From 100 Ma to present the expected model predicts around 10°C of cooling to surface temperatures, which is consistent with the low levels of denudation throughout the Cenozoic predicted by geomorphological studies [START_REF] Beauvais | Modes, tempo, and spatial variability of Cenozoic cratonic denudation: The West African example[END_REF][START_REF] Grimaud | Cenozoic sediment budget of West Africa and the Niger delta[END_REF].

The GN09, GN10, GM12 thermal history reproduces the data well and shows cooling at a moderate rate between 100 and 40°C from 170 Ma to 120 Ma before low levels of cooling from 120 Ma to present-day. The GM13 single sample (Fig. 4a) also shows low amounts of cooling through low temperatures (i.e., <40°C) from c. 120 Ma but has experienced more rapid cooling at c. 140 Ma.

Furthest from the coast, GM14, GN11 and GN12 yield a model where the top sample in the profile cools from 100°C at 170 Ma to 60°C at 135 Ma (Fig. 4a). The ExTH then predicts protracted cooling to bring the samples to present-day surface temperatures. The ExTH reproduces the data reasonably well and the path is consistent with the geomorphological constraints in the region.

Guinea South

All models across the Guinea South dataset (Fig. 2) show a similarly protracted cooling history (Fig. 4b). All samples cool through the base of the AFT PAZ during the Late Triassic to Late Jurassic, with the exception of GN18, which cooled through this boundary during the mid-Permian. Although slow cooling is predominant, there is some variation in the cooling rate through the Cretaceous (c. 150-70 Ma) with some samples showing slower cooling and cooling through lower temperatures than others. For example, GN18 experiences cooling from 61 to 42°C at a rate of 0.24°C/Myr between 150 and 100 Ma, compared to the joint model of GN31

and GN33 cooling from 90 to 42°C at a rate of 0.6°C/Myr.

Ivory Coast

Two samples, SA14 and SA3, were collected from the coastline and when modelled together produce an ExTH that shows rapid cooling in the Early Cretaceous between 140 and 110 Ma (Fig. 4c), followed by minimal cooling until present day. The wide 95% credible intervals in the late Cretaceous imply that the initial cooling may have been greater than the expected model shows, but it is also possible that these samples experienced some minor burial related reheating (<60°C). The DL03BisA and SOK1A ExTH model shows cooling through the base of the PAZ at c. 200 Ma, followed by a moderate cooling rate of c. 0.6-0.7°C/Myr through the Late Jurassic and Early Cretaceous, before the cooling rate drops quickly over the last 110 Ma (Fig. 4c).

Further inland, BEZ-VAA yields an ExTH model showing relatively moderate cooling through the PAZ during the Cretaceous with only minor, slow cooling through the Cenozoic (Fig. 4c). The MK64-SEN ExTH model shows cooling through the PAZ during the Late Triassic and Jurassic (c. 215-150 Ma), and residence at near surface temperatures from the Cretaceous to the present-day (Fig. 4c). The ExTH model for the furthest sample from the coast, KG33, shows rapid cooling between 150 and 140 Ma followed by low-temperature (<40°C) protracted cooling from the Cretaceous to present-day.

Two samples (YAL and YOR) sit west of the main transect (Fig. 2) on the older palaeoto meso-Archean basement. YOR resided at the foot of an escarpment, which marks a jump in elevation from c. 300 m to up to 1000 m. YAL is sampled inland of these maximum elevations, c. 30 km northwest of YOR, from an elevation of 511 m. The individual models for these samples show an initially rapid phase of cooling at two distinct times. For YAL, cooling through the PAZ is inferred to occur between 170 and 155 Ma, while YOR infers rapid cooling from temperatures hotter from the base of the PAZ to c. 70°C at 120 Ma (Fig. 4c). As these samples are collected from a faulted region marking the boundary between the older Archean basement and the younger Archean and Proterozoic basement (Fig. 1), fault reactivation may have caused the distinct thermal histories.

Discussion

Spatial patterns of palaeo-temperature and denudation

Using the ExTH for the models presented above, we can infer the palaeotemperature pattern across the West Africa margin during specific time periods (Fig. 5). Assuming a geothermal gradient over time then allows us to estimate the amount of denudation or burial that may have occurred over time (Fig. 6). In Figure 6 we highlight our denudation/burial estimates assuming a standard geotherm of 25°C/km and use values derived from this geotherm in our discussion. However, we acknowledge that the geothermal gradient was likely elevated during continental breakup along passive and transform margins [START_REF] Balázs | Contrasting transform and passive margin subsidence history and heat flow evolution: insights from 3D thermo-mechanical modelling[END_REF] and the time of CAMP emplacement. We also acknowledge that the geothermal gradient more generally can vary over time and space due to deep thermal processes and due to the addition and removal of material with different thermal properties (e.g., Łuszczak et al., 2017) with implications for the temporal and spatial pattern of denudation. To reflect this, in Figure 6 we show a range for denudation/burial estimates for geothermal gradients between 20-70°C/km.

The main observation that can be made using this range of values is that if geothermal gradients were higher and decayed over time, the total amount of material removed in the late Jurassic -Cretaceous would have been significantly less that that predicted using a gradient of 25°485 C/km. However, it is not clear how far an elevated geotherm would have reached in the context of our sample locations relative to the main rift zone, and how this would have evolved in tandem with post-rift surface, tectonic and thermal processes. Acquiring thermochronometric data from different depths from borehole profiles would potentially yield greater insights into the evolution of the geotherm.

At 200 Ma, MK64 and SEN reside at temperatures cooler that 120°C (Fig. 5). All other samples were either residing at temperatures hotter than 100 -120°C at 200 Ma or were brought to temperatures hotter than 120°C after 200 Ma but prior to cooling through the PAZ during the Cretaceous (Fig. 5). If these palaeotemperatures were solely due to burial, then this would equate to a total of 4 to 5 km of overburden that has been removed since the late Jurassic (for a geothermal gradient of (25°C/km)). Most samples show a progressive decrease in the amount of cooling and consequently in denudation from the Early Cretaceous to the Cenozoic (Fig. 6).

Some samples show slight increases in cooling in the Late Cenozoic, however, due to the low temperature nature of this cooling it is uncertain.

The emplacement of CAMP rocks may have caused the temperature of upper crustal rocks to increase by direct heating due to their proximity to intrusions, associated regional hydrothermal activity, or burial of under thick layers of extrusive rocks, and/or changes in the geothermal gradient during magmatism and changes in thermal conductivity of the rock blanketing (i.e., the extrusive volcanic rocks) the previously exposed basement. CAMP lavas would have covered the 'pre-rift' landscape, which likely had some Palaeozoic cover [START_REF] Hubbard | The Phanerozoic cover sequences preserved as xenoliths in the kimberlite of eastern Sierra Leone[END_REF]. Given the location of the Guinean samples within the CAMP it is possible that this event reset the AFT data at c. 200 Ma.

With the exception of samples YAL and YOR, the Ivory Coast samples lie out with the preserved CAMP intrusion field (Fig. 1b). If these samples were reset at c. 200 Ma with the Guinean samples, then it is likely it was due to burial under CAMP extrusive rocks or more regional hydrothermal activity at the time. MK64 and SEN have not been completely reset at this time, which may be due to their residence at a higher elevation and consequently not being covered by as thick a cover of volcanic rocks or less thick cover of Palaeozoic sediments.

However, as the AFT data do not record when the data were reset, only when they cooled below the closure temperature, burial of the samples entirely beneath a sedimentary cover during the Palaeozoic cannot be excluded as a possibility.

For both transects, the samples nearest to the coast reside at the highest temperatures at 150 Ma (Ivory Coast samples are hotter than the base of the PAZ) and show episodes of rapid cooling (Fig. 5). In Guinea North, this rapid cooling occurs during the 110-90 Ma interval, after an initially moderate phase of cooling, whereas in Ivory Coast rapid cooling occurs during 130-110 Ma. Both of these events involve cooling of c. 30-40°C at a rate of 1.5-2°C/Myr. For a geothermal gradient of 25°C/km, this cooling equates to a denudation thickness of 1.2-1.6 km (Fig. 6) at a rate of 60 to 80 m/Myr.

In Guinea North, the interior samples cool at a low to moderate rate through the Early to mid-Cretaceous (~0.5 -0.6°C/Myr over 150-90 Ma) before the cooling rate and the total amount of cooling rapidly decreases through the Late Cretaceous and into the Cenozoic (~0.1 -0.2°C/Myr over 70-0 Ma). Similar trends in cooling rate and total amount of cooling through the Cretaceous and the Cenozoic are also observed along the interior of the Ivory Coast transect.

Some of the initial cooling in the Early Cretaceous may be attributable to thermal relaxation following the CAMP event but it is also possible that a significant proportion of the cooling through the Early to mid-Cretaceous is attributable to exhumation driven by moderate erosion rates in response to continental break-up. This style of landscape evolution is comparable to that observed along passive continental margins (e.g., [START_REF] Gallagher | Denudation and uplift at passive margins: the record on the Atlantic Margin of southern Africa[END_REF][START_REF] Persano | Apatite (U-Th)/He age constraints on the development of the Great Escarpment on the southeastern Australian passive margin[END_REF][START_REF] Spotila | Origin of the Blue Ridge escarpment along the passive margin of Eastern North America[END_REF][START_REF] Balestrieri | Geomorphic development of the escarpment of the Eritrean margin, southern Red Sea from combined apatite fission-track and (U-Th)/He thermochronometry[END_REF][START_REF] Campanile | Sedimentation record in the Konkan-Kerala Basin: implications for the evolution of the Western Ghats and the Western Indian passive margin[END_REF][START_REF] Peulvast | Low post-Cenomanian denudation depths across the Brazilian Northeast: implications for long-term landscape evolution at a transform continental margin[END_REF].

Due to the temperature sensitivity of the AFT thermochronometer, we are not able to resolve any detail on thermal events colder than 60°C. However, our data remain informative as they require denudation to be less than c. 1-2 km (Fig. 6). In the absence of data or any additional information driving inflections (e.g., episodes of rapid cooling or reheating) in the low-temperature thermal history during the Cenozoic our models infer simple monotonic cooling. Our data and models are therefore consistent with observations by [START_REF] Beauvais | Modes, tempo, and spatial variability of Cenozoic cratonic denudation: The West African example[END_REF] and [START_REF] Grimaud | Very long-term incision dynamics of big rivers[END_REF][START_REF] Grimaud | Cenozoic sediment budget of West Africa and the Niger delta[END_REF] that denudation was limited to c. 2-10 m/Myr across this part of West Africa. Application of a lower temperature thermochronometer such as apatite (U-Th)/He would potentially provide additional insights on the timing and rate of Cenozoic cooling.

5.2

Thermal and tectonic processes along the West African margin.

Our new AFT data is an important contribution to a growing thermochronological record with coverage from the Guinean margin at along the Central Atlantic to the Benin margin along the Equatorial Atlantic (Fig. 7, Fig. 8). This regional dataset helps to better understand the spatial and temporal record of exhumation and the driving geological processes. All data from the West African transform margin are presented in an AFT age-MTL plot (Fig. 7a) and show a partial 'boomerang'-style relationship. 'Boomerang' relationships in AFT datasets are formed by a population of 'young' AFT ages and one of 'old' AFT ages, both with long MTLs, and a zone between these age populations where the AFT-MTL relationship follows a U-shaped trend.

The data can be interpreted as preserving an early (old) rapid cooling event and then recording a later (young) rapid cooling event, with intervening samples exhibiting the effects of partial resetting [START_REF] Green | On the thermo-tectonic evolution of Northern England: Evidence from fission track analysis[END_REF]. In the West African data set the younger peak is evident corresponding to long MTLs (>12.5 μm) with Late Jurassic and Early-to mid-Cretaceous (145 -90 Ma) AFT ages, signifying the importance of this period for the exhumation and cooling of West Africa.

However, the older peak is absent. This suggests that all samples across this region of West Africa experienced some thermal annealing after the Cambrian/Early Ordovician.

The data defining the 'younger' peak are scattered due to variation in the AFT ages and

MTLs. The long period of time over which samples yield long MTLs (Fig. 7a) may be indicative of temporal variations for the onset of cooling and variations in cooling rate across

West Africa in the Mesozoic. The variation of the MTLs may be directly related to cooling, but it should be noted that Fig. 7a presents lengths that have not been corrected for their orientation to the crystallographic c-axis (e.g., [START_REF] Ketcham | Improved modeling of fission-track annealing in apatite[END_REF], which may improve consistency in the length measurements [START_REF] Ketcham | Reproducibility of thermal history reconstruction from apatite fission-track and (U-Th)/He data[END_REF]. [START_REF] Gunnell | Radiometric ages of laterites and constraints on long-term denudation rates in West Africa[END_REF], [START_REF] Clift | Apatite fission track analysis of Sites 959 and 960 on the transform continental margin of Ghana, West Africa[END_REF], [START_REF] Bigot-Cormier | Denudation of the Côte d'Ivoire-Ghana transform continental margin from apatite fission tracks[END_REF] and [START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF] do not state that a c-axis correction was made to the length measurements. The absence of this additional data is common in older studies as annealing anisotropy and appropriate numerical corrections were not typically considered.

Some studies (e.g., [START_REF] De Grave | The thermo-tectonic history of the Song-Kul plateau, Kyrgyz Tien Shan: Constraints by apatite and titanite thermochronometry and zircon U/Pb dating[END_REF][START_REF] Van Ranst | New insights from low-temperature thermochronology into the tectonic and geomorphologic evolution of the south-eastern Brazilian highlands and passive margin[END_REF] question the extent of anisotropy in fossil fission tracks and whether it is appropriate to make a c-axis correction, and this consideration may have been adopted by [START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF]. [START_REF] Clift | Apatite fission track analysis of Sites 959 and 960 on the transform continental margin of Ghana, West Africa[END_REF] and [START_REF] Bigot-Cormier | Denudation of the Côte d'Ivoire-Ghana transform continental margin from apatite fission tracks[END_REF] also do not include compositional information (e.g., Dpar or Cl wt. %),

which is now commonplace in AFT studies. The absence of these data may help to explain some of the conflicts in timing in the thermal histories observed along the offshore marginal ridge by [START_REF] Clift | Apatite fission track analysis of Sites 959 and 960 on the transform continental margin of Ghana, West Africa[END_REF] and [START_REF] Bigot-Cormier | Denudation of the Côte d'Ivoire-Ghana transform continental margin from apatite fission tracks[END_REF].

An AFT age-elevation relationship is most apparent in the Guinean samples, where a clear positive correlation is observed with increasing AFT age over c. 1200 m of elevation gain (Fig. 7b). Despite AFT ages reaching up to c. 400 Ma in Benin, no age-elevation relationship is apparent, with samples residing in a fairly narrow range of elevations c. 200-400 m., suggesting prevailing lateral variations in erosion. It is also apparent that AFT ages greater than 200 Ma (i.e., CAMP emplacement) are only found in datasets in the eastern part of the West African margin (e.g., Benin and Ghana). This lends further support to the thermal influence of the CAMP on the Guinean and Ivory Coast datasets and subsequent cooling driven by progressive exhumation of the Guinea and Ivory Coast regions of the West African margin (Fig. 9a,b).

The absence of the CAMP influence on the eastern datasets has made it possible to better observe the spatial variation in response of surface processes to major tectonic processes (Fig. 9). The lack of CAMP influence is inferred due to the preservation of relatively older (i.e., pre-CAMP) AFT ages in Benin [START_REF] Wildman | Long-term evolution of the West African transform margin: estimates of denudation from Benin using apatite thermochronology[END_REF], Burkina Faso [START_REF] Gunnell | Radiometric ages of laterites and constraints on long-term denudation rates in West Africa[END_REF] and Ghana [START_REF] Lisker | Denudation and uplift across the Ghana transform margin as indicated by new apatite fission track data[END_REF][START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF] The jump in AFT ages in the range 110-130 Ma along a coastal strip in Benin to ages greater than 300 Ma further inland coincident with the location of the present-day continental drainage-divide at c. 375 km inland (Fig. 8), is comparable to models of passive margin landscape evolution, where erosion during escarpment retreat or plateau downwearing drives exhumation [START_REF] Wildman | Long-term evolution of the West African transform margin: estimates of denudation from Benin using apatite thermochronology[END_REF]. Due to the flexural isostatic response to erosion, denudation magnitudes can be several kilometres and is greatest in the region between the coast and the continental divide (Fig. 9) [START_REF] Gallagher | Denudation and uplift at passive margins: the record on the Atlantic Margin of southern Africa[END_REF][START_REF] Braun | A review of numerical modeling studies of passive margin escarpments leading to a new analytical expression for the rate of escarpment migration velocity[END_REF][START_REF] Wildman | Long-term evolution of the West African transform margin: estimates of denudation from Benin using apatite thermochronology[END_REF]. While the increase in AFT age is not as dramatic along the Guinean and Ivory Coast margins a similar scenario may be invoked. Along the Guinea North transect, the trend of increasing AFT age with distance from the coast, and with increasing elevation, may also mark the influence of a drainage divide at c. 200 km inland, where the oldest AFT ages are observed.

In Ivory Coast, AFT ages progressively increase from the coast to c. 350 km inland, where MK64 and SEN-1 preserve older AFT ages.

The patterns of long wavelength deformation, driving the exhumation patterns across the Central and Equatorial Atlantic, are consistent with the model of transform margin evolution described by [START_REF] Ye | Paleogeographic and structural evolution of northwestern Africa and its Atlantic margins since the early Mesozoic[END_REF] involving the formation and persistence of flexural margin upwarps (e.g., [START_REF] Gilchrist | Differential denudation and flexural isostasy in formation of rifted-margin upwarps[END_REF][START_REF] Gallagher | The onshore record of passive margin evolution[END_REF]) (Fig. 9). Within this model, the divergent Central Atlantic further North than our study area exhibits a longer wavelength, coast parallel, up-warp compared to the transform Equatorial Atlantic. The difference in up-warp wavelength is attributed to the age of the respective margins, with wavelength increasing for the older central Atlantic margin due to greater thermal relaxation and lithospheric strengthening [START_REF] Ye | Paleogeographic and structural evolution of northwestern Africa and its Atlantic margins since the early Mesozoic[END_REF]. The main pattern of denudation following rifting in the Equatorial Atlantic is extends c. 350 km inland from the coast from Guinea along to Benin. However, in Guinea we have no data >350 km from the coast along the Guinea North transect.

Across Guinea, Ivory Coast and Benin, samples inland of the divide infer Mesozoic cooling suggesting erosion of interior flanks of the marginal up-warp beyond the drainage divide (Fig. 9) that may have fed sediment into interior basins (see [START_REF] Ye | Paleogeographic and structural evolution of northwestern Africa and its Atlantic margins since the early Mesozoic[END_REF]. Ye et al.

(2017) attribute the long wavelength of the Ivory Coast upwarp to processes in the asthenosphere and highlight the emplacement timing (150 and 135 Ma) of Leo-Man shield kimberlitic province in support for chemical and physical changes in the lithospheric mantle, which drives vertical motions (e.g., [START_REF] Ault | Phanerozoic surface history of the Slave craton[END_REF][START_REF] Stanley | Erosion patterns and mantle sources of topographic change across the southern African Plateau derived from the shallow and deep records of kimberlites[END_REF]. However, the inference from our thermal histories that crustal cooling across the Ivory Coast was already underway in the late Jurassic-Early Cretaceous suggests a closer link to the tectonic and or thermal processes forming the up-warp along the divergent Guinean margin.

Thermal histories from samples from the coastline along Guinea North and Ivory Coast also support more rapid cooling events (Fig. 4a, Fig. 4c). These rapid cooling events may be attributed to more rapid exhumation due to erosion of short-wavelength uplift rift flanks. Along the Ivory Coast equatorial margin, the timing of this rapid exhumation is coeval with the main phase of transform faulting in the Equatorial Atlantic suggesting a causal link (Fig. 9c). Data from Benin [START_REF] Wildman | Long-term evolution of the West African transform margin: estimates of denudation from Benin using apatite thermochronology[END_REF] and Ghana [START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF][START_REF] Lisker | Denudation and uplift across the Ghana transform margin as indicated by new apatite fission track data[END_REF] causing differential exhumation across the structure (Fig. 9) [START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF]. Possibly fault reactivation or hydrothermal circulation can explain the offset ages along the Ghanian coast (Fig. 8) [START_REF] Lisker | Denudation and uplift across the Ghana transform margin as indicated by new apatite fission track data[END_REF].

The rapid exhumation along the Guinea coastline is approximately 80 Ma after the onset of seafloor spreading in the Central Atlantic and overlaps with the end of transform faulting in the Equatorial Atlantic (Fig. 9d), possibly suggesting some regional tectonic uplift in response to continental break-up in the Equatorial Atlantic. 9d,e). Reactivation of brittle faults during the syn-and post-rift phase has also been proposed along several type-example passive continental margins (e.g., [START_REF] Ksienzyk | From orogen to passive margin: constraints from fission track and (U-Th)/He analyses on Mesozoic uplift and fault reactivation in SW Norway[END_REF][START_REF] Cogné | Post-rift reactivation of the onshore margin of southeast Brazil: Evidence from apatite (U-Th)/He and fission-track data[END_REF][START_REF] Wildman | The chronology and tectonic style of landscape evolution along the elevated Atlantic continental margin of South Africa resolved by joint apatite fission track and (U-Th-Sm)/He thermochronology[END_REF]. Moreover, inversion during the Late Albian has been observed in structures along the Guinea plateau [START_REF] Ye | Paleogeographic and structural evolution of northwestern Africa and its Atlantic margins since the early Mesozoic[END_REF] and so fault reactivation along the Guinean coast, driving the rapid cooling between 110 and 90 Ma cannot be ruled out. Uplift and tectonic reactivation across the Central and Equatorial margins following continental break-up in the Equatorial Atlantic and through the active-transform phase is therefore possible but still poorly resolved by the AFT datasets.

The limitations on the temperature sensitivity of AFT thermochronology and the low magnitude of Mesozoic cooling make resolving the thermal history and the chronology of denudation and burial through the Cenozoic challenging. Along the coastal regions of the Central and Equatorial Atlantic there are remnants of Late Cretaceous sedimentary basins (i.e., embayments). However, the past extent of these basins, in terms of both their regional coverage and thickness, is uncertain. Our data from Guinea and Ivory Coast and that presented by [START_REF] Wildman | Long-term evolution of the West African transform margin: estimates of denudation from Benin using apatite thermochronology[END_REF] from Benin do not exclude burial of the coast (e.g., Fig. 9e) and along the fringes of interior basins (e.g., Fig. 9d) as a possibility but they do infer that any burial was not sufficient to cause significant track annealing and so the total amount of heating related to this burial was <60°C.

More generally across West Africa, the chronology of landscape evolution has been resolved using dated relict lateritic landforms [START_REF] Beauvais | Modes, tempo, and spatial variability of Cenozoic cratonic denudation: The West African example[END_REF][START_REF] Grimaud | Very long-term incision dynamics of big rivers[END_REF][START_REF] Chardon | Stabilization of large drainage basins over geological time scales: Cenozoic West Africa, hot spot swell growth, and the Niger River[END_REF]. These datasets all appear to converge on the conclusion that in West Africa there was greater tectonic stability and lower magnitudes of erosion during the Cenozoic and certainly since the Eocene. Surface uplift during the Cenozoic, driven by mantle upwellings, has been proposed as a mechanism to form the 'basin-and-swell' topography described across Africa (e.g., [START_REF] Burke | The African plate[END_REF]. The growth of the Hoggar Swell, to the northeast of our study area, since the mid-Eocene caused significant denudation and increase in sediment flux from the Niger-Benue catchment. Elsewhere in West Africa any minor surface uplift contribution of mantle driven upwelling did not trigger km-scale erosion. Instead, erosion would have been limited to enhanced dissection of successive weathering profiles.

The regional drainage network response to tectonic and isostatically driven surface deformation will determine how the eroded sediments are partitioned into offshore and interior basins (e.g., [START_REF] Ye | Paleogeographic and structural evolution of northwestern Africa and its Atlantic margins since the early Mesozoic[END_REF][START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF]. A source-to-sink analysis by [START_REF] Ye | Evolution topographique, tectonique et sédimentaire syn-à post-rift de la marge transformante ouest Africaine[END_REF] observe that in the Ivory Coast-Ghana basin the largest accumulation of siliciclastic sediments occurs during the Campanian (85-72 Ma). Over the mid-Late Cretaceous (104 -66 Ma), accumulation rates increase from 4.3 to 13.3 x 10 3 km 3 /Myr. At this time, most samples across the Guinea and Ivory Coast margins, except for the coastal samples in Guinea, show very low cooling rates, or rapidly decreasing cooling rates, reflecting slow, low-magnitude exhumation. These low cooling rates continue throughout the Cenozoic, while offshore accumulation rates first decrease in the Eocene to 3.3 x 10 3 km 3 /Myr before increasing to 12.3 x 10 3 km 3 /Myr in the Plio-Pleistocene, with a short-lived peak of 18.3 x 10 3 km 3 /Myr in the Early Oligocene (34-31 Ma). The general trend observed in the Ivory Coast -Ghana basin of high Late Cretaceous accumulation rates, low rates in the early Cenozoic and high rates in the Late Cenozoic is also observed in other basins in the Equatorial Atlantic suggesting a discordance between the timing of uplift driven erosion and peaks in offshore accumulation [START_REF] Ye | Evolution topographique, tectonique et sédimentaire syn-à post-rift de la marge transformante ouest Africaine[END_REF]. [START_REF] Ye | Evolution topographique, tectonique et sédimentaire syn-à post-rift de la marge transformante ouest Africaine[END_REF] conclude that the disagreement between the timing and magnitude of onshore denudation and the record of accumulation can be attributed to drainage reorganization and changes in size of paleo-catchments and/or post-depositional sedimentary redistribution on sediment supply.

Conclusions

The landscape of the West African continental margin has largely evolved in response to thermal, tectonic, and surface processes since the emplacement of the Central Atlantic Magmatic Province in the Late Jurassic. This event marked the beginning of a series of tectonic events including rifting and continental break-up in the Central Atlantic followed by oblique movement and break-up in the Equatorial Atlantic. We infer the following from our new AFT data, thermal histories and their association with previously published AFT data. The AFT data and thermal histories infer low cooling rates and total amount of cooling through the Late Cretaceous and Cenozoic suggesting low erosion rates prevailed, which is consistent with Cenozoic denudation estimates of 2-10 m/Myr observed from Geomorphology studies. However, due to the lower temperature limit of the AFT thermochronometer, sub-km scale denudation cannot be confidently resolved. Similarly, burial of the coast and along the fringes of interior basins during the Late Cretaceous is a possibility but the total amount of heating related to this burial was <60°C. Continental break-up in the Equatorial Atlantic, 5: 'Active-transform' phase [START_REF] Ye | Post-rift stratigraphic architectures along the African margin of the Equatorial Atlantic: Part I the influence of extension obliquity[END_REF] in Equatorial Atlantic. (Bigot-Cromier et al., 2005). Burial inferred inland in Benin along the margins of the Iullemmenden Basin [START_REF] Wildman | Long-term evolution of the West African transform margin: estimates of denudation from Benin using apatite thermochronology[END_REF] and possibly the northern slopes of the offshore marginal ridge [START_REF] Clift | Apatite fission track analysis of Sites 959 and 960 on the transform continental margin of Ghana, West Africa[END_REF]. (e) Denudation magnitudes are inferred to be low from the Late Cretaceous until present-day. Some interior regions continue to have low-moderate rates of denudation and another phase of intracontinental fault inversion may have occurred. (Galbraith, 2005). AFT ages are central ages calculated with 1σ standard error. 'GN' samples were etched using 5M

HNO3 after [START_REF] Gleadow | Thermal history of granitic rocks from western Victoria: A fission-track dating study[END_REF] and ages were calculated using a ζ = 338±38 with a standard CN-5 glass. All other samples were etched using 5.5M HNO3 after [START_REF] Donelick | Apatite fission-track analysis[END_REF] and ages calculated with ζ = 303.1±9.6 using a standard IRMM540 glass. Dispersion is the standard deviation of the single-grain ages as a percentage of their central age. Mean track lengths are corrected for their orientation to the c-axis after [START_REF] Ketcham | Improved modeling of fission-track annealing in apatite[END_REF]. Sample location and rock type information are found in Supplementary Data, Table S1.

Supplementary Data

Data S1: Further details on the thermal history modelling methodology

  with thermal histories inferring Late Triassic -Early Jurassic cooling attributed to post-CAMP thermal relaxation, early to mid-Cretaceous cooling (145 -90 Ma) due to rifting related exhumation, and Late Cretaceous-Cenozoic cooling due to exhumation driven by post-rift tectonic reactivation along the Bole-Nangodi shear zone (Fig. 1a). Further east, along the Benin sector of the margin Wildman et al. (2019) present AFT ages ranging from 106 ± 5 to 401 ± 45 Ma with single-grain apatite (U-Th)/He ages (corrected for alpha-ejection) ranging from 26 ± 4 to 500 ± 95. These data and associated thermal history modelling are used to advocate for relatively rapid exhumation driven cooling in the Early Cretaceous along the coastal margin and the interior Iullemmeden Basin (Fig. 1a) and very slow cooling of the region between the coast Iullemmeden Basin through the Palaeozoic. The thermochronometery study in Benin does not preclude heating related to burial along the coastal basins and interior basins during the mid-Cretaceous (110 -90 Ma) and Late Cretaceous to Early Cenozoic, respectively. Offshore, Clift et al. (1998) and Bigot-Cormier et al. (2005) present AFT data from sedimentary samples collected from Ocean Drilling Program borehole sites through and surface samples along the marginal ridge (Fig. 1a). Clift et al. (1998) present central AFT ages from the borehole sites ranging from 88 ± 4 Ma to 113 ± 4 Ma. Central ages form samples from the southern continental slope of the ridge, reported by Bigot-Cormier et al. (2005) have a comparable but slightly larger range from 64.2 ± 4.7 to 125.3 ± 10.6 Ma, with the majority of AFT ages less than 100 Ma. Thermal history modelling of the data presented by Clift et al. (1998) infers pre-depositional cooling, which represents evolution of the source region (i.e., the evolving continental margin), between 130 and 110 Ma in three samples and between 110 and 100 Ma in another. Minor post-depositional heating and fission-track annealing is inferred following deposition. Clift et al. (1998) suggest the first cooling event was driven by uplift and erosion during transform deformation and a later cooling event by a post-transform period of inversion, respectively. As the majority of the AFT ages presented by Bigot-Cromier et al.

  (Fig. 1; Fig. 8). Samples with older (> 200 Ma) AFT ages and relatively short MTLs in these studies reflect regions of the intracontinental craton that have experienced protracted cooling and slow rates and/or low magnitudes of erosion throughout the Palaeozoic and Mesozoic. However, Fernie et al. (2018) do present thermal histories that cool through the PAZ between 200 and 150 Ma, which they attribute to post-CAMP thermal relaxation.

  ) suggest cooling driven by exhumation related to rifting and associated processes occurred along the entire Equatorial Atlantic margin and further into the continental interior (e.g., Benue Trough in Ghana, Gao rift in Benin). The abrupt transition in AFT ages in Ghana across the Bole-Nangodi shear zone (Fig. 1a, Fig 8) invoke intracontinental reactivation, during the rift period,

  [START_REF] Fernie | Thermochronological insights into reactivation of a continental shear zone in response to Equatorial Atlantic rifting (northern Ghana)[END_REF] present central AFT ages as young as 65 ± 11 Ma and tentatively propose a Late-Cretaceous -early Cenozoic, minor cooling event, attributed to inversion during NNW-SSE shortening within western Africa driving a phase of enhanced exhumation. This timing would also be consistent with the timing of exhumation of the offshore marginal ridge proposed by[START_REF] Bigot-Cormier | Denudation of the Côte d'Ivoire-Ghana transform continental margin from apatite fission tracks[END_REF] (Fig.

  of the CAMP, either by direct heating or by burial under thick magmatic rocks, was significant across Guinea and Ivory Coast and sufficient to heat the presently exposed crust to temperatures >110°C. Further east, the absence of the CAMP facilitates preservation of older pre-rift AFT ages. (ii) Patterns of AFT data with distance from the coast are comparable to those expected for traditional models of passive margin evolution, where most erosion occurs coastward of a continental divide during the main phase of rifting. Coastward limbs of long-wavelength flexural upwarps (as described by Ye et al., 2017) in Guinea and Ivory Coast were eroded, at moderate rates, in in response to the opening of the Central Atlantic in the Late Jurassic-Early Cretaceous and in response to the onset of transform movement in the Equatorial Atlantic during the early to mid-Cretaceous (145 -90 Ma). Landward limbs inland of the upwarp crest also experienced suggest significant interior erosion and fed interior basins. (iii) Fast cooling along the Ivory Coast coastline between c. 130 and 110 Ma, is attributed to rapid erosion of a short-wavelength rift shoulder formed during the onset of transform faulting in the Equatorial Atlantic. The onset of fast cooling along the Guinean coast at 110-90 Ma is coeval with the timing of continental break-up in the Equatorial Atlantic. While some of the regional West African AFT data, particularly in Ghana, suggest fault reactivation of basement structures during the Cretaceous and potentially in the Early Cenozoic, it is poorly resolved elsewhere but remains possible (iv)

Figure 1 :

 1 Figure 1: (a) Tectonic map of West Africa showing major sedimentary basins, deformation belts,

Figure 2 :

 2 Figure 2: Topographic map of the Guinea-Ivory Coast region of West Africa showing major

Figure 3 :

 3 Figure 3: AFT data plots: (a) c-axis projected mean track length vs. Central AFT age, (b) c-axis

Figure 4 :

 4 Figure 4: Thermal history models for samples from (a) Guinea North (b) Guinea South, and (c)

Figure 5 :

 5 Figure 5: Topographic profiles for the Guinea North and Ivory Coast transects with sample

Figure 6 :

 6 Figure 6: Topographic profiles for the Guinea North and Ivory Coast transects with estimates

Figure 7 :

 7 Figure 7: West African AFT data plots: (a) Mean track length vs. Central AFT age. Measured

Fig. 8 :

 8 Fig. 8: Maps of West African published AFT data coloured to indicate AFT age. Sample source

Fig. 9 :

 9 Fig. 9: Maps of West Africa summarising tectonic, thermal, and surface evolution over the Late

Fig. S1 :

 S1 Fig. S1: Radial plots for single-grain AFT ages.

Fig. S2 :

 S2 Fig. S2: Apatite fission-track length distributions

Fig. S6 :

 S6 Fig. S6: Detailed models and data predictions for Guinea North inland samples.

  201 Ma followed by ongoing activity until c. 192 Ma based on 40 Ar/ 39 Ar and zircon 206 Pb/ 238 U ages (see Marzoli et al. 2018 for a review).

	Large
	volumes of mafic intrusive material were emplaced during this magmatism in the form of long,
	dense dyke swarms and voluminous sills

Table 1 :

 1 Apatite Fission Track data. ρs, ρi, ρd are track density of induced, spontaneous, dosimeter tracks. P(χ2) is p-value of the chi-sq age homogeneity test

Table S1 :

 S1 Location and rock type information for samples

Table S2 :

 S2 Summary of QTQt model set up

Acknowledgements

This study was funded by Total Energies R&D through the Transform Source-to-Sink Project (TS2P). We would like to thank David Mark Webster for assisting with sample collection and sample preparation during his PhD research and Robert McDonald for technical assistance at the University of Glasgow. We are indebted to the Direction Nationale de la Géologie (DNG) and particularly the late Nassirou Bah for facilitating fieldwork and sampling in Guinea. Frank Lisker and an anonymous reviewer are thanked for detailed comments on this work.