Substrate Entrainment, Depositional Relief, and Sediment Capture: Impact of a Submarine Landslide on Flow Process and Sediment Supply

Submarine landslides can generate complicated patterns of seafloor relief that influence subsequent flow behaviour and sediment dispersal patterns. In subsurface studies, the term mass transport deposits (MTDs) is commonly used and covers a range of processes and resultant deposits. While the large-scale morphology of submarine landslide deposits can be resolved in seismic reflection data, the nature of their upper surface and its impact on both facies distributions and stratal architecture of overlying deposits is rarely resolvable. However, field-based studies often allow a more detailed characterisation of the deposit. The early post-rift Middle Jurassic deep-water succession of the Los Molles Formation is exceptionally well-exposed along a dip-orientated WSW-ENE outcrop belt in the Chacay Melehue depocentre, Neuquén Basin, Argentina. We correlate 27 sedimentary logs constrained by marker beds to document the sedimentology and architecture of a >47 m thick and at least 9.6 km long debrite, which contains two different types of megaclasts. The debrite overlies ramps and steps, indicating erosion and substrate entrainment. Two distinct sandstone-dominated units overlie the debrite. The lower sandstone unit is characterised by: 1) abrupt thickness changes, wedging and progressive rotation of laminae in sandstone beds associated with growth strata; and 2) detached sandstone load balls within the underlying debrite. The combination of these features suggests syn-sedimentary foundering processes due to density instabilities at the top of the fluid-saturated mud-rich debrite. The debrite relief controlled the spatial distribution of foundered sandstones. The upper sandstone unit is characterised by thin-bedded deposits, locally overlain by medium-to thick-bedded lobe axis/off-axis deposits. The thin-beds show local thinning and onlapping onto the debrite, where it develops its highest relief. Facies distributions and stacking patterns record the progradation of submarine lobes and their complex interaction with long-lived debrite-related topography. The emplacement of a kilometre-scale debrite in an otherwise mud-rich basinal setting and accumulation of overlying sand-rich deposits suggests a genetic link between the mass-wasting event and transient coarse clastic sediment supply to an otherwise sand-starved part of the basin. Therefore, submarine landslides demonstrably impact the routing and behaviour of subsequent sediment gravity flows, which must be considered when predicting facies distributions and palaeoenvironments above MTDs in subsurface datasets.

Submarine landslides can generate complicated patterns of seafloor relief that influence subsequent flow behaviour and sediment dispersal patterns. In subsurface studies, the term mass transport deposits (MTDs) is commonly used and covers a range of processes and resultant deposits. While the large-scale morphology of submarine landslide deposits can be resolved in seismic reflection data, the nature of their upper surface and its impact on both facies distributions and stratal architecture of overlying deposits is rarely resolvable. However, field-based studies often allow a more detailed characterisation of the deposit. The early post-rift Middle Jurassic deep-water succession of the Los Molles Formation is exceptionally well-exposed along a dip-orientated WSW-ENE outcrop belt in the Chacay Melehue depocentre, Neuquén Basin, Argentina. We correlate 27 sedimentary logs constrained by marker beds to document the sedimentology and architecture of a >47 m thick and at least 9.6 km long debrite, which contains two different types of megaclasts. The debrite overlies ramps and steps, indicating erosion and substrate entrainment. Two distinct sandstone-dominated units overlie the debrite. The lower sandstone unit is characterised by: 1) abrupt thickness changes, wedging and progressive rotation of laminae in sandstone beds associated with growth strata; and 2) detached sandstone load balls within the underlying debrite. The combination of these features suggests syn-sedimentary foundering processes due to density instabilities at the top of the fluid-saturated mud-rich debrite. The debrite relief controlled the spatial distribution of foundered sandstones. The upper sandstone unit is characterised by thin-bedded deposits, locally overlain by medium-to thick-bedded lobe axis/off-axis deposits. The thin-beds show local thinning and onlapping onto the debrite, where it develops its highest relief. Facies distributions and stacking patterns record the progradation of submarine lobes and their complex interaction with long-lived debriterelated topography. The emplacement of a kilometre-scale debrite in an otherwise mudrich basinal setting and accumulation of overlying sand-rich deposits suggests a genetic link between the mass-wasting event and transient coarse clastic sediment supply to an otherwise sand-starved part of the basin. Therefore, submarine landslides demonstrably INTRODUCTION Submarine landslide deposits, olistostromes (Flores, 1955), or Mass Transport Deposits (MTDs) (Nardin et al., 1979), are sedimentary bodies that have been translated downslope from high to low gradient slopes as a result of mass failure and gravitational processes (Hampton et al., 1995;Moscardelli and Wood, 2008;Ogata et al., 2012;Festa et al., 2016;Kneller et al., 2016). The typically cohesive nature of the flows enables the transportation of large clasts (>4.1m diameter; herein named megaclasts, sensu Blair and McPherson, 1999) (Labaume et al., 1987;Payros et al., 1999;McGilvery and Cook, 2003;Lee et al., 2004;Jackson, 2011;Ogata et al., 2012;Hodgson et al., 2019;Nwoko et al., 2020a). Megaclasts within MTDs are sourced either from headwall areas or entrained from the substrate (Festa et al., 2016;Ogata et al., 2019). These features, accompanied by synand post-depositional faulting (Dykstra, 2005;Dykstra et al., 2011), generate the topographically irregular upper surfaces of MTDs (Moscardelli et al., 2006;Bull et al., 2009).
Deep-water sediment gravity flows interact with the rugose topography of MTDs, which influences flow behaviour, deceleration and steadiness (Lowe and Guy, 2000;Armitage et al., 2009;Jackson and Johnson, 2009;Fairweather, 2014;Ortiz-Karpf et al., 2015Steventon et al., 2021), and therefore dispersal patterns and depositional architecture . MTD surface relief has been shown to affect facies distribution and associated sedimentary architecture; this has been reported from both outcrop (Pickering and Corregidor, 2005;Armitage et al., 2009;Dykstra et al., 2011;Fallgatter et al., 2017;Brooks et al., 2018;Valdez et al., 2019) and subsurface studies (Ortiz-Karpf et al., 2017;Nwoko et al., 2020b). However, MTDs may continue to deform after initial emplacement through creeping processes (e.g. Butler and McCaffrey, 2010) or secondary mass movements (Sobiesiak et al., 2016). Furthermore, the high water content within newly deposited MTDs promotes active dewatering at their upper surface (Mulder and Alexander, 2001;Talling et al., 2012;Browne et al., 2020) associated with local instabilities and movement (Iverson, 1997;Major and Iverson, 1999;Van der Merwe et al., 2009). Fluids can also generate overpressure along with the interface between MTDs and its sediment cover, exploiting pathways created by internal MTD deformation (Ogata et al., 2012;Migeon et al., 2014). Therefore, the interaction between the initial topographic relief of MTDs, dewatering processes, post-depositional deformation and subsequent sediment gravity flows (and their deposits) is highly dynamic and inherently complex (e.g. Alves, 2015). A better understanding of sedimentary processes above MTDs can help subsurface predictions of facies distributions, which might have been overlooked due to variable seismic resolution and core coverage. Therefore, detailed field-based studies can help to bridge the resolution gap.
Here, we aim to understand an exceptionally well-exposed debrite and overlying sand-rich strata in the Bathonian Los Molles Formation, which were physically correlated over 9.6 km along a depositional dip transect in the Chacay Melehue depocenter (Neuquén Basin, Argentina). The objectives of this study are to 1) document the anatomy and stratigraphic architecture of the debrite, 2) investigate the impact of the dynamic upper relief on the overlying heterolithic and sand-rich strata, and 3) discuss the role that mass-wasting processes may have played as a trigger for subsequent sandrich sediment supply.

