Mass Balance of Cenozoic Andes-Amazon Source to Sink System — Marañ ó n Basin , Peru

We investigate the mass balance of the Cenozoic Andes-Amazon source to sink system using rock uplift proxies and solid sedimentation of the Marañón Basin in Peru. The evolution of sedimentation rates is calibrated with regional structural restored cross-section. The quantification of eroded sediments from reliefs to sedimentary basin is achieved with ×10 Myr resolution and compared to present day proxies from the HYBAM (HYdrologie et Biogéochimie du Bassin Amazonien) Critical Zone Observatory. Erosion of the early Andean landforms started during the Upper Mesozoic period, but sediment rates significantly increase during the Neogene. This is in agreement with the calibrated increase of rock uplift in the Andean orogenic belt.


Introduction
Quantification of sedimentary basin is one of the numerous steps to understand the interaction and resulting product of the solid earth with the external envelopes of the Earth (Sediment Routing Systems; e.g., [1]).Foreland basins represent one of the repositories of the mountain building history and its evolution [2].Vertical movements and horizontal shortening of the sedimentary wedge [3] are the main topics of interest of foreland basin studies, whereas the movements of mass and flux associated with sedimentary signal represent of least studied aspects of these geological structures [4].Due to the focus of interest on the uplifted structures at the transition between the main relief and sedimentary basin, efforts have been made towards the qualification and quantification of erosion in the hinterland (e.g., [2,5]).The depositional record of the foreland sedimentation and erosion in a growing mountain range is associated with unconformities that make the deconvolution of geological processes difficult to establish [6].
The northern Central Andes of Peru are characterized by the absence of Altiplano, the flat slab subduction of the Nazca oceanic plate [7][8][9] and an extensive deformation of the eastern Andean orogenic wedge with an amount of shortening greater than 140 km [10].As a result, this region shows the easternmost propagation of the Subandean orogenic front [11] and record the foreland basin system [12] related to the rise of the Andean orogenic wedge and its erosion.
Geosciences 2018, 8, x FOR PEER REVIEW 2 of 17 Total suspended sediment yield exported from the current Andean orogenic wedge to the Amazon retro-foreland basin has been estimated by the HYBAM (HYdrologie et Biogéochimie du Bassin Amazonien) international network (http://www.ore-hybam.org/) at different latitudes (Figure 1A) [13][14][15][16].In northern Bolivia and Brazil, the foreland basin system is wide and approximately 30-50% of the exported mass of sediments is deposited in the foredeep depozone [13,16].The average production for the entire foreland basin catchment is estimated at 640 Mt•year −1 (±30%) by the latest measures of Vauchel et al. [16].In contrast, in the Amazonian Oriente foreland basin of Ecuador (Figure 1B), sediment budget along the Napo River shows an unusual increase in the suspended sediment flux from 12 Mt•year −1 to 42 Mt•year −1 [14].With a 9 years dataset, Armijos et al. [18,19] further confirmed this increase (8 to 42 Mt•year −1 ).This phenomenon is associated with active tectonics, uplift and erosion of the entire foreland basin systems.Sediment budget studies in northern Bolivia and Ecuador show, therefore, an extreme latitudinal variability in the Andean-Amazonian foreland basin systems, with subzones experimenting either uplift or subsidence.
South of the uplifted and eroded Oriente basin, the long term sedimentary record of the Marañón Basin is preserved through important sedimentary archives [17] (Figure 1C,D).From the Cretaceous to the Eocene period (Figure 1C) the Marañón Basin shows a depocenter located between the present day eastern slope of the Andes and the cratonic domains to the east.The Eocene to recent isopach map illustrates the change of dynamics of the sedimentation, with an important amount of sediments (>3000 m) focused in the foredeep depozone (Figure 1D).The sedimentation extends further out of the study area, as expressed by the evolution of the Amazonian Basin and its transcontinentalization since the Miocene [20].
Despite important work on present day/recent flux of sediments from the Andes to the Amazonian Basin, limited studies have been capturing the geological time scale sedimentation rates related to the growth and erosion of the Andes [21,22].Within the Central Andean foreland during the Miocene period to recent times [21], the sedimentation rate has significantly increased ~7.9-6 Ma (Figure 2).In northern Bolivia and Brazil, the foreland basin system is wide and approximately 30-50% of the exported mass of sediments is deposited in the foredeep depozone [13,16].The average production for the entire foreland basin catchment is estimated at 640 Mt•year −1 (±30%) by the latest measures of Vauchel et al. [16].In contrast, in the Amazonian Oriente foreland basin of Ecuador (Figure 1B), sediment budget along the Napo River shows an unusual increase in the suspended sediment flux from 12 Mt•year −1 to 42 Mt•year −1 [14].With a 9 years dataset, Armijos et al. [18,19] further confirmed this increase (8 to 42 Mt•year −1 ).This phenomenon is associated with active tectonics, uplift and erosion of the entire foreland basin systems.Sediment budget studies in northern Bolivia and Ecuador show, therefore, an extreme latitudinal variability in the Andean-Amazonian foreland basin systems, with subzones experimenting either uplift or subsidence.
South of the uplifted and eroded Oriente basin, the long term sedimentary record of the Marañón Basin is preserved through important sedimentary archives [17] (Figure 1C,D).From the Cretaceous to the Eocene period (Figure 1C) the Marañón Basin shows a depocenter located between the present day eastern slope of the Andes and the cratonic domains to the east.The Eocene to recent isopach map illustrates the change of dynamics of the sedimentation, with an important amount of sediments (>3000 m) focused in the foredeep depozone (Figure 1D).The sedimentation extends further out of the study area, as expressed by the evolution of the Amazonian Basin and its transcontinentalization since the Miocene [20].
Despite important work on present day/recent flux of sediments from the Andes to the Amazonian Basin, limited studies have been capturing the geological time scale sedimentation rates related to the growth and erosion of the Andes [21,22].Within the Central Andean foreland during the Miocene period to recent times [21], the sedimentation rate has significantly increased ~7.9-6 Ma (Figure 2).[8,9].F1 to F4 appoint long term foreland basin sequences defined in this study.Andean tectonics phases are sourced from [23,24].The sedimentation rates published nearby the study area [21,22].BFU: Basal foreland unconformity; In2: Inca phase 2; Q1 and Q2:3: Quechua phase.This increase is associated with coarse sediments of a fluvial megafan deposited under climate change and monsoon-like intensification.On the western side of the Andes, the sedimentation rate in the preserved record of the Lima Basin [22] show a longer sequence with cyclic variations of sedimentation and erosion (Figure 2).Three peaks are recorded at ~23 Ma, ~12 Ma and at the initiation of the Pliocene.Enhanced erosion in the Oligocene is documented in the Lima Basin but no documentation of this has been published in the foreland Amazonian foreland domain.We aim to fill this gap by our work.
In northern Peru, the latitudinal variability in the Andean-Amazonian foreland basin systems is well illustrated by the Marañón foreland systems architecture [12], which passes from a wedge-top depozone in the south to a foredeep depozone in the north (Figure 3).The orogenic wedge of the Marañón foreland basin systems is formed by the Eastern Cordillera and the Santiago, Moyabamba and Huallaga Subandean basins, and developed since the late Oligocene [9].[8,9].F1 to F4 appoint long term foreland basin sequences defined in this study.Andean tectonics phases are sourced from [23,24].The sedimentation rates published nearby the study area [21,22].BFU: Basal foreland unconformity; In2: Inca phase 2; Q1 and Q2:3: Quechua phase.
This increase is associated with coarse sediments of a fluvial megafan deposited under climate change and monsoon-like intensification.On the western side of the Andes, the sedimentation rate in the preserved record of the Lima Basin [22] show a longer sequence with cyclic variations of sedimentation and erosion (Figure 2).Three peaks are recorded at ~23 Ma, ~12 Ma and at the initiation of the Pliocene.Enhanced erosion in the Oligocene is documented in the Lima Basin but no documentation of this has been published in the foreland Amazonian foreland domain.We aim to fill this gap by our work.
In northern Peru, the latitudinal variability in the Andean-Amazonian foreland basin systems is well illustrated by the Marañón foreland systems architecture [12], which passes from a wedge-top depozone in the south to a foredeep depozone in the north (Figure 3).The orogenic wedge of the Marañón foreland basin systems is formed by the Eastern Cordillera and the Santiago, Moyabamba and Huallaga Subandean basins, and developed since the late Oligocene [9].[17] are projected on present-day cross-section (see location on the map of Figure 1).Earthquake location and depth are sourced from Devlin et al. [25].
Thrust deformation was controlled by interference of thick-skinned and thin-skinned tectonics related to a regional late Permian salt décollement.In its southern part, the Marañón foreland basin systems is well constrained by surface and sub-surface data, which facilitates constructing a threestep restoration calibrated from thermochronology [8] (Figure 3).The first stage corresponds to the middle Eocene restoration, which corresponds to the end of a period of strong erosion sealing an early Eocene compressive weak deformation.A considerable increase of the deformation has been evidenced since the middle Miocene.The total amount of shortening has been estimated at 70 km in the current stage.
The aim of this study is (1) to obtain a regional knowledge of the Marañón Basin sedimentary infill in terms of foreland basin evolution (sink), (2) constrain the Upper Mesozoic and Cenozoic chronostratigraphic framework to establish this evolution, (3) estimate the sedimentation rates for each long term sequence in the foreland basin system along with the exhumation history (source) to outline the mass balance of this portion of the Andean source to sink system during the Cenozoic.

Materials and Methods
Lithology and physical properties have been compiled from 21 exploration boreholes drilled within the Marañón Basin (Figure 1B).Biostratigraphic studies have been used to build the stratigraphic framework over the study area [26][27][28].We consider in this study the Albian, late Cretaceous and Cenozoic series deposited on the Andean Basal Foreland Unconformity (BFU; Figure 2).Wireline sonic logs have been used to depth convert the seismic framework and to compute the porosity depth profiles [29,30] from 15 wells of the Marañón Basin (Figure 4).Three main lithological classes  [17] are projected on present-day cross-section (see location on the map of Figure 1).Earthquake location and depth are sourced from Devlin et al. [25].
Thrust deformation was controlled by interference of thick-skinned and thin-skinned tectonics related to a regional late Permian salt décollement.In its southern part, the Marañón foreland basin systems is well constrained by surface and sub-surface data, which facilitates constructing a three-step restoration calibrated from thermochronology [8] (Figure 3).The first stage corresponds to the middle Eocene restoration, which corresponds to the end of a period of strong erosion sealing an early Eocene compressive weak deformation.A considerable increase of the deformation has been evidenced since the middle Miocene.The total amount of shortening has been estimated at 70 km in the current stage.
The aim of this study is (1) to obtain a regional knowledge of the Marañón Basin sedimentary infill in terms of foreland basin evolution (sink), (2) constrain the Upper Mesozoic and Cenozoic chronostratigraphic framework to establish this evolution, (3) estimate the sedimentation rates for each long term sequence in the foreland basin system along with the exhumation history (source) to outline the mass balance of this portion of the Andean source to sink system during the Cenozoic.

Materials and Methods
Lithology and physical properties have been compiled from 21 exploration boreholes drilled within the Marañón Basin (Figure 1B).Biostratigraphic studies have been used to build the stratigraphic framework over the study area [26][27][28].We consider in this study the Albian, late Cretaceous and Cenozoic series deposited on the Andean Basal Foreland Unconformity (BFU; Figure 2).Wireline sonic logs have been used to depth convert the seismic framework and to compute the porosity depth profiles [29,30] from 15 wells of the Marañón Basin (Figure 4).Three main lithological classes (clastic, carbonate, organic matter) have been quantified from boreholes drilled across the regional cross-section R1, R3 and R5.
For each structural cross-section and each foreland basin sequence, areas have been computed and clastic sedimentation rates have been calculated [32], and then applied in 3D in the Marañón Basin.Results of quantification are summarized in Table 1.Detailed aspects of the methods are within Appendix A.

Marañón Basin Foreland Sequences
Overall, we define and quantify 4 long term sequences based on biostratigraphy, lithological and provenance change:

•
Foreland sequence 1 (F1, Early Cretaceous (Albian) to late Cretaceous, 110-65 Ma) is the first long term sequence on the Cushabatay Basal Foreland Unconformity (Figure 2); it includes the Cushabatay-Raya, Agua Caliente-Chonta and Vivian-Cachiyacu-Huchpayacu fluvial to shallow marine cyclic transgressive sequences and is capped by the basal regional unconformity of the Paleocene Casa Blanca sandstones; Albian and Maastrichtian sandstones of this first sequence have a cratonic provenance as shown by the presence of Precambrian-inherited zircon grains [49].

•
Foreland sequence 2 (F2: Paleocene-early Eocene, 65-45 Ma) is the second long term sequence and is characterized by a continental environment (Figure 2); it recorded the first siliciclastic supply from the Andean uplift, and yield the first Cretaceous zircon grains [33]; it probably corresponds to the distal foreland deposits of the Marañón fold and thrust belt (presently preserved in the Western Cordillera uplift [23]).

•
Foreland sequence 3 (F3: Middle Eocene-Oligocene, 45-23 Ma) is the third long term sequence, which overlies the Basal Pozo regional erosional unconformity (Figure 2), interpreted as an unloading orogenic period [50].This stage of the foreland basin systems is represented in the sequential restoration in Figure 3.The detrital zircon grains are sourced from the Cenozoic volcanism and the uplifting of the Western Cordillera [49].

•
Foreland sequence four (F4: Neogene, 23-0 Ma) is the last long term sequence recording an acceleration of the deformation and sedimentation (Figure 3), and is contemporaneous with the uplift of the Eastern Cordillera as confirmed by zircon grains provenance [49].This period during the Miocene is characterized by lacustrine conditions alternating with episodes of fluvial drainage and marginal marine influence [51,52], and finally by a complex deltaic to fluvial environment.The present-day surface of the basin is marked by an important alluvial/fluvial network that conveys sediment from the Andes to the Amazon River, and finally to the Atlantic Ocean.

Marañón Basin Foreland Thickness
Three structural regional cross-sections (R1, R3 and R5), representative of the latitudinal evolution of the Marañón foreland basin, have been used to compute the 2D sedimentation rates.The four long term sequences have been delineated using the boundaries defined with seismic and wells.The cross-section R1 through the Huallaga-Marañón foreland shows that both basins are subjected to strong erosion.This amount of eroded thickness has been taken into account in our sediment budget.The cross-sections R3 and R5 represent the northern foredeep depozone of the Marañón Basin, where the foreland sequences are practically preserved.
The sediment thickness variation of the Marañón Basin (Figure 5) is illustrated with three SW-NE regional sections (R1, R3, R5).To the north, R1 section cross the Huallaga and Marañón wedge-top domains (Figure 1B, Figure 5A).The thickest deposits associated with sequence F1-F4 are within syncline structures with up to 7.12 km of sediments.Towards the NE, for the sequences F1 to F4, the Marañón Basin records a maximum thickness of 3.84 km of sediments that thins eastward to 2.56 km.Along R3 section (Figure 5B), the sequence F1-F4 is recorded at the front of the subandean wedge-top within the Marañón foredeep with 6.23 km of sediments.It is thinning towards the east with 2.56 km where the thin F1-F2 sequence (<60 m) at the Iquitos forebulge area is overlapped by a thicker F3-F4 sequence.Finally to the south of the study area, the R5 section shows the thickest sedimentary infill with 8.29 km for the sequences F1-F4 (Figure 5C).It thins gradually towards the NE to reach about 1.89 km, of which 1.67 km are associated with the F3-F4 sequences.

Physical Properties of the Marañón Basin Sediments
The physical properties and the computed results are summarized in Table 1.Porosity of all lithology decreases with depth with near surface porosity values from 25% to 60% and minimum value at the bottom of the drilled part of the basin ranging from 1 to <30% (Figure 4).The four foreland sequences have been sampled to correct for the porosity when computing the solid sediment estimates.The foreland sequence F1 median porosity is 18.9% (8% standard deviation).This sequence is dominated by clastic sediments with carbonate contribution restricted to the Chonta Formation.The foreland sequence F2 shows a slight increase of median porosity with 20% (10% s.d.).Its main lithology is clastics with a low fraction of organic matter.The foreland sequence F3 median porosity is 33% (4% s.d.) and it has the same lithological content as sequence F2.The foreland sequence F4 median porosity is 41% (5% s.d.) and its lithological content is dominated by clastics with the presence of 10% of carbonates and minor amounts of organic matter.

Quantification of Sedimentation Rates in the Foreland Basin Systems
The sedimentation rates along the cross-sections R1, R3 and R5 (Figure 5) and the whole Marañón Basin (Table 1) have been computed using the physical properties extracted from the subsurface data exposed above.
The sedimentation rate since the early Cretaceous (Albian, 110 Ma) based on the sections R1, and R3 show a decrease at the transition with the Paleogene (Figure 6), followed by an increase until recent times.On section R5 the sedimentation rates have increased since the early Cretaceous (Figure 6).Clastic deposition represents a higher fraction of sedimentation, with the figure never being lower than 75%.The carbonate or organics (coal, organic matter) do not represent a significant portion of the sediments that are preserved in the studied intervals (Table 1 and Appendix A) The sedimentation rate since the early Cretaceous (Albian, 110 Ma) based on the sections R1, and R3 show a decrease at the transition with the Paleogene (Figure 6), followed by an increase until recent times.On section R5 the sedimentation rates have increased since the early Cretaceous (Figure 6).Clastic deposition represents a higher fraction of sedimentation, with the figure never being lower than 75%.The carbonate or organics (coal, organic matter) do not represent a significant portion of the sediments that are preserved in the studied intervals (Table 1 and Appendix A) To the South of the Marañón Basin (R1, Figure 6), the computed sedimentation rates show a slight decrease less than two times from the late Cretaceous to the Paleogene.This is followed by an increase higher than twice from the middle Eocene to the late Oligocene periods for solid sediment and clastic rates of sedimentation.The highest figure of sedimentation, higher by two to three folds, To the South of the Marañón Basin (R1, Figure 6), the computed sedimentation rates show a slight decrease less than two times from the late Cretaceous to the Paleogene.This is followed by an increase higher than twice from the middle Eocene to the late Oligocene periods for solid sediment and clastic rates of sedimentation.The highest figure of sedimentation, higher by two to three folds, is occurring from the late Oligocene to the Neogene to recent sequence.This, despite the erosion observed in the section R1 (See Figure 5).Section R3 shows a stronger sedimentation rate during the Neogene.In sections R1 and R3, the sedimentation rate calculation of F2 is probably biased as a result of the Basal Pozo (base of F3) erosion that increases to the north (Figure 2).This may explain the computed sedimentation rates trend variation from R1 to R5 observed for the Paleogene.
Once these results on regional cross section are transferred to the whole Marañón Basin, taking into account the isopach maps computed regionally (Table 1), we can document the basin-wide evolution of sediments flux and the Qs computed (Figure 7A).
As a modern suspended sediment load figure, we refer to the work of Filizola and Guyot [53] based on sediment concentration monitoring of the Amazon River catchment at Obidos (Brazil) (Figure 1A), coupled with water discharge measurements using Doppler technology.The modern sediment discharge assessment is ranging from 880 to 1200 Mt•year −1 .In the Ucayali piggyback Basin, the reference represents the sediment flux deposited/preserved in the floodplain (Figure 7B), estimated at 220 Mt•year −1 by Santini et al. [15].
The figures of sediment supply (Qs) since the Albian have been a lower order of magnitude than the modern estimated of sediment discharge (Figure 7B).The Qs during the Cenozoic has been increasing by steps with two to three times the value of the previous sequence.The Neogene to recent show values 50 to 100 times lower than the present day observed sediment discharge of the River Amazon at Obidos, and are 18 times lower than the modern sedimentation rates in the neighboring Ucayali Basin.The Amazon Submarine fan volume shows an initiation during the Neogene of a paleodischarge of suspended material of 900-1100 Mt•year −1 [54] (Figure 7B).
Neogene.In sections R1 and R3, the sedimentation rate calculation of F2 is probably biased as a result of the Basal Pozo (base of F3) erosion that increases to the north (Figure 2).This may explain the computed sedimentation rates trend variation from R1 to R5 observed for the Paleogene.
Once these results on regional cross section are transferred to the whole Marañón Basin, taking into account the isopach maps computed regionally (Table 1), we can document the basin-wide evolution of sediments flux and the Qs computed (Figure 7A).S1).The red dotted vertical lines are events synthesized by Hoorn et al. [51] of the exhumation history of the Andes based on Apatite Fission Track Analysis (AFTA) age cluster analysis.S1).The red dotted vertical lines are events synthesized by Hoorn et al. [51] of the exhumation history of the Andes based on Apatite Fission Track Analysis (AFTA) age cluster analysis.

