Tectonostratigraphic analysis of the Late Triassic-Early Jurassic syn-rift sequence of the Neuquén Basin in the Sañicó depocentre , Neuquén Province , Argentina

The initial syn-rift infill (Upper Triassic-Lower Jurassic) of the Neuquén Basin involved a complex stacking of lava, pyroclastic and sedimentary units, in which volcanic rocks compose up to 72% of the infill. The high-quality outcrops of the southern Sañicó depocentre, which was uplifted and exposed during the Andean orogeny, were studied in order to determine the lithofacies arrangements, the tectonostratigraphic framework and the main controls over the syn-rift succession. The tectonostratigraphic study made it possible to determine three sections characterized by the predominance of effusive, pyroclastic or sedimentary deposits. These sections indicate three different evolutionary stages of the volcanic rifting. The onset of this rift stage was dominated by effusive rocks with minor ignimbrite and alluvial volcaniclastic units, which indicates the occurrence of composite volcanoes at the boundaries of the depocentre. The mid-rift stage was characterized by the development of fault-controlled, large ignimbrite deposits intebedded with volcaniclastic alluvial deposits and lava units, which suggests a volcanotectonic depression related to asymmetrical and incremental subsidence. The end of the volcanic rifting event was marked by a decrease in the volume of the volcanic rocks and the occurrence of low-relief effusive aggradational volcanic edifices and lava fields in close relation with alluvial, deltaic and lacustrine units. The analysis of the synrift sequence in the Sañicó depocentre showed the complex relationships between tectonic and volcanic controls which occurred in different magnitudes and times during the rifting evolution.


Introduction
The Neuquén Basin is located in the Andean margin between 32 and 40°S latitude (Fig. 1). It constitutes one of the most distinctive features which originated during the Upper Triassic-Lower Jurassic on the western edge of Gondwana by continental extension. Syn-rift depocentres were developed during the initial stages of the basin as isolated troughs following a multi-directional pattern, linked with profuse magmatic activity (Legarreta and Uliana, 1996;Franzese and Spalletti, 2001;Howell et al., 2005). The initial depocentres, characterized as a series of long, narrow half-grabens, were filled by a complex array of volcanic and volcaniclastic materials which are known as the Precuyano Cycle (Gulisano et al., 1984). In particular, it is important to highlight the huge volumes of lava flows and pyroclastic materials that constitute the syn-rift sequence and their strong imprint on the sedimentary systems. Although there have been a number of contributions on the evolution and fill of the basin (Howell et al., 2005, and papers therein), studies of the early depocentres of the Neuquén Basin have concentrated on either stratigraphic analyses (Gulisano, 1981;Gulisano and Pando, 1981;Gulisano et al., 1984;Riccardi and Gulisano, 1990;Leanza, 1990;Gulisano and Gutiérrez Plieming, 1994;Llambías et al., 2007;Schiuma and Llambías, 2008) or on tectonostratigraphic analyses, instead of focusing on the occurrence of extension and volcanism, and the relationships between the latter and its imprint over sedimentary environments (Álvarez and Ramos, 1999;Franzese et al., 2006;Franzese et al., 2007;Muravchik et al., 2008;Giambiagi et al., 2008a;Giambiagi et al., 2008b;Giambiagi et al., 2008c;Pángaro et al., 2009). The southern part of the Neuquén Basin is characterized by reactivated normal faults and reverse faults related to the evolution of the Andean chain that led to wellexposed Mesozoic successions of the initial stage of the basin. The aims of this study are to carry out a tectonostratigraphic analysis of the volcanic syn-rift deposits cropping out in the Sañicó depocentre -to the south of the Neuquén Province, Argentina-focusing both on the depositional processes and on the relationship between tectonic and volcanic controls in the development of the sedimentary sequence, and to devise an evolution scheme of the depocentre. The study of the syn-rift succession was carried out through detailed geological mapping and measuring of stratigraphic and sedimentary logs. Petrographic analyses were made to establish the compositional characteristics of lava and pyroclastic rocks and the nature of the provenance of the sedimentary facies.

