Hybrid coastal edges in the Neuquén Basin ( Allen Formation , Upper Cretaceous , Argentina )

The Allen Formation records the first Ingression Atlantic to the Neuquén Basin during the Late Cretaceous. The definition of lithofacies and facies associations interpretation for stratigraphic sections in Paso Córdoba and Salitral Moreno area, Río Negro, Argentina allowed to establish the depositional system that characterized this transgression in the northeastern edge of the Basin. In this paper we present sedimentological analysis of conglomeratic, sandstone, heterolithic and pelitic facies, which allowed the interpretation of tidal channels (TC), intertidal flats (ITF), stormsinfluenced tidal flat (SITF), subtidal flat (STF) and shoreface deposits (SF) parts of the depositional environment. These deposits represent a sedimentary records preserved example of hybrid systems, in which, the base of the sequence has greater tidal influence, while the upper portion is dominated by wave action. The paleocurrent data indicating a NNWSSE direction to the shore and correlations and spatial distribution of facies associations propose paleogeographic and paleoenvironmental interpretations to Malargüe Group base. Then in this paper the relationship of this coastal environment presents with wind systems previously defined in this area for the Allen Formation.


Introduction
The correct definition of the depositional environments is an essential component of basin analysis, especially in oil basins where it is of great interest to have detailed sedimentological studies that contribute to regional correlation and characterization of the various units that comprise it, in particular in sections of the stratigraphic column characterized by transitional environments where continuity is highly variable and difficult to predict. The Neuquén Basin, located in west-central Argentina (Fig. 1A, B, and C) between 34°-41°S and 66°-71°W, is the leading oil and gas producing sector in the country. The Allen Formation of the Malargüe Group (Fig. 1D) is one of the units that forms the basin's stratigraphic sequence and, recently, has provided a wide variety of fossils which have been of great interest to paleontological studies. However, the knowledge we have regarding the sedimentary paleoenvironments that characterized this lithostratigraphic unit is still limited in many areas of the basin.
The definition to date is that the Allen Formation records a continental to marine sedimentation coinciding with the first Atlantic transgression during the Late Cretaceous. The variety and arrangement of fossils found suggest some complexity in estimating the layout of the coastline, although the depositional environments recognized so far in this formation, such as estuaries and tidal flats (Andreis et al., 1974;Barrio, 1990;Armas and Sánchez, 2011), brackish lakes  and aeolian fields (Armas and Sánchez, 2013) in part allow the corresponding paleogeography to be established. This work contributes to this last aspect and presents the interpretation of the sedimentary system linked to a transitional environment with variations in the influence of tides and the action of waves.
The differentiations in coastal depositional environments caused by waves and tides have been defined since the early 70s by Galloway (1975), Wright (1977) and Hayes (1979) with their classification diagrams, among others. Dalrymple and Zaitlin (1992) suggest the definition of tidal flats as environments dominated by tides with a gentle slope and little influence of waves, unlike beach or shoreface environments, in which the slopes are greater and are subject to strong wave action. Although existing models of shoreface depositional systems rarely consider the effects of tides, this oc-curs in a wide range of coastal systems and presents great variability controlled by relative changes in sea level, the influence of storm and fair-weather waves and tides (Dashtgard et al., 2012). What is certain is that among these coastal environments there are considerable mixtures of energies which are given very little study. However, some tools such as the identification of muddy heterolithic facies, including flasser, wavy and lenticular stratification or the deposits generated by storm waves, including hummocky cross stratification are essential to distinguish the dominant power (wave or tidal, respectively) of the sedimentary paleoenvironment (Yang et al., 2005). A term which has emerged in recent times is the reference to hybrid depositional systems, which occur where the influence of waves, seas and rivers overlap (Boyd et al., 1992;Dalrymple and Zaitlin, 1992;Basilici et al., 2012).
The aim of this paper is to present detailed sedimentological study conducted for the Allen Formation by showing the types of paleoenvironments that occurred with the onset of the Atlantic ingression on the eastern edge of the basin, and thus contribute to the definition of the paleogeography of the Upper Cretaceous.
The research was carried out in Río Negro's northeastern area (Fig. 1C), by conducting a detailed survey of continuous stratigraphic sections of up to 2 km along the north side of the Salitral Moreno (pSM Fig. 2) and in the vicinity of the town of Paso Córdoba (pPC1 and pPC2 Fig. 2). In addition to preparing sedimentological profiles and photomosaics mapping of sedimentary bodies, research includes the identification and interpretation of lithofacies-essential units that allow the characterization of sedimentation paleoenvironments. The definition of these lithofacies involves combining observations made about spatial relationships and internal characteristics (lithology and sedimentary structures) with comparative data from other studies of stratigraphic units, particularly from studies of modern sedimentary environments (Walker, 1990). Furthermore, in this work, the interpreted lithofacies were grouped into facies associations, considering genetically related features and sedimentary features including geometry, continuity and shape of the lithologic units (Miall, 1984). This definition of facies associations in conjunction with the statistical analysis of paleocurrent data make it possible to create a model of deposition for the Allen Formation.