Andean Relief Exhumation History from AFTA
Over the last 110 Myr the exhumation history of the Andean orogenic domain surrounding the Marañón Basin shows three stages of evolution (Figure 7C).The first stage following the unconformity associated with the BFU lasted until the end of the Eocene, with little to no significant rate of exhumation.This stage is followed by the first significant exhumation event from the end of the Eocene to the Base of the Neogene where the number of sampled age of exhumation decreases.
The following stage marks the most important exhumation event of the Andean domain up to the Plio-Pleistocene.Denudation estimates from the AFTA modeling path that range from Northern Bolivia lead to an estimate of erosion since the Eocene of 0.16 ± 0.08 mm/year [34].Between 10 • S and 0 • , mass estimates and excess erosion of the Andean reliefs show a value ranging from 1 to 1.8 mm/year [55].These effects on mass balance of this discrepancy in denudation/erosion estimates will be discussed below and compared with the preserved volumetric in the sink.

Andean Foreland Sequence and Source to Sink Evolution
The Marañón FBS deposits can be divided into four long term sequences (F1 to F4) defined from previous studies, wells stratigraphic correlations and regional unconformities identified on seismic data.The first foreland sequence F1 is the Albian and late Cretaceous in age (Cushabatay Formation) and overlies the Andean basal foreland unconformity (BFU).It consists of three fluvial to shallow marine cyclic transgressive sequences that are probably deposited in a backbulge depozone.Albian and Maastrichtian sandstones have cratonic provenance [49].
The Paleocene-early Eocene long term sequence F2 recorded the first siliciclastic supply from the Andean uplift, and represents a continental distal foredeep to backbulge depozone (see R1 sequential restoration of Figure 3).The thickness of F2 decreases towards the south, which is probably due to a strongest erosion of the Basal Pozo erosional surface.This middle Eocene regional erosional surface sealed the deformation and configuration of the first stage of the Huallaga-Marañón FBS.It recorded a quiescence of thrust tectonics and an unloading orogenic period also described in Ecuador and Colombia [50].
The middle Eocene-Oligocene F3 long term sequence recorded the beginning of the second stage of the Huallaga-Marañón FBS, and the uplifting of the Western Cordillera [48] (Figure 1B).The Marañón Basin evolved in a complex environment alternating between marginal marine and fluvial environments.
During the Neogene, the F4 sequence recorded the uplift of the Eastern Cordillera and an evolution of the Marañón Basin from distal to proximal foredeep depozone with increasing rate of sedimentation (Figure 6), contemporaneous with the formation of the Iquitos forebulge [56,57].The early and middle Miocene period is still characterized by a complex environment alternating between marginal marine and fluvial environments [51,52], due probably to global high sea level and fast subsidence.The transcontinental Amazon River drainage occurred in the late Miocene [20,58] and provoked a radical change in the Marañón foredeep sedimentation that became exclusively fluvial [59].During the Pliocene, the eastward propagation of thrust tectonics combined the Fitzcarrald Arch uplift [7] led to the uplift and erosion of the southern part of the Marañón Basin.The strongest Neogene sedimentation rate is represented in the foredeep depozone of the section R3.

