Audrey Gailler

Philippe Charvis

Ernst R Flueh

Segmentation of the Nazca and South American plates along the Ecuador subduction zone from wide angle seismic profiles

We describe the deep structure of the south Colombian -northern Ecuador convergent margin using travel time inversion of wide-angle seismic data recently collected offshore. The margin appears segmented into three contrasting zones. In the North Zone, affected by four great subduction earthquakes during the 20 th century, normal oceanic crust subducts beneath the oceanic Cretaceous substratum of the margin underlined by seismic velocities as high as 6.0-6.5-km/s. In the Central Zone the subducting oceanic crust is over-thickened beneath the Carnegie Ridge. A steeper slope and a welldeveloped, high velocity, Cretaceous oceanic basement characterizes the margin wedge. This area coincides with a gap in significant subduction earthquake activity. In the South Zone, the subducting oceanic crust is normal. The fore-arc is characterized by large sedimentary basins suggesting significant subsidence. Velocities in the margin wedge are significantly lower and denote a different nature or a higher degree of fracturing. Even if the distance between the three profiles exceeds 150 km, the structural segmentation obtained along the Ecuadorian margin correlates well with the distribution of seismic activity and the neotectonic zonation.

Introduction and tectonic setting

The Ecuador margin (3.0°N-3.5°S) represents a region of intense crustal deformation, related to the subduction of the oceanic Nazca plate beneath the South-American continental plate, along a N90°E direction and at a rate of ~58-mm/yr [START_REF] Trenkamp | Wide plate margin deformation, southern Central America and northwestern South America, CASA GPS observations[END_REF]. The margin undergoes shortening perpendicular to the trench axis, involving the north Andean block expulsion at ~6 mm/y [START_REF] Trenkamp | Wide plate margin deformation, southern Central America and northwestern South America, CASA GPS observations[END_REF] to the Northeast, parallel to the margin (Fig. 1). The oceanic Nazca plate in this area has an age to less than 26 My [START_REF] Lonsdale | Creation of the Cocos and Nazca plates by fission of the Farallon plate[END_REF], and carries the Carnegie Ridge (CR), which testifies to the interaction between the Galapagos hotspot (GHS) and the Cocos-Nazca spreading centre since ~20 My (e.g., [START_REF] Sallares | Crustal thickness constraints on the geodynamic evolution of the Galapagos Volcanic Province[END_REF]). Its subduction beneath the Ecuadorian coast induces lateral variations of seismicity, arc magmatism, deformations, vertical motion and sediment distribution along the margin [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF][START_REF] Gutsher | Tectonic segmentation of the North Andean margin: impact of the Carnegie Ridge collision[END_REF]. The Nazca plate is also marked by huge scarps such as the Yaquina Graben, interpreted as a transform fault associated with the remnant East Malpelo spreading axis [START_REF] Lonsdale | Creation of the Cocos and Nazca plates by fission of the Farallon plate[END_REF][START_REF] Hey | Tectonic evolution of the Cocos-Nazca spreading center[END_REF][START_REF] Lonsdale | Structure and tectonic history of the eastern Panama Basin[END_REF] (Fig. 1). The Ecuador margin consists of the oceanic North-Andean block which was accreted to the Brazilian continental craton during the Cretaceous [START_REF] Jaillard | Palaeogene deformations of the forearc zone of south Ecuador in relation to the geodynamic evolution[END_REF]. These oceanic terranes forming the Ecuadorian margin substratum seem to extend southward into northern Peru (i.e., Banco Peru; Fig. 1) [START_REF] Shepherd | Coastal Structure of the Continental-Margin, Northwest Peru and Southwest Ecuador[END_REF]. They consist of tholeitic/andesitic basalts and pillow lavas associated with siliceous sediments, defined as the Piñon Formation. In Ecuador, the Piñon Formation is overlain by an intra-oceanic volcanic arcs series, such as the Cayo formation (e.g., [START_REF] Jaillard | Basin development in an accretionnary, oceanic-floored fore-arc setting : southern coastal Ecuador during late Cretaceaous -late Eocene time[END_REF][START_REF] Reynaud | Oceanic plateau and island arcs of southwestern Ecuador: their place in the geodynamic evolution of northwestern South America[END_REF]). To the north (south Colombia, Gorgona Island), outcrops of Cretaceous oceanic terranes (mafic to ultramafic rocks) are associated with the Diabasico Formation [START_REF] Goosens | Geochemistry of tholeiites of the basic igneous complex of Northwestern South America[END_REF][START_REF] Mamberti | The Early Cretaceous San Juan Plutonic Suite, Ecuador: A magma chamber in an oceanic plateau?[END_REF][START_REF] Marriner | The Petrology and Geochemistry of Cretaceous to Recent Volcanism in Colombia -the Magmatic History of an Accretionary Plate Margin[END_REF].

Several studies on deformation of Andean coastal basins (e.g. [START_REF] Deniaud | Enregistrement sédimentaire et structural de l'évolution géodynamique des Andes équatoriennes au cours du Néogène : étude des bassins d'avant-arc et bilans de masse[END_REF]) have shown that the margin has undergone a generally compressive regime since accretion of the oceanic terranes to the Ecuadorian paleo-margin (~90 My ago). This compressive regime appears largely disturbed in the Gulf of Guayaquil (GG) area, where the tectonic structures testify to a N-S trending opening of the GG basin. The development of this deepest sedimentary basin of the Ecuadorian fore-arc during the last 5 Ma is believed to have resulted from the escape of the north Andean block to the N-E (Fig. 1) [START_REF] Witt | Development of the Gulf of Guayaquil (Ecuador) during the Quaternary as an effect of the North Andean block tectonic escape[END_REF]. Recent coastal uplift facing CR is indicated by Pliocene marine terraces (Tablazos) exposed at 200-300-m elevations [START_REF] Vries | The geology of late Cenozoic marine terraces (tablazos) in northwestern Peru[END_REF][START_REF] Pedoja | Plio-Quaternary uplift of the Manta Peninsula and La Plata Island and the subduction of the Carnegie Ridge, central coast of Ecuador[END_REF]. To the north, the morphology of the margin, the diversion of drainage networks and increases in coastal accretion support an active subsidence in relation to the seismic activity of the margin [START_REF] Dumont | Fan beheading and drainage diversion as evidence of a 3200-2800 BP earthquake event in the Esmeraldas-Tumaco seismic zone: A case study for the effects of great subduction earthquakes[END_REF]. To the south, the GG area appears also as a globally subsiding zone (e.g. [START_REF] Witt | Development of the Gulf of Guayaquil (Ecuador) during the Quaternary as an effect of the North Andean block tectonic escape[END_REF]). Our group conducted two major marine seismic experiments along the south Colombia-Ecuador convergent margin: the SISTEUR cruise in 2000 [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF], and the SALIERI cruise in 2001 [START_REF] Flueh | Cruise Report SO159 -Salieri[END_REF]. These surveys provided a unique set of high-quality wide-angle reflection/refraction and multichannel seismic data (MCS) that enabled us to image the structure of both plates as well as the geometry and properties of the inter-plate zone. This paper, based on the tomographic inversion of 3 wide-angle seismic profiles shot across the margin (Fig. 1), provides the first overview of the structure of the crust and upper mantle of the Nazca and South American plates at the South Colombia -Ecuador subduction zone. We will then tentatively correlate long wavelength structural variations observed along the margin with subduction processes and large subduction earthquakes.

Data acquisition