Geological environment
The Allen Formation is part of the stratigraphic sequence that fills the Neuquén Basin. From the Triassic period to the present, the active subduction and development of the magmatic arc on the western edge of Gondwana controlled the geological evolu-tion of the basin. During the first phase of syn-rift, the low rate of subduction of the proto-Pacific plate generated a regional continental intraplate extension. At this stage the basin was filled with continental and volcaniclastic deposits (Vergani et al., 1995;Howell et al., 2005). Later, in the post-rift phase, the development of an active subduction zone associated  with the evolution of a magmatic arc which generated backarc subsidence produced deposits of marine and continental sediments (Franzese et al., 2003;Howell et al., 2005). In the earliest stages of the development of the foreland basin in late Cretaceous is linked to a decrease in the angle of slab subduction produced by compression and flexural subsidence (Vergani et al., 1995;Ramos, 1999) and the uplift of the thrust belt. This stage is marked by the deposit of the Neuquén Group (Tunik et al., 2010;Aguirre Urreta et al., 2011) with a continental sedimentation and a transitional environment (Armas and Sánchez, 2011) linked to the first Atlantic transgression into the basin, of which the Malargüe Group offers the greatest evidence (Fig. 1D). The Andean tectonic finally caused the complete folding and uplift of the Mesozoic succession, exposing the wide variety of depositional environments (Franzese and Spalletti, 2001;Howell et al., 2005).

Stratigraphy of the area
The Allen Formation is the lower unit of the Malargüe Group (Fig. 1D) and is considered to be late Campanian to early Maastrichtian based on foraminifera (Ballent, 1980) and magnetostratigraphic (Dingus et al., 2000) studies. The stratotype of this formation was defined by Uliana and Dellapé (1981) in eastern area of the Bajo de Añelo, where the relation between base and top is clearly exposed. The deposits are mostly clastic interbedded with banks of limestone and layers of anhydrite (Uliana and Dellapé, 1981) for which Barrio (1990) defined continental and shallow marine facies associated with semiarid conditions. The interpreted sedimentary paleoenvironments range from purely continental such as ephemeral lacustrine, aeolian and fluvial systems to coastal marine paleoenvironments with development of estuaries and tidal flats Sánchez, 2011, 2013), followed by a lagoon sedimentary stage from marsh to sea with carbonate precipitation in an area protected from waves, ending with a retraction leading to the accumulation of evaporites (Page et al., 1999;Barrio, 1990). The fossil record shows that this formation includes vertebrate skeletal remains, dinosaur eggshells, theropod teeth, remains of birds, scutes from turtles and crocodiles, lungfish teeth, indeterminate plant impressions, freshwater ostracods and fresh and salt water pelecypods (Uliana and Dellapé, 1981;Hugo and Leanza, 2001;Coria et al., 2007;Salgado et al., 2007).
Allen Formation deposits record the first Atlantic ingression (Uliana, 1979), and this transgressive process was accompanied by a gradual subsidence of the basin with the development of a low-gradient intercontinental coastal plain or barrier beach (Andreis et al., 1974;Uliana and Biddle, 1988).