GEOLOGICAL SETTING
The Neuquén Basin is located in central-western Argentina and central-eastern Chile, covering an area of 160,000 km 2 ( Figure 1A). The basin is bounded to the north-east by the Sierra Pintada, to the south by the North Patagonian Massif, and since the Early Jurassic, by the early Andean magmatic arc to the west (Legarreta and Gulisano, 1989;Suárez and de la Cruz, 1997;Franzese and Spalletti, 2001;Howell et al., 2005). The Neuquén Basin contains a >6 km-thick sedimentary succession that spans the Mesozoic to the Late Cenozoic and records several unconformities related to tectonic phases (Vergani et al., 1995;Legarreta and Uliana, 1996;Howell et al., 2005). Three key tectonic phases are recognised (Vergani et al., 1995;Franzese and Spalletti, 2001;Franzese et al., 2003): 1) Triassic-to-Early Jurassic rifting and the onset of subsidence; 2) Early Jurassic-to-Early Cretaceous post-rift thermal subsidence associated with the development of the Andean magmatic arc and back-arc basin; and 3) Late Cretaceous-to-Early Cenozoic Andean compression and foreland basin development. In the western sector of the Central Neuquén Basin, the deep-to shallow-marine deposits of the early post-rift Cuyo Group (Lower-to-Middle Jurassic) (Gulisano et al., 1984) unconformably overlie the continental syn-rift volcano-sedimentary deposits of the Precuyano Group (Gulisano et al., 1984;Gulisano and Gutiérrez Pleimling, 1995;Legarreta and Uliana, 1996;Pángaro et al., 2009;Leanza et al., 2013) or the Palaeozoic basement of the Choiyoi Group (Llambías et al., 2003(Llambías et al., , 2007 (Figure 2A).
Our investigation focuses on the Early Bathonian stratigraphy of the Upper Los Molles Formation, which forms a ∼70 m thick interval characterised by ammonite-rich black shales and heterolithic successions comprising tuff layers with an intervening MTD and sandstone deposits ( Figure 1B).

METHODOLOGY
The sedimentology and stratigraphic architecture of a 70 m thick interval ( Figures 2C, 3) within the Upper Los Molles Formation were investigated along a 9.6 km long and WSW-ENE orientated outcrop belt ( Figures 1B,C). The succession dips 10-20°to the SE, with minimal structural overprint from the later tectonic inversion. Twenty-seven (27) sedimentary logs were measured at 1:25 to 1:40 scale along this transect (CML-0 to CML-27 from SW to NE) to document the broad depositional architecture of 4 different units (Unit 1, 2, 3, 4A and B) ( Figure 1C). Ten detailed logs were measured at a 1:2 scale at specific locations to document fine-scale thickness and facies changes. Four marker beds were used to build a robust physical correlation between sedimentary logs (Figures 2, 3). The marker beds are 1) Datum A, or the "Burro" marker bed, a light-grey indurated graded siltstone at the base of the study interval (Figures 2A, 3, 4A); 2) a gravelly thin- Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 757617 6 bed ( Figure 4F) and 3) a tuff layer ( Figure 4G), both within one of the studied units (Unit 4A); and 4) Datum B, a tuff layer overlying the study interval (100-150 m above) across the study area ( Figures 2D, 3). Uncrewed Aerial Vehicle (UAV) photogrammetry (Figures 2D,5) was used in conjunction with standard field techniques, such as mapping and logging, to capture the micro-and macro-scale features of the investigated stratigraphic units. Fifty-eight (58) palaeocurrent measurements were collected from ripples, cross-bedding, flame structure and convolute lamination vergence from bed tops of sandstones, and plotted in rose diagrams.

RESULTS: SUBDIVISION AND CHARACTERISATION OF SEDIMENTARY UNITS 1-4
The study interval is subdivided informally into four different units ( Figure 2C), based on their distinctive facies ( Table 1) and stratal relationships.
The central and eastern sectors contain a discrete stratigraphic interval that exhibits deformed bedding ( Figure 4B). This interval is thickest (at least 10 m; Figure 3) in the central sector, where, internally, it exhibits an array of imbricated decametre-scale east-verging thrusts (offset <2 m) and associated drag folds. The thrusts originate from a bed parallel surface, leaving the underlying bedding undeformed ( Figure 4B). In the eastern sector, a thin (∼5 m thick) interval of intense deformation is characterised by open folds and minor thrusts (offset <1 m) ( Figure 4B). Unit 1 stratigraphy in the western sector lacks any deformation.