Mass Balance
In the northern Marañón Basin, sedimentation rates computed from the regional cross-section R5 show a progressive increase from the Albian to the present that seems to represent a classic evolution of the FBS [11].To the south (sections R3 and R1), the sedimentation rate shows a decrease during the Paleocene to middle Eocene sequence due to the middle Eocene erosional basal surface of F3 (Basal Pozo) that recorded a period of orogenic unloading interrupting the foreland basin sedimentation.A second Marañón foreland basin system started in the late Eocene-Oligocene.Sedimentation rates significantly increase during the Neogene, driven by a fast subsidence due to the orogenic loading of the uplifting Eastern Cordillera.The computed potential result of erosion of the source area related to the study area, with estimate a surface of 155 × 10 3 km 2 to a maximum of 210 × 10 3 km 2 , and the published erosion figures (minimum 0.08 mm/year and maximum 1.8 mm/year) [34,55], the volume corresponding to the Eocene to recent times (45 Ma) is ranging from 558 × 10 3 km 3 up to 17,010 × 10 3 km 3 .Compared to the preserved volume of sediments, ~420 ± 40 × 10 3 km 3 in the Marañón Basin within sequences F3 and F4 (this study), these figures represent 1.3 to 40 times more sediments potentially produced .The higher figure will mean that a significant volume of sediments will have been exported out of the study area.Yet no such anomalous volume of sediments exists east of the Andes.Even the most recent (~10 Ma Miocene) significant sedimentary feature of the south American continent, i.e., the Amazon submarine Fan represent only 700 × 10 3 km 3 [54].
The sedimentation rates calculated for the Neogene is 18 times lower than the modern sedimentation rates measured in the Ucayali Basin (Figure 7B).This difference of values is probably due to the complex evolution of the Marañón FBS that we have not completely deciphered.The identified potential erosion and unconformities in the Marañón Basin show that numerous events have occurred with non-preserved sedimentary archive at the location of study area.Part of this record could have been displaced outside or recycled within the study area, this is one of the most challenging aspects of the source to sink system and mass balance studies [60].It shows also that it remains important periods of erosion, hiatuses and no sedimentation to be found.