Facies analysis
The assessment of the 10 lithofacies presented in Table 1 is the result of the facies analysis, and the code used for their designation was based on the criteria outlined by Miall (1996) for fluvial facies, with the addition of a third capital-letter 'T' in reference to tidal environment. Thus, the first letter refers to the particle size, the second to the sedimentary structure and the third to the type of agent. The definition of lithofacies makes it possible to group them in "hierarchical units", by determining the different facies associations.

Facies associations: Tidal channels (TC)
This association is composed of the GtT, StT and ShT lithofacies (Table 1) comprising lenticular bodies with maximum thicknesses varying between 0.40 and 1.2 m, limited by gently undulating and planar surfaces. GtT and StT facies associated with three-dimensional gravel and sandy bedforms alternate at levels of 10 to 15 cm and present planar, tangential to the base and/or sigmoidal cross stratification (Figs. 3A, B). The foresets of cross stratification vary between 0.15 and 0.70 m with dip angles from 8° to 20°. Thin mud drapes are preserved on these with variations in thickness (from 0.01 m to 0.05 m) and continuity, alternating with the sandy sets (Fig. 3B).
The paleocurrents are bipolar, but with a predominant direction (pPC2- Fig. 2), and reactivation surfaces have a very low angle (Fig. 3D). The heterolithic of climbing ripple lithofacies (HrT) characterize the base of stratification and are affected by syn-sedimentary or slightly post-depositional deformation, covering thicknesses ranging between 0.05 and 0.80 m. The strata forming this association are interbedded with HfT ( Fig. 3C), HlT, HrT and FlT lithofacies corresponding to intertidal flat deposits defined below (Fig. 3D).
Interpretation: The lithofacies of this association are primarily assigned to a variety of hydrodynamic regimes. The general characteristics match those described as the product of the migration of sinuous and straight-crested sand ripples under conditions of bidirectional currents, with the presence of a dominant flow and a subordinate one (Mowbray and Visser, 1984). The sigmoidal sets in particular are generated by vortices in the flow separation area related to the acceleration and subsequent deceleration of the tidal flow. The alternating layers of sandstone and claystone reflect ebb and flow cycles in tidally dominated environments (Plink-Björklund, 2005). The sandy sets are deposited in the strong flow of the dominating current, while the mud layers correspond to the period of slack water where the flow velocity decreases and allows the deposition of suspended materials (Mellere and Steel, 1995). The presence of mud layers, the bipolarity of the paleocurrents and the type of fill allow interpreting this association as tidal channels with multiepisodic fill due to medium-scale migrating dunes where the HhT lithofacies represents the final event of filling up (Reading, 1996;Neuwerth et al., 2006). It is suggested that these channels were dominated by high-energy currents and bimodal stratification with a predominant direction reflects the direction of a dominant current and a subordinate one (Plink-Björklund, 2005). In the regional paleogeographic reconstructions for the area (Sellwood and Valdés, 2006), the first can be associated with flow events, followed by periods of slack water that allows the deposition of the suspended load parallel with the foreset surface that generates the mud drapes (Yokokawa et al., 1995;Ghosh et al., 2004). During the ebb tibe, generated structures are exposed, and according to their intensity can erode or not the mud drapes. The changes in the dip angles recorded by the foresets are related to changes in the transport of the bed load for short periods, related to the semidiurnal/diurnal variations (Mowbray and Visser, 1984;Richards, 1994).

Lithofacies Descriptions Interpretations
GtT Fine grained, regular sorted clast to matrix supported conglomerate with tangential cross-bedding at the base. Bodies of tabular and in wedge geometry up to 0.4 m thick with net wavy and erosive base.

ShcsT
Coarse to medium grained, well sorted sandstone with hummocky cross-stratification of large scale and swalley cross-bedded, with wavelengths from 2 to 5 m. Bodies up to 1.2 m thick and lateral extension up to 30 m. Internally strata show planar lamination set from 0.007 m to 0.02 m thick and very fine mud drapes.
Combined action of high energy f lows and storm events.

ShT
Fine grained sandstone and siltstones with parting lineation, low angle planar cross-bedded (5°) and parallel stratification to quasiplanar lamination with set from 0.01 m to 0.05 m thick. Bodies tabular up to 1 m thick with wavy base.
Upper flow regime, transition from dune to planar layer.