Granular-to medium-grained sandstones
Normally-graded, well-sorted thinbeds. Foreset heights range from 5 to 7 cm, and angles vary between 10°a nd 35°. Erosional bases are common-sharp contacts, with planar base and undulatory top 5-10 cm Deposition and tractional reworking by turbidity currents (Tinterri, 2011) F7: Medium-bedded banded sandstones Banded sandstones with sharp alternation between darker and lighter bands. Lighter bands are grain-supported, while darker bands are matrix-supported and lack mudstone clasts Sandstones comprising alternation between matrix-poor light bands and matrix-rich dark bands (0.2-2 cm thick). Similar grain size (fine to medium) along with different bands. Heterolithic bedforms and pinch-andswell geometries can be developed. The bed bases can be structureless 10-50 cm Deposits beneath mud-rich transitional plug flow formed by tractional reworking within the upper stage plane bed flow regime (Baas et al., 2009(Baas et al., , 2011(Baas et al., , 2016Stevenson et al., 2020) F8: Thick-bedded banded sandstones with mudstone clast Banded sandstones with diffuse alternation between darker and lighter bands. Lighter bands are grain-supported, while darker bands are matrix-supported, with abundant mudstone clasts Sandstones comprising banding between matrix-poor light bands and matrix-rich mudstone clast (millimetric scale) bearing dark bands (0.5-2 cm).
Banding is diffuse and can be developed throughout the bed or from the middle to the top parts of a bed, commonly overlaid by convolute lamination. Laminae show local tilting and increasing spacing between laminae 50-150 cm Rapid aggradation and episodic damping of near bed turbulence due to clay flocs disaggregation (Lowe and Guy, 2000). Increasing spacing between laminae is attributed to growth strata due to foundering processes Erosion and incorporation of intrabasinal clasts. The entrained substrate was rapidly disaggregated within the flow resulting in clast-rich and clay-rich divisions at the bed top. The flows increased in concentration but had not developed stable density stratification (Haughton et al., 2003;Davis et al., 2009;Hodgson, 2009 turbidity currents (Allen, 1982;Trabucho-Alexandre et al., 2012;Könitzer et al., 2014), whereas the medium-bedded sandstones represent the deposits of medium-to high-density turbidity current . The laterally extensive character, mudstone dominance, and overall eastward (downdip) fining and thinning trend of thin-bedded sandstones of Unit 1 suggest deposition from low-energy sediment gravity flows in distal areas (e.g., Mutti, 1977), with possible distal lobe fringe deposits (Boulesteix et al., 2020). The discrete intervals of deformed bedding found in the upper parts of Unit 1 represent a post-depositional sheared zone linked to the overlying Unit 2.

Unit 2
Description: Unit 2 has an unconformable basal contact that truncates Unit 1 in the central sector ( Figure 3). The relief of the basal contact is characterised by down-and up-stepping segments (ramps, >2°) linked by bedding-parallel segments (flats). The average thickness of Unit 2 is 20-30 m but can locally reach >47 m in the central sector and abruptly thins to <8 m towards the eastern and western sectors ( Figure 3). This increase in thickness coincides with deeper erosion on the basal surface. Unit 2 is characterised by a matrix-supported medium-grained muddy sandstone to sandy mudstone and is very poorly sorted throughout, ungraded, and with a chaotic distribution of outsized clasts (F14; Figure 5). Clasts range in character and size from granular quartz grains and rounded volcanic epiclasts to much larger megaclasts (>4.1-140 m long)  of either conglomeratic or heterolithic lithology (Figures 2C,D,4F,5). The chaotic distribution of polymictic clast encased into a muddy sand matrix is responsible for the block-in-matrix fabric (e.g. Ogata et al., 2012). Typically, conglomeratic megaclasts are rounded, elongated and weakly deformed (Figures 2D, 4F,G) and are clast-supported, with well-rounded to sub-angular clasts (0.03-1 m diameter) and fragments of thick-shelled bivalves (oysters; Figure 4D). These oyster-bearing conglomeratic megaclasts are preferentially located near the base of Unit 2 ( Figures 2C,D, 3B). In contrast, heterolithic megaclasts are angular and characterised by internally folded packages of planar laminated and normally graded thin-bedded material ( Figures 2C, 5) and preferentially distributed toward the top of Unit 2 (Figures 3, 5). This heterogeneity promotes a homogeneous matrix-rich texture in the middle division.
Apart from the irregular basal contact of Unit 2, thickness changes within the unit are strongly controlled by the rugose upper surface. Kilometre-scale wavelength (1-3 km) and metrescale amplitudes (0.5-8 m) are responsible for a complex supra debrite topography.
Interpretation: The sedimentary characteristics of this unit, such as the chaotic distribution of (mega)clast floating onto a

Lithofacies
Lithology Description Thickness Process interpretation rations and erosional beds when lying above fine-grained intervals. Banding is locally developed at bed tops, alternating between matrix-poor light bands and matrix-rich dark bands (0.2-2 cm thick). Similar grain size (fine to medium) along with different bands. Heterolithic bedforms and pinch-andswell geometries can be developed high aggradation rates (Kneller and Branney, 1995;Sumner et al., 2008;Talling et al., 2012). The banding represents planar lamination (Bouma Tb division) (Stevenson et al., 2020) F12: Thick-bedded structureless matrix-poor sandstones with normallygraded mudstone clast Structureless sandstones with a mudstone clast at the base Structureless thick-bedded, mediumto coarse-grained, crudely normallygraded sandstones, with low-matrix content. They contain some mudstone clasts (0.1-1 m) with diffuse boundaries preferentially located at the base, which show coarse tail grading. Mudstone diapirs along the basal interface are common 0.5-2 m Deposition under high-density turbidity currents (sensu Lowe, 1982), formed by incremental layer-by-layer deposition with very high aggradation rates (Kneller and Branney, 1995;Sumner et al., 2008;Talling et al., 2012). Mudstone clast is entrained due to erosion of an unconsolidated debrite (Unit 2) and synsedimentary buoyancy product of density instabilities (Owen, 1987(Owen, , 2003 muddy sandstone to sandy mudstone matrix, suggest nearinstantaneous deposition from a flow with high yield strength and buoyant support that could transport clasts up to 140 m long (Stow and Johansson, 2000;Mulder and Alexander, 2001). Using Datum A, the debrite formed a long-wavelength mounded top (Figure 3), attributed to the parental flow's cohesive nature and en-masse freezing. We interpret Unit 2 as a cohesive debris flow deposit . The ramp and flat geometry at the base of Unit 2 indicate the debris flow's erosive nature (e.g., Lucente and Pini, 2003). The similarity in composition between the heterolithic megaclasts and underlying Unit 1 suggests entrainment of deep-marine substrate blocks due to the shear stress exerted by the overriding debris flow (Van der Merwe et al., 2009;Watt et al., 2012;Ogata et al., 2014a;Hodgson et al., 2019). In contrast, the sand content in the matrix and the oyster-bearing conglomerate megaclasts suggest a shallow marine origin of the mass failure (e.g., Ogata et al., 2012). Alternatively, the megaclast bearing shells could come from the remobilization of older slope strata, including shallow-marine deposits (La Primavera Fm; Suárez and de la Cruz, 1997;Llambías and Leanza, 2005). The Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 757617 two distinct megaclast sources suggest long-distance transport of clasts and flow bulking through local substrate entrainment (e.g., Sobiesiak et al., 2016).

Unit 3
Description: Unit 3 (0-4 m thick) is composed of thick sandstone beds (0.5-2 m) with sharp, irregular concave-up bases and abrupt pinchout terminations, which result in a disconnected distribution of packages of wedge-shaped sandstone bodies ( Figure 6) (see architecture section). Unit 3 is only present where Unit 2 is relatively thin (in the eastern-central and western sectors) and is absent in the central region where Unit 2 shows its maximum thickness (Figure 3). Where Unit 3 is absent, Unit 4 overlies Unit 2 ( Figure 3). Unit 3 comprises bed types characterised by two main amalgamated divisions (lower and upper divisions) with some grain size breaks lacking any mudstone-and siltstone-rich bounding intervals ( Figures 6A, 7). The basal interface of these sandstone bodies shows centimetrescale undulations characterised by abundant load casts, semidetached ball structures, and mudstone intrusions (diapirs and injectites) originating from Unit 2 ( Figure 6G). Two different types of thick-bedded amalgamated sandstone facies dominate the lower division, which varies along the transect. Grain-size breaks define amalgamation surfaces within sandstones. In the western sector, and more rarely in the eastern sector, the lower divisions are characterised by thickbedded (0.5-2 m thick), structureless, weakly normally-graded, moderately-to poorly-sorted sandstones (F12). At bed bases,  Table 1 and Panels 4, 6 for more detail. See Panel 8 for the lateral variability of each facies association.
The lower division of Unit 3 sandstone bodies are overlain by an upper division (up to 2 m thick), which comprises coarse to very fine-grained, normally graded, moderately-to poorly-sorted sandstones (0.5-1.7 m) ( Figure 6A). Banding can be developed throughout the bed or overlying a structureless division ( Figure 6A). The banding is characterised by an alternation between lighter matrix-poor bands and darker matrix-rich bands that comprise bedding parallel millimetric mudstone clast with sharp boundaries (mudstone clast type B) (F8; Figure 6C). Contacts between bands are diffuse ( Figures  6B,F). The spacing between the individual bands (0.5-2 cm) increases from the margin to central parts of the sandstone body (Figure 8), commonly showing rotation ( Figure 6E). These sandstones develop symmetrical and asymmetrical convolute lamination at bed tops (predominant vergence towards NE; Figure 6D). Decimetre-scale long and centimetre-scale thick mudstone injections can be observed within this division ( Figure 6B).
Interpretation: The wedge-shaped and deformed concave-up basal contacts of the sandstone bodies beds in Unit 3 are interpreted to reflect the interaction with the rugose upper surface and syn-sedimentary foundering of sand into the underlying mud-rich debrite (Figures 6, 9). Foundering is driven by instability due to the density contrast between the sand deposited above a less dense debrite (density loading) and lateral changes in sediment load (uneven loading) (Owen, 1987(Owen, , 2003 produced by the short-wavelength rugosity of the upper surface. The lack of sedimentary structures in the lower divisions of bed types recognised in Unit 3 is interpreted as a product of hindered settling from highly-concentrated gravity-flows, FIGURE 8 | Correlation panel focusing on the strata (Unit 3 and 4) overlying the debrite (Unit 2). Note the colour bar next to each log representing the facies associations (see Panel 5). Unit 3 is only present in western and eastern sectors whilst absent in the debrite high (CML-12). The lower part of Unit 4A thins and onlaps the debrite, while the upper one shows a larger lateral extent. Note the gravelly and tuff marker beds (red and white dashed lines, respectively). Unit 4A consists of an alternation between fine-grained lobe fringes and coarser healing lobe fringes in the western sector. Coarse-grained healing fringes pinch out, developing finegrained lobe fringes in the eastern sector. Unit 4B consists of amalgamated thick-bedded sandstone of lobe axis, thinning into medium-bedded dominate lobe off-axis environment in the western sector. The sand-rich lobe thins and fines towards the east, pinching out in the central sector and interfingering with the lobe fringe deposits of Unit 4A. Note that the pinch-out terminations are developed where the debrite relief is highest. The rose diagram shows the details of ripples (green), vergent convolute lamination (orange), Flame structures (yellow) and cross-bedding (red). Mean vectors of each type are shown, all suggesting a NE trend, except in Unit 4B, where the ripples suggest an E-directed palaeoflow indicating deflection processes. The perimeter of the rose diagrams corresponds to 100% of the value.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 757617 The shape of the thick-bedded sandstone bodies depends on the size and geometry of the thin-bedded megaclast. See Panel 5 for more detail. (H) Sketch of the architecture and evolutionary model of Type 2 architecture: While foundering, the sandstone might be protruded by the megaclast due to its higher competence than the surrounding debrite matrix.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 757617 resulting in turbulence damping and rapid deposition , inhibiting any period of traction (Sumner et al., 2008). The normally-graded lower divisions were produced by incremental layer-by-layer deposition from high concentration gravity flows, such as high-density turbidity currents (sensu Lowe, 1982). In contrast, the thick-bedded argillaceous sandstones with ungraded mudstone clasts observed in the distal areas (eastern sector) are interpreted as moderatestrength cohesive debrites (sensu Talling et al., 2012). The decimetre-scale mudstone clasts (type A) were transported due to the matrix strength of the debris flows and their positive buoyancy with respect to the encasing matrix. The lateral facies transition from high-density turbidites to moderate strength cohesive debrite suggest a flow transformation due to the entrainment of cohesive material from the underlying debrite (e.g., Kane and Pontén, 2012;Baker et al., 2017). The unconsolidated state of the debrite might have enhanced the substrate entrainment of decimetre-scale mudstone clasts (type A) and disaggregation (as indicated by the diffuse boundaries: Figure 6G), increasing the amount of mud and, therefore, the cohesiveness of the flow. Based on facies juxtaposition, the foundered sandstones can be subdivided into two different facies associations: 1) proximal and 2) distal, foundered sandstones facies associations (Figure 7). Both high-density turbidites and moderate strength cohesive debrites are characterised by rapid deposition (incremental deposition and en masse freezing, respectively), triggering the liquefaction of the fluid-saturated and unconsolidated upper surface of the debrite and foundering of sand ( Figure 9). The undulations of the concave-up basal interface reflect complex interactions with the substrate: as the denser sand sank into the fluid-saturated muddy substrate, the buoyancy of mud promoted the syn-to postdepositional intrusion (mud diapirs and injectites) of the substrate into sandstones. The most advanced stage of foundering is observed when detached sand-balls develop ( Figure 5; Owen, 2003;Tinterri et al., 2016).
In contrast, the banded sandstone characteristic of the upper divisions is interpreted to be formed under episodic near-bed turbulence damping at high rates of deposition (Lowe and Guy, 2000). The juxtaposition of the banded sandstones over the mudstone-clast bearing sandstones of the lower divisions suggests highly stratified flows, mixing and upwards transfer of centimetre-scale mudstone clast (type B) and the cohesive material from the disaggregation of the entrained decimeter-scale mudstone clast (type A). This enrichment in cohesive clayey material triggered the periodic suppression of turbulence and, therefore, banding development. The banding passes into convolute laminations towards the top, indicating moderate rates of deposition. The vergence of convoluted laminations suggests a syn-sedimentary shear-stress exerted by the overriding flow (McClelland et al., 2011;Butler et al., 2016) and flow-rebound produced by the underlying debrite relief (e.g., Tinterri et al., 2016).

Unit 4
Description: Unit 4 (10-27.3 m thick) has a sharp and concordant contact with the underlying Unit 2 and Unit 3 ( Figure 3). It comprises two subunits: a lower heterolithic interval (Unit 4A) and an upper sandstone-prone interval (Unit 4B; Figure 2C).
Unit 4A is thin-to medium-bedded (0.01-0.5 m; Figure 4H) heterolithic succession (F1, F2, F3 and F5) with a maximum thickness of 22 m, thinning to 8 m in the central sector above where Unit 2 is thickest (Figure 3). Most of the thin-beds (0.01-0.1 m thick) are fine-to medium-grained normallygraded sandstones, matrix-poor, moderately well-sorted, and structureless with common planar-parallel lamination and/or starved-ripple lamination near bed tops (F3; Figure 4C). Palaeocurrent measurement shows a consistent flow trend towards the NE (Figure 9). Unit 4A also comprises coarse-to granular normal-graded sandstones, relatively low matrix content with common grain size-breaks (0.07-0.2 m thick), erosive bases and sharp-planar tops (F4; Figure 4D) and two medium-bedded matrix-supported conglomerates with sandstone clasts (F14; 0.25 and 0.35 m thick, respectively) that pinch out towards the central sector ( Figure 8). One of these thin gravelly beds, which lacks any lateral thinning or fining trend ( Figure 4I), was traced across the exposure (gravelly marker bed; dashed red line in Figure 8). In addition, a 0.15 m thick tuff layer ( Figure 4J) was also used for correlation purposes (tuff marker bed; dashed white line in Figure 3). The medium-bedded sandstones (0.1-0.5 m thick) are structureless, ungraded, with planar-parallel and convolute lamination at bed tops, except one that shows cross-bedding (F6; Figure 4D). These sandstones have sharp bed bases and tops and lack mudstone clasts. In the eastern sector, Unit 4A is dominated by thin-to medium-bedded heterolithic succession that lacks any gravelly (F5) or matrix-supported conglomerate beds (F14).
Unit 4B (5.7 m thick in the western sector) thins eastwards along a 4.3 km transect until it pinches out, where Unit 2 is thickest (Figure 3). In the western sector, it dominantly comprises medium-(F5; Figure 4K) to thick-bedded sandstones (F11; Figure 4K), with less common "bipartite" sandstone beds (F9 and F10) composed of a matrix-poor lower division and a matrix-rich upper division with mudstone clasts. Conformable bases and sharp tops characterise the thick-bedded sandstones (0.5-1.2 m). Where the thick-bedded sandstones are not amalgamated and are intercalated centimetre-thick beds of fine-grained material (F1 and F2), bed bases are loaded locally. The thick-bedded sandstones are normally graded from medium to fine sand, well-sorted with rare centimetric mudstone clasts at the bed top. Soft-sediment deformation structures, such as centimetrescale flames with NE vergence, are also common at bed bases and along amalgamation surfaces ( Figure 4L). Banded sandstones are medium-bedded (0.1-0.5 m), fine-to medium-grained, and characterised by alternating between light-and dark-coloured bands, ranging from 0.2 to 2 cm thick (F7; Figure 4L). Both band types show a similar maximum grain size, although the darker bands are matrix-rich, and light bands are matrix-poor. Banding is generally sub-parallel to bedding. Although banded sandstones are more commonly associated with thick-bedded structureless sandstones, the banded sandstones can be individual event beds, with banding above the structureless basal division. The mediumbedded bipartite sandstone beds (0.1-0.5 m) consist of a medium-grained, matrix-poor and structureless lower division, which is overlain by a fine-grained matrix-rich upper division characterised by poor sorting and abundant mudstone clasts (0.05-0.3 m) with low sphericity and variable roundness (F9 and F10; Figure 4M). The lower and upper divisions show a gradual upwards increase in matrix content rather than across a sharp boundary. When Unit 4B is traced eastwards towards the central sector, the sandstone package transitions into a few thinbeds (0.1 m thick) of weakly graded, very poorly-sorted matrixrich sandstone, lacking the mudstone clasts observed in western areas. Unit 4B is absent in the eastern sector.
Interpretation: In Unit 4A, the thin sandstone beds showing planar and cross ripple laminations support an interpretation as low-to medium-density turbidites . The starved-ripple lamination observed in thin-bedded sandstones is interpreted as the reworking of sand deposited by dilute flows with low sedimentation rates (Talling et al., 2007;Jobe et al., 2012). The intercalation of thin-bedded sandstones with finergrained deposits suggests a lobe fringe environment (Lobe fringe facies association, Figure 5) (Prélat et al., 2009;Spychala et al., 2017b). The abundant coarse-grained to gravelly thin-bedded sandstones in the western sector record intermittent energetic coarse-grained flows, suggesting sporadic sediment bypass processes (Stevenson et al., 2015). However, the low matrix content within the granular beds suggests a sediment source area where only coarse-to granular grain size was available. The intercalation of such different facies suggests the juxtaposition of depositional environments of contrasting energy and/or different sediment sources. Either scenario could be possible given the complex sediment routing patterns and multiple transverse or axial sources available in the Neuquén Basin during the early post-rift setting (Vicente, 2005;Privat et al., 2021) and by analogy to other post-rift settings (e.g., Lien, 2005;Fugelli and Olsen, 2007;Hansen et al., 2021). The mass failure would trigger a new coarse-grained source due to slide scar position and geometry, promoting intermittent sand supply to an otherwise sand-starved environment (see "Origin and role of the Mass-wasting process as a trigger for turbidite systems development" in discussion). The downdip variability in the thickness of Unit 4A (from 22 to 8.5 m thick), reduction in gravelly sandstone content and the stratigraphic thinning between the granular marker bed (red dashed line in correlation) and the top debrite (Unit 2) reveals the existence of subtle relief on the debrite surface ( Figure 3B). Furthermore, the two poorly sorted ungraded muddy sandstones, which are interpreted to be debrites due to their chaotic distribution of clast within the argillaceous matrix, also pinch out towards the central sector. The ripples and convolute laminae with SW vergence (Figure 3) contrast with the consistent NE paleoflow, suggesting local flow deflection (cf. Tinterri et al., 2016) in the central sector, where the debrite relief is highest, indicating the interaction between sediment gravity-flows and the upper surface of the debrite.
Massive medium-to thick-bedded deposits of Unit 4B are interpreted as high-density turbidites formed by incremental layer-by-layer deposition with high aggradation rates (Kneller and Branney, 1995), interpreted to represent proximal lobe axis environments (lobe axis facies association; Figure 5) (e.g., Prélat et al., 2009;Kane et al., 2017). The location of these facies in the westernmost sector, and the palaeoflow measurements, suggests that the western sector was relatively proximal. Banded sandstones represent the deposits of mud-rich transitional flows formed by tractional reworking (Stevenson et al., 2020). The bipartite beds consisting of a basal structureless to planar laminated sandstone division, overlain by a linked mudstone clast-rich upper division are interpreted as hybrid event beds (HEBs), formed from transitional flows deposited under highdeceleration rates Hodgson, 2009;Kane and Pontén, 2012) in more distal environments than the banded sandstones (Stevenson et al., 2020). The gradual and diffuse boundary between the basal turbidite and the upper debrite suggest vertical segregation of particles within the cohesive flow . The facies evolution of Unit 4B from proximal (western sector) to distal (eastern sector) of thickbedded sandstones into hybrid event beds likely represents the downdip transition from lobe-axis/off-axis environments (lobe axis facies association: Figure 5) (sensu Prélat et al., 2009) into lobe-fringe environments (lobe fringe facies association; Figure 5) (e.g., Spychala et al., 2017a;Kane et al., 2017), persisting until the frontal/oblique pinchout (e.g. Hansen et al., 2019).

Large Scale Architecture: Debrite Relief
Using Datum A ("Burro" marker bed), the upper surface of the >9.6 km long debrite forms a broad convex-up relief that reaches a maximum in the central section coincident with the deepest incision (at least 22.5 m of erosional relief; Figure 3B). The spatial association of the thickest part of the debrite with the deepest incision support a genetic link between the geometry of the flatramp-flat shaped basal shear zone and the mounded top. The morphology of the basal surface can buttress material translated downslope and develop positive topographic features, such as pressure ridges (Moscardelli et al., 2006;Bull et al., 2009). Bed-bybed correlation within Unit 4A shows that where the upper surface of the debrite develops the highest relief (∼8 m of positive relief with respect to the western sector), Unit 3 is absent, and it is overlain by Unit 4A (Figure 3), showing a laterally continuous stratigraphic interval with metre-scale thickness variations (Figure 8). Unit 4A thins from 22 m (CML-1) and 13 m (CML-2) to 6 m (CML-12) across the highest part of the debrite ( Figure 3B). The lower part of Unit 4A pinches out in the central sector, developing onlaps of individual beds, and supporting the existence of a gentle relief (Bakke et al., 2013;Soutter et al., 2019). In contrast, the upper part of Unit 4A shows tabular architecture with a lateral continuity of over 7 km. Unit 4A is overlain by Unit 4B, which shows a progressive thinning of the submarine lobe from the western to the central sector over 5.6 km, from 5.7 m (CML-3) to 1.7 m (CML-12) and 1 m (CML-14) with a mean thinning rate of 0.9 m/km. The submarine lobe pinches out between CML-14 and CML-22 (<2 km), interfingered with unit 4A ( Figure 3B). The lack of onlap geometries against underlying deposits and subtle thinning rates consistent with unconfined settings (e.g., Prélat et al., 2009) suggests a lack of a pronounced pre-existing relief. However, the coincidence of lobe pinch-out in the area where the debrite relief is highest and where the underlying Unit 4A is thinnest might reflect subtle residual relief.

Small-Scale Architecture: Foundered Sandstones
The steeply-dipping unconformable base, internal deformation and abrupt thickness changes of Unit 3 sandstones contrast with their flat and conformable tops (see stereoplots in Figure 6). These sandstone bodies can be subdivided into three different types by their architecture.

Type 1
Description: The thinner foundered sandstone bodies range between 0.5 and 2 m thick and are only formed by the banded sandstones ( Figures 9A,B). They are characterised by 5-25 m wide lenticular shapes, with thickness/width ratios varying from 1:5 to 1:18. These sandstone bodies show relatively constant thinning rates (∼0.25 cm/m) towards their pinch outs. They are characterised by: 1) advancing onlap terminations onto Unit 2 at the base, and a vertical change into; 2) progressive rotation of laminae and the wedging of the sandstones ( Figure 9C).
Interpretation: The onlap termination indicates the interaction between the parental sediment gravity flow and pre-existing debrite-related relief (e.g., Bakke et al., 2013), and the sediment load was insufficient to trigger the soft-sediment deformation along the upper surface. In contrast, the overlying rotation and wedging represent growth strata associated with the syn-sedimentary foundering. This juxtaposition of terminations indicates that foundering did not start since the onset of deposition of sand due to insufficient stress to trigger the softsediment deformation. This supports an incremental layer-bylayer deposition of these sandstones rather than the freezing of the parental flow.

Type 2
Description: Thick-bedded foundered sandstones (up to 4 m thick) are characterised by irregular stepped bases and abrupt thickness variations (up to 2 m thinning over 1 m laterally) ( Figure 9D). They are composed by the juxtaposition of two different divisions: 1) lower and 2) upper divisions ( Figure 9E). The lower divisions comprise structureless sandstones with poorly-developed amalgamation surfaces (F12 or F13). They rarely exceed 10 m laterally and 3 m in thickness (thickness/ width 1:2 to 2:1) and are characterised by both abrupt onlaps terminations and wedging. In contrast, the upper divisions are characterised by banded sandstones (F8), which are more laterally extensive than the underlying division, with a maximum length of 50 m and rarely exceed 1 m in thickness, and thin laterally towards margins (thickness/width 1:10).
Interpretation: The coexistence of onlaps terminations and wedging indicate that the foundering began at the onset of deposition and the existence of pre-existing topography along the upper surface. The sediment load was enough to trigger the foundering because the debrite relief strongly influenced the initial high-concentration flows, which promoted a loss in flow capacity and deposition under high aggradation rates. The rapid deposition and foundering are responsible for the poorlydeveloped onlap terminations and amalgamation surfaces. The deposition of lower division deposits promoted a reduction in debrite rugosity due to the infilling of topographic lows, enabling the deposition of laterally more extensive deposits. The rotation and wedging in the banding of the upper division is less well developed than in Type 1 sandstone bodies. This suggests a progressive reduction in syn-sedimentary deformation and an increase in seafloor stability (e.g., Owen, 1987Owen, , 2003.

Type 3
Description: These sandstone bodies show similar facies juxtaposition as in Type 2. In this case, the sandstones terminate against heterolithic megaclasts due to their preferential location towards the top of the debrite (Figures 5,9F). In these cases, the geometries of the foundered sandstones diverge from the concave-up geometry, dependent on the shape of the megaclast. Some megaclasts disconnect bodies laterally, whereas others only impact the base, with the top part of the sandstones undisturbed ( Figure 9G).
Interpretation: The fluid-saturated matrix and rigid megaclast respond differently to the shear stress exerted by the deposition of sand creating. This differential compaction (Ogata et al., 2014b) controls the architecture of the foundered sandstones, creating complex bodies ( Figure 9H).

Basal Shear Zone and Impact on the Substrate
As submarine landslides travel across the seafloor, they exert shear stress on the substrate, coupled with significant overpressure (Bull et al., 2009;Hodgson et al., 2019;Payros and Pujalte, 2019). This leads to substrate entrainment (Eggenhuisen et al., 2011;Hodgson et al., 2019) and/or deformation (Butler and McCaffrey, 2010;Watt et al., 2012;Dakin et al., 2013;Ogata et al., 2014b). The debris flow (Unit 2) incised at least 22.5 m into the substrate (Unit 1; Figure 3). In the central sector, the basal shear surface forms ramps (up to 800 m long, >2°) and flats (up to 1,550 km long; Figure 3B; between the logs CML-9 and CML-10) (see Lucente and Pini, 2003;Frey Martinez et al., 2005 for flatramp-flat geometry). The stress applied to the substrate during the emplacement is accommodated by both stratigraphic intervals consisting of deformed packages (basal shear-zone) and interfaces consisting on a plane (basal shear-surface) (Alves and Lourenço, 2010), such as the discrete basal shear zone located in upper Unit 1. The absence of contractional features in the deposits underlying Unit 1 supports the deformation as the result of the shear stress produced by the Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 757617 emplacement of the debrite rather than tectonism. The basal shear zone has variable thickness and deformation styles. It is absent in the western sector, whereas erosion and deformed intervals record a high degree of basal shear stress in the central sector ( Figure 3B). In the central sector, the deformed package (up to 10 m thick) is characterised by decametre-scale thrusts with metre-scale offsets and drag folding ( Figure 4B). The predominance of imbricate thrusting over folding, and lack of internal disaggregation within the package, indicate competent substrate rheology (e.g., Van Der Merwe et al., 2011). The eastward vergence of the compressional structures ( Figure 4B) indicates an eastward emplacement direction for the debris flow (Twiss and Moores, 1992), consistent with the palaeoflow indicators in the bounding strata. The thrusting is attributed to bulldozing by the entrenched debris flow (e.g., Jackson, 2011;Hodgson et al., 2019;Payros and Pujalte, 2019), representing the initial stage of substrate entrainment. Entrainment of megaclasts into a debris flow has been reported in other systems in the subsurface (Moscardelli et al., 2006;Alves and Cartwright, 2009;Sawyer et al., 2009;Dakin et al., 2013;Ortiz-Karpf et al., 2015;Soutter et al., 2018;Nwoko et al., 2020a) and more rarely at outcrop (Ogata et al., 2014b;Sobiesiak et al., 2016;Hodgson et al., 2019;Cumberpatch et al., 2021). The progressive increase in thickness and degree of strain along the basal shear zone of Unit 2 and the enrichment in rafted heterolithic megaclasts ( Figures 3B, 5) suggest downdip evolution of the debris flow, which might have affected the parental debris flow rheology (e.g., Hodgson et al., 2019;Payros and Pujalte, 2019) and bulking of the flow (Gee et al., 2006;Alves and Cartwright, 2009;Butler and McCaffrey, 2010;Hodgson et al., 2019). The preferential location of heterolithic megaclast towards the top of Unit 2 might be related to internal granular convection cells created along with the debris flow, enhanced by the buoyancy of less dense rafted megaclast compared to debrite matrix  and kinetic sieving (Legros, 2002). In contrast, conglomerate megaclasts are always found at the base of the debrite due to their higher density than the surrounding debrite matrix.

Dynamic Debrite Topography and Impact on Overlying Strata
The absence of Unit 3 sandstones over the thickest part of the debrite suggests that the sediment gravity flows were strongly stratified and influenced by the debrite relief ( Figure 8). The sandstone bodies are also disconnected at finer scales, revealing short wavelength (metre-scale) and amplitude (decimetre-scale) rugosity on the debrite surface. The existence of simultaneous short wavelength and amplitude rugosity superimposed on a large-scale wavelength relief on the upper surface of an MTD has also been reported by Armitage et al. (2009), defined as "surface-topography hierarchy," in the Cretaceous Tres Pasos Formation at the Sierra Contreras (Chile) and by Fairweather (2014) in Carboniferous Paganzo Basin at Cerro Bola (Argentina). In this study, the deposition of sand in preexisting lows filled the short-wavelength rugosity and triggered the loading of individual sandstone bodies onto the mud-rich debrite (See "Small-scale architecture: Foundered sandstones" sections), leaving the large-scale relief underfilled (Figures 3,  8). The foundering process is evidence of substrate liquefaction and highlights the dynamic interface between the debrite and subsequent flows and their deposits. A similar scenario was proposed by Van der Merwe et al. (2009), Van der Merwe et al. (2011) in the Vischkuil Formation in the Laingsburg depocentre (Karoo Basin).
The ability of supra MTD rugosity to pond turbidity currents travelling across their upper surface is a well-known phenomenon . However, the presence of Unit 3 foundered sandstones up-dip and down-dip of the debrite high ( Figure 3B), and its consistent NE paleocurrent trend, suggest connected sediment transport routes across the debrite with no evidence of flow ponding or stripping (e.g. Armitage et al., 2009;Fairweather, 2014). The highly-stratified grown-hugging parental flows of Unit 3 would have been ponded in proximal parts (western sector) if a fully enclosing topography existed given their reduced ability to surmount obstacles (Al-Ja' Aidi et al., 2004;Bakke et al., 2013), resulting in sand starvation over the debrite in distal settings (Sinclair and Tomasso, 2002;Kneller et al., 2016). The overlying Unit 4A can be traced laterally across the study area, with metre-scale thinning where the debrite is thickest (see CML-12; Figure 3). Apart from this, the advancing onlap geometries of the thin beds and the divergence in the overall NE-orientated paleocurrents (rose diagram Unit 4; Figure 3) indicate the progressive healing of the large-scale wavelength debrite relief, with some deflection of turbidity currents ( Figure 10).
The thin sandstone beds of the upper part of Unit 4A healed the debrite high. However, the gravelly beds thin and fine from proximal to distal areas (western to eastern sectors), and the two debrites pinch out in proximal areas (western sector), suggesting subtle remnant topography (Figure 8). The different lateral continuity of individual beds is explained by different rheologies of individual sediment gravity flows, which affect the flow efficiency (Al Ja' Aidi et al., 2004). Cohesive debris flows are more influenced by irregular relief, while low-density turbidity currents are less affected by seafloor topography (Bakke et al., 2013;Soutter et al., 2019). This suggests that laterally continuous thick accumulations of lobe fringe deposits can develop on gentle topographies, while the submarine lobes' axial parts were restricted to lower relief areas. The interaction of thin-bedded turbidites successions with gentle topography has also been reported in other deep-water settings (i.e., "aggradational lobe fringes"; Spychala et al., 2017b). The deposition of lobe fringe successions reduced confinement, which enabled the deposition of the Unit 4B submarine lobe.
The submarine lobe is characterised by a progressive thinning and fining, developing pinch-out geometries and interfingering with Unit 4A in the area where the relief of the debrite is highest. The development of pinch-out geometries over the areas where the debrite shows a mounded relief and where the Unit 4A onlaps and thins suggests that the relief was not wholly healed with the deposition of Unit 4A and affected the parental flows of Unit 4B.
One explanation is that the exposure exhumes the Unit 4B lobe obliquely, with the medium-to thick-bedded lobe axis deposits in the westernmost sector (CML-1 to CML-4) transitioning into an HEB-dominated fringe, being highly impacted by gentle seafloor topography (Soutter et al., 2019;Privat et al., 2021). Alternatively, the seafloor relief could have promoted the modification of flow pathways and deflection of flows, thus changing the downdip orientation ( Figure 8). All these scenarios suggest a confined and uncontained (see Southern et al., 2015) lobe-type depositional system. The precise dispersal pattern of the flows remains unknown due to the outcrop limitations. Nonetheless, the documented stratigraphic evolution reveals that long-lived debrite relief and progressive healing by deposition of aggradational lobe fringes enabled the progradation of sand-rich submarine lobe, albeit with changes in flow rheology the bed style and element-scale pinchout (Figure 8).
Another explanation is that the debrite relief in the central sector might have been rejuvenated through volume changes in the debrite due to differential compaction by loading the lobe itself in the proximal sector. However, it seems unlikely that the deposition of a 5.7 m thick lobe could promote a volume loss in an 8-47 m thick debrite, given that both units are separated by 8-22 m thick thin-bedded interval. In contrast, another two hypotheses are deemed more plausible; fluid loss and fault-controlled mechanical subsidence. Fluid loss-controlled evacuation could have promoted differential subsidence of the upper surface of the debrite and overlying units (e.g., Browne et al., 2020). Alternatively, given the early post-rift setting, mechanical subsidence by an east-facing and N-S striking fault (Manceda and Figueroa, 1995) could have generated more accommodation in the western part of the study area (see Cristallini et al., 2006). However, this implies a very localised and rapid reactivation, and there is no other evidence for postrift tectonism identified.

Origin and Role of the Mass-Wasting Process as a Trigger for Turbidite Systems Development
The emplacement of the >9.6 km long and erosional debrite in the Chacay Melehue depocentre reflects an abrupt change in sedimentation patterns, which were previously dominated by dilute mud-rich flows (Unit 1). The first significant sand influxes in the depocentre for ∼6 Myr since the Aalenian (interval II of Gulisano and Gutiérrez Pleimling, 1995) are recorded by the sand-rich deposits (Unit 3) immediately overlying the debrite (Figure 10). The juxtaposition of sandrich turbidites over debrites (metres to hundreds of meters thick) have been reported in other systems (Kleverlaan, 1987;Labaume et al., 1987;Payros et al., 1999;Fallgatter et al., 2017). These authors suggest that the debris flow underwent a period of mixing with ambient water, leading to the generation of an overriding cogenetic turbidity current. The foundering phenomenon reported here reveals a close spatiotemporal relationship between the debrite emplacement (Unit 2) and overlying sandstone deposition (Unit 3). An alternative mechanism is that the mass-failure event altered the basin margin physiography such that a sand source was captured. Mass-wasting processes responsible for the evacuation of material from shelf edge and upper slope areas alter the bathymetric configuration of basin FIGURE 10 | Down-dip oriented schematic diagram illustrating the relief created by the debrite and the impact on younger sand-rich units. Foundered sandstones fill the small-scale rugosity, leaving the kilometre-scale accommodation underfilled. The submarine gravity flows are deflected by long-lived subtle debrite-related relief (right block). Partial healing and drapping of the debrite with the progradation of submarine lobes, which are gently impacted by the long-lived inherited relief. Note the white dashed line representing the correlation shown in Panels 3, 8.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 757617 margins and promote the funnelling of sediment stored in shallow marine environments through slide scars (e.g., Moscardelli and Wood, 2008;Kneller et al., 2016;Steventon et al., 2020) (Figure 11). The role played as a trigger mechanism for sand delivery into deep-water setting by the mass-wasting event is a plausible scenario given the sandstarvation recorded in the coeval deposits of Los Molles Fm along the eastern margin of the Chacay Melehue depocentre (Veiga et al., 2013). Given the palaeoflow and kinematic indicators, the thickness patterns of the studied units ( Figure 3B), and previous studies on sediment supply from the volcanic arc (Vicente, 2005), we propose that the mass failure originated to the west of Chacay Melehue, where a major syn-rift fault is located close to the volcanic arc (<30 km; Manceda and Figueroa, 1995;De La Cruz and Suarez, 1997;Vicente, 2005). The role of the western volcanic arc as a source area for the early post-rift sediment supply in the Chacay Melehue depocentre is supported by the southeastwards directed paleocurrents measured in the Aalenian turbidite system at the base of the Los Molles Formation (Vicente, 2005; Figure 2) and the abundance of pyroclastic deposits within Los Molles Formation stratigraphy (Zöllner and Amos, 1973;Rosenfeld and Volldaeimer, 1980;Gulisano and Gutiérrez Pleimling, 1995;De La Cruz and Suarez, 1997;Llambías and Leanza, 2005). The oyster-bearing conglomerate megaclast and well-rounded volcanic epiclast within the matrix of the debrite reflect long-lived reworking in shallow-marine settings prior to the mass failure, suggesting a shallow-water origin or remobilization of older slope strata, including shallow-marine deposits. This could represent the downslope transfer of sand following the collapse of reworked volcaniclastic deposits along the magmatic arc ( Figure 9). The evolution from the initial mass-wasting sediment supply responsible for erratically distributed foundered sandstone bodies (Unit 3) to a more mature system with the subtle distribution and diversity of lobe architectural elements (Unit 4) reflects the evolution to a more organised sediment supply system ( Figure 10). This was abruptly superseded by a return to sand-starved conditions dominated by dilute mud-rich flows and hemipelagic deposition until the end of the Lower Callovian (Gulisano and Gutiérrez Pleimling, 1995).

CONCLUSION
We document the anatomy and architecture of a >9.6 km long exhumed debrite and show how its twofold short-and longwavelength relief and composition provided a likely input route for the subsequent sand-rich deep-water system and influenced flow behaviour depositional patterns. The basal surface of the debrite formed ramps and steps, indicating deep incision and entrainment of the substrate that included megaclasts. The foundering of overlying sands, their resultant geometry and spatial distribution, and the down-dip increase in mud content, indicate a dynamic and rugose upper surface to the debrite and complex flow-deposit interactions. The spatial distribution of the foundered sandstones indicates groundhugging flows and the existence of debrite relief, which was progressively but not entirely healed by the submarine lobe. However, the architecture and facies distribution of the submarine lobe and their parental flows were still impacted by the long-lived, possibly rejuvenated, debrite-related topography.
The debrite emplacement coincided with an abrupt change in sediment supply to the Chacay Melehue depocentre from longterm mud-rich sedimentation to a transient sand-rich system. This change in depositional character is interpreted to have resulted from the funnelling of sediment stored in shallow marine environments to the west through a slide scar created by the debris flow, thus reconfiguring the sediment delivery pathway. Therefore, this study highlights that basin margin mass failures and their deposits play a key role in sediment dispersal patterns into deep-water settings, as well as the behaviour of subsequent sediment gravity flows travelling across their upper surface. Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 757617