During the SISTEUR cruise, wide-angle seismic data were acquired using a seismic source composed of an array of eight 16-l airguns, shot at 125 m intervals. During the SALIERI cruise, the seismic source was composed of up to three 32-l airguns; in this case, the shot interval was set at a time of 60 s leading to a ∼140 m shot spacing. 24 UTIG-IRD OBS [START_REF] Nakamura | Development of upgraded ocean bottom seismograph[END_REF] and up to 10 land-seismometers were used during the SISTEUR cruise, whereas 26 UTIG-IRD OBS and 24 GEOMAR OBH and OBS [START_REF] Flueh | Cruise Report SO159 -Salieri[END_REF] were used during the SALIERI cruise. Most of the wide-angle data were of good quality, showing clear intracrustal refracted phases (Pg), deeper arrivals refracted in the upper mantle (Pn), and reflections at the boundary between the subducting oceanic crust and the mantle (PmP) (Fig. 234, table [START_REF] Trenkamp | Wide plate margin deformation, southern Central America and northwestern South America, CASA GPS observations[END_REF][START_REF] Lonsdale | Creation of the Cocos and Nazca plates by fission of the Farallon plate[END_REF][START_REF] Sallares | Crustal thickness constraints on the geodynamic evolution of the Galapagos Volcanic Province[END_REF]. Pg is here used as a generic term to indicate all the first arrivals that are not refracted in the mantle (sediments and basement). Significant variations were observed between record sections for instruments located in deep water, shallow-water and on-land instruments. This is probably related to the strong change in the bathymetry (more than 3000 m) across the margin but also to the strong lateral variation in the crustal structure along each profile.

Record sections for OBS located seaward from the trench

OBS located on top of the Nazca plate, seaward from the trench, show Pg phases with apparent velocities ranging from 2.5 to ~6.0 km/s. South of the CR, Pg phases rapidly reach apparent velocities higher than 5 km/s at near offset (Fig. 4, OBH 45 and Fig. 3, OBS 04). In contrast, the apparent velocity of Pg phases is less than 4 km/s north of the CR (OBH 116, Fig. 2). At larger offsets, apparent velocities of first arrivals generally exceed 7.0-7.5 km/s on most record sections. This has been interpreted to correspond to refraction in the upper mantle (Pn), except westward of OBS 04 (Fig. 3) located on top of the CR along profile SIS-4 where the Pn phase is lacking. In this section (Fig. 3) the apparent velocity, up to 40 km offset, is 7 km/s and probably corresponds to refraction within the subducting oceanic crust. Good quality PmP arrivals are also observed.

Record sections for OBS located landward from the trench

For OBS located on the continental slope and shelf Pg arrivals are characterized by low apparent velocities (from 2 to 4 km/s) up to 50 km along profile SAL-2 (OBS 57, Fig. 4). A shadow zone is observable on many instruments located along profiles SAL-6 and SIS-4 (Fig. 2 and3; OBS 121 and 10 respectively). In both cases, this shadow zone occurs between Pg phases with apparent velocities of ~4 km/s and ~6 km/s. Along profile SAL-6, the shadow zone is ~0.5 s thick. At larger offsets (>50-60 km), Pn arrivals are identified as first arrivals on most record sections along the three profiles. PmP arrivals are observed to be variable in quality (e.g., SAL-6; Fig. 2, OBS 121).

Travel time picking

Travel times of refracted (Pg, Pn) and reflected (PmP) phases identified in the record sections were picked manually (table 1-3). Most of the time, the PmP cannot be identified on the shallow-water instruments. The estimated picking error takes into account the quality of the phase and a possible systematic shift in the arrival identification in order to weight the different picks during the inversion. A 40-50 ms and 60-70 ms picking error was assigned respectively to near offsets and far offset Pg arrivals. A 80-90 ms error was assumed for Pn phases, and 100 ms for PmP phases.

Travel Time Tomography

Method

A joint refraction and reflection travel time tomography method was used to build a two-dimensional velocity model along the three seismic lines [START_REF] Korenaga | Crustal structure of the southeast Greenland margin from joint refraction and reflection seismic tomography[END_REF]. The method allows the simultaneous inversion of first arrival travel times (refracted phases) and a reflection phase, in order to simultaneously determine a 2-D velocity field and the geometry of a floating reflector (not associated with a change in velocity in the model). In this study, we choose to invert the reflector corresponding to the crust-mantle boundary (PmP phase). Detailed descriptions of this method as well as parameters used for the regularization of the inversion can be found in [START_REF] Korenaga | Crustal structure of the southeast Greenland margin from joint refraction and reflection seismic tomography[END_REF]. Using this inversion procedure means that we have to take into account three major issues: The first issue is due to the inherent non-uniqueness of the inversion solution. The final tomographic model depends on the starting velocity model chosen for the inversion, on the method of linearization and on the geometry of the experiment. The second issue concerns the estimate of uncertainties in the final velocity model, i.e. how well-resolved are the different parts of the model. This is done by estimating the degree of dependence of the obtained solution on the initial starting model, as well as to estimate uncertainties, via Monte Carlo analysis. The method involves a number of inversions with a variety of initial models using travel time picks with random errors applied. If all the Monte Carlo realizations have the same probability and if the initial models cover the full region of non-null probability within the space of parameters, the standard deviation of the obtained solutions can be interpreted as a measure of the final model parameters uncertainty. The third issue concerns evaluating the trade-off between depth and velocity parameters linked to the reflected phase travel time inversion and depending on the source-receiver geometry. A conventional method to estimate the degree of velocity-depth ambiguity is to perform sensitivity tests, such as checkerboard tests analysis. This procedure consists of using synthetic anomalies of various sizes at specific places in the models to evaluate at what point these anomalies are well-resolved by the travel time inversion.

Model construction and parameterization

In order to take into account the issues detailed above, the inversion is conducted using a three step approach on each of the three profiles: (1) Estimation of the 1D reference velocity model; (2) generation of a set of 100 1D Monte Carlo starting models; and (3) iterative inversion of all the starting models. This approach is similar to that in [START_REF] Korenaga | Crustal structure of the southeast Greenland margin from joint refraction and reflection seismic tomography[END_REF] (See Fig. 5 in [START_REF] Sallares | Seismic structure of Cocos and Malpelo Volcanic Ridges and implications for hot spot-ridge interaction[END_REF] for more details). The main difference between our and Sallarès et al. [START_REF] Sallares | Seismic structure of Cocos and Malpelo Volcanic Ridges and implications for hot spot-ridge interaction[END_REF] approach is that we adopt as final velocity model the best result of all the Monte Carlo realizations, characterized by the best RMS-χ 2 couple (Fig. 5), rather than the average of all the inversion models. The velocity features modelled along each profile are similar in both the final best-fitting models (Fig. 5) and the mean final models derived from the Monte Carlo approach (Fig. 6-a1, 6-b1, 6-c1) but they are smoothed in the latter.

Resolution and accuracy

Model parameter uncertainties

The derivative weight sum (DWS, Fig. 6) [START_REF] Toomey | Tomographic Inversion of Local Earthquake Data from the Hengill-Grensdalur Central Volcano Complex, Iceland[END_REF], allows a concise representation of the ray coverage in the best final model (Fig. 5). However, no quantitative measurement concerning the resolution and accuracy of the model can be associated with this type of representation. The mean deviations of velocity and depth parameters from each Monte Carlo ensemble averaged can be interpreted as a measure of the models uncertainties and therefore used to estimate the robustness of our iterative scheme (Fig. 6), which correlates well with the ray coverage imprint (DWS). Along profile SIS-4 (Fig. 6-b2 and 6-b3), the > 10 km OBS spacing seaward from the trench can be associated with a lack of short offset refractions and consequently to a larger uncertainty in the upper part of the oceanic crust (down to 7 km depth). Beyond 100-km distance, only refractions sample the lower part of the model so it cannot be constrained properly between 15 and 22 km depth. Along profile SAL-2 (Fig. 6-c2 and 6-c3), the amplitude of arrivals recorded by OBS deployed on the continental margin becomes rapidly weak due to a thick sedimentary cover, so no clear PmP or deep refracted phase can be identified (Fig. 4, OBS 57). Rays poorly sample the model, from 90 km to the landward end of the line, and consequently the Moho depth and seismic velocities are poorly constrained in this area beneath 14 km depth. Along profile SAL-6 (Fig. 6-a2 and 6-a3), the part of the model located between 50 and 70 km distance and including the lower part of the continental margin (below 10 km depth) and the subducting oceanic plate (up to 15 km depth) shows the largest velocity uncertainties (0.2-0.25 km/s). Absolute velocities are not well constrained there. In this zone the DWS shows a lack in ray coverage associated with the shadow zone observed on several record sections (Fig. 2, OBH 116 and OBS 121). This shadow zone appears as a thick low velocity zone that encloses the lower part of the margin and the inter-plate contact zone in both the best fitting model (Fig. 5a) and the final averaged Monte Carlo model (Fig. 6-a1). The 100 2D final models set linked to the Monte Carlo analysis shows that the low velocity zone occurs systematically and independently of the initial velocity model considered. Checkerboard tests will help to estimate the resolvability of this low velocity zone.

Checkerboard Tests

A number of checkerboard tests (Fig. 7) has been performed in order to determine the size and amplitude of synthetic velocity anomalies resolved in different parts of the velocity models. Along all profiles, the shape of the velocity anomalies is well recovered beneath the OBS network (≤10 km spacing) and only a few small artefacts (< 0.2 km/s) near the crust-mantle boundary interface are observed. The differences are especially significant in profile SIS-4, in agreement with the sparser acquisition geometry and the less clear PmP and Pn phases recorded along this line, which affect the resolvability of the deep parts of the model (Fig. 7-b). Beneath zones where OBS are more than 15 km spaced, the pattern of velocity anomalies spreads in the direction of the mean ray-paths (Fig. 7-c2), and velocity artefacts are larger and higher (0.2-0.4 km/s; Fig. 7-c2). Along profile SAL-6, synthetic tests enable checking of the resolvability of the low velocity zone observed in both the best fitting model and the mean final model. The location and amplitude of velocity anomalies are reasonably well recovered within model uncertainty limits (Fig. 7-a2), indicating that the experiment geometry and data set allow a distinction between positive and negative variations at 6 to 17 km depth (Fig. 7-a2). The sensitivity of the data is therefore considered sufficient to resolve structures with size and amplitudes similar to the low velocity zone at this depth range implying that the low-velocity zone observed in the final models is not an artefact linked to the inversion procedure; even if the absolute value of velocity in this area is poorly constrained by the data set.

Results of tomographic inversion

The tomographic inversion of travel times provides three 2D velocity models across the Ecuador convergent margin and the geometry of the Moho of the Nazca plate (Fig. 5). Since there is no clear reflector separating the overriding margin from the oceanic subducting plate in the wide-angle seismic data, except along profile SIS-4 [START_REF] Graindorge | Seismic evidences of Carnegie Ridge subduction and the Ecuadorian margin structure at 1.4°S[END_REF], we have digitized and converted to depth the inter-plate contact reflectors observed along coincident available MCS sections [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF][START_REF] Calahorrano | Structure de la marge du Golfe de Guayaquil (Equateur) et propriétés physiques du chenal de subduction à partir de données de sismique marine réflexion et réfraction[END_REF][START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF][START_REF] Marcaillou | Régimes tectoniques et thermiques de la marge Nord Equateur -Sud Colombie (0°-3.5°N) -Implications sur la sismogenèse[END_REF]. The top of the subducting plate and the decollement, which mark respectively the base and the top of the subduction channel, are displayed as well as several major faults affecting the margin (Fig. 5). Only long-wavelength velocity anomalies are determined from travel time inversion and so it is likely that velocity in the subduction channel is not properly calculated. Consequently the thickness of the subduction channel observed in the tomographic models must be considered as an upperbound. Nevertheless, the comparison of wide-angle seismic models and MCS data enables the inter-plate contact in the velocity model to be located, helping to correlate and interpret velocity structures. The uncertainty on the Moho depth is ±1 km in the better-constrained central part of the model (Fig. 6). The crust-mantle transition is marked by a transition from crustal velocities (< 7 km/s) to mantle velocities (7.5 to 8.0 km/s). The agreement between the velocity transition and the Moho location is quite good in the central part of the model, where both parameters are well constrained according to the outcome of the uncertainty analysis. The dip of the subducting plate is estimated from the dip of the Moho along the different profiles, which is consistent with the dip of the top of the oceanic crust derived from MCS data (Fig. 5,table 4).

Interpretation

The tomographic models are not split into discrete layers. Nevertheless, to support the geological interpretation of the models, several model bodies or layers are differentiated based mainly on the velocities. Along profile SAL-6 the base of the sedimentary cover inferred from the coincident MCS profile SIS-44 (Fig. 8) correspond to velocities ranging from 3.5 to 4.0 km/s. At 4 km depth seismic velocities reach 3.6 km/s for sands (30% porosity) or marine shales [START_REF] Japsen | Constraints on velocity-depth trends from rock physics models[END_REF]. In this paper, velocities lower than 3.5 km/s are assumed to represent sedimentary rocks. Based on this assumption it is likely that the thickness of the sedimentary cover is underestimated because at depth velocities in highly compacted sediments are likely to exceed 3.5 km/s. The margin basement displays velocities ranging from 3.5 to 5.5 km/s along profile SAL-2 but velocities reach 6.2 km/s at depth along profiles SAL-6 and SIS-4. These differences will be discussed in the next sections. The structure of the oceanic crust is discussed using a classical 3 layer scheme: oceanic layer 1 consisting of sediments with velocities lower than 3.5 km/s; oceanic layer 2 (ol2) characterized by velocities ranging from 3.5 to 6.5 km/s and a high vertical velocity gradient and oceanic layer 3 (ol3) characterized by velocities higher than 6.5 km/s, its base defined by the Moho discontinuity. According to the velocity models, we find considerable contrasts in structure, morphology, and sedimentary cover for both plates. We propose a zonation of the South Colombia-Ecuador convergent margin into three contrasted areas (table 4, Fig. 1): the north zone (2°N-1°N), the central zone (1°N-2.5°S) and the south zone (2.5°-3.5°S).

The North Zone : Profile SAL-6 (North Ecuador, 2°N-1°N)

 Along profile SAL-6, the velocity structure of the oceanic crust is well constrained from 10 to 40 km distance because of the 7 OBS located west of the trench (Fig. 5a). A ~1.5 to 2.0 km thick sedimentary cover (V p < 3.5 km/s) overlays oceanic layer 2 characterized by a large vertical velocity gradient, ranging from 3.5 km/s at ~5 km depth to 6.5 km/s at ~7 km depth. The small velocity gradient in the lower part of the oceanic crust (V p from 6.5 to 7.2 km/s) is in good agreement with average oceanic layer 3 velocities [START_REF] Mutter | Variations in thickness of layer 3 dominate oceanic crustal structure[END_REF]. The Moho is 10 to 11 km deep, leading to a 6.0 to 7.0 km thick igneous oceanic crust with an ol2/ol3 thickness ratio of ~0.59.

The velocity structure is consistent with "normal" oceanic crust showing no evidence of hotspot influence (e.g. [START_REF] White | Oceanic crustal thickness from seismic measurements and rare earth element inversions[END_REF]). The hypothesis suggesting a continuation of the Cocos-Malpelo Ridge east of the Yaquina Graben, between 2°N and 4°N [START_REF] Sallares | Crustal thickness constraints on the geodynamic evolution of the Galapagos Volcanic Province[END_REF], is not supported by our dataset. East of the trench (distance 50-90 km, Fig. 5a), the Moho dips at 8.0° ±0.5°. The base of the crust can be followed landward for a distance of 80-90 km with a good accuracy (Fig. 5a). Seismic velocities immediately above the Moho remain constant at around 7.0-7.5 km/s even if they are likely overestimated due the tomographic approach used. On the contrary, the velocity increases rapidly beneath the MCS reflector marking the top of the subducting oceanic crust (black dots, Fig. 5a). Its velocity increases from 5.0 km/s at 55-km to 6.0 km/s at 95-km distance in the down-dip direction, which possibly indicates modifications of the upper part of the oceanic plate during its burial.  The sedimentary layers observed along MCS line SIS-44 (Fig. 8) are consistent with velocities lower than 3.5 km/s modelled along the coincident wide-angle seismic line SAL-6 (Fig. 5a). The trench shows a 2 to 3 km thick sedimentary infill along the north Ecuador margin. This estimate is consistent with previous studies [START_REF] Calahorrano | Structure de la marge du Golfe de Guayaquil (Equateur) et propriétés physiques du chenal de subduction à partir de données de sismique marine réflexion et réfraction[END_REF][START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF][START_REF] Marcaillou | Régimes tectoniques et thermiques de la marge Nord Equateur -Sud Colombie (0°-3.5°N) -Implications sur la sismogenèse[END_REF] along the south Colombia margin, indicating that the trench infill is probably due to recent and massive turbiditic intakes, via canyons, in association with hemipelagic sedimentation. The continental slope immediately east of the trench is covered by 2.0-3.5 km thick sedimentary deposit. This correlates with a small, probably ancient, accretionary prism which has been described in this area using MCS data (Fig. 8) [START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF][START_REF] Marcaillou | Régimes tectoniques et thermiques de la marge Nord Equateur -Sud Colombie (0°-3.5°N) -Implications sur la sismogenèse[END_REF].

On the continental margin, the thickness of the sedimentary cover is quite uniform around ~2.5 km on average. Two sedimentary basins lie on the fore-arc (55-90-km and 90-125-km distance, Fig. 5a), and the sedimentary basin becomes thinner at the front of the margin. The continental platform is characterized by a ~600 m deep and ~90 km wide morphological reentrant limited landward by a subsiding coastal area [START_REF] Dumont | Fan beheading and drainage diversion as evidence of a 3200-2800 BP earthquake event in the Esmeraldas-Tumaco seismic zone: A case study for the effects of great subduction earthquakes[END_REF]. The subsidence of the continental platform is supported by the presence of the 2-3 km thick sedimentary basins covering the margin basement. Correlations with on-land drills made in the fore-arc basins on top of the margin [START_REF] Mountney | Quantitative analysis of Miocene to Recent forearc basin evolution along the Colombian convergent margin[END_REF] suggest a Palaeogene to Neogene calcareous, clay and sand origin for this sedimentary infill (e.g. [START_REF] Deniaud | Enregistrement sédimentaire et structural de l'évolution géodynamique des Andes équatoriennes au cours du Néogène : étude des bassins d'avant-arc et bilans de masse[END_REF]).  The margin basement is characterized by velocities ranging from 5.4 to 6.2 km/s (Fig. 9) that occupies a large part of the deep margin, and can be followed from ~ 35 km east of the trench to the coast at ~ 7 km depth (see table 5). The highest basement velocities (~6.2 km/s) are reached at 50-km distance from the trench, at ~10 km depth.  The reflector interpreted from the MCS data (Fig. 8) as a major splay fault dipping to the east [START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF] marks the top of a large low velocity zone inferred from our dataset (Fig. 5a). It clearly soles out at the inter-plate boundary at 13-14 km depth (45 km from the trench), and bounds the westward limit of the "high velocity" (V p ≥ 6.0 km/s) upper plate basement. The seaward section of this system forms in fact a large low velocity zone and appears significantly disturbed at the deformation front (Fig. 8).

The Central Zone: profile SIS-4 (Carnegie Ridge, 1°N-2.5°S)

The final velocity model we obtain along profile SIS-4, located on the southern flank of the CR, is in agreement with the layered model based on the forward modelling of the same dataset [START_REF] Graindorge | Seismic evidences of Carnegie Ridge subduction and the Ecuadorian margin structure at 1.4°S[END_REF] but exhibits some differences.  The oceanic crust velocity structure is followed up to 90-km westward from the trench (Fig. 5b). The total thickness of the crust is ~15 km here, based on the Moho location obtained from PmP travel times inversion (15 to 18 km depth, Fig. 6-b1). ol2 is ~3.5 km thick, similar to that observed along SAL-6 (~3 km) whereas ol3 is 10 km thick, two fold the thickness of ol3 beneath line SAL-6. A travel time tomography analysis along a wide-angle line extending from the top of CR to the oceanic basin across its northern flank [START_REF] Sallares | Seismic structure of the Carnegie ridge and the nature of the Galapagos hotspot[END_REF] (see location of refraction line in Fig. 1) reveals that the maximum crustal thickness beneath the crest of the ridge is 19 km (~100 km north of profile SIS-4) with a similar velocity structure. If we compare this profile with profile SIS-4, it appears that the thickness of ol2 remains uniform through CR whereas the thickness of ol3 increases with increasing total thickness. The crustal thickening is dominated by ol3 thickening (ol2/ol3 thickness ratio <0.3), which confirms that the CR has originated from the interaction between the GHS and the Galapagos spreading centre [START_REF] Sallares | Crustal thickness constraints on the geodynamic evolution of the Galapagos Volcanic Province[END_REF][START_REF] Barckhausen | Revised tectonic boundaries in the Cocos Plate off Costa Rica: Implications for the segmentation of the convergent margin and for plate tectonic models[END_REF][START_REF] Hey | Recent Plate Motions in Galapagos Area[END_REF][START_REF] Wilson | History of Rift Propagation and Magnetization Intensity for the Cocos-Nazca Spreading Center[END_REF], generating increased melt production within the mantle source. Due to the scarcity of data, the model is not well constrained east of 130-km distance. Reflection from the top of the oceanic crust along coincident MCS lines and modelling of intra-crustal wide-angle reflections provide a dip ranging from 4° beneath the slope to ~ 10° dip beneath the Ecuadorian coast at ~15 km depth [START_REF] Graindorge | Seismic evidences of Carnegie Ridge subduction and the Ecuadorian margin structure at 1.4°S[END_REF]. Deeper and farther inland the Wadati-Benioff zone dip varies between 35° and 25° between 1°N and 1.5°S, reaching at least a ~150 km depth in the north and ~ 200 km depth in the south [START_REF] Guillier | Seismological evidence on the geometry of the orogenic system in central-northern Ecuador (South America)[END_REF].  On the CR and in the oceanic basin south of it, the sedimentary cover and trench fill are as thin as 0.5 to 1.0 km assuming a velocity lower than 3.5 km/s. The lack of sediments on top of the CR and on its southern flank could be explained by the prominent topographic feature of the CR, which acts as a barrier and influences marine current trajectories, sedimentary flux, deposition [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF][START_REF] Gutsher | Tectonic segmentation of the North Andean margin: impact of the Carnegie Ridge collision[END_REF][START_REF] Michaud | Fields of multi-kilometer scale sub-circular depressions in the Carnegie ridge sedimentary blanket : effect of underwater carbonate dissolution ?[END_REF][START_REF]Southeast Pacific Paleoceanographic Transects Sites 1232 -1242[END_REF], and the erosive power of strong marine currents. Similarly, the continental shelf shows a very thin sedimentary cover (less than ~1 km thick), which confirms that no subduction-accretion occurs along this segment [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF][START_REF] Gutsher | Tectonic segmentation of the North Andean margin: impact of the Carnegie Ridge collision[END_REF][START_REF] Michaud | Fields of multi-kilometer scale sub-circular depressions in the Carnegie ridge sedimentary blanket : effect of underwater carbonate dissolution ?[END_REF][START_REF]Southeast Pacific Paleoceanographic Transects Sites 1232 -1242[END_REF]. The coastal area facing the CR has been uplifted as indicated by Quaternary marine terraces exposed at 200-300 m elevation [START_REF] Vries | The geology of late Cenozoic marine terraces (tablazos) in northwestern Peru[END_REF]. As a consequence, no rivers coming from the Andes cross the Coastal Cordillera and precludes sediment influx towards the trench [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF][START_REF] Gutsher | Tectonic segmentation of the North Andean margin: impact of the Carnegie Ridge collision[END_REF][START_REF] Michaud | Fields of multi-kilometer scale sub-circular depressions in the Carnegie ridge sedimentary blanket : effect of underwater carbonate dissolution ?[END_REF][START_REF]Southeast Pacific Paleoceanographic Transects Sites 1232 -1242[END_REF].  At the southern flank of the CR, the margin wedge consist of a body with velocities that reach 5.7 km/s at 30 km east of the trench and 6.2 km/s at 40 km east of the trench, at ~5 km depth (table 5). These velocities are higher than velocities described in similar context and will be discussed in the next section.

The South Zone: Profile SAL-2 (Gulf of Guayaquil, 2.5°S-3.5°S)

 Along profile SAL-2 the structure of the oceanic plate is followed up to ~ 47 km westward from the trench. Seaward from the trench, the Moho is located at ~11.5 km depth (Fig. 5c). The crust of the oceanic plate is ~ 8 km thick, only slightly thicker than along the northern line (SAL-6). Major differences are observable in the very thin sedimentary cover which has a maximum thickness of ~ 0.5 km (for velocities < 3.5 km/s), a thinner ol2 (~ 2.5 km thick) with velocities from ~ 3.5 to 6.0 km/s, and a 5 km thick ol3 (V p rising from 6.2 to 7.0 km/s). Despite a velocity structure close to normal oceanic crust, the low ol2/ol3 thickness ratio (~0.40) and the slight thickening of ol3 (~1 km thicker) suggest the crust underwent some influence of the GHS activity during its emplacement [START_REF] Mutter | Variations in thickness of layer 3 dominate oceanic crustal structure[END_REF]. East of the trench, the oceanic Moho dips at 5.0° ± 0.5°. The thickness of the subducting oceanic crust deduced from depth-migrated MCS reflectors is 6.5-7.0 km. The velocity structure of ol2 and ol3 can be followed from the trench to ~95-km distance and is characterized by isovelocity curves parallel to the Moho at the base of the crust. The velocity at the top of the oceanic crust increases slightly during the down-warping process (from 4 km/s at 4 km depth to 5.5 km/s at 13 km depth), similar to our observations along SAL-6. The tomographic model shows the margin sedimentary cover (V p ≤ 3.5 km/s) progressively thickening from ~ 1 km at margin front to up to 7 km in the inner part of the GG (between 50-km to 135-km distance). Because of the compaction processes that likely occur at depth, it is possible that velocities higher than 3.5 km/s characterize the lower part of the sedimentary pile. Consequently, the sedimentary thickness in the basins might be underestimated but this will not affect our interpretation. Interpretation of industrial multichannel seismic profiles show that the continental slope (distance 50-115 km, Fig. 5) is affected by a series of N-S trending listric normal faults controlled by the subduction processes [START_REF] Witt | Development of the Gulf of Guayaquil (Ecuador) during the Quaternary as an effect of the North Andean block tectonic escape[END_REF][START_REF] Witt | Constraints on the tectonic evolution of the North Andean Block trailing edge: evolution of the Gulf of Guayaquil-Tumbes Basin and the intermountain basins of the Central Ecuadorian Andes[END_REF]. On the contrary, the continental platform (east of 115-km distance along profile SAL-2) is mainly controlled by the northward tectonic escape of the North Andean Block since the Pleistocene. But both segments seem to have been subsiding zones since the Miocene [START_REF] Witt | Development of the Gulf of Guayaquil (Ecuador) during the Quaternary as an effect of the North Andean block tectonic escape[END_REF][START_REF] Witt | Constraints on the tectonic evolution of the North Andean Block trailing edge: evolution of the Gulf of Guayaquil-Tumbes Basin and the intermountain basins of the Central Ecuadorian Andes[END_REF]. According to this interpretation, the up to 7-km thick eastern basin located on the platform (distance 105-145-km) may be related both to subduction-erosion processes (i.e. tectonic erosion of the upper plate's underside) since the Miocene, and to the N-S extensive regime affecting the GG since the Pleistocene. Whereas the sedimentary basin lying on the upper part of the slope (distance 80-110-km, Fig. 5c) is more likely to be related to tectonic erosion at the base of the margin [START_REF] Witt | Development of the Gulf of Guayaquil (Ecuador) during the Quaternary as an effect of the North Andean block tectonic escape[END_REF][START_REF] Witt | Constraints on the tectonic evolution of the North Andean Block trailing edge: evolution of the Gulf of Guayaquil-Tumbes Basin and the intermountain basins of the Central Ecuadorian Andes[END_REF]. Westward, a thin sedimentary cover and infill characterize the deformation front and the trench area. Thus, the important sediment influx coming from the Andes does not seem to reach the trench, but is predominantly trapped in the large sedimentary basins lying on the continental platform and on the slope, in the region of greater subsidence. At the deformation front, a thinner sedimentary cover is associated with a progressive uplift of the margin basement possibly indicating weak compressive deformation expressed at the wedge [START_REF] Calahorrano | Structure de la marge du Golfe de Guayaquil (Equateur) et propriétés physiques du chenal de subduction à partir de données de sismique marine réflexion et réfraction[END_REF]. The subsidence of the margin fore-arc in front of the GG, made evident by the thick sedimentary basin, supports a concept of large scale tectonic erosion. In such a context, the strong dip of the margin front appears as a consequence of frontal erosion. Moreover, several MCS sections acquired in this area show basal erosion associated with the subduction of down going seamounts carried by the oceanic plate [START_REF] Sage | Interplate patchiness and subduction-erosion mechanisms: Evidence from depthmigrated seismic images at the central Ecuador convergent margin[END_REF].  In this zone, seismic velocities ranging from ~3.5 to 5.2 km/s characterize the basement of the margin. Here, velocities in the margin wedge are lower than the ones described along northern and central Ecuador profiles (Fig. 5, table 5). From the trench to 80-km distance the continental shelf shows a progressive uplift of this layer up to 1 km below the seafloor. At 100-km distance another high point of this layer reaches up to 2.5 km beneath the seafloor. It separates the upper slope from the eastern basin. Velocities obtained for the margin substratum, 60 km east of the trench, are as low as 5.2 km/s at 10 km depth, with a constant vertical gradient of ~0.2 km/s/km (Fig. 9). The GG roughly separates the northern fore-arc units (North Andean Block) made of oceanic terranes from the southern fore-arc units of continental composition [START_REF] Jaillard | Preface -Andean geodynamics: main issues and contributions from the 4th ISAG[END_REF][START_REF] Mourier | An Accreted Continental Terrane in Northwestern Peru[END_REF] but the limit between oceanic and continental basement is not precisely known. The Piñon formation outcrops on-land in northern Peru and possibly extends off-shore beneath Banco Peru (Fig. 1) [START_REF] Shepherd | Coastal Structure of the Continental-Margin, Northwest Peru and Southwest Ecuador[END_REF]. Shadow zones appear along several wide-angle seismic sections (e.g., Fig. 4). They cannot be modelled from inversion of first arrival travel times alone but possibly denote a low velocity subduction channel that supports the basal erosion of the margin's underside. In fact, the resulting debris dragged along the inter-plate contact zone generally create a small velocity inversion which gradually disappears with the down-going process(e.g., [START_REF] Sage | Interplate patchiness and subduction-erosion mechanisms: Evidence from depthmigrated seismic images at the central Ecuador convergent margin[END_REF]).

Discussion

 In this study, we corroborate that the Moho beneath the CR presents a strong dissymmetry with a sharp transition from the thick crust associated with volcanic ridges to normal oceanic crust associated with the oceanic basin beneath the northern flank [START_REF] Sallares | Seismic structure of the Carnegie ridge and the nature of the Galapagos hotspot[END_REF] (whereas the crust of the ridge thins progressively southward [START_REF] Sallares | Crustal thickness constraints on the geodynamic evolution of the Galapagos Volcanic Province[END_REF][START_REF] Sallares | Seismic structure of the Carnegie ridge and the nature of the Galapagos hotspot[END_REF]). We show that the oceanic crust south of the CR southern flank (profile SAL-2) is slightly thicker than "normal" oceanic crust probably because of the influence of the GHS, that is outlined by numerous seamounts described from multibeam bathymetric data (Fig. 1) [START_REF] Sage | Interplate patchiness and subduction-erosion mechanisms: Evidence from depthmigrated seismic images at the central Ecuador convergent margin[END_REF]. This is consistent with the present-day thermal anomaly associated with the GHS that influences the generation of new oceanic crust within at least 350 km radius [START_REF] Sallares | Seismic structure of the Carnegie ridge and the nature of the Galapagos hotspot[END_REF][START_REF] Canales | Crustal thickness along the western Galapagos Spreading Center and the compensation of the Galapagos hotspot swell[END_REF]. To the north (profile SAL-6), we find "normal" oceanic crust thickness which precludes the continuation of the Cocos-Malpelo Ridge east of the Yaquina Graben [START_REF] Sallares | Seismic structure of the Carnegie ridge and the nature of the Galapagos hotspot[END_REF].  The basement wedge offshore central Ecuador near La Plata Island exhibits higher velocities at similar distance from the trench than any other well-constrained margin wedge (table 5). Margins with similar velocities are the Nicoya Peninsula [START_REF] Christeson | Structure of the Costa Rica convergent margin, offshore Nicoya Peninsula[END_REF] consisting of a well-known ophiolitic complex, the Chile margin offshore Antofagasta where the basement consists of basaltic to andesitic rocks [START_REF] Sallares | Structure and tectonics of the erosional convergent margin off Antofagasta, north Chile (23 degrees 30 ' S)[END_REF], the Shumagin Islands margin wedge which consists of interstratified volcanic and sedimentary rocks [START_REF] Lewis | Structural Development of an Accretionary Prism by Thrust and Strike-Slip Faulting -Shumagin Region, Aleutian Trench[END_REF] and the Cascadia margin offshore Oregon associated with the Siletz terrane of oceanic origin [START_REF] Gerdom | The continental margin off Oregon from seismic investigations[END_REF][START_REF] Graindorge | Crustal structure beneath the Strait of Juan de Fuca and Southern Vancouver Island from seismic gravity analyses[END_REF][START_REF] Trehu | Crustal Architecture of the Cascadia Fore-Arc[END_REF]. The margin of Ecuador is part of the North Andean Block and mainly consists of oceanic cretaceous rocks outcropping extensively on-land and on La Plata Island immediately north of profile SIS-4 (Fig. 1). Then, the margin wedge consists undoubtedly of the same type of oceanic rocks. Furthermore, the velocities encounted in the margin basement are similar to the one determined in the upper crust of numerous oceanic plateaus like the Kerguelen Plateau [START_REF] Charvis | Deep structure of the Kerguelen-Heard Plateau and hot spot-related activity[END_REF], the Ontong Java Plateau [START_REF] Gladczenko | Crustal structure of the Ontong Java Plateau: Modeling of new gravity and existing seismic data[END_REF], the Caribbean Plateau [START_REF] Sallares | Seismic velocity structure across the middle American landbridge in northern Costa Rica[END_REF] or the Cocos and Carnegie ridges [START_REF] Sallares | Seismic structure of Cocos and Malpelo Volcanic Ridges and implications for hot spot-ridge interaction[END_REF][START_REF] Sallares | Seismic structure of the Carnegie ridge and the nature of the Galapagos hotspot[END_REF] (refer to the compilation in [START_REF] Christeson | Structure of the Costa Rica convergent margin, offshore Nicoya Peninsula[END_REF] for more details).  Velocities in the margin wedge offshore North Ecuador -South Colombia are comparable to those determined in the Cascadia margin and Eastern Nankai Trough, which include accreted oceanic fragment [START_REF] Dessa | Deep seismic imaging of the eastern Nankai trough, Japan, from multifold ocean bottom seismometer data by combined travel time tomography and prestack depth migration[END_REF][START_REF] Nakanishi | Detailed subduction structure across the eastern Nankai Trough obtained from ocean bottom seismographic profiles[END_REF]. They are not as high as the velocities observed in front of the Carnegie ridge at the toe of the margin. Nevertheless, outcrops of oceanic rocks on-land at the border between Colombia and Ecuador, as well as on Gorgone Island less than 140 km north of profile SAL-6 [START_REF] Goosens | Geochemistry of tholeiites of the basic igneous complex of Northwestern South America[END_REF][START_REF] Mamberti | The Early Cretaceous San Juan Plutonic Suite, Ecuador: A magma chamber in an oceanic plateau?[END_REF][START_REF] Marriner | The Petrology and Geochemistry of Cretaceous to Recent Volcanism in Colombia -the Magmatic History of an Accretionary Plate Margin[END_REF][START_REF] Nivia | Geochemistry and origin of the amaime and volcanic sequences southwestern Colombia[END_REF], prove a similar oceanic nature for the basement wedge which is part of the North Andean Block. Along the Peru margin, despite rather high seismic velocities, the margin wedge consists of Palaeozoic or Precambrian foliated metamorphic rocks of continental origin [START_REF] Krabbenhoft | Crustal structure of the Peruvian continental margin from wide-angle seismic studies[END_REF]. Velocities in the margin wedge beneath the Gulf of Guayaquil are lower than those observed along the central and north Ecuador margin and could denote either a different nature or a higher degree of fracturing.  A progressive increase of seismic velocities in the backstop from 3.0-3.5 km/s near the trench to 4.5-6.2 km/s over a horizontal distance of 20 km is observed in all of the well-resolved model based on forward modelling or on tomographic inversion (see table 5). Velocities lower than 3.5 km/s are most likely to represent sedimentary rocks either accreted at the margin or deriving from the continent. Low velocities observed at the margin toe are usually interpreted as associated to fracturing, or alteration processes that affect the rocks of the margin basement and increase the porosity of these rocks [START_REF] Christeson | Structure of the Costa Rica convergent margin, offshore Nicoya Peninsula[END_REF][START_REF] Sallares | Structure and tectonics of the erosional convergent margin off Antofagasta, north Chile (23 degrees 30 ' S)[END_REF][START_REF] Krabbenhoft | Crustal structure of the Peruvian continental margin from wide-angle seismic studies[END_REF][START_REF] Ranero | Subduction erosion along the Middle America convergent margin[END_REF]. These processes decrease landward where the margin wedge is thicker and more mechanically resistant.  The seismic velocities obtained along profile SAL-2 for the margin basement of the GG cannot be immediately linked to accreted oceanic terranes as the GG probably represents the southern limit of the North Andean Block [START_REF] Witt | Constraints on the tectonic evolution of the North Andean Block trailing edge: evolution of the Gulf of Guayaquil-Tumbes Basin and the intermountain basins of the Central Ecuadorian Andes[END_REF]. Velocities in the margin basement range reach 4.5 to 5.0 km/s at 13 km depth. They are similar to velocities in the OSA Peninsula (Costa Rica) where the rocks outcropping are a middle mid-Eocene/Miocene melange dominated by basalt, chert, and limestone resulting from accretion of seamounts [START_REF] Vannucchi | From seamount accretion to tectonic erosion: Formation of Osa Melange and the effects of Cocos Ridge subduction in southern Costa Rica[END_REF]. It is also comparable to the velocity structure described along the Western Nankai Trough, where the margin basement is interpreted as an old metamorphosed accretionary complex, the Shimanto Belt [START_REF] Kodaira | Western Nankai Trough seismogenic zone: Results from a wide-angle ocean bottom seismic survey[END_REF][START_REF] Takahashi | Seismic structure of western end of the Nankai trough seismogenic zone[END_REF]. So, it is not possible from seismic velocity alone to determine the nature of the basement beneath the GG. The low velocities obtained along profile SAL-2 are lower than velocities described northward in front of the Carnegie Ridge (this study) interpreted to be related to the oceanic North Andean Block. But they are also lower than velocities described in the basement of the margin off northern Peru [START_REF] Krabbenhoft | Crustal structure of the Peruvian continental margin from wide-angle seismic studies[END_REF] which are related to the metamorphic basement. We assume that these low velocities in the GG basement should be mainly related to a high degree of fracturing associated with the NS extension in the Gulf itself and with the erosion of the outer part of the margin [START_REF] Witt | Constraints on the tectonic evolution of the North Andean Block trailing edge: evolution of the Gulf of Guayaquil-Tumbes Basin and the intermountain basins of the Central Ecuadorian Andes[END_REF].  Tectonic erosion is favoured in regions where convergence rates exceed ~60-mm/yr and where the sedimentary cover is <1 km [START_REF] Clift | Controls on tectonic accretion versus erosion in subduction zones: Implications for the origin and recycling of the continental crust[END_REF]. Rates of subduction-erosion can be computed efficiently from the offshore subsidence record and the geometry of the margin [START_REF] Kukowski | Subduction erosion -the "normal" mode of fore-arc material transfer along the Chilean margin?[END_REF]. Along the margin, from Colombia to the GG, the slope of the margin is 15° to 20° (Fig. 1 and5). Subsidence can be inferred by the thickness of low velocities and non-compacted sediments. The steepness of the slope, the lack of an accretionary prism except near the Colombian border, subsidence of the middle slope and of the forearc basin suggest significant erosion at the base of the overriding plate [START_REF] Kukowski | Subduction erosion -the "normal" mode of fore-arc material transfer along the Chilean margin?[END_REF]. Subduction-erosion can be envisaged as being linked to strong coupling between the overriding and subducting plate, although processes like hydro-fracturing at the base of the overriding fore-arc crust are also identified as a convincing explanation for basal tectonic erosion. Fluid overpressure, finally, is thought to generate fractures in the basal region of the upper plate [START_REF] Kukowski | Subduction erosion -the "normal" mode of fore-arc material transfer along the Chilean margin?[END_REF][START_REF] Lallemand | Coulomb Theory Applied to Accretionary and Nonaccretionary Wedges -Possible Causes for Tectonic Erosion and or Frontal Accretion[END_REF][START_REF] Huene | Generic model of subduction erosion[END_REF]. This loosens fragments that are then dragged into the subduction channel where they may undergo further downward transport [START_REF] Kukowski | Subduction erosion -the "normal" mode of fore-arc material transfer along the Chilean margin?[END_REF][START_REF] Lallemand | Coulomb Theory Applied to Accretionary and Nonaccretionary Wedges -Possible Causes for Tectonic Erosion and or Frontal Accretion[END_REF][START_REF] Huene | Generic model of subduction erosion[END_REF]. In north Ecuador, a small accretionary prism observed at the front of the margin suggests accretion whereas the strong subsidence of the margin supports basal erosion at the deep part of the basement. The margin appears to be at the limit between subduction-accretion which is more developed northward (towards Colombia) and subduction-erosion as it has been suggested from previous studies [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF][START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF][START_REF] Marcaillou | Régimes tectoniques et thermiques de la marge Nord Equateur -Sud Colombie (0°-3.5°N) -Implications sur la sismogenèse[END_REF][START_REF] Gutscher | Tectonic segmentation of the North Andean margin: impact of the Carnegie Ridge collision[END_REF].

At the CR, the increased buoyancy of the Nazca plate is likely to cause the observed general uplift of the coast in this area [START_REF] Pedoja | Plio-Quaternary uplift of the Manta Peninsula and La Plata Island and the subduction of the Carnegie Ridge, central coast of Ecuador[END_REF]. There is no evidence of subduction-accretion along this segment of the margin on the MCS data [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF]. The high velocity oceanic basement, which extends westward only 30 km from the trench (Fig. 5b) and seems to extend beneath the whole margin, as well as the lack of fore-arc basin, do not support a basal erosion scenario. Furthermore, the MCS data reveal basal erosion related to the presence of seamounts on the subducting plate, which affects the margin south of latitude 1°40' S [START_REF] Sage | Interplate patchiness and subduction-erosion mechanisms: Evidence from depthmigrated seismic images at the central Ecuador convergent margin[END_REF]. South of the CR, the strong subsidence of the margin, the low seismic velocity of the possibly highly fractured basement and evidence for low velocity material in the subduction channel despite the lack of sediments in the trench support an intense erosion of the margin associated with a north-south extension of the continental platform. In this area, we suggest basal erosion is related to high-pressure fluids rather than a strong coupling.  The segmentation of the margin is primarily related to the subduction of the CR that controls the sediment distribution and increases significantly the interplate coupling in the central zone. The north zone presents a well-defined sequence of 4 earthquakes rupturing the seismogenic interface during the 20 th century (Fig. 1) [START_REF] Beck | The Rupture Process of the Great 1979 Colombia Earthquake -Evidence for the Asperity Model[END_REF][START_REF] Kanamori | Variable Rupture Mode of the Subduction Zone Along the Ecuador -Colombia Coast[END_REF][START_REF] Mendoza | Seismicity Associated with the Great Colombia-Ecuador Earthquakes of 1942, 1958, and 1979 -Implications for Barrier Models of Earthquake Rupture[END_REF]. The regular shallow seismicity along the margin, based on worldwide seismic networks (Fig. 10), is concentrated immediately East of the trench mainly with thrust focal mechanisms. This supports that the subduction is partially locked north of 2°S as inferred from geodetic measurements and modelling [START_REF] Trenkamp | Wide plate margin deformation, southern Central America and northwestern South America, CASA GPS observations[END_REF]. The oceanic nature of the basement of the overriding plate in the Central and North Zone contributes to increase the interplate coupling because of its higher density. The contact between the top of the downgoing oceanic crust with strong material (velocities > 5 km/s) was interpreted as an estimate of the region of seismic onset [START_REF] Pacheco | Nature of Seismic Coupling Along Simple Plate Boundaries of the Subduction Type[END_REF]. In our case, the lack of shallow seismicity (<20 km) in the GG (Fig. 10) is consistent with the change in the velocity structure of the overriding plate and also possibly with a different process of basal erosion associated with fluid overpressure. We suggest this region is characterized by a low seismic coupling. The two large earthquakes that occurred in 1901 and 1953 are located near the northern and southern border of the GG (Fig. 10) and suggest a rapid change in seismic coupling to the north and to the south of the gulf.

Conclusions

Wide angle seismic models obtained by travel time inversion at the Ecuadorian convergent margin reveal variations in the crustal structure of the subducting Nazca plate and the overriding South American plate that strongly control the deformation of the margin and the seismic cycle :

(1) The North Zone exhibits normal oceanic crust with no evidence for the influence of the GHS. The overriding plate is characterized by a well-developed Cretaceous oceanic substratum, described onland as the Diabasico and Piñon formations, underlain by high seismic velocities (5.5-6.2 km/s). The structure of the margin suggests it was affected by transient phases of erosion and accretion. The inner part of the margin shows a large low velocity zone. Its upper limit correlates with a major splay fault that originates on the inter-plate contact zone and could have favoured the propagation of the 1958 great earthquake rupture [START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF].

(2) The Central Zone reveals an over-thickened oceanic crust (15 to 19 km) beneath the CR. The Cretaceous oceanic basement is well-developed and is likely to have been uplifted beneath the margin. The buoyancy of the Nazca plate and the high density of the margin basement both contribute to an increased interplate coupling. The fore-arc shows neither evidence of subsidence nor of frontal accretion.

(3) The South Zone is characterized by a 7 km thick oceanic crust, which was slightly influenced by the GHS during its emplacement 20 My ago. The continental slope and the platform show clues of subsidence, possibly related to intense tectonic erosion of the upper plate's base. We suggest that the low seismic velocities described for the margin basement could be related to an intense fracturing associated with the presence of fluids and low seismic coupling. Seafloor topography is derived from satellite altimetry and ship depth soundings [START_REF] Smith | Global sea floor topography from satellite altimetry and ship depth soundings[END_REF] completed by available swath mapping [START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF][START_REF] Flueh | Cruise Report SO159 -Salieri[END_REF]. It is contoured at 500 m intervals between -5000 and -500 m, and 25 m intervals between -250 and the coastline. Heavy solid black lines show locations of the three wide-angle seismic profiles used in this work. White circles: locations of OBSs along each profile. Bold circles: same as white circles for sections shown in Fig. 2 to 5. White triangles: land station locations. White line is wide-angle seismic refraction line published in [START_REF] Sallares | Seismic structure of the Carnegie ridge and the nature of the Galapagos hotspot[END_REF]. Stars: epicentres of the six great subduction earthquakes that occurred along the margin since 1901, including four in the northern part (3.5°N-0.5°S). The most important (1906, Mw=8.8) [START_REF] Kanamori | Variable Rupture Mode of the Subduction Zone Along the Ecuador -Colombia Coast[END_REF] has partially been reactivated three times from south to north: 1942 (Mw=7.9), 1958 (Mw=7.8) and 1979 (Mw=8.2) [START_REF] Beck | The Rupture Process of the Great 1979 Colombia Earthquake -Evidence for the Asperity Model[END_REF][START_REF] Mendoza | Seismicity Associated with the Great Colombia-Ecuador Earthquakes of 1942, 1958, and 1979 -Implications for Barrier Models of Earthquake Rupture[END_REF]. In contrast, no major earthquakes have been reported in the southern part of the study area (0.5°S-3°S), except two poorly documented events in 1901 and 1953. Dashed contours: rupture zone of the 1958 and 1979 earthquakes [START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF]. Bold arrow: relative convergence between Nazca Plate and North Andean Block from GPS data. DGM: Dolores-Guayaquil Megashear. Location of Cayo and Piñon formation outcrops and their equivalent in Colombia (Diabasico formation) are in red and orange (after the geologic maps of Ecuador and Colombia). (2) Corresponding velocity and Moho depth uncertainties (contours interval is set at 0.05 km/s). The velocity uncertainties correspond to the average of the absolute difference between the average velocity derived from the 100 final models at each node of the model and the velocity obtained for each of the 100 final models at the same node.

Figure and table captions

(3) Derivative weight sum (DWS) provides a quantitative representation of the ray coverage (in white) calculated in the best final velocity models (Fig. 5). Uncertainties are much larger in the lower part of the models (imaged by PmP and Pn phases only), than in the upper crust (~ 0.15 to 0.20 km/s) along all profiles. In profile SIS-4, the number of Pn and PmP picked was too small to allow a regular sampling of the crust-mantle boundary (as confirmed by the DWS) so the lower crust is poorly resolved. Uncertainty of Moho depths is usually lower than 5 km, the maximum is reached beneath the margin where the Moho is deeper and velocity uncertainties higher than 0.15 km/s. In the shallower part of the floating reflector (seaward each profile) the depth uncertainty is generally smaller than 3 km. The synthetic models consist of the addition of two sinusoidal velocity anomalies (8*5 km in the upper part of the model; and 25*12 km in the lower part) featuring periodic positive and negative amplitudes (up to ± 5%) to the best fitting velocity models (Fig. 5). The same source receiver geometry is used to compute synthetic travel times but random noise is added. Subsequently, these synthetic data sets are inverted, starting from the corresponding unperturbed models. (1) True model perturbation showing the amplitude and geometry of velocity anomalies; (2) recovery of anomalies obtained after four iterations. The Moho geometry is successfully recovered with an uncertainty of less than 1 km. Thus, we can suggest that the velocity-depth ambiguity raised by the tomography is insignificant in our study and that the data sets have enough spatial resolution to resolve velocity anomalies with the size and amplitude of the ones considered here. [START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF]. This MCS line is coincident with profile SAL-6 (see location in fig. 1). (top) Line drawing of line SIS-44 showing the splay fault, the decollement and the top of oceanic crust as interpreted from [START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF] and used in this study (Fig. 5). (bottom) Close-up showing the splay fault that is connected to a summit graben associated with a gentle rise of seafloor. The splay fault appears to connect to the decollement at depth. [START_REF] Dziewonski | Determination of earthquake source parameters from waveform data for studies of global and regional seismicity[END_REF]. On line data available at http://www.globalcmt.org/ Table 1: Profiles characteristics and model parameterization, showing the number of instruments, the number of picked arrivals, the model dimensions and the grid spacing chosen for the tomographic inversion. OBS/OBH spacing ranges from 5 km (SAL-6 and SAL-2) to 10-12 km (SIS-4).

Table 2:

Inversion parameters used to obtain the final tomographic models displayed in Fig. 5. Different combinations of correlation length values were tested, and we have noted that varying these parameters by 50% does not significantly affect the overall smoothness in the resulting final models. In order to force convergence, we choose to lower the smoothing and damping constraints values during iterations, dividing the inversion into two to three steps depending on the model considered. Several combinations of values were tested, but the best results were obtained with the procedure depicted in this table. The depth kernel weighting is unvarying (w=1) allowing an equal weight to both refracted and reflected phases during inversion.

Table 3: Quantification of the fit between observed and computed travel times for the best fitting (Fig. 5) and the mean models (Fig. 6).

Table 4: Synthetic view on the results obtained along the three profiles.

Table 5:

Worldwide margin wedge velocities from wide-angle surveys (the highest velocity 40 km from the trench is used to sort the data). This table is after the compilation of [START_REF] Christeson | Structure of the Costa Rica convergent margin, offshore Nicoya Peninsula[END_REF] with some additional recent studies. The low (< 3.0 km/s) surficial velocities associated with a slope apron are not included, for tomographic inversion this lower limit was arbitrarily taken at 3.5 km/s (velocities in parenthesis)

Figure 1 :

 1 Figure 1: General map of the eastern part of the Nazca plate and of the southern Colombia and Ecuador margin.Seafloor topography is derived from satellite altimetry and ship depth soundings[START_REF] Smith | Global sea floor topography from satellite altimetry and ship depth soundings[END_REF] completed by available swath mapping[START_REF] Collot | the Sisteur scientific party, Exploring the Ecuador-Colombia active margin and inter-plate seismogenic zone[END_REF][START_REF] Flueh | Cruise Report SO159 -Salieri[END_REF]. It is contoured at 500 m intervals between -5000 and -500 m, and 25 m intervals between -250 and the coastline. Heavy solid black lines show locations of the three wide-angle seismic profiles used in this work. White circles: locations of OBSs along each profile. Bold circles: same as white circles for sections shown in Fig.2to 5. White triangles: land station locations. White line is wide-angle seismic refraction line published in[START_REF] Sallares | Seismic structure of the Carnegie ridge and the nature of the Galapagos hotspot[END_REF]. Stars: epicentres of the six great subduction earthquakes that occurred along the margin since 1901, including four in the northern part (3.5°N-0.5°S). The most important (1906, Mw=8.8)[START_REF] Kanamori | Variable Rupture Mode of the Subduction Zone Along the Ecuador -Colombia Coast[END_REF] has partially been reactivated three times from south to north: 1942 (Mw=7.9), 1958 (Mw=7.8) and 1979 (Mw=8.2)[START_REF] Beck | The Rupture Process of the Great 1979 Colombia Earthquake -Evidence for the Asperity Model[END_REF][START_REF] Mendoza | Seismicity Associated with the Great Colombia-Ecuador Earthquakes of 1942, 1958, and 1979 -Implications for Barrier Models of Earthquake Rupture[END_REF]. In contrast, no major earthquakes have been reported in the southern part of the study area (0.5°S-3°S), except two poorly documented events in 1901 and 1953. Dashed contours: rupture zone of the 1958 and 1979 earthquakes[START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF]. Bold arrow: relative convergence between Nazca Plate and North Andean Block from GPS

Figure 2 :

 2 Figure 2: Examples of wide-angle seismic data, picked and calculated travel times and corresponding ray paths for the best fitting velocity model (Fig. 5a) along profile SAL-6 (North Ecuador) for OBH 116 (left) and OBS 121 (right). (top) Seismic record section for vertical geophone after data preprocessing (deconvolution, 5-15 Hz Butterworth filter, equalization); (middle) picked (solid circles with error bars) and predicted (open circles) travel times for first arrival phases (Pg and Pn) and refraction from the Moho (PmP); (bottom) Corresponding ray tracing for all identified seismic phases. The black arrows point at the trench location.

Figure 3 :

 3 Figure 3: same as Fig. 2 for profile SIS-4 (Carnegie Ridge), OBS 04 (left) and OBS 10 (right) for the best fitting velocity model (Fig. 5b).

Figure 4 :

 4 Figure 4: same as Fig. 2 for profile SAL-2 (Gulf of Guayaquil), OBH 45 (left) and OBS 57 (right) except hydrophone component is shown for OBH 45. Best fitting velocity model (Fig. 5c).

Figure 5 :

 5 Figure 5: Results of tomographic inversion of Pg, Pn and PmP arrivals. Iteration steps and RMS travel time misfits are displayed in table 2. Models are displayed from north to south : (a) profile 6, Northern Ecuador; (b) profile 4, Carnegie Ridge; (c) profile 2, Gulf of Guayaquil. White circles: OBS and land stations locations; solid black lines: location of the Moho interface inverted from PmP arrivals; white dashed line: base of sedimentary cover derived from coincident MCS data (Fig.8); black dots: top of the oceanic crust from coincident MCS data converted to depth; whites dots: top of the subduction channel (decollement) derived from coincident MCS data converted to depth; red dots: interpreted splay fault derived from coincident MCS data converted into depth; the vertical grey line correlate the trench location from one profile to another.

Figure 6 :

 6 Figure 6: Results of the Monte-Carlo analysis and ray tracing representation for (a) profile SAL-6; (b) profile SIS-4; (c) profile SAL-2. (1) Final velocity model derived from averaging all Monte-Carlo ensembles, the average models are computed by taking the mean values of the 100 final models at each point of the grid; the two dashed lines indicates the variation in Moho location for the 100 Monte-Carlo inversions.(2) Corresponding velocity and Moho depth uncertainties (contours interval is set at 0.05 km/s). The velocity uncertainties correspond to the average of the absolute difference between the average velocity derived from the 100 final models at each node of the model and the velocity obtained for each of the 100 final models at the same node.(3) Derivative weight sum (DWS) provides a quantitative representation of the ray coverage (in white) calculated in the best final velocity models (Fig.5). Uncertainties are much larger in the lower part of the models (imaged by PmP and Pn phases only), than in the upper crust (~ 0.15 to 0.20 km/s) along all profiles. In profile SIS-4, the number of Pn and PmP picked was too small to allow a regular sampling of the crust-mantle boundary (as confirmed by the DWS) so the lower crust is poorly resolved. Uncertainty of Moho depths is usually lower than 5 km, the maximum is reached beneath the margin where the Moho is deeper and velocity uncertainties higher than 0.15 km/s. In the shallower part of the floating reflector (seaward each profile) the depth uncertainty is generally smaller than 3 km.

Figure 7 :

 7 Figure 7: Results of the checkerboard tests for (a) profile SAL-6; (b) profile SIS-4; (c) profile SAL-2.The synthetic models consist of the addition of two sinusoidal velocity anomalies (8*5 km in the upper part of the model; and 25*12 km in the lower part) featuring periodic positive and negative amplitudes (up to ± 5%) to the best fitting velocity models (Fig.5). The same source receiver geometry is used to compute synthetic travel times but random noise is added. Subsequently, these synthetic data sets are inverted, starting from the corresponding unperturbed models. (1) True model perturbation showing the amplitude and geometry of velocity anomalies; (2) recovery of anomalies obtained after four iterations. The Moho geometry is successfully recovered with an uncertainty of less than 1 km. Thus, we can suggest that the velocity-depth ambiguity raised by the tomography is insignificant in our study and that the data sets have enough spatial resolution to resolve velocity anomalies with the size and amplitude of the ones considered here.

Figure 8 :

 8 Figure 8: Multichannel seismic reflection line SIS-44 across the margin wedge[START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF]. This MCS line is coincident with profile SAL-6 (see location in fig.1). (top) Line drawing of line SIS-44 showing the splay fault, the decollement and the top of oceanic crust as interpreted from[START_REF] Collot | Are rupture zone limits of great subduction earthquakes controlled by upper plate structures? Evidences from Multichannel seismic reflection data acquired across the northern Ecuador-southwest Colombia margin[END_REF] and used in this study (Fig.5). (bottom) Close-up showing the splay fault that is connected to a summit graben associated with a gentle rise of seafloor. The splay fault appears to connect to the decollement at depth.

Figure 9 :

 9 Figure 9: Velocity logs showing vertical velocity variation below the top of basement (assumed to start at 4 km/s) versus depth across the margin 25 km east of the trench (left) and 60 km east of the trench (right) for profile SAL-6 (red), profile SIS-4 (green) and profile SAL-2 (blue). The thickness of the curves represents the uncertainty of the velocity inferred from the previous Monte-Carlo analysis. (dashed lines): location of the interplate contact zone.

Figure 10 :

 10 Figure 10: Seismicity of the Ecuador-Colombia margin. (a) Hypocenters from 1973 to 2006 are from the USGS National Earthquakes Information Center. On line data available at http://neic.usgs.gov/neis/epic/epic.html. (b) Earthquake fault plane solutions from the Global Centroid Moment Tensor (CMTS) database from 1976 to 2006[START_REF] Dziewonski | Determination of earthquake source parameters from waveform data for studies of global and regional seismicity[END_REF]. On line data available at http://www.globalcmt.org/

Acknowledgements

The French scientific party was supported by Institut de Recherche pour le Développement (IRD), Institut National des Sciences de l'Univers (INSU) and Institut Français pour l'Exploitation de la Mer (IFREMER), which provided ship time during the SISTEUR experiment. The SALIERI-2001 experiment was funded by the German Ministry of Education, Research, Science and Technology (BMBF) under project no. 03G0159A. We are also indebted to H. Andersen, master of the R/V Sonne, and its crew for their support during the SALIERI experiment and to Captain R. Davila and the crew of R/V Orion (INOCAR) from which OBS were operated during the SISTEUR experiment. During the Sisteur experiment the R/V Atalante (chief scientist Jean-Yves Collot) shot over the OBS networks. Also special thanks to Y. Hello, A. Anglade and B. Yates who successfully operated the OBS at sea and to the scientific parties of both campaigns. Finally we thank Tim Minshull for his detailed review which helped us to improve the initial version of this manuscript.

Chile, Valparaiso 3.0-4.0 3.0-4.2 [START_REF] Ye | Crustal structure along the EDGE transect beneath the Kodiak shelf off Alaska derived from OBH seismic refraction data[END_REF] Taiwan, Ryukyu trench 3.5-3.7 3.5-4.3 [START_REF] Flueh | Seismic investigation of the continental margin off-and onshore Valparaiso, Chile[END_REF] Cascadia, Washington 3.5-4.2 3.2-4.4 [START_REF] Mcintosh | Crustal structure beneath the Nanao forearc basin from TAICRUST MCS/OBS Line 14[END_REF] Western Nankai Trough 3.5-4.4 3.5-4.5 [START_REF] Vannucchi | From seamount accretion to tectonic erosion: Formation of Osa Melange and the effects of Cocos Ridge subduction in southern Costa Rica[END_REF] Ecuador, Gulf of Guayaquil Cascadia, Oregon 3.5-4.0 3.5-4.5 [START_REF] Lewis | Structural Development of an Accretionary Prism by Thrust and Strike-Slip Faulting -Shumagin Region, Aleutian Trench[END_REF] Alaska, Shumagin