Figure 1 .
Figure 1.(A) Andean and Amazonian source to sink profile on the Napo and Madeira River basin, the flux is sourced from [13-16] (SEB, San Sebastian St.; FDO, Francisco de Orellana St.; ROC, Nuevo Rocafuerte St. and BEL, Bellavista St.); (B) Digital Elevation Map (DEM) map with stratigraphic and structural cross sections constructed for this study.Dotted lines represent the outline of the sedimentary basins: Marañón, Santiago, Huallaga and Oriente.Biostratigraphic studies from previous reports are showed in white filled circles; (C) Isopach maps of the Cretaceous-Early Paleogene periods and (D) Late Paleogene-recent [17].

Figure 1 .
Figure 1.(A) Andean and Amazonian source to sink profile on the Napo and Madeira River basin, the flux is sourced from [13-16] (SEB, San Sebastian St.; FDO, Francisco de Orellana St.; ROC, Nuevo Rocafuerte St. and BEL, Bellavista St.); (B) Digital Elevation Map (DEM) map with stratigraphic and structural cross sections constructed for this study.Dotted lines represent the outline of the sedimentary basins: Marañón, Santiago, Huallaga and Oriente.Biostratigraphic studies from previous reports are showed in white filled circles; (C) Isopach maps of the Cretaceous-Early Paleogene periods and (D) Late Paleogene-recent [17].