StT
Medium to fine grained, well sorted sandstone with tangential/ sigmoidal cross-bedding and planar cross-bedding with mud drapes. Bodies of tabular or in wedge geometry up to 3.5 m thick with net and erosive base. Paleocurrents directions show bipolar distribution.
Migration of 3D and 2D bedforms under bidirectional current conditions with the presence of a dominant flow and subordinate f low.

SpT
Coarse to medium grained, well sorted sandstone with planar and asymptotic base cross-bedding, and planar lamination. Bodies of wedge geometry up to 0.5 m thick and lateral extension up to 8 m.
Migration of straight crest bars in shoreface.

SrT
Coarse to medium grained, well sorted sandstone with climbing, asymmetric and symmetric in phase ripples. The ripples presents distinctive internal structure characterized by chevron-like laminations and rounded crest. The wavelengths vary from 0.03m to 0.15 m. Bodies tabular up to 0.35 m thick with net base.

HfT
Very fine grained sandstone and mud with flasser and lenticular bedding. Bodies of tabular and in wedge geometry up to 0.7 m thick.
Tidal currents with fluctuation in flow rates.

HrT
Very fine grained sandstone and mud with climbing ripples lamination gradually varying from type I to type II (Allen 1984). Bodies of tabular geometry up to 0.5 m thick.
Fluctuation in aggradations/ downstream migration rates bedforms, from subcritical flows tabular (type I and II) to supercritical flow (type III).

HlT
Very fine grained sandstone and mud with parallel lamination. Bodies of tabular geometry up to 0.7 m thick.
The sandstone deposited during peak tidal flow and the mud drapes deposited during slackwater periods.

FlT
Siltstone and mud with parallel lamination. Bodies of tabular geometry up to 0.5 m thick with net base.

Facies association: Intertidal flats (ITF)
This association is characterized by the predominant presence of heterolithic facies, including HfT, HrT, HlT and FlT lithofacies (Table 1). These are arranged with a recurring alternation, presenting few bioturbations, root molds, and continuous levels or cemented or carbonate nodules; also, common are violet and ocher colored patches or sheets, as well as millimetric levels of coal (Fig. 3C). This association forms tabular bodies ranging from 0.20 to 1.2 m in thickness with a wavy net basis interspersed with the association interpreted as tidal channels (Fig. 3D) best recorded in the southwest of the study area.
Interpretation: The characters described in this association are similar to those of the intertidal flats (Reading, 1996). Heterolytic stratification is generated by the modified flows caused by frequent beating of the sea and finely laminated sets record quadrature periods while the grouping of the thick slabs, the sicigia periods (Fan and Congxian, 2002;Eriksson et al., 2006). This preserved heterolithic stratification then suggests the action of flows with high variability at speeds that allow the transport and deposition of sands and mud, recurrently (Leckie and Singh, 1991).
The presence of millimetric carbonate levels make it possible to interpret events of reduced clastic input under conditions of restricted water environments associated with marshes, shallow marine or tidal flats (Holz, 2003;Munir et al., 2005). In particular marsh deposits are characterized by the recognition of siltypelitic sediments rich in organic matter with lesser amounts of sandstone, and levels of calcareous levels and abundant root molds in compact layers (Lander et al., 1990). This is due to the fact that marshes are confined to areas above the average elevation of water (Dalrymple et al., 1990) and poor drainage on these elevated areas within the flats favors the preservation of organic material. Because they are reducing environments, wetlands are conducive to the development of soils with gley horizons, which is why mottling is common and abundant (Aslan and Autin, 1998).
Although the low sedimentation rate favors the action of organisms, it is possible that the conditions of salinity and low nutrient availability have been the cause of the low density of bioturbations (Weimer et al., 1981). The stacking of multiple units is interpreted as vertical accretion that characterizes these flats, but it is common that such provision be affected by the installation of tidal channels (Fig. 3D).