Geological Setting
The Neuquén Basin is located on the eastern side of the Andes in central Argentina, and it constitutes one of the most important hydrocarbon basins in South America. It is limited by a volcanic arc to the west and by two mesocratonic areas, the Sierra Pintada System and the North Patagonian Massif, to the north-east and south-east respectively (Fig. 1). The basin has a multiphase tectonic history involving an early stage of extension (Upper Triassic-Lower Jurassic) followed by a post-rift phase superimposed with the development of the magmatic arc (Ramos and Folguera, 2005;Ramos, 2009). Subsequent stages of inversion and extension, which occurred during Mesozoic and Cenozoic times, relate to the evolution of the Andean orogeny (Vergani et al., 1995;Ramos and Folguera, 2005;Howell et al., 2005;Cristallini et al., 2006;Mpodozis and Ramos, 2008). The basin recorded at least two episodes of rift during its origin. In the first syn-rift stage (Rhaetian-Pliensbachian) the initial configuration of the basin was characterized by the development of isolated deep depressions bounded by normal faults and filled with volcano-sedimentary successions (Franzese and Spalletti, 2001). Subsurface studies show them as a series of half-grabens with variable polarity, intersected by en-echelon transfer faults (Vergani et al., 1995;Cristallini et al., 2006). The stratigraphy of the initial syn-rift of the basin is complex, involving a continental volcanosedimentary succession known as the Precuyano Cycle (Gulisano et al., 1984;Gulisano, 1993), which is constituted by different lithostratigraphic units (Franzese and Spalletti, 2001). The second stage of rifting (Pliensbachian-Toarcian) involved mainly marine deposits known as the Cuyano Cycle, which overlapped with the onset of the transition to the post-rift stage (Vergani et al., 1995).
During the Andean orogenic cycle, different deformation periods caused the inversion and uplifting of the basin, initially during the Late Jurassic and subsequently since the Late Cretaceous Ramos and Folguera, 2005;Morabito et al., 2011). The southern part of the basin records the typical areas which have usually been taken as a case study for the analysis of the initial syn-rift of the basin (Gulisano and Pando, 1981;Leanza, 1990;DʼElia and Franzese, 2005;Muravchik and Franzese, 2005;Franzese et al., 2006;Franzese et al., 2007;DʼElia, 2008;Muravchik et al., 2008;Bilmes et al., 2008;Spalletti et al., 2010;DʼElia, 2010;Muravchik and D ̕ Elia, 2010;Muravchik et al., 2011;D'Elia et al., 2012). The study area is located on the southern boundary of the Neuquén Basin, in the limit between the Neuquén and Río Negro provinces, Patagonia, Argentina (Fig. 1). In this area, the initial depocentres have recently been redefined, among them the Sañicó depocentre D'Elia, 2010;D'Elia et al., 2012). Reactivated normal faults and reverse faults led to well-developed successions of the basin ( Fig. 2A, C), which were correlated by different authors with pre-rift, syn-rift and transitional syn-rift to post-rift units (Fig. 2B;Franzese and Spalletti, 2001;Muravchik et al., 2008). Locally, Cenozoic volcanic and sedimentary rocks overlie the initial fill of the basin.
The syn-rift fill lies unconformably over the pre-rift units. It is constituted by the Sañicó For- Argentina, and in the Coastal Cordillera, Chile. Location of the study area is highlighted in a black square. Modified from Franzese and Spalletti (2001). mation (Fig. 2C;Stipanicic, 1967;Stipanicic et al., 1968). It is important to point out that some authors include the siliciclastic deposits of the Upper Triassic Paso Flores Formation (Frenguelli, 1948;Morel and Ganuza, 2002), that crops out in the vicinities of the study area, to the Precuyano Cycle (Franzese and Spalletti, 2001). Nevertheless, most of the Paso Flores Formation lies outside of the Sañicó Depocentre. The significant hiatus between this unit and the Sañicó Formation and the strong difference in lateral and vertical distribution of the two units suggest that the Paso Flores Formation should not be incorporated to the syn-rift evolu-tion of the Sañicó Depocentre (Fig. 2). Therefore, in this paper we will concentrate in the analysis of the Sañicó Formation. The Sañicó Formation is composed by lava, pyroclastic and siliciclastic deposits (Lambert and Galli, 1950;Galli, 1953;Galli, 1969) constrained to a Hettangian-Sinemurian age on the basis of stratigraphic analysis (Stipanicic, 1967;Stipanicic et al., 1968). Recent data taken from the underlying Piedra del Águila Formation close to the study area, constrained this unit to the Sinemurian (Spalletti et al., 2010). The transitional syn-rift to post-rift deposits unconformably overlies both basement and syn-rift units ( Fig. 2A) across a major discontinuity that is related to the widespread marine transgression of the Neuquén Basin (Gulisano et al., 1984). These deposits were grouped into the Piedra Pintada Formation, which is composed of coeval siliciclastic and carbonate deposits ( Fig. 2B; Herbst, 1966;Stipanicic et al., 1968;Stipanicic, 1969) of Pliensbachian age (Damborenea et al., 1975;Damborenea and Manceñido, 1993), which were interpreted as shallow-marine siliciclastic and carbonate deposits that change laterally and upward into off-shore deposits (Gulisano and Pando, 1981;Gulisano et al., 1984).

The Sañicó depocentre
The tectonostratigraphic framework of the Sañicó depocentre was established on the basis of structural and geological mapping, as well as on the sections logged (Figs. 2A; 7A and 8B). It was mainly filled by a volcano-sedimentary continental succession (the Sañicó Formation), and secondarily by marine deposits of the Cuyano Cycle ( Fig.  2A, C). The faulted margins of the Sañicó depocentre are constituted at present by high degree to partially inverted fault systems (i.e., the Cerro Michal and the Cerro Corona fault systems. Both boundaries of the depocentre are shown as thickskinned emergent anticlines forming NE-oriented folds with minor ENE-and NNW-trending structures ( Fig. 2A, C). Within the depocentre, minor partially reactivated normal faults parallel and oblique to the margin structures were registered ( Fig. 2A). Several features may be considered as indicators of fault reactivation and inversion which are mainly related to abrupt changes in topography and stratigraphy associated with the major tectonic boundaries. The Cerro Michal fault system ( Fig.  2A, C) separates the syn-rift succession in a thin sequence of lava flows to the NW of the faults from a thick syn-rift succession with more than a thousand metres of lava, pyroclastic and sedimentary rocks to the SE (Fig. 7A). To the SE of the Cerro Michal fault system, the syn-rift fill of the depocentre constitutes a large syncline which can be traced for 25 km to the Cerro Corona fault system, where the syn-rift succession is absent or reaches its minimum thickness (Figs. 2A,2C and 7A). This suggests that the Cerro Michal and Cerro Corona fault systems acted as the boundaries of the depocentre. Reverse faults located on the western side of the study area caused the regional tilting of the syn-rift sequence to the east. Perpendicular to the bounding fault systems the syn-rift fill may be followed for approximately 5 km ( Fig. 2A). Beyond this point it is unconformably covered by deposits corresponding to the Cuyano Cycle. Although the Sañicó depocentre was affected by the Andean deformation, the geometries and dimensions recognized herein are comparable with those of the syn-rift depocentres documented in subsurface (Vergani et al., 1995;Cristallini et al., 2006;Pángaro et al., 2009).

Stratigraphy of the Sañicó depocentre
According to the different rock-forming and fragment-forming processes, three major groups of facies were distinguished: lava, pyroclastic and sedimentary. The first two are associated with primary volcanic processes, while the latter is related to secondary surface processes (i.e., resedimentation, reworking, weathering, mass wasting and erosion). In order to avoid the possible descriptive and genetic imbalances that frequently exist between volcanorelated, pyroclastic-related and sedimentary-related facies, in this work a non-genetic descriptive scheme is adopted to describe the facies (McPhie et al., 1993;Branney and Kokelaar, 2002;Walker, 2006, respectively). Subsequently, each facies was interpreted in terms of its rheological and hydrodynamic conditions (Fig. 3).

Lava lithofacies
They constitute up to 38% of the total thickness of the volcanic syn-rift succession. Five lithofacies have been defined according to their composition, on the basis of phenocryst assemblages, and their internal structures (Fig. 3). The coherent andesite lithofacies (cA; Fig. 4A The lava facies of the syn-rift sequence reveal a multi-compositional trend with predominance of andesite lithofacies. The range of textural arrangements and differences in phenocryst abundance indicate FIG. 3. Facies description and interpretation. The codes for lava facies were created for this work. The codes for pyroclastic facies (PDC=Pyroclastic Density Currents) were taken from Branney and Kokelaar (2002), while the codes for sedimentary facies were modified from Miall (1978) and Smith (1986). polygenetic textures, suggesting a two-stage cooling history (i.e., in the magma chamber or near surface and surface conditions) and complex feeder systems (Wilson, 1989;Marsh, 2000;Best and Christiansen, 2001). The holocrystaline to porphyritic glassy texture indicated by spherulites (high temperature devitrification of glass) suggests that liquids cooled below their normal point of crystallization or supercooled (McArthur et al., 1998;Best and Christiansen, 2001). Most of the structural patterns imply that lava facies were deposited as subaerial or free-water shallow intrusion, although in some cases hyaloclastic structures suggest damp conditions (McPhie et al., 1993). The coherent facies involve areas of slow cooling or gentle shear rate (Kilburn, 2000;Stewart and McPhie, 2003). By contrast, the occurrence of jigsaw breccias implies that the brecciation was developed by brittle fractures which occurred in response to contraction during the cooling, whereas the stratified autobreccias with rotated clasts, suggest that the main fragmentation mechanism was brittle fracture in response to transport (McPhie et al., 1993;Stewart and McPhie, 2003;Stewart and McPhie, 2006). On the basis of textural, structural and geometrical parameters different volcanic bodies may be defined, such as andesitic lava flows, lava domes and dikes composed of coherent or both coherent and brecciated lava facies (Figs. 3 and 4A-D); dacitic and rhyodacitic necks and lava domes generally constituted by coherent facies (Figs. 3 and 4E-G); and rhoylitic lava flows formed by coherent or autobrecciated facies (Cas and Wright, 1987;McPhie et al., 1993;Németh and Ulrike, 2007).

Pyroclastic lithofacies
They are the second most important lithology and constitute up to 34% of the volcanic syn-rift succession (Figs. 7A and 8A). Seven lithofacies have been defined (Fig. 3) according to combinations of texture, structure and composition determined on the basis of cognate crystal fragments. Except the massive lithic breccia lithofacies, which presents an andesitic composition, the rest of the lithofacies have a dacitic to rhyodacitic composition with minor amounts of rhyolites. Evidence of high temperature welding (i.e., welded shards, matrix spherulites, rheomorphic deformation) is generally absent. In most of the cases the pyroclastic facies show a low-grade of welding and diagenetic alteration, vapour phase crystallization or silicification. The pyroclastic facies of the volcanic syn-rift succession show different eruptive mechanisms, as well as different processes of transport and deposition. The occurrence of non-to poorly vesiculated porphyritic glassy juvenile clasts with an andesitic composition in the massive breccia facies (mBr; Fig. 5A) suggests that this was originated by gravitational collapse of lava flows or lava domes, causing short-lived pyroclastic flows, which formed block and ash flow deposits (Cas and Wright, 1987;McPhie et al., 1993). Pumiceous lapilli and cuspate shards in acid pyroclastic deposits indicate that they were caused by a magmatic eruption with an effective degree of fragmentation (Cas and Wright, 1987). Most of the acid pyroclastic facies (mlBr, mLT, dsLT, //sT and xsT; Fig. 3) were formed by progressive aggradation near the lower depositional flow boundary of a pyroclastic density current (PDC) (Branney and Kokelaar, 2002), whereas the pmL and sT facies were deposited by fallout from eruption plumes, forming fall layers. The massive and diffusestratified lithic breccia facies (mlBr-dslBr; Fig. 5B) may be interpreted as a coarse facies of ignimbrite, deposited through the lower flow boundary of a PDC in which turbulence and traction were suppressed (Branney and Kokelaar, 2002). The absence of tractional structures or diffusive stratification with linear fabric indicates rapid progressive aggradation from the fluid escape-dominated to grain flow-dominated depositional flow boundary of a PDC (Branney and Kokelaar, 2002;Kokelaar et al., 2007). The massive lapilli-tuff facies (mLT; Fig. 5C, D) indicates rapid aggradation from the high-concentration fluid escape-dominated depositional flow boundary of a PDC, whereas the massive lapilli-tuff facies with a directional fabric (mLTf) records deposition from a flow-boundary zone dominated by fluid escape processes, but with a component of granular flow (Branney and Kokelaar, 2002). In this case study the massive lapilli-tuff facies with eutaxitic fabric (emLT) commonly represents cold-burial compaction associated with diagenesis (McPhie et al., 1993;Gifkins et al., 2005) and minority welding deformation during post-emplacement stages (McArthur et al., 1998). The diffuse-stratified lapilli-tuff facies (dsLT; Fig. 5E) is deposited when flow-boundary zone conditions are intermediate between a fluid escape-dominated and a traction-dominated flow-boundary zone, causing the granular flow-dominated boundary of a PDC (Branney and Kokelaar, 2002;Kokelaar et al., 2007). Stratified and cross bedded tuff facies (//sT, xsT; Fig. 5) indicate deposition from the traction-dominated flow boundaries of a PDC (Branney and Kokelaar, 2002). The massive pumice lapillistone facies (pmL; Fig. 5G) is interpreted as proximal deposits by fallout from eruption plumes. The presence of pumice fall layers within a sheet of mLT records temporary cessations of flow (i.e., Plinian or sub-Plinian pumice fall units) (Cas and Wright, 1987). The thin stratified, laminated or masive tuff facies (sT) includes distal fall deposits from eruption plumes or from ash-rich pyroclastic density currents (Cas and Wright, 1987;Branney and Kokelaar, 2002).

Sedimentary lithofacies
The sedimentary facies constitute up to 28% of the volcanic syn-rift succession and are composed by siliciclastic and carbonate rocks. Clastic sedimentary facies (volcano-related and non-volcanorelated) were considered as epiclastic facies (sensu Cas and Wright, 1987). For the purposes of this contribution, the use of the term 'volcaniclastic' is restricted to the original definition (Fisher, 1961 in Fisher andSchmincke, 1984): deposits composed predominantly of volcanic particles. To define the relationship between sedimentary environments and the volcanic setting, a provenance analysis was carried out, determining the monomictic or polymictic composition and clast-forming processes of the deposits (i.e., monomictic ash-rich, monomictic lava-rich, polymictic volcaniclastic or polymictic epiclastic facies). According to these concepts, eleven siliciclastic lithofacies and two carbonate facies have been defined (Fig. 3). Siliciclastic lithofacies are the most common sedimentary facies within the syn-rift successions (Figs. 7A and 8A). Most of the siliciclastic facies within the studied succession show an absence of tractional structures. Therefore, this feature reveals the prevalence of transport and deposition from sediment gravity flow processes occurring in the syn-rift sequence (Fig. 3). The massive matrix-supported to clast-supported breccia facies (Brmm), with large megablocks within an unsorted monomictic andesitic matrix with different degrees of fragmentation (Fig. 6A, C), has characteristic features of volcanic avalanche debris flow deposits (Schneider and Fisher, 1998;Clavero et al., 2002;Bret et al., 2003;Bernard et al., 2009). The massive matrix-supported conglomerate facies (Gmm) can be differentiated from avalanche deposits by their better sorting, their lobate or tabular pattern and the significantly smaller height and width of their outcrops (Fig. 6B). These characteristics suggest debris flow deposits (Coussot and Meunier, 1996;Bernard et al., 2009). The massive clast-supported breccia facies (Brm; Fig. 6G) and the massive to horizontally stratified clast-supported conglomerate facies (Gm; Fig. 6D), which show a polymodal matrix and an absence of tractional structures, as well as the horizontally stratified pebbly sandstone/ sandstone facies (SGh/Sh(a); Fig. 6F) and the massive sandstone facies (Sm), which present massive or horizontal stratification with major transitional boundaries, suggest the accumulation of gravelly and sandy hyperconcentrated sheet flow deposits as a result of rapid deposition from sediment-laden currents (Smith, 1987;Smith and Lowe, 1991). The large-scale cross stratified clast-supported conglomerates facies (Gp(b)) is composed at the same conglomerates facies that constituted the Gm lithofacies, developed in tabular bodies, with irregular bases which show large-scale foresets of 2 m to 3 m thick (Fig. 6E). The medium-scale cross-stratified clast-supported conglomerate facies (Gp(a)), and the cross-stratified sandstone facies (Sx) form lenticular bodies with erosive, concave bases. These bodies represent channel-fill deposits (Miall, 2006), associated with bedload transport and deposition from stream flows. The horizontally laminated sandstone facies (Sh(b)) with tabular bodies suggest sheet flood deposits in upper-flow regime conditions (Miall, 2006). The horizontally laminated siltstone facies (Fl; Fig. 6H) may represent settling from suspensions in low-energy water bodies (Miall, 2006). On the basis of the provenance analysis, the facies discussed above showed different compositional features. The Brmm, Brm and Sm facies present only a monomictic volcaniclastic composition (Brmm, monomictic andesitic composition; Brm and Sm, monomictic pyroclastic composition, i.e., pumiceous breccias and tuffaceous sandstones). The Gmm, Gm, Gp(a), SGh/Sh(a) and Sx facies present three different compositional types: monomictic andesitic composition, polymictic volcaniclastic composition or polymictic epiclastic composition (volcanic and basement lithoclasts), while the Gp(b) and Sh facies are constituted by polymictic epiclastic conglomerates.
The compositional features of the epiclastic facies reveal different degrees of relationship between the sediments fed to the sedimentary systems and the characteristics of the volcanic landscapes. The monomictic compositions generally suggest both a rapid resedimentation in syn-eruptive (effusive or explosive) stages (McPhie et al., 1993) and poorly integrated drainage networks, whereas a polymictic volcaniclastic composition represents a reworking and deposition during either minor syn-eruptive or inter-eruptive periods with poorly integrated drainage networks. The polymictic epiclastic compositions are constituted by an admixture of volcanic and igneousmetamorphic fragments originated by the erosion of basement rocks as well as the ones that constitute the syn-rift succession. This feature indicates integrated drainage networks and a decreasing volcanic supply to the depocentre.
The carbonate lithofacies (Fig. 3) are more restricted and localized than the other lithofacies, being developed towards the top of the syn-rift sequence (Fig. 7A). The carbonates are always extremely silicified and show a very low preservation of their primary sedimentary features, which generally occur interbedded with the Fl lithofacies. The horizontally laminated limestone facies (Ll(a); Fig. 6I) and the microbial laminated limestone facies (Ll(b)) may be interpreted as the non-marine precipitation of carbonates in shallow water bodies (Riding, 2000). The microbial laminated limestone facies (Ll(b); Fig. 6H, I) represents stromatolitic deposits in lowenergy condition with little or no sediment supply (Talbot and Allen, 2002).

Spatial distribution of the lithofacies within the syn-rift sequence in the Sañicó depocentre
The stratigraphic and facies analyses reveal the large volume of effusive and explosive volcanic products within the Sañicó Formation, which constitute up to 72% of the syn-rift succession. Furthermore, different mechanisms of transport and deposition, as well as complex interplays between volcanic materials and sedimentary facies, were established (Fig. 3). Systematic vertical and lateral variations in the distribution of the volcanic and sedimentary lithofacies (Figs. 7A and 8A) indicate changes in the tectonic and volcanic controls during the syn-rift evolution. The tectonostratigraphic evolution can be recognized on the basis of the prevailing volcanic rocks in close relation with the four different compositional types of epiclastic facies documented throughout the syn-rift sequence. As a result, three informal sections can be devised: a lower section, middle section and upper section (Figs. 7A, B and 9). Volcanic rocks are volumetrically the most important in the lower and middle sections, while they are less frequent in the upper section of the syn-rift succession (Fig. 7A).
The lower section is mainly characterized by andesitic lava rocks and andesitic-dominated volcaniclastic rocks (monomictic andesitic composition and polymictic volcaniclastic composition). To a lesser extent, rhyolitic/rhyodacitic lavas and pyroclastic current deposits occur ( Fig. 7A and 8A). This section was deposited in a highly compartmentalized depocentre, as depicted by large thickness variation and a high degree of preservation both beyond the boundaries of the depocentre and within it (Fig. 7A). Transverse to the faulted margins of the depocentre in the footwalls of the Cerro Michal and Cerro Corona fault systems, this section is composed of a thin succession of coherent and autobrecciated andesitic lava flows and lava domes (cA+aA facies) and dike swarms (cA) oriented parallel to the boundary fault systems. Towards the hanging walls, the section becomes thicker, reaching its maximum thickness close to the Cerro Michal fault system (1,100 m). In the hanging walls, the sections are constituted by large successions of coherent and autobrecciated andesitic lava flows (cA+aA facies) that prevail towards the base, and pyroclastic current deposits and gravity flow-dominated deposits towards the top (Fig. 7A). In some cases, andesitic lava domes and dikes (cA±aA facies), rhyodacitic lava domes (cRD facies) and dacitic necks (cD facies) associated with interior faults occur. The pyroclastic current deposits (mLT+dsLT±mlBr±xT facies) appear as huge lenticular bodies with their axes perpendicular to the boundary fault systems (Fig. 8A, B). The sedimentary lithofacies appear associated in different types of units. Volcanic debris avalanche deposits and debris flow deposits (Brmm+Gmm) are the dominating sedimentary lithofacies which occur as successions several decametres thick and hundreds of metres wide (Fig. 8A, B). To a lesser extent, major lenticular bodies several metres thick and hundreds of metres wide-mainly dominated by sandy and gravelly hyperconcentrated flow deposits (Gm+SGh/Sh(a))-were recorded (Fig. 8A). Inside them, tractional facies such as Gp(a) and Sx and debris flow deposits (Gmm) may occasionally appear. It is important to highlight that parallel to the boundary fault systems this section shows a change in thickness, reaching their minimum towards the tip point of the fault, in which the different facies onlap the basement rocks (Fig. 8A, B).
The middle section are mostly composed of thick rhyolitic/rhyodacite pyroclastic current deposits which appear intebedded with pyroclastic-monomictic volcaniclastic deposits and andesitic-rhyolitic lava flows and lava domes related to extensional faults (Figs. 7A and 9). The architecture of this section is very different from the lower section, prevailing only inside the depocentre (Fig. 9) and a more regular geometry and lateral continuity of the lithofacies are common. Pyroclastic current deposits crop out as laterally continuous units up to 350 m thick, constituting more than 85% of the volume of the middle section (Fig. 7A). At least three units of pyroclastic current deposits were recognized and they appear intercalated with either thin successions (<30 m thick) of pumice/tuff-rich hyperconcentrated flow deposits (Brm+Sm) or andesitic lava flows/lava domes (cA+aA) related to fractures (Fig. 7A). Pyroclastic units are lapilli-tuff-dominated (i.e., mLT+dsLT±// sT±xsT facies) and in some cases show pmL facies at the base and mlBr facies at the base or in the middle part (Fig. 7A). Each of the pyroclastic units show large change in thickness related to both internal and boundary faults of the depocentre (Fig. 7A, B) during the eruptive events. These structural features not only controlled the accommodation for these units, but also the distribution of internal lithofacies, causing the deposition of mlBr facies close to the boundary fault systems of the depocentre (Fig. 7A). Towards the top of this section, scarce rhyolitic lava flows (cR±aR) related to fractures were recorded (Fig. 7A).
The upper section of the syn-rift sequence is characterized by a decrease in the volume of the volcanic rocks (<55%), as well as a reduction of the accommodation space, presenting a more limited distribution and less thickness in the depocentre (Figs. 7A, B and 9). It is composed of andesitic lava flows with epiclastic (i.e., polymictic volcaniclastic and polymictic epiclastic facies) and carbonate sedimentary rocks. These lithofacies are asymmetrically distributed both laterally and vertically within the depocentre (Fig. 7A, B). Close to the hanging wall of the Cerro Corona boundary this section is constituted by coherent, and to a lesser extent autobrecciated, andesitic lava flow (cA±aA) successions up to 250 m thick. On the other side of the depocentre, near the Cerro Michal fault system hanging wall, this section is constituted by coherent andesitic lava flows (cA) intercalated with epiclastic and carbonate deposits (Fig. 7A). Towards the base, the epiclastic facies composed of debris flow deposits with a polymictic volcaniclastic provenance have a large distribution. Towards the top, polymictic epiclastic and carbonate facies interbebbed with thin coherent andesitic lava flows (cA) occur. Close to the Cerro Michal boundary, polymictic epiclastic rocks are developed as tabular bodies composed of large-scale crossstratified conglomerates (Gp(b)) with Gm and Sh(b) lithofacies (Fig. 7A). These bodies form sequences up to 5 m thick with coarsening-upward grain-size trends that can be related to gradual changes in the depositional processes. Towards the inner depocentre, these facies turn into a thin succession of limestones ((Ll(a)+Ll(b)) and fine-grained siliciclastic deposits (Ll) with minor Sx lithofacies.

Discussion: evolution and tectonostratigraphic controls over the volcanic syn-rift sequence
The fault systems of the Sañicó depocentre controlled the accumulation of a thick syn-rift succession, in the hanging wall of the main boundary fault of the depocentre (Cerro Michal fault system; Figs. 7A, B and 9). Although the stratigraphic analysis constrained the syn-rift to the Hettangian-Sinemurian (i.e., Sañicó Formation; Stipanicic, 1967;Stipanicic et al., 1968), a recent geochronological study developed close to the studied area, constrained this volcanic syn-rift sequence to the Sinemurian age (i.e., 191,7±2,8 Ma to Pliensbachian age; Spalletti et al., 2010), representing a total period of approximately 4 Ma.
Stratigraphic and facies analysis determined the main component of the syn-rift sequence with their hydrodynamic and rheological mechanisms of transport and deposition (Fig. 3). In order to establish the tectonostratigraphic framework of the Sañicó depocentre, three sections were defined on the basis of the prevailing volcanic rocks in close relation with the four different compositional types of epiclastic facies documented (Figs. 7A and 9). Thus, the relationships and architectural patterns between the lava, pyroclastic and sedimentary components that constitute the syn-rift sequence were determined (Figs. 7B and 8B). The change in the amount of volcanic versus sedimentary materials, the variations in the composition and mechanisms of transport and deposition of the sedimentary facies, and different tectonic and volcanic controls over the fill are the key elements that marked the signature of each section. Therefore, the tectonostratigraphic evolution of the syn-rift succession in the Sañicó depocentre may be carried out by means of the analysis of each section. These sections in turn represent three major evolutionary stages of the volcanic syn-rift sequence: the early-rift stage, mid-rift stage and late-rift stage.
The early-rift stage (lower section) corresponds to the onset of the syn-tectonic infill of the Sañicó depocentre, which is characterized by a volcano-sedimentary wedge dominated by effusive volcanic products and volcaniclastic rocks (Figs. 7A, 8A and 9) caused by the resedimentation and reworking of the former. During this stage the syn-rift units are wedge-shaped and highly compartmentalized with large thickness variations indicating a half-graben geometry (Fig. 7). This condition is common in the early evolution of the rift basin (Gawthorpe and Leeder, 2000;Morley, 1999a). Volcanic rocks compose up to 61% of the section, showing a discernible polarity with the sedimentary deposits (Fig. 7A). In the footwalls of the boundary fault systems of the depocentre, this section is composed of a thin succession of lava flows, lava domes and dike swarms (cA+aA facies), whereas towards the hanging walls it is formed by a large succession of lava flows (cA+aA facies), pyroclastic current deposits and gravity flow-dominated deposits (Fig. 7A). To a lesser extent, dikes, lava domes and necks of diverse composition are present. The huge lenticular bodies of pyroclastic current deposits (mLT+dsLT±mlBr±xT facies) represent the filling of paleovalleys (Fig. 8A) from short-lived granular fluid-based pyroclastic density currents caused by simple-pulse collapse from pyroclastic column (Branney and Kokelaar,FIG. 9. Lateral and vertical distribution of sections defined for the syn-rift sequence in the Sañicó depocentre (see figs. 3 and 7 for the location and details of the logged sections).
2002). The sedimentary lithofacies with andesitic dominated provenance, that is, the volcanic debris avalanche deposits and debris flow deposits (i.e., Brmm+Gmm) and the major lenticular bodies dominated by hyperconcentrated flow deposits (i.e., Gm+SGh/Sh(a)±Gp(a)±Sx±Gmm), suggest shortand long-lived degradational stages of the volcanic landscapes (Davidson and De Silva, 2000) due to alluvial processes (sensu Blair and McPherson, 1994) in either syn-eruption or inter-eruption periods with less integrated drainage networks (Palmer et al., 1993;Zanchetta et al., 2004). Most of the features of the volcanic and sedimentary units defined herein were described for many composite volcanoes in the world, such as the compositional range from andesitic to rhyodacitic products dominated by thick andesitic lava successions (Cas and Wright, 1987), the pyroclastic current deposits that filled valleys (Davidson and De Silva, 2000), major degradational events composed of volcanic debris avalanches (Schneider and Fisher, 1998;Reubi and Hernández, 2000;Belousov et al., 1999;Clavero et al., 2002;Shea et al., 2008), large channelized volcaniclastic alluvial units (Palmer et al., 1993;Zanchetta et al., 2004) and minor acid lava domes and necks (cRD, cD facies) that worked as satellite vents (Davidson and De Silva, 2000). The characteristics of this units and their polarity within the depocentre suggest that during the early-rift stage composite volcanoes occurred on the boundary of the Sañicó half-graben (Fig. 10). The evolution of the syn-rift sequence was followed by the mid-rift stage (middle section), which is pyroclastic-dominated (Fig. 7A). The syn-rift deposits of this stage are recorded only inside the depocentre and, although it presents thickness changes related to internal faults, more regular geometries are common (Figs. 3, 7 and 9). The middle section is mainly composed of thick acid pyroclastic current deposits intebedded with thin pumiceous and ash-rich volcaniclastic deposits and andesiticrhyolitic lava units related with extensional faults. The internal and boundary faults of the depocentre (Figs. 3 and 7) controlled the thickness changes and facies distribution of the pyroclastic units during the eruptive events. This reveals that previous faults worked as volcano-tectonic extensional faults during the eruptive events. At least three volcano-tectonic periods were recorded on the basis of the occurrence of large-pyroclastic current deposits that choked the accommodation space of the depocentre. The distribution and internal facies arrangement of the lapilli-tuff-dominated pyroclastic current deposits (i.e., mLT+dsLT±//sT±xsT±pmL±mlBr facies) suggest that they originated as a result of long-lived, sustained, quasi-steady, granular fluid-based pyroclastic density currents, due to complex-pulse collapse from pyroclastic column Branney and Kokelaar, 2002). The ignimbrite deposits are interbedded with andesitic lava flows/lava domes (cA+aA) mainly related to fractures and thin successions of pumiceous and ash-rich hyperconcentrated flow deposits (Brm+Sm). This assemblage may indicative of rapid resedimentation from alluvial processes, common in syn-eruptive stages and poorly integrated drainage networks (Smith, 1986;Smith, 1991;Smith and Lowe, 1991;Muravchik et al., 2011). Both units indicate that the volcano-tectonic events were not only related with asymmetric subsidence, but also with episodic subsidence. Some authors classify this kind volcanic edifice as volcano-tectonic depressions (Moore and Kokelaar, 1997;Moore and Kokelaar, 1998;Lipman, 2000). These types of calderas were recorded in different extensional or transtensional settings related to magmatism in the world (the Scafell caldera, Branney and Kokelaar, 1994; the Glencoe Graben, Moore and Kokelaar, 1997;Moore and Kokelaar, 1998; the Long Valley caldera, ; the Taupo and Okataina calderas, Spinks et al., 2005; the Sierra Madre Occidental caldera, Aguirre-Díaz et al., 2008; the Cerro Aguas Calientes caldera; Petrinovic et al., 2010). Therefore, a volcanotectonic depression related with asymmetrical and incremental subsidence was superimposed during the mid-rift stage of the Sañicó depocentre (Fig. 10).
The late-rift stage (upper section) records a decrease in the volume of volcanic rocks and a reduction of the accommodation space (Figs. 7,9 and 10). During this period, andesitic lava units with epiclastic and carbonate sedimentary rocks were asymmetrically distributed both laterally and vertically within the depocentre. Close to the hanging wall of the Cerro Corona boundary (Fig. 7), low-relief effusive aggradational volcanic edifices and lava fields may be interpreted on the basis of a coherent, and to a lesser extent autobrecciated, andesitic lava flow (cA±aA) succession (Brown and Bell, 2007). On the other side of the depocentre (Fig. 7), at the base of this section, near the Cerro Michal fault system hanging wall, coherent andesitic lava flows (cA) interbedded with epiclastic polymictic volcaniclastic debris flow deposits indicate low-relief effusive aggradational volcanic edifices and lava fields (Brown and Bell, 2007) with highly integrated drainage networks as a result of the erosion and reworking of the different volcanic rocks of the syn-rift. Towards the top of the section (Fig. 7), thin coherent andesitic lava flows with epiclastic and carbonate facies are developed. Tabular bodies characterized by large-scale crossstratified conglomerates (Gp(b)+Gm+Sh(b) facies) suggest a coarse-grained delta front and alluvial facies (Reading and Collinson, 2002;Blair and McPherson, 2008). The polymictic epiclastic compositon of these facies indicates highly integrated drainage networks and less influence of active volcanism. Towards the inner depocentre, the sedimentary environment turns into a thin succession of limestones and fine-grained siliciclastic deposits interpreted as shallow lacustrine units (i.e., Ll(a)+Ll(b)+Ll; Talbot and Allen, 2002) with underflow deposits (Sx) in a prodelta to off-shore lacustrine position (Reading and Collinson, 2002;Blair and McPherson, 2008;Rohais et al., 2008).
Most of the studies of rift-basins were made for non-volcanic extensional depocentres (e.g., Schlische, 1991;Schlische, 1992;Morley, 1995;Schlische and Anders, 1996;Gawthorpe and Leeder, 2000;Gawthorpe et al., 2003;Young et al., 2003;Morley et al., 2004;Jackson et al., 2005), which is why they do not discuss the interrelationships between the tectonic structures and volcanism, the large volume of volcanic and volcaniclastic materials provided to the depocentres, the aggradation above the geomorphologic threshold, and the lateral and vertical variations of environments in these places. The tectonostratigraphic framework given above showed how these tectonic and volcanic variables were present in different magnitudes and times during the syn-rift evolution. The tectonic controls over the syn-rift sequence are shown at different scales. The Sañicó depocentre may be classified as a volcanic half-graben with a NW-SE orientation and oblique internal tectonic features (Figs. 3 and 7). According to the pervasive and isolated fabrics suggested by some authors (Varela et al., 1991), the external and internal structural configuration of the depocentre seems to be inherited from pre-existing fabrics of the basement. This condition has been recognized for the Neuquén Basin (Franzese and Spalletti, 2001) and other Gondwanic basins (e.g., Cuyo Basin; Ramos and Kay, 1991), as well as for other extensional basins in the world (Morley, 1999b;Moustafa, 2002;Morley et al., 2004). Although recent studies have shown an almost unidirectional large-scale pattern (NNW-SSE) for the initial depocentres (Cristallini et al., 2006;Giambiagi et al., 2008a;Giambiagi et al., 2008b;Giambiagi et al., 2008c, Cristallini et al., 2009, related to NNE regional extension (Bechis et al., 2010), a similar orientation to those of the Sañicó depocentre were indicated for the early stage of the Neuquén Basin (e.g., Sierra de Chacaico; Franzese et al., 2007; Sierra de Reyes and Sierra de Cara Cura; Giambiagi et al., 2008a;Giambiagi et al., 2008b). On a smaller scale, the extensional faults not only controlled the accommodation space of the depocentre, but also provided magma pathways to the surface and a focus for dike emplacement (e.g., volcanic necks and dike swarms), which could have contributed to the rift extension (Abebe et al., 2007;Spinks et al., 2005;Cassidy et al., 2009). Thus, the tectonic control was a major feature at the location of the volcanic building (i.e., composite volcanoes on the boundary of the depocentre), taking over at the same time the polarity of the sedimentary facies in both the volcano-dominated stages of the syn-rift (e.g., lower section) and a sedimentary-dominated stage of the syn-rift (i.e., upper section).
The main volcanic controls documented in the syn-rift sequence may be summarized as the superimposition of magmatism with tectonic features and the influence of such controls over the volcanic landscapes and sedimentary environments. The superimposition of magmatism with tectonic fea-tures is indicated by the repeated volcano-tectonic episodes recorded. They suggest the recurrent plumbing of magmas along the Sañicó depocentre and the magma-chamber roof collapse causing the formation of a polygonal caldera (i.e., mid-rift stage; Fig. 10), usually identified in extensional settings of this kind . In order to establish the influence of the volcanic controls on the effect of the ancient landscapes, the interrelations between the magmatic systems with the upper part of the crust is an essential aspect to be considered (Cas and Wright, 1987;Németh and Ulrike, 2007). As examples, constructional volcanic edifices interpreted during the early-rift stage (i.e., composite volcanoes) controlled the distribution of the volcanic materials beyond the limit of the depocentre (Fig. 10), whereas an excavational volcanic building identified for the mid-rift stage (i.e., volcano-tectonic depressions) controlled the preservation of pyroclastic deposits within the depocentre (Fig. 10). Finally, the degree of volcanic activity and the types of eruption processes (i.e., effusive or explosive) are also important elements to be considered in order to determine their effects over the sedimentary paleoenvironments of the syn-rift sequence. During high volcanic supply, the overfilled condition of the depocentre caused disordered stacking of the sedimentary successions (e.g., lower and middle sections). Besides, the effusive-or explosive-dominated eruptive processes within the syn-rift sequence controlled the grain size of the volcaniclastic debris provided to the depocentre (e.g., coarse-grained during effusive periods, and tuffaceous fine-grained during the explosive stages). On the other hand, the decrease in volcanic supply identified in the upper section caused the development of well-organized stacking patterns in the sedimentary sequences (e.g., coarsening-upward succession of delta front deposits), a better arrangement of the sedimentary environment, and highly integrated drainage networks.
An integrated model may be suggested for the Sañicó depocentre which developed around 4 Ma during Sinemurian age. The onset of the rift sequence began with the development of a highly compartmentalized half-graben depocentre simultaneously with a growth of polygenetic constructional volcanic edifices in boundary position (i.e., composite volcanoes defined in the lower section; Fig. 10). Subsequently, volcano-tectonic subsidence events were superimposed, originating polygonal calderas (Fig. 10). Finally, minor volcanic edifices with an organized sedimentary environment were developed during the end of the volcanism and the extensional event (Fig. 10). It is important to highlight that on the basis of similar stratigraphic patterns recorded in other depocentres of the basin (Pángaro et al., 2002;Franzese et al., 2006;Franzese et al., 2007;Llambías et al., 2007;D'Elia et al., 2012), as well as in other volcanic extensional basins in the world (e.g., the volcanic Taupo zone; Price et al., 2005;Rowland et al., 2010), it can be stated that this model could be useful to characterize the volcanic-rift depocentres at the onset of the Neuquén Basin.
Other remarkable issue is the evolution of the depocentre in relation to the large-scale Upper Triassic-Lower Jurassic extensional event and how classical stratigraphical units of this age are grouped within the rifting event. In this case it is possible that the Upper Triassic Paso Flores Formation is related to the first-order extensional event that originated de Neuquén basin, but certainly is not related with the tectonostratigraphic evolution of the Sañicó depocentre. This indicate that, at least, there are two future topics that must pay attention: a. accurate evolution time of the isolated depocentres of the basin-in this case approximately 4 Ma-and b. the possibility of total reorganization of the accommodation space between different syn-rift event during the rifting.

c onclusions
The Sañicó depocentre is a extensional half-graben with a NW-SE orientation and oblique internal tectonic features that were active during the early rifting of the Neuquén Basin. The graben border system of the Sañicó half-graben was controlled by the Cerro Michal fault system and its associated structures. The syn-rift fill in the Sañicó area constitutes a complex arrangement of units very diverse in nature in which volcanic rocks compose up to 72% of the infill. The stratigraphic and facies analysis established the different components of the syn-rift sequence, i.e., lava, pyroclastic and sedimentary deposits. The predominance of these was identified in different parts of the fill. The main mechanisms of transport and deposition identified are lava flows, pyroclastic currents and sediment gravity flows. The tectonostratigraphic analysis revealed that the evolution of the syn-rift sequence in the Sañicó depocentre may be considered in three sections (the lower, middle, and upper sections) characterized by the predominance of effusive, pyroclastic or sedimentary rocks and distinct architectural patterns. The signature of each stage of the syn-rift evolution mainly depends on the change in the amount of volcanic versus sedimentary products in close relationship with the structural features of the depocentre. The onset of the syn-tectonic infill of the Sañicó depocentre (lower section) was characterized by andesitic-rhyodacitic lava flows with the predominance of thick andesitic lava sequences, pyroclastic current deposits that filled valleys, major degradational events composed of volcanic debris avalanches, large channelized volcaniclastic alluvial units and minor acid lava units that worked as satellite vents. The arrangement and distribution of these units may be interpreted as composite volcanoes on the boundary of the Sañicó depocentre. The evolution of the syn-rift sequence continued with the development of fault-controlled, large ignimbrite deposits with an acid composition, intebedded with volcaniclastic alluvial deposits and andesitic-rhyolitic lava units related to extensional faults (middle section), which suggests a volcanotectonic depression related with asymmetrical and incremental subsidence. The end of this syn-rift event was marked by a decrease in the volume of volcanic rocks, a reduction of the accommodation space and the development of highly integrated drainage networks. During this period (upper section), lowrelief effusive aggradational volcanic edifices and lava fields in close relation with alluvial, deltaic and lacustrine units occurred.
The analysis of the syn-rift sequence in the Sañicó depocentre determined the complex relationships between tectonic and volcanic variables, occurring in different magnitudes and times during the syn-rift evolution. The tectonic and volcanic controls over the syn-rift sequence were established at different scales, affecting the accommodation space, the magma-feeding systems, the types and location of the volcanic buildings, as well as the polarity, types and distribution of syn-rift deposits.
On the basis of facies and tectonostratigraphic analyses, a conceptual model was devised to understand the relationships between volcanism and extensional tectonics, which could be useful to characterize the volcanic-rift depocentres of the Neuquén Basin.