Figure 3 .
Figure 3. Huallaga-Marañón balanced cross-section (R1) and sequential restoration.Initial section (Middle Eocene) was obtained by flattening the base of the Pozo Formation (±45 Ma), which sealed the previous deformation and configuration of the Foreland Basin System (FBS).The intermediary stage corresponds to the Late Early Miocene stage and is calibrated by Apatite Fission Track AFT dating in the Chazuta Thrust and Biabo Anticline.AFT samples[17] are projected on present-day cross-section (see location on the map of Figure1).Earthquake location and depth are sourced from Devlin et al.[25].

Figure 3 .
Figure 3. Huallaga-Marañón balanced cross-section (R1) and sequential restoration.Initial section (Middle Eocene) was obtained by flattening the base of the Pozo Formation (±45 Ma), which sealed the previous deformation and configuration of the Foreland Basin System (FBS).The intermediary stage corresponds to the Late Early Miocene stage and is calibrated by Apatite Fission Track AFT dating in the Chazuta Thrust and Biabo Anticline.AFT samples[17] are projected on present-day cross-section (see location on the map of Figure1).Earthquake location and depth are sourced from Devlin et al.[25].

Figure 4 .
Figure 4. Physical properties of sediments for the whole drilled sections of the Marañón Basin, and clastic sediments porosities frequency of sequences F1 to F4.The sandstone and shale compaction curves are sourced from Sclater and Christie [31].