Facies association: Subtidal flats (STF)
This association is composed of the interbedded of 10-30 cm units constituted by ShT, SrT and HfT lithofacies ( Table 1). The main feature is the domain of upper-regime sedimentary structures such as sandstone with parallel and quasiplanar lamination (Fig.  4A) and parting lineation. These units toward the top show the reworking of the wave ripples (Fig. 4B), and the increase in pelitic material and layers of gypsum more than 5 cm thick accompanied by the development of flaser and lenticular stratification (HfT). Regarding the SrT lithofacies not only does it occur at the top of the units but also in strata of up to 30 cm where the passage of symmetric and in-phase ripples to asymmetric ripples is preserved ( Fig. 4A and C). The geometry of this association is tabular, with flat base and up to 3 m thick, and has an aerial extent of more than 6 km. It is mostly preserved southern part of the area, in the vicinity of Salitral Moreno and is closely associated with the facies associations interpreted below as storminfluenced tidal flats and shoreface deposits.
Interpretation: Tractive sedimentary structures whose characteristics suggest deposits of upper flow regime sand flats dominate this association (Plink-Björklund, 2005). The presence of parallel lamination in this association demonstrates the migration of low-amplitude waves on flat surfaces by currents that reached speeds of 2 m/sec and depths less than 2-3 m or the combination of storm events and tidal action (Dalrymple et al., 1990;Fielding, 2006). The varied bedforms and the stability field of these structures is a function of changes in the speed and size of grain for a constant water depth and temperature (Fielding, 2006). The presence of parallel stratification and parting lineation characterize the upper plane bed phase (Richards, 1994;Fielding, 2006). The latter indicates the action of micro-vortexes under conditions of high-energy fluxes, which select and deposit sand grains (Allen, 1984). The flat low-angle cross stratification corresponds to long-wavelength symmetrical dunes (Saunderson and Locket, 1983) generated in stage between an upper plane bed phase and anti-dune plane (Fielding, 2006).
Two types of processes can be the cause of the presence of lenticular and flasser stratification of this association. Tide action and the variations in the velocity of flow and ebb currents contribute to the formation of such structures and to begin to form based on the migration of ripple during the decrease in the velocity of the current of tide flows. After this phase, the relative increase in the amount of mud in suspension induces 'decantation deposition' to cover the topography of the bedform at a late stage. During the period of slack water, the deposition of mud will reach a maximum thickness and if the accumulation is thick enough, the layer will be preserved, since its granulometry will require high current power for its mobilization (Reineck and Wunderlich, 1968;Hawley, 1981). When the threshold speed for the remobilization of pelitic material is reached, the removal will affect the stoss side of the ripple, and again will go into suspension while the bedform migrates and the remnants located in the foresets and troughs are preserved (Fan and Congxian, 2002). This process causes the progressive burial of the mud patches, generating flasser stratification. If the rate after the slack water period did not increase, a continuous layer of mud could be preserved. Considering that for flasser stratification to take place a high supply of suspended matter is needed, another major cause could be attributed to storm action (McCave, 1970(McCave, , 1971) since storms are the only ones capable of producing a considerable increase in suspended pelitic material (Hawley, 1981). The close relationship with ShT lithofacies make it possible to consider this as the most probable cause under conditions of high flow rate (ShT). The origin of this association is interpreted as a sand tidal flat in a subtidal environment (Dalrymple et al., 1990). (Table 1). The units are limited by undulating planar surfaces which, in some cases, reveal impressions of raindrops (Fig. 4F). The ShcsM lithofacies has a quasiplanar parallel stratification internally and wave ripple (SrT) toward the top, predominantly ripples in phase (Fig. 4G). While wavelength values of hummocky cross stratification are around 2 m, these structures also occur in large scales reaching wavelengths of up to 20 m. The basal contact is gently undulating, following the morphology of the underlying strata. thick (Fig. 4H). Internally the layers show planar lamination with a set from 0.007 to 0.02 m and very thin levels of mud in some cases. The geometry of this association is tabular or wedge-shaped up to 10 m thick with an aerial extension of more than 10 km. It has a lateral and vertical relationship with the facies associations interpreted above as subtidal flat and intratidal flat.

Facies association: Storm-influenced tidal f lats (SITF) This association includes ShcsT (Figs. 4D and E) as the dominant lithofacies, interbedded with the SrT and HfT lithofacies
Interpretation: This association suggests deposition by oscillatory periodic flows of high energy and high net sedimentation rate. The contact type shows that the depositional flow had low erodibility given the presence of mud covering the underlying strata, whereby the formation of the hummocky cross-stratification indicates accretion processes. In this process, sedimentation begins with the deposition of parallel plates in flat-bed phase, which initially adapt to the roof of the underlying strata and toward the top convert into weakly asymmetric shapes. The last deposited phase of sandstone correspond to two-dimensional symmetrical ripples, indicating oscillatory flow (Basilici et al., 2012). The presence of impressions of raindrops present in this facies association indicates its emergence during the post-storm phase. SrT development at the top of the macroforms suggests the action of waves and a high rate of sedimentation (Reineck and Singh, 1980) given the presence ripples in phase. The characteristics of this association which dominate the oscillatory flows and considering its relationship with other associations make it possible to interpret it as tidal flat deposits with strong control of storm events (Basilici et al., 2012).

Facies association: Shoreface (SF)
This association is composed of sandstone lithofacies with planar stratification (SpT) and low-angle cross stratification (ShT), and trough cross stratification (StT) to a lesser extent (Fig. 5A and B) but with little development of mud drape and strongly flattened tops. Paleocurrent directions show a bimodal and oblique arrangement (pPC1 and pPC2- Fig. 2), although conformably unimodal stacked sets predominate with direction toward the W and SW limited by low angle surfaces. The units with opposite paleoflows are stacked with a clear wedge-shaped geometry limited by undulating planar surfaces that truncate the strata. The development of wave ripplemark (Fig. 5 A, B, and C) up to 5 cm in wavelength (SrT) towards the top of each unit is characteristic. The maximum thickness is 3 m to 3 km of lateral continuity, limited by a very low angle surface with dipping toward the SE.
Interpretation: The wave-reworked, wedgeshaped lithofacies in this association can be interpreted as deposits of a shoreface environment. The bipolarity of the paleocurrents suggests alternating unidirectional currents perpendicular to the shoreline, the dominant headed toward the W and SW (towards the coast) and subordinate current to the SE. The obliquity these directions present and their bipolarity can be assigned to alternating episodes of transport toward the coast by asymmetric waves and transport toward the sea by rip currents, plus the presence of parallel lamination suggests wave wash in the shoreface zone (Colquhoun, 1995). Based on these characteristics, considering the dominance of sediment transport toward the coast, it is possible to interpret berm deposits (Hine, 1979). In these environments the low-angle cross stratification (ShT) suggests wash processes (Clifton, 1971), and the gentle slope with deposits is interpreted as a smooth surface dipping toward the sea, characteristic of shoreface. In those cases where the lamination has an opposite dipping (SpT), it is interpreted as suggesting backshore deposits, suggesting the washing of berm under storm conditions. In environments with tidal influence, significant differences in the range of these, between syzygy and quadrature, causes growth of the beach due to the aggradation of small berms on the straight section of the beach during syzygy tide. During quadrature tide in periods of fair-weather waves, the tide does not exceed the berm crest and hence the sands and gravels are deposited in the shoreface along the high-tide swash line forming a stratification consistent with the shoreface. During syzygy tide due to the increase in the tidal range, these sediments are transported toward the berm reworking the top and producing its lateral and vertical growth (Hine, 1979).

Depositional system
These facies associations that make up this depositional system hinder its definition, since there are tidal and wave influences, both controls of great importance. However, it is understood that it consists of tidal channels, tidal flats affected or not by storms and shoreface deposits.
The main development of this system's lower section is near Paso Córdoba, where tidal influence dominates, and to a lesser extent wave action, with deposits of intertidal (ITF) and subtidal flats (STF) and tidal channels (TC) with east-west and northeastsoutheast directions (Fig. 6A), with evidence of the dominance of flow currents according to paleocurrent data (pPC2- Fig. 2). In sectors, the channels' lateral and vertical relationships with the tidal flat make it possible to infer sinuous conditions with lateral accretion. As for the upper section of the system, the shoreface associations (SF) of the storm-influenced tidal flat (SITF) and subtidal flat (STF) are the dominant ones (Fig. 6B). The shoreface presents an arrangement of bimodal paleocurrents and further development of unimodal sets toward the W and SW, constituting the record of backshore washing. Considering that the bipolarity of the paleoflows suggests alternating unidirectional currents perpendicular to the coastline, it is estimated that it had a NNW-SSE direction. More precisely in the northern area (pPC1- Fig. 2) there is evidence of berm deposits with clear records of washing due to wave action as well as tides. To the southwest, near Salitral Moreno (pSM- Fig. 2) associated with beach deposits, there are extensive sand cords with a lateral continuity of more than 12 km and integrated by large hummocky structures whose characteristics described above make it possible to link them to storm-influenced tidal flats (SITF) in shallow marine environments, where tides and waves alike represented important controls on the dynamics. Moreover, there is the wave action evidenced by such structures and ripples, with paleocurrent data involving NNW and SE directions for the paleoflows. With respect to the tidal action, this is associated due to the presence of heterolithic cross-stratification, the mud drapes that cyclically alternate with climbing ripples, the flasser and wavy stratification, and bipolarity in the paleocurrent. The characteristics of this system as a whole indicate hybrid depositional conditions, in some cases showing subaerial exposure, suggesting an intertidal environment. Therefore, a tidal flat, which toward the top dominates the influence of the waves over the tide and the development of shoreface, is defined.

Discussion
Although the sedimentary record is evident in the joint action of tides and waves, there is always a predominance of control over the other. Some sedimentological features are characteristic of these and as a result it is important that they be found, interpreted and preserved. For example, muddy deposits associated with tidal flats and periods of low wave energy are key tools given that the erosional remnant serves to distinguish such an environment influenced by tides. In this paper, that tool made it possible to define facies associations with strong tidal influence. The deposits may have low preservation, and occur in drape forming rhythmites, or concentrated at the base of strata as intraclasts (Plink-Björklund, 2005;Yang et al., 2005). The presence of mud drape in the sedimentary record also contributes to infer periods of subaerial exposure of the plain because intense solar insolation and high temperatures promote the consolidation of muds during exposure and thus increases the potential for preservation (Anderson and Howell, 1984;Krogel and Fleming, 1998). Moreover, as mentioned in the introduction of this article, hummocky cross stratification indicates the domain of waves in the sedimentary record. It is known that the wavelength of this structure correlates with the wavelength of the wave, and tidal plains are shorter (<2m) than typical shorefaces (Amos et al., 1996). The smallest size is attributed to the reduction in the size of waves that occurs when waves propagate in the intertidal zone (Clifton, 1971). However, as posed by Yang et al. (2005), the presence of hummocky cross stratification with small wavelengths in itself would not be indicative of an intertidal origin because they have registered these structures lacustrine shorefaces where waves are small (Greenwood and Sherman, 1986).
These kinds of mixed-energy coastal environments with tidal and wave action are more common than you may think, and one of the best known is called open-coast tidal flats. In this type of environment, meso to macrotidal ranges are required, with highenergy waves; they characterized by concave upward profiles and development of tidal beaches with slopes between 0.15° and 3° (Yang et al., 2005). At times it can be difficult to estimate slope variations in the sedimentary record, but some clues are preserved. The increased thickness of the bar deposits suggest an increasing slope, while the presence of mud-rich flats is more likely to develop along coasts with smoother gradients. The slope determines the variation in wave energy, while the sediment grain size is also an important factor in the gradient of the coast, showing a directly proportional relationship between both controls (Yang et al., 2005). With respect to the gradient of the coast, in the case of the deposits used in this paper, besides the presence of muddy and heterolithic facies, and the wavelength values of the ShcsM, numerous authors from Windhausen (1914) and Wichmann (1927) to Uliana and Dellapé (1981);  Barrio (1990); Aguirre-Urreta et al., (2011) have indicated the existence of a shallow epicontinental seas during the Upper Cretaceous from the Atlantic transgression during the Maastrichtian, associated with the opening of the South Atlantic Ocean. The Malargüe Group has always been interpreted as the recording of the first Atlantic ingression to the basin, although recent investigations show the influence of tides on the Anacleto Formation (top of the Neuquén Group, Fig. 1D) which along with the Allen Formation constitute the transgressive systems tract (Armas and Sánchez, 2011) in the area of Cinco Saltos ( Fig. 7A and B). As mentioned above, it is difficult to define the beginning of this ingression in northwestern Río Negro due to the wide variation of fossils (continental and marine) supplied by the Allen Formation in the area. However, the interpretation of depositional systems based on the facies analysis and considering the paleogeographic reconstruction of the Upper Cretaceous according to Sellwood and Valdés (2006), it is possible to estimate some aspects of paleogeography in the central-eastern edge of the Neuquén Basin (Fig. 7). The geographical distribution of outcrops and defined facies associations propose the development of a coastal and tidal system with a NNE-SSW coastline in line with the paleocurrent analysis. The HCS is one of the most recurrent lithofacies in the system, and is associated with tidal deposits indicating deposition in shallow waters (Li et al., 2000;Yang et al., 2005). The distribution of the associations and the field relationships suggest that the best exhibitions of shoreface deposits are found to the north ( Fig. 7C and D), while the further development of tidal channels and flats occur south of the area (Fig. 7C and E). For the latter it is important to mention that Armas (2013) postulates the development of a coastal lagoon to the northwest of the studied sections at the top of the Anacleto Formation, so it is estimated that such tidal channels are part of an inlet connecting the lagoon with the sea to the east. We could also discuss whether the basis of the sections surveyed in this paper corresponds to the Allen Formation or actually forms part of the top of the Anacleto Formation. The discontinuity of the outcrops and field relationships in the study area suggests future research for adjusting the formational boundaries, although in transitional environments like these, this goal is highly complex. A very important feature of these coastal deposits is that they are laterally associated with large aeolian systems (Fig. 7C) identified by Armas and Sánchez (2013) in the Allen Formation in the area of Paso Córdoba. In this area, the authors interpreted a complex aeolian fluvial system and complex dunes with sinuous crests that form a 'draas', where aeolian system of transverse dunes with wet interdunes is installed, associated with a decrease in the availability of sediment (Armas and Sánchez, 2013). This change in type of aeolian system and decreased availability are associated with the increase in flooded areas to the north, due to the installation of the estuarine system in the area of Cinco Saltos (Fig. 7B) for the lower member of the Allen Formation (Armas and Sánchez, 2011) and the development of lagoons registered for the middle member in that area (Andreis et al., 1974). These aeolian systems are presented in concordance with the deposits studied in this paper, interpreted as tidal channel and tidal flat facies with strong wave influence, records of storm events and shoreface deposits, whose architectures are typical of epicontinental seas.
According to the tectonic history of the Neuquén Basin and sedimentological studies available to date, the deposition of Malargüe Group began at the moment Anacleto and Allen Formations came into contact (Uliana and Dellapé, 1981;Armas and Sánchez, 2011). This stage was characterized by a continued rise in eustatic level which favored the increase in accommodation and the complete installation of epeiric sea that covered much of the basin during the first Atlantic transgression. It is important to highlight that the study area for the stratigraphic period analyzed was located in the distal part of the foreland basin according to the model of Posamentier and Allen (1993), characterized by a low subsidence rate; therefore the accommodation pattern at the time of deposition was primarily controlled by the eustatic factor (Armas and Sánchez, 2011). This factor as well as an increase in the gradient of the floor of the basin toward the NE during the Late Cretaceous (Uliana and Dellapé, 1981) were the main controls on the progress of the marine systems toward continent with a subordinate influence of tectonics and subsidence on the development of the stratigraphic sequence.

Conclusions
In stratigraphic sections of the Allen Formation of the study area, ten lithofacies (one conglomeratic, five of sandstones, three heterolithics and one pelitic) were identified, which allowed the definition of four facies associations, which were interpreted as tidal channels, intertidal flat, subtidal flat, storminfluenced tidal flat and shoreface deposits.
This detailed sedimentological study contributes in part to the reconstruction of the paleogeography of the eastern edge of the basin for the Upper Cretaceous, revealing the arrangement of the coastline and the complexity of the transition environment during the Atlantic ingression. A hybrid coastal system of tidal flats with a large storm influence in some areas linked to aeolian systems previously defined in the study area by Armas and Sánchez (2013) is established for the Allen Formation.