Figure 4 .
Figure 4. Physical properties of sediments for the whole drilled sections of the Marañón Basin, and clastic sediments porosities frequency of sequences F1 to F4.The sandstone and shale compaction curves are sourced from Sclater and Christie [31].

Figure 4 .
Figure 4. Physical properties of sediments for the whole drilled sections of the Marañón Basin, and clastic sediments porosities frequency of sequences F1 to F4.The sandstone and shale compaction curves are sourced from Sclater and Christie [31].

Figure 7 .
Figure 7. (A) Sediment supply and flux computed for the Marañón Basin over the last 110 Ma; (B) Modern sediment supply corresponds to the flux measured at Obidos [53].In the Ucayali Basin, the modern sedimentation rates represent 220 Mt•year −1 [15].Note the logarithmic scale to allow for "comparison" of the modern (high) and geological (low) rates.Amazon submarine fan flux are sourced from Wetzel [54] and references within; (C) The relative frequency exhumation history is based on a compilation of apatite fission track ages published in the vicinity of the study area (see Supplementary Materials TableS1).The red dotted vertical lines are events synthesized by Hoorn et al.[51] of the exhumation history of the Andes based on Apatite Fission Track Analysis (AFTA) age cluster analysis.

Figure 7 .
Figure 7. (A) Sediment supply and flux computed for the Marañón Basin over the last 110 Ma; (B) Modern sediment supply corresponds to the flux measured at Obidos [53].In the Ucayali Basin, the modern sedimentation rates represent 220 Mt•year −1 [15].Note the logarithmic scale to allow for "comparison" of the modern (high) and geological (low) rates.Amazon submarine fan flux are sourced from Wetzel [54] and references within; (C) The relative frequency exhumation history is based on a compilation of apatite fission track ages published in the vicinity of the study area (see Supplementary Materials TableS1).The red dotted vertical lines are events synthesized by Hoorn et al.[51] of the exhumation history of the Andes based on Apatite Fission Track Analysis (AFTA) age cluster analysis.

Table 1 .
Marañón Basin volumetric, physical properties of sediments and sedimentation rates computed from regional cross sections R1, R3 and R5 and isopach maps.