¹Smithsonian Tropical Research Institute, Balboa, Ancon, Panamá; Corporación Geológica ARES, Bogotá, Colombia e-mail: gbayona@cgares.org
²Ecopetrol S.A.- Instituto Colombiano del Petróleo, A.A. 4185, Bucaramanga, Santander, Colombia

A foreland basin is a dynamic system whose depositional systems migrate in response to changes in tectonic uplift patterns, sedimentary filling processes and isostatic rebound of the lithosphere. The Paleocene-middle Miocene foreland system of the Llanos foothills and Llanos basin of Colombia includes regional unconformities, abrupt changes in lithology/stacking patterns and flooding surfaces bounding reservoir and seal units. Here we integrate a systematic biostratigraphic study, stratal architecture and tectonic subsidence analyses, regional seismic profiles, and provenance data to define the diachronism of such surfaces and to document the direction of migration of foreland depozones.

In a flexurally-deformed basin, sandstone composition, rates of accommodation and sediment supply vary across and along the basin. We show how a coeval depositional profile in the Llanos foothills-Llanos foreland basin consists of litharenites interbedded with mudstones (seal rock, supplied from the orogenic front to the west) that correlate cratonward with organic-rich mudstones and coal (source rock), and to amalgamated fluvial-estuarine quartzarenites (reservoir rock, supplied from the craton to the east) adjacent to a subaerial forebulge (unconformity). This system migrated northward and eastward during the Paleocene, westward during the early-middle Eocene, and eastward during the Oligocene. In the lower-middle Miocene succession of the Llanos basin, identification of flooding events indicates a westward encroaching of a shallow-water lacustrine system that covered an eastward-directed fluvial-deltaic system. A similar process has been documented in other basins in Venezuela and Bolivia, indicating the regional extent of such flooding event may be related to the onset of Andean-scale mountain-building processes.

Keywords: nonmarine foreland, sequence stratigraphy, Colombia, Andes, tectonics, depositional system, biostratigraphy analysis, subsidence, paleocene, miocene, Mirador formation, Barco formation, Llanos Orientales.

En una cuenca flexural, la composición de las areniscas, los patrones de acomodación y aporte de sedimentos varían dentro de la cuenca. Este artículo presenta un perfil deposicional entre el Piedemonte hasta la cuenca de los Llanos el cual incluye litoareniscas interestratificadas con lodolitas (roca sello derivada del frente de deformación) que correlacionan hacia el Este con lodolitas carbonosas y carbón (roca fuente). Siguiendo hacia el Este continúan cuarzo areniscas fluvio-estuarino (roca reservorio derivada de áreas cratónicas) adyacentes a la zona de levantamiento flexural, con exposición subaérea (discordancia). Este sistema migró hacia el Norte y Este durante el Paleoceno, hacia el Oeste durante el Eoceno temprano-medio, y hacia el Este durante el Oligoceno. En la sucesión del Mioceno inferior-medio de los Llanos, la identificación de eventos de inundación indica un avance hacia el Oeste de sistemas lacustre someros, cubriendo el sistema fluvio-deltaico que avanzaba hacia el Este desde el Piedemonte. Un proceso similar ha sido documentado en otras cuencas en Venezuela y Bolivia, indicando la extensión regional de este evento de inundación, el cual puede estar relacionado con el inicio del levantamiento de los Andes.

Palabras clave: cuenca de antepaís continental, estratigrafía de secuencias, Colombia, Andes, tectónica, sistemas deposicionales, análisis bioestratigráficos, subsidencia, tectónica, Paleoceno, Mioceno, formación Mirador, formación Barco, Llanos Orientales.

INTRODUCTION

Sequence stratigraphy models have emphasized in the utility of correlating key stratigraphic surfaces with chronostratigraphic value, such as unconformities and marine flooding surfaces (Galloway, 1989; Van Wagoner, Mitchum, Campion, & Rahmanain, 1990), rather than using lithostratigraphic units. The sequence stratigraphic approach has proved to be very powerful along marine to marginal passive-margin basins for defining reservoir units filling incised valleys and for predicting the position of reservoir units in offshore settings (Vail, Mitchum, & Thompson, 1977; Van Wagoner et al., 1990). In tectonically-active continental basins, however, upstream river profiles are more vulnerable to tectonic and climate fluctuations, whereas downstream river profiles are more susceptible to sea-level fluctuations (see review in Holbrook, Scott, & Oboh-Ikuenobe, 2006). This paper concentrates on correlation of the Paleocene-middle Miocene succession that includes two reservoir units in the Llanos foothills, the Barco and Mirador formations. These fluvial units accumulated in upstream settings and at the migrating boundary of the forebulge-foredeep depozones of the Llanos foothills and Llanos foreland basin of Colombia (Bayona et al., 2006).

The flexural geometry of a foreland basin, combined with the lateral changes in rates of accommodation space/sediment supply, makes it possible to correlate the unconformity formed by forebulge uplift with areas where the rate of accommodation space/sediment supply is the highest within the foredeep. The merging of different key stratigraphic surfaces and migration of foreland depozones is well documented by empirical data (Gómez, Jordan, Allmendinger, & Cardozo, 2005a), as well as conceptual (DeCelles & Giles, 1996; Currie, 1997) and numerical (Flemings & Jordan, 1990) models.

In one section or well, depositional sequences in a migrating nonmarine foreland system are recognized on the basis of regional unconformities and changes of facies assemblages with distinctive architectural elements (e.g., low-sinuosity versus high-sinuosity fluvial systems; Currie, 1997). When the foredeep depozone migrates cratonward, the rapid increase of rate of accommodation space and decrease of sediment supply produces the vertical juxtaposition of fine-grained on coarse-grained lithologies due to a rise of base-level position. When foredeep migration occurs in the opposite direction, coarse-grained lithologies unconformably overlie the former depositional sequence due to a base-level drop. Those abrupt changes in lithology/stacking patterns can be identified in outcrop (Céspedes & Peña, 1995; Jaramillo, 1999; Gómez et al., 2005a, Bayona et al., 2006) or wells (Fajardo, Rojas, Cristancho, & Consorcio G&G Going System, 2000). Spatial and temporal variations of rates of accommodation space/sediment supply in the nonmarine Paleogene Magdalena-Llanos foreland basin of Colombia have been associated with migration of a flexural wave (Gómez et al., 2005a, Bayona et al., 2006).

This paper determines whether those abrupt changes in lithology/stacking patterns (or "base-level turnarounds" in the terminology of Ramon and Cross, 1997) can be used for correlation of reservoir units between the Llanos foothills and Llanos foreland basin. We focus on compositional and stratal architecture of reservoir units, as well as in abrupt changes in lithology/stacking patterns in units overlying reservoir units. The diachronism of such changes is documented by rigorous palynological analysis of fine-grained strata overlying reservoir beds. Provenance is also addressed by comparing units above and below base-level turnarounds because a compositional variation of framework sand grains may indicate whether tectonic uplift had exposed new source areas in the orogen (Rogers, 1994).

Regional setting and study areas

Three major orogenic belts are the result of the complex interaction of the Nazca, Caribbean and South America plates: the Western Cordillera (WC), the Central Cordillera (CC) and the Eastern Cordillera (EC) (Figure 1). The EC is a doubly vergent thrust belt placing Cretaceous and Paleogene rocks over a thick Cenozoic succession in the Magdalena Valley to the west and the Llanos basin to the east. Paleozoic and other basement rocks are exposed at the axial zone of the Eastern Cordillera (Santander and Floresta massifs) and along the southern Llanos foothills segment (Quetame and Garzon massifs). A summary of geological evolution of the Colombian Andes (Western, Central and Eastern Cordilleras) can be found in Cediel, Shaw, and Cáceres (2003).

Several stages of deformation within the EC have been recognized, including deformation in the latest Cretaceous (Branquet, Laumonier, Cheilletz, & Giuliani, 1999; Bayona et al., 2006), early Paleogene (Sarmiento-Rojas, 2001; Restrepo-Pace, Colmenares, Higuera, & Mayorga, 2004) and Oligocene (Toro, Roure, Bordas-Le Flonch, Le Cornec-Lance, & Sassi, 2004; Gómez et al., 2005a). Previous studies of the Cenozoic evolution of the EC have agreed that: (1) uplift of the CC predated that of the EC (Cooper et al., 1995; Gómez et al., 2005a), (2) the EC represents the inversion of Mesozoic graben-bounding faults (see a review and details in Mora et al., 2006), and (3) the major pulse of rock and surface uplift of the EC occurred in the Late Miocene-Pliocene (Cooper et al., 1995; Toro et al., 2004; Gómez et al., 2005a).

The Paleocene-Miocene succession in the Llanos foothills and Llanos basin record the filling of a complex nonmarine foreland system that developed since latest Maastrichtian time (Sarmiento-Rojas, 2001; Bayona et al., 2006; Cortés et al., 2006). The depositional environment changed from dominantly marine during the latest Cretaceous to marginal and continental during the Paleocene. The primary mechanisms driving this shift is the increasing rate of sediment supply associated with exhumation and denudation of the Central Cordillera, according to Gómez et al. (2005a). Uplifts or bypassing areas disrupted the Paleocene foreland configuration along the axial zone and eastern borders of the basin (Fabre, 1981, 1987; Sarmiento-Rojas, 2001; Bayona, Cortés, Jaramillo and Llinás, 2003; Restrepo-Pace et al., 2004, Cortes, 2004; Bayona et al., 2006). Uplifts along the western margin of the basin in Early Eocene time are recorded by an angular unconformity underlying Eocene strata in the Magdalena Valley. For the Late Eocene-Oligocene, Toro et al. (2004), Gómez et al. (2005a) and Bayona et al. (2006) used the evidence of growth-strata patterns, flexural subsidence and thermochronology to document the presence of a western (Magdalena Valley) and an eastern (axial EC-Llanos) depocenters filled with continental and marginal synorogenic strata. There is no record of Miocene strata in the axial zone or on the flanks of the EC, but continental synorogenic deposition continued both in the Llanos (Leon Formation) and Magdalena basins.

Stratigraphic correlation of the paleocene-middle miocene succession

The Paleocene - middle Miocene foreland succession of the Llanos foothills and Llanos basin includes, in stratigraphic order, the Barco, Cuervos, Mirador, Carbonera and Leon formation. At least eight flooding and unconformities surfaces (Cooper et al., 1995) or 12 base-level turnarounds (Reyes et al., 1991; Fajardo et al., 2000) have been identified in the eastward-thinning clastic wedge of the Llanos basin (Cooper et al., 1995). These studies proposed that the Paleocene Barco Formation and the Eocene Mirador Formation, the reservoir units in the Llanos foothills, are bounded at the base by an unconformity (or an abrupt decrease of accommodation space/sediment supply rates) and at the top by a flooding surface (or an abrupt increase of accommodation space/sediment supply rates). In the Llanos basin, sandstones of the lowermost Carbonera Formation overlie Paleocene strata or the basal Cenozoic unconformity in distal settings (Jaramillo et al., 2007).

We studied the Paleocene-middle Miocene succession along three Stratigraphic Lines of Correlation (SLC), named the northern, central and southern SLC (Figure 1), in order to characterize the regional extent and stratigraphic significance of those base-level turnarounds in the Llanos foothills and Llanos basin. Key sedimentological, vertical stacking patterns, framework sandstone composition, and paleontological (pollen-spore-dinoflagellate, foraminifera, macrofossils content) data were analyzed for a total of 21 wells using electric logs, cores and cuttings. Three composite stratigraphic sections on the Llanos foothills and axial zone of the Eastern Cordillera were constructed to link depositional systems between the Llanos foothills and Llanos basin. One thousand three hundred fifty two palynological samples from wells and outcrops were analyzed to support stratal correlations (Table 1; see details for construction of the Cenozoic palynological zonation in Jaramillo & Rueda 2004; Jaramillo, Muñoz, Cogollo, & Parra, 2005). All the above data, in conjunction with the vertical arrangement of lithofacies, provide the support for the interpretation of depositional environments.

Major unconformities and Miocene flooding surfaces were selected mainly as stratigraphic datums for stratigraphic correlations presented in Figures 2 to 4. Stratigraphic descriptions focus on intervals of change in lithology/stacking patterns. From base to top these intervals are: the Barco-Cuervos contact, the Cuervos-Mirador contact, the Mirador-Carbonera contact, lower and upper contacts of basal Carbonera sandstones in the Llanos basin, and changes across Lower and Middle Miocene flooding surfaces. In order to understand the diachronism of changes in lithology/stacking patterns within the Paleocene-middle Miocene succession, we used our palynological control as a guide to the placement Cenozoic epoch boundaries (e.g. lower Paleocene/upper Paleocene; Oligocene/Miocene), as precisely as age resolution allowed.

Barco-Cuervos and Cuervos-Mirador contacts

Distribution and age. The Barco and Cuervos formations (or Lower Socha-Upper Socha formations in the axial zone of the Eastern Cordillera) comprise the lowermost nonmarine foreland succession in the Llanos foothills and Llanos basin. The Barco Formation unconformably overlies fine-grained strata of the Colon-Mito Juan Formation (or Guaduas Formation in the axial zone of the Eastern Cordillera). The Barco-Cuervos succession is recorded in the northern and central SLC, but coeval strata are truncated 10 - 30 km north of the southern SLC. The Barco-Cuervos succession thins eastward and forms an eastward-stepping stratal package onlapping the basal Cenozoic unconformity in the western Llanos basin. Strata beneath the unconformity are older toward the craton, whereas strata above the unconformity are younger cratonward. The age of the Barco Formation is lower Paleocene in the central SLC and mostly upper Paleocene in the northern SLC. The Upper Paleocene Cuervos Formation correlates eastward in the Llanos basin with the Barco Formation (Figures 2 and 3). Lower Mirador beds accumulated earlier in the Llanos foothills than equivalent sandstones beds in the axial zone of the Eastern Cordillera (Figures 2 and 3).

Lithology and stacking patterns. The Barco Formation rests upon a regional disconformity (Figures 2 and 3). This unit consists of amalgamated quartzose sandstones with sedimentary structures that change upsection from (1) cross-bedded fine- to coarse-grained fining-upward sandstones that are locally conglomeratic in Cocuy (south of northern SLC in Figure 1; Fabre 1981) to (2) interbeds of fining-upward sandstone and mudstone beds, and to (3) bidirectional cross-bedded and bioturbated heterolithic laminated sandstones at the top with organic rich mudstones and coal beds in the Llanos foothills. In the Llanos basin, the Barco Formation changes from aggraded fine-to-coarse-grained sandstones to fining-upward sandstones at the top.

The contact between the Barco and Cuervos formations in the northern and central SLC is marked by an increase of fine-grained lithologies with a dominant aggradational stacking pattern. Lower beds of the Cuervos Formation consist mainly of laminated to locally bioturbated dark-gray organic-rich claystones and mudstones with thin coal seams and pass upsection to fining- and coarsening-upward successions with cross-bedded, ripple and wavy laminated sandstones.

Sandstone composition. An upward increase in polycrystalline quartz and unstable lithic fragments has been documented from quartzarenites to sublitharenites in the Barco - Lower Socha units, and to litharenites in the Cuervos - Upper Socha formations (Fabre, 1981; Jaramillo, Roa, & Torres, 1993; Céspedes & Peña, 1995; Vergara & Rodríguez, 1996; Mesa, 1997, 2003, 2004). Lithic fragments comprise up to the 33% of the framework grains and include mainly sedimentary fragments (siltstones, chert) and a minor component of metamorphic (gneiss, schist, phyllite), and igneous rock fragments. Quartzarenite composition dominates in lower Mirador Formation sandstone beds (Mesa, 1997, 2004), but Cretaceous foraminifera fragments occur locally in the matrix of conglomerates (Céspedes & Peña, 1995) and chert, argillaceous, gneiss, schist and igneous fragments (Mesa, 1997, 2004) have been reported.

Palynological association. In the axial zone of the EC, the Lower Socha Formation includes continental palynofloras (angiosperms, mostly palms) (Pardo, 2004). In the Llanos foothills, palynofloras in the Barco Formation are dominated by very few species, with Nypa pollen (a brackish-water palm) as the most common element. Few dinoflagellates grains were reported in the Llanos foothills of the central SLC near the top of the Barco Formation (well C&BA in Figure 3). Lower beds of Cuervos and Upper Socha formations are dominated by a low-diversity assemblage composed mainly of palms, Araceaea, Proteaceae, Bombacacae and legumes suggesting a forest with tropical affinities. At the Lower Socha-Upper Socha contact, the pollen content increases significantly (Pardo, 2004), whereas morphospecies diversity increases at the Cuervos-Mirador contact (Jaramillo et al., 2006).

Depositional environment. The vertical arrangement of lithofacies and palynological association of the Barco-Cuervos-Mirador succession has been interpreted as product of deposition in fluvial systems. The uppermost Barco sandstones show more fluvial influence in the axial Eastern Cordillera than in the Llanos foothills; in the latter area, bi-directional cross beds, heterolithic lamination (Cazier et al., 1995; Reyes, 1996) and few dinoflagellate grains suggest more tidally-influenced deposition. The Cuervos-Upper Socha accumulated in coastal plains with more estuarine influence toward the Llanos foothills. The excellent pollen recovery in lower beds of the Cuervos and Upper Socha units suggest a very humid environment with small fluctuations of the water table that was near the surface (Pardo, 2004).

Summary. The Barco-Cuervos succession comprises a conformable fluvial-coastal plain package bounded by two disconformities (or abrupt decrease of accommodation/sediment supply rate). This succession is an eastward-stepping stratal package onlapping the basal Cenozoic unconformity. Amalgamated quartzose sandstone beds overlie the lower (Barco Formation) and upper (Mirador Formation) disconformities, but Barco sandstones are younger eastward, whereas Mirador sandstones are younger westward. The Barco-Cuervos contact is younger cratonward and records (1) an abrupt increase of accommodation space, and (2) input of sedimentary and metamorphic lithic fragments to the coastal plain system. The Barco-Cuervos contact records the lateral migration of depositional environments, whereas the Cuervos-Mirador contact records an abrupt decrease of accommodation/sediment supply rate within a continental setting.

Mirador-Carbonera contact in the Llanos foothills; lower and upper contacts of the basal Carbonera sandstones in the Llanos basin

Distribution and age. Strata of lower and middle Eocene age (Mirador and Picacho formations) are reported only to the west of the Llanos foothills (Jaramillo et al., 2007) (Figure 3). Upper Eocene strata include the contact between the Mirador-Carbonera formations in the Llanos foothills (or the Picacho-Concentracion formations in the axial Eastern Cordillera). In the Llanos basin, either a basal sandstone or fine-grained successions of the Carbonera Formation unconformably cover Paleocene and Campanian beds in the northern and central SLC (Figures 2 and 3), and Upper Cretaceous and Paleozoic beds in the southern SLC (Figure 4). The age of lower Carbonera strata is younger northward and eastward and varies from late Eocene to mid-Oligocene. Oligocene foreland deposits of the Carbonera Formation formed a wedge-like basin geometry; the thickness of Oligocene strata in the Llanos basin is at least half of the thickness preserved in the Llanos foothills (Figures 3 and 4). Upper Oligocene strata completely covered the Llanos basin with exception of distal uplifted blocks in the southern Llanos basin (Figure 4).

Lithology and stacking patterns. Lithofacies associations of the upper Mirador Formation vary along the Llanos foothills. In the northern SLC, the upper Mirador sandstones show complete fining-upward successions and coal interbeds (Reyes, 2004). In the Llanos foothills of the central SLC, fine to medium-grained quartzarenites have more interbeds of organic mudstones with plant remains that show evidence of deposition in a brackish-water environment, including: (1) a diverse ichnofacies association (Ophiomorpha, Thalassionoides, Psilonichnus and Diplocraterion; Pulham, Mitchell, MacDonald, & Daly, 1997), (2) couplets in foreset laminations, and (3) wavy and flaser lamination. Overlying Carbonera strata in the Llanos foothills of the northern and central SLC consist dominantly of 250 m-thick dark green and gray laminated and poorly bioturbated mudstones (Jaramillo, 1999; Jaramillo & Dilcher, 2001; Mora & Parra, 2004) grading to coarsening-upward successions (Figure 3).

In the Llanos basin, the basal sandstone of the Carbonera Formation consists of amalgamated to fining-upward successions of fine- to coarse-grained quartzarenites, locally conglomeratic. In the central SLC, both the thickness of this basal sandstone and the thickness of the sandstone-mudstone transition increase eastward (Figure 3). In the southern SLC, the thickness increases eastward up to >100 m (wells F1 and Ca-1 in Figure 4), and decreases farther east. The thickness changes are attributed to the filling of incised valleys or onlap over structural highs (e.g., Las Brujas high). Dark-colored mudstones, shales and thin coal interbeds diachronously overlie these basal sandstones. At the western segment of the southern SLC, this fine-grained interval is followed by a fining-upward succession, which is not present in any other area.

Sandstone composition. Sandstone of the Mirador Formation and the basal sandstones of the Carbonera are quartzarenites (Mesa, 1997, 2004). The lower half of the Carbonera Formation is comprised of quartz-arenites and sublitharenites (Fabre 1981; Moreno & Velásquez 1993; Cardona & Gutiérrrez 1995). In the axial zone of the Eastern Cordillera, conglomerate beds of the Concentracion Formation contain chert and fossiliferous limestones of Cretaceous age (Reyes & Valentino 1976).

Depositional Environment. Upper Mirador sandstones have been interpreted as fluvial channels and mouth-bar sands deposited on a coastal plain in the central Llanos foothills (Cazier et al., 1995; Fajardo, 1995; Warren & Pulham, 2001), with more continental influence in the northern Llanos foothills (Reyes, 2004). Lithological associations of lowermost Carbonera beds in the Llanos basin have been interpreted as a change from channel-fill processes in fluvial systems to tidally-influenced coastal plains and delta bays to the top of the basal sandstones (Fajardo et al., 2000).

Summary. Quartzose sandstones of the Mirador Formation in the Llanos foothills and the basal sandstones of the Carbonera Formation in the Llanos basin overlie two different stratigraphic surfaces. The former unit disconformably overlies Paleocene strata, whereas the latter is an eastward-stepping stratal package that onlap the basal Cenozoic unconformity in the Llanos basin. The nature of the upper contacts of these sandstone units also differ. The Carbonera-Mirador contact in the Llanos foothills area shows evidence of a late Eocene marine transgression, as documented by ichnofossils and paleontological data. In contrast, the upper contact of the basal sandstones in the Llanos basin is more transitional, and the overlying fine-grained strata only show evidence of eastward and northward migration of coastal-plain environments.

Lower and Middle Miocene flooding surfaces of the upper Carbonera-Leon formations

Distribution and age. Strata deposited during early and middle Miocene time correspond to the upper Carbonera and Leon formations in the Llanos foothills and Llanos basin. Sandstones of the upper Carbonera record the last pulses of pre-Andean deformation, whereas shales of the Leon Formation record a basin-wide flooding event during the Middle Miocene. The Carbonera-Leon contact is defined by the onset of fine-grained deposition that occurred earlier in the eastern Llanos basin than in the Llanos foothills (Figure 4).

Lithology and stacking patterns. In the Llanos foothills of the northern and central SLC, lower Miocene Carbonera strata consist of coarsening-upward successions that begin with laminated mudstones grading to tabular and wavy laminated, locally bioturbated, fine-grained sandstones (Fajardo et al., 2000). In the central Llanos foothills, these successions include coal interbeds, feldspar-bearing fine-grained muddy sandstones and locally conglomeratic massive to cross-bedded sandstones (Mora & Parra, 2004).

In the northern and central Llanos basin, fining- and coarsening-upward successions constitute the Carbonera Formation. This unit grades eastward to a finer-grained succession. In the southern SLC, upper Carbonera strata change from coarsening-upward successions, which include calcareous mudstones to the south and east, to variable grain-size trends toward the western margin (Figure 4).

Lower strata of the Leon Formation consist of wavy laminated, bioturbated, and varicolored mudstones interbedded with tabular-bedded, bioturbated quartzarenites (Geoestratos-Dunia, 2003) in the Llanos foothills. In the Llanos basin, lower Leon strata consist of a thick interval of dark-colored laminated mudstones and shales. Sandstone interbeds in the Leon Formation increase northward and westward (Cooper et al., 1995; Fajardo et al., 2000).

Sandstone composition. Compositional analysis of well cuttings indicates an increment of feldspar and lithic fragments, and the sandstones at the western margin of the Llanos basin are compositionally more immature than sandstones at the eastern margin of the basin.

Palynological association. Poaceae and Mauritia (monostand palms commonly found of riparian forests) dominate the assemblage in upper Oligocene strata, suggesting development of savannas. Lower-middle Miocene palynofloras have a lower diversity than Eocene-lower Oligocene palynofloras. In some fine-grained intervals in the upper Carbonera and at the Carbonera-Leon contact, fresh-water (lacustrine algae), but also brackish-water palynological associations, dinoflagellate cysts and foraminifera test linings are observed (Figure 4), and contributed to the evidence for defining flooding events in the lower-middle Miocene strata. Samples from these levels have been analysed for foraminifera and thecamoebians content, showing strong evidences of alternating fresh to brackish-water conditions (Flavia Fiorini, STRI post-doctoral fellow, personal communication, 2007). Therefore, flooding surfaces were identified in wells on the basis of biostratigraphic data and the shale signature at the base of coarsening-upward successions in gamma-ray profiles.

Summary. Fluvial-deltaic depositional systems of the lower Miocene Carbonera formation advanced eastward, transporting synorogenic detritus toward shallow lacustrine environments. Fluvial systems draining the Guyana craton also supplied detritus into the lacustrine setting. Westward-stepping and encroaching of fine-grained lacustrine deposits indicate a decreasing influx of synorogenic detritus and an increase of accommodation space in the Llanos basin.

SEISMIC IMAGES OF STRATIGRAPHIC SURFACES

Seismic facies, unconformity surfaces, and onlap/truncations of reflectors were identified on 2D seismic reflection profiles. Unconformities, abrupt changes in lithology/stacking patterns and flooding events identified in wells were tied to seismic reflectors and traced regionally along composite seismic lines. Upper Paleocene Barco-Cuervos strata in the northern SLC are defined as discontinuous seismic reflectors onlapping onto an unconformity toward the eastward trend (Figure 5). The unconformity underlying Carbonera strata in the Llanos basin is identified by truncation of seismic reflectors (Figure 5).

Strong seismic reflectors in the southern SLC can be traced for several tens of kilometers, providing opportunities to decipher onlap of seismic reflectors and lateral continuity of seismic facies. The eastward onlap of reflectors of the diachronous basal sandstone and mudstones of the lower Carbonera Formation can be traced laterally for several tens of kilometers (Figure 6). In the eastern side of Figure 6, reflectors corresponding to lower Miocene strata tend to be more continuous. The uppermost continuous reflector corresponds to the Carbonera-Leon contact that separates a seismic stratigraphy facies of several continuous reflectors from a facies with fewer and more discontinuous reflectors. The former represents the sandstone-mudstone successions of the Carbonera Formation and the latter represents the dominance of fine-grained lithologies of similar acoustic impedance of the Leon Formation. However, when the top of the Carbonera Formation reflector is traced westward, the equivalent reflector on the western side of Figure 6 corresponds to reflectors within the Carbonera Formation.

ONE DIMENSIONAL SUBSIDENCE ANALYSIS AND EUSTASY

Backstripping techniques were used to decompact the measured stratigraphic thickness of each section to define a relation between changes in lithology/stacking pattern and slope configurations of tectonic subsidence curves. We followed the methods specified in Allen and Allen (1992) for construction of tectonic subsidence curves, and we used first-order sea-level curve of Haq, Hardenbol, and Vail (1987) for correction of eustasy.

The results of our one-dimensional backstripping indicate that in most cases the tectonic subsidence curves are not straight lines, suggesting changes in their tectonic subsidence history. The Barco-Cuervos and Mirador-Carbonera contacts do not necessarily coincide with the point of an abrupt increase of tectonic subsidence (Figure 7). The increase of tectonic subsidence across the Barco-Cuervos contact occurred earlier in the Llanos foothills and axial Eastern Cordillera of the central SLC (59 - 62 my), then advanced to the northern SLC (58 - 60 my), and finally migrated to the Llanos basin at around 57 - 59 my.

The Late Eocene increase in tectonic subsidence and the onset of fine-grained Carbonera/Concentracion deposition become younger basinward. However, fine-grained deposition occurred slightly earlier in the Llanos foothills (Middle-Late Eocene; 41 - 47 my) than in the axial Eastern Cordillera and in the western side of the southern SLC (Late Eocene 36 - 42 my). In the central segment of the southern SLC and in the Llanos basin of the central SLC, fine-grained deposition began in early Oligocene (30 - 33 my), and reached structural highs on the northern and southern Llanos basins at the Oligocene-Miocene boundary (23 - 25 my). A slight increase in the rate of tectonic subsidence occurs nearly at the top of the Carbonera-Leon formations in the Llanos foothills and Llanos basin (aprox. 14 my), and the slope gently increases after 10 my.

DISCUSSION

The Barco and Mirador formations, the two main reservoir units in the Llanos foothills, are bounded at the base by an unconformity and at the top by abrupt changes in lithology/stacking patterns (or base-level turnarounds). Even though these surfaces can be recognized in outcrops, wells and seismic lines, our palynological data indicate that strata overlying those surfaces are diachronous. In addition, changes in stratal architecture and sandstone composition across these surfaces, and the eastward onlap of seismic reflectors on the unconformity may be explained by migration of a flexural basin.

The Paleocene Barco-Cuervos succession and Upper Eocene-Oligocene strata document the cratonward migration of depositional environments as the foredeep depozone migrated eastward. In both cases the onset of sandstone deposition is diachronous and becomes younger eastward (Figure 8). The onlap of seismic reflector on the unconformity and the quartzose composition of sandstone beds indicate accumulation adjacent to the forebulge. At the same time and toward the west, a conformable change from amalgamation of fluvial-estuarine sand bars (reservoir units) to deposition of fine-grained siliciclastic sediments on coastal and flood plains represent the abrupt change in lithology/stacking patterns of fine-grained strata (e.g., Cuervos Formation) overlying those quartzose units (e.g., Barco Formation). These changes are not flooding surfaces, since there is not a marine transgression or accumulation below a deep water column (i.e. a lake). Toward the western side of the foredeep depozone, sedimentary and metamorphic rock fragments in upper Paleocene Cuervos sandstones, and feldspar-bearing sandstones in Oligocene Carbonera strata accumulated adjacent to uplifted hinterland blocks (Bayona et al., 2006). Therefore, the Barco-Cuervos succession and the lower beds of the Carbonera Formation record two episodes of eastward migration of the foreland basin (Figure 8), which are also coincident with events of increasing rates of tectonic subsidence.

The diachronous onset of Mirador and Carbonera deposition in the Llanos foothills and Eastern Cordillera indicate a period of westward migration of the foreland basin. Amalgamated quartzose sandstones of the Mirador and Picacho formations in the Llanos foothills and axial zone of the Eastern Cordillera accumulated in a period of very slow tectonic subsidence during early and middle Eocene time. The thick record of Eocene deposition in the axial zone of the Eastern Cordillera and the absence of Eocene strata in the Llanos basin support the westward migration of the flexural wave (Figure 8).

The westward encroaching of shallow-water deposition during early to middle Miocene time caused the flooding of the fluvial-deltaic system in the western part of the Llanos basin. Flooding events are well recognized across the Llanos basin, allowing a chronostratigraphic correlation in the foreland basin (Figure 6 and 8). This pattern of deposition suggests a regional transgression of a broad shallow-water lacustrine system, as reported in other South American basins. Miocene lacustrine and lagoonal deposits have been also documented in the Colombian Amazonas basin (Hoorn, Guerrero, Sarmiento, & Lorente, 1995), in southern Andean foreland basins (the Pebas sea of Marshall & Lundberg, 1996; Uba, Heubeck, & Hulka, 2005), and in the northern Magdalena Valley (Gómez et al., 2005a) as well as in sedimentary basins with other tectonic regimes such as in the 7 km-thick extensional Urumaco basin in northern Venezuela (Hambalek et al., 1994).

CONCLUSIONS

• The Paleocene-Oligocene succession between the axial zone of the Eastern Cordillera and the Llanos foreland basin of Colombia records three migration events of the foreland depozones that controlled the present geometry of reservoir units. Diachronous deposition of the Barco-Cuervos formation during the Paleocene recorded the eastward onlap of foreland deposits and the northward advance of synorogenic deposition. Later in Eocene time, westward migration of the foredeep depozone favored deposition of amalgamated quartzose sandstones of the Mirador Formation and bypass on the Llanos basin. In late Eocene to Oligocene time, the foredeep depozone advanced rapidly eastward, favoring thick fine-grained deposition in the Llanos foothills and diachronous onset of fluvial sandstone deposition on the Llanos basin.

• Unconformities that underlie amalgamated fluvial sandstones of the Barco and Mirador formations in the Llanos foothills should be considered only as the base of synorogenic successions, whereas the conformable contact between fluvial-estuarine sandstone and marginal mudstones should be used in a first instance to constrain the diachronous migration of the flexural wave.

• During the early and middle Miocene, the foreland basin recorded a westward-stepping transgression of shallow-water systems onto fluvial-deltaic systems. Lower to middle Miocene flooding surfaces are the most appropriate for strata correlation in the Llanos foreland basin because of the basinwide extension, seismic signature, and lack of significant diachronism of strata bounding those surfaces. Andean-scale mountain-building processes and rise of sea level in the Miocene interacted to generate the conditions for rapid and basinwide encroaching of a lacustrine-lagoonal system over a nonmarine foreland basin.

• Although lithology and vertical stacking patterns are fundamental for the understanding of stratal architecture, the success of correlation in a continental nonmarine foreland basin depends upon the integration of data from other techniques. The correlation relies primarily on recognition of the diachronism of potential key stratigraphic surfaces within a chronostratigraphic framework (Figure 8). Secondly, tectonic subsidence and provenance analyses supply important evidence to determine whether deposition took place during an episode of increasing tectonic subsidence and unroofing of new source areas. The tracing of key stratigraphic surfaces along seismic lines allows consideration lateral changes of lithofacies and strata stacking pattern within the basin. The integration of all these elements is fundamental to the determination of depositional sequences within a tectonically-active continental foreland basin.

ACKNOWLEGMENTS

This research was supported by the Instituto Colombiano del Petroleo (ICP), Ecopetrol S.A., the Smithsonian Paleobiology Endowment Fund, and the Unrestricted Endowments SI Grants. Thanks to the Biostratigraphic Team at the Instituto Colombiano del Petroleo (ICP) for its continuous support. Discussions with to J. Aristizabal, E. Cardozo, J. Rubiano, N. Gamba and P. Villamarín contributed to understand different interpretations of the Paleogene stratigraphy in the Llanos basin, and comments of two anonymous reviewers contributed to improve the flow of the English. Constructive comments of Tom Becker, Sanjeev Gupta, Guy Plint and John Holbrook helped to shape the content of this manuscript.

REFERENCES

Allen, P., & Allen, J. (1992). Basin analysis, principles and applications. London, Blackwell Scientific Publications.        [ Links ]

Audemard, F. E., & Audemard, F.A. (2002). Structure of the Merida Andes, Venezuela: relations with the South America-Caribbean geodynamic interaction. Tectonophysics, 345 (1-4): 299-327.        [ Links ]

Baby, P., Rochat, P., Mascle, G., & Herail, G. (1997). Neogene shortening contribution to crustal thickening in the back arc of the Central Andes. Geology, 25 (10): 883-886.        [ Links ]

Bayona, G., Cortés, M., Aristizabal, J.J., Jaramillo, C., Ojeda, G., Reyes-Harker, A., Rueda, M., & Villamarín, P. (2006). Pinch out of upper cretaceous-oligocene reservoir units in the Llanos Basin of Colombia: a result of flexural deformation in a broken foreland. America Association of Petroleum Geologists, Annual Meeting, Abstract Volume, 8, Houston, USA.        [ Links ]

Bayona, G., Cortés, M., Jaramillo, C., & Llinás, R.D., (2003). The Fusagasugá succession: a record of the complex Latest Cretaceous-pre-Miocene deformation between the Magdalena Valley and Sabana de Bogotá areas. VIII Simposio Bolivariano de Exploración en las Cuencas Subandinas, 180-193.         [ Links ]

Branquet, Y., Laumonier, B., Cheilletz, A., & Giuliani, G. (1999). Emeralds in the Eastern Cordillera of Colombia: two tectonic settings for one mineralization. Geology, 27 (7): 597-600.        [ Links ]

Cardona, P. J., & Gutiérrez, G. Z. (1995). Estratigrafía y ambientes de depósito de la formación Carbonera en un área al Noroeste de Yopal - Casanare (Colombia). Tesis de pregrado Fac.Ciencias, Universidad Nacional de Colombia, Bogota, 83 pp.        [ Links ]

Cazier, E.C., Hayward, A.B., Espinosa, G., Velandia, J., Mugniot, J.F., & Leel, W. C. (1995). Petroleum geology of the Cusiana field, llanos basin foothills, Colombia. America Association of Petroleum Geologists Bulletin, 79 (10): 1444-1463.        [ Links ]

Cediel, F., Shaw, R., & Cáceres, C. (2003). Tectonic assembly of the northern Andean block, in Bartolini, C., Buffler, R., & Blickwede, J., eds., The Circum-Gulf of Mexico and the Caribbean: hydrocarbon habitats, basin formation and plate tectonics, Colombia. America Association of Petroleum Geologists Memoir, 79: 815-848.        [ Links ]

Céspedes, S., & Peña, L. (1995). Relaciones estratigráficas y ambientes de depósito de las formaciones del terciario inferior aflorante entre Tunja y Paz del Río, Boyacá. Tesis de pregrado Fac. de Ciencias, Universidad Nacional de Colombia, Bogotá, 50 pp.        [ Links ]

Coates, A. G., Collins, L. S., Aubry, M., & Berggren, W.A. (2004). The geology of the Darien, Panamá, and the late miocene-pliocene collision of the Panamá arc with northwestern South America. Geological Society of America, Bulletin, 116 (11-12): 1327–1344; doi: 10.1130/B25275.1.        [ Links ]

Cooper, M. A., Addison, F.T., Alvarez, R., Coral, M., Graham, R.H., Hayward, A. B., Howe, S., Martínez, J., Naar, J., Peñas. R., Pulham, A. J., & Taborda, A. (1995). Basin development and tectonic history of the Llanos basin, eastern cordillera, and middle Magdalena Valley, Colombia. America Association of Petroleum Geologists, Bulletin, 79 (10): 1421-1443.        [ Links ]

Cortés, M., Bayona, G., Aristizabal, J., Ojeda, G., Reyes-Harker, A. & Gamban, N. (2006). Structure and kinematics of the eastern foothills of the eastern cordillera of colombia from balanced cross-sections and forward modelling. Asociación colombiana de geólogos y geofísicos del petróleo. Memorias del IX Simposio Bolivariano de Cuencas Subandinas, 14, in CD.        [ Links ]

Currie, B. (1997). Sequence stratigraphy of nonmarine Jurassic-Cretaceous rocks, Central Cordilleran foreland-basin system. Geological Society of America Bulletin. 109 (9): 1206-1222.        [ Links ]

DeCelles, P., & Giles, K. (1996). Foreland basin systems. Basin Research, 8 (2): 105-123.        [ Links ]

Fabre, A. (1981). Estratigrafía de la Sierra Nevada del Cocuy, Boyacá y Arauca, Cordillera Oriental (Colombia). Geología Norandina, 4: 3-12.        [ Links ]

Fabre, A. (1987). Tectonique et géneration d'hydrocarbures: un modèle de l'evolution de la Cordillère Orientale de Colombie et dubBassin de Llanos pendant le crétacé et le tertiaire. Archives des Sciences Genève, 40: 145-190.        [ Links ]

Fajardo, A. (1995). 4-D Stratigraphic architecture and 3-D reservoir fluid-flow model of the Mirador formation, Cusiana field, foothills area of the Cordillera Oriental, Colombia. Master’s Thesis, Colorado School of Mines, USA.        [ Links ]

Fajardo, A., Rojas, E., Cristancho, J, & Consorcio G&G Going System, L. (2000). Definición del Modelo estratigráfico en el intervalo Cretáceo tardío a mioceno medio en la cuenca Llanos Orientales y Piedemonte llanero. Informe final, Ecopetrol S.A. - Instituto Colombiano del Petróleo (ICP).        [ Links ]

Flemings, P.B., & Jordan, T.E. (1990). Stratigraphic modelling of foreland basins: Interpreting thrust deformation and lithosphere rheology. Geology, 18 (5): 430-434.        [ Links ]

Galloway, W. E. (1989). Genetic stratigraphic sequences in basin analysis in: architecture and genesis of flooding-surface bounded depositional units. America Association of Petroleum Geologists Bulletin, 73 (2): 125-142.        [ Links ]

Geoestratos-Dunia (2003). Control geológico de campo del bloque Sirirí y áreas aledañas. Informe final, Ecopetrol S.A. - Instituto Colombiano del Petroleo (ICP), 181.        [ Links ]

Gómez, E., Jordan, T., Allmendinger, R.W., & Cardozo, N. (2005a). Development of the Colombian foreland-basin system as a consequence of diachronous exhumation of the northern Andes. Geological Society of America Bulletin, 117 (9-10): 1272-1292; doi: 10.1130/B25456.1.        [ Links ]

Gómez, E., Jordan, T., Allmendinger, R. W., Hegarty, K., & Kelley, S. (2005b). Syntectonic cenozoic sedimentation in the northern middle Magdalena valley basin of Colombia and implications for exhumation of the northern Andes. Geological Society of America Bulletin, 117 (5-6): 547–569.        [ Links ]

Guerrero, J., & Sarmiento, G. (1996). Estratigrafía física, palinológica, sedimentológica y secuencial del cretácico superior y paleoceno del Piedemonte Llanero, implicaciones en exploración petrolera. Geología Colombiana, 20: 3-66.        [ Links ]

Hambalek, N., Rull, V., De Digiacomo, E., & Diaz de Gamero, M.L. (1994). Evolución paleoecológica y paleoambiental de la secuencia del Neogeno en el surco de Urumaco: estudio palinológico y litológico. Boletín de la Sociedad Venezolana de Geología, 191: 7-19.        [ Links ]

Haq, B.U., Hardenbol, J., & Vail, P. (1987). Chronology of fluctuating sea levels since the triassic. Science, 235:1156-1166.        [ Links ]

Holbrook, J., Scott, R.W., & Oboh-Ikuenobe, F. (2006). Base-level buffers and butresses: a model for upstream versus downstream control on fluvial geometry and architecture within sequences. Journal of Sedimentary Research, 76 (1): 162-174.        [ Links ]

Hoorn, C., Guerrero, J., Sarmiento, G.A., & Lorente, M.A. (1995). Andean tectonics as a cause for changing drainage patterns in Miocene northern South America. Geology, 23 (3): 237–240.        [ Links ]

Jaramillo, C. (1999). Middle paleogene palynology of Colombia, South America: Biostratigraphic, sequence stratigraphic and diversity implications. Ph. D. Thesis, University of Florida, USA.        [ Links ]

Jaramillo, C. (2002). Response of tropical vegetation to paleogene warming. Paleobiology, 28: 222-243.        [ Links ]

Jaramillo, C., & Dilcher, D.L. (2001). Middle paleogene palynology of central Colombia, South America: A study of pollen and spores from tropical latitudes. Palaeontographica B, 258: 87-213.        [ Links ]

Jaramillo, C., Muñoz, F., Cogollo, M., & Parra, F. (2005). Quantitative biostratigraphy for the Cuervos formation (paleocene) of the llanos foothills, Colombia: Improving palynological resolution for oil exploration, in Powell, A.J., & Riding, J.B., eds., recent developments in Applied Biostratigraphy. The Micropalaeontological Society, Special Publications, 145–159.        [ Links ]

Jaramillo, C., Rueda, M., & Mora, G. (2006). Cenozoic plant diversity in the neotropics. Science, 311 (5769):1893-1896.         [ Links ]

Jaramillo, C., & Rueda, M. (2004). Impact of biostratigraphy on oil exploration. Memorias de la tercera convencion técnica de la Asociación Colombiana de Geólogos y Geofísicos del Petróleo, 7, Bogotá, Colombia.        [ Links ]

Jaramillo, C., Rueda, M., Bayona, G., Santos, C., Flórez, & Parra, F. (2007). Biostratigraphy breaking paradigms. The absence of the eocene Mirador formation in the eastern llanos of Colombia. SEPM Special issue. Biostratigraphic Applications.        [ Links ]

Jaramillo, L., Roa, E., & Torres, M. (1993). Relaciones estratigráficas entre las unidades arenosas del paleógeno (paleoceno) del Piedemonte llanero y la parte media de la cordillera oriental. Tesis de Pregrado Fac.Ciencias,, Universidad Nacional de Colombia, Bogotá, 40 pp.        [ Links ]

Marshall, L.G., & Lundberg, J.G. (1996). Miocene deposits in the Amazonian foreland basin (technical comments). Science, 273 (5271): 123–124.        [ Links ]

Mesa, A. (1997). Diagenesis and reservoir quality of the Guadalupe, Barco and Mirador formations (Campanean to Eocene), Llanos basin, Colombia. Ph. D. Thesis, Johannes Gutenberg Universitat, Mainz, Germany, 141 pp.        [ Links ]

Mesa, A. (2003). Descripción petrológica de muestras de zanja del pozo Gibraltar 1, ST-3 y ST-5 durante las actividades de perforación. Informe final, Ecopetrol S.A. - Instituto Colombiano del Petróleo (ICP).        [ Links ]

Mesa, A. (2004). Descripción detallada petrográfica del pozo Gibraltar 2. Informe final, Ecopetrol S.A. - Instituto Colombiano del Petróleo (ICP).        [ Links ]

Mora, A., & Parra, M. (2004). Secciones estratigráficas de las formaciones Guadalupe, Barco y Carbonera, anticlinal del Guavio, Piedecuesta. Informe final, Instituto Colombiano del Petróleo, 28 pp.        [ Links ]

Mora, A., Parra, M., Strecker, M.R., Kammer, A., Dimate, C., & Rodríguez, F. (2006). Cenozoic contractional reactivation of Mesozoic extensional structures in the Eastern Cordillera of Colombia. Tectonics, 25. TC2010, doi:10.1029/2005TC001854.        [ Links ]

Moreno, J., & Velásquez. M. (1993). Estratigrafía y tectónica en los alrededores del municipio de Nuchia departamento de Casanare, Colombia. Tesis de Pregrado Fac. Ciencias, Universidad Nacional de Colombia, Bogotá.        [ Links ]

Ojeda, G. Bayona, G., Pinilla, J., Cortés, M., & Gamba, N. (2006). Subsidence and geodynamic analysis of The Llanos Basin. Linking mountain building and basin filling processes. Asociación Colombiana de Geólogos y Geofísicos del Petróleo, Memorias del IX Simposio Bolivariano de Cuencas Subandinas, in CD.        [ Links ]

Pardo, A. (2004). Paleocene-eocene palynology and palynofacies from northeastern Colombia and western Venezuela. Ph. D. Thesis, Universite de Liege, 103 pp.        [ Links ]

Pulham, A. J., Mitchell, A., MacDonald, D., & Daly, C. (1997). Reservoir modeling in the Cusiana field, llanos foothills, Eastern Cordillera: characterization of a deeply-buried, low-porosity reservoir. VII Simposio Bolivariano, 198-216.        [ Links ]

Ramón, J.C., & Cross, T.A. (1997). Characterization and prediction of reservoir architecture and petrophysical properties in fluvial channel sandstones, Middle Magdalena Basin, Colombia. CT&F – Ciencia, Technología y Futuro, 1 (3): 19-46.        [ Links ]

Restrepo-Pace, P., Colmenares, F., Higuera, C., & Moyorga, M. (2004). A fold-and-thrust belt along the western flank of the Eastern Cordillera of Colombia - Style, kinematics, and timing constraints derived from seismic data and detailed surface mapping, in McClay, K., ed., Thrust tectonics and hydrocarbon systems, Colombia. America Association of Petroleum Geologists Memoir, 82: 598-613.        [ Links ]

Reyes, A. (1996). Sedimentology and stratigraphy and its control on porosity and permeability Barco formation, Cusiana field, Colombia. Master’s thesis, University of Reading, Postgraduate Research Institute of Sedimentology, 70 pp.        [ Links ]

Reyes, A. (2004). Taller corazón Gibraltar-2: Bucaramanga. Reporte Interno, Ecopetrol S.A.        [ Links ]

Reyes, I., & Valentino, M.T. (1976). Geología del yacimiento y variabilidad de las características geoquímicas del mineral de hierro en la región de Paz Vieja (municipio de Paz del Río, departamento de Boyacá). Primer Congreso Colombiano de Geología, Memorias, Bogotá, Colombia, 267-324.        [ Links ]

Reyes, J., Silva, M., Munar, F., Lasso, A., Bohórquez, J., Valderrama, J., Cadena, A., Velasco, J., Rendón, A., & Blanco, N. (1991). Objetivos Estratigráficos en la subcuenca Apiay-Ariari. Memorias del IV Simposio Bolivariano, Bogotá, Colombia, 27.        [ Links ]

Rogers, R.R. (1994). Nature and origin of through-going dicontinuities in nonmarine foreland basin strata, upper cretaceous, Montana. Implications for aequence Analysis. Geology, 22: 1119-1122.        [ Links ]

Sarmiento-Rojas, L.F. (2001). Mesozoic rifting and cenozoic basin inversion history of the Eastern Cordillera, Colombian Andes, Inferences from tectonic models. Bogotá, Ecopetrol-Netherlands research, School of Sedimentary Geology, 295.        [ Links ]

Toro, J., Roure, F., Bordas-Le Flonch, N., Le Cornec-Lance, S., & Sassu, W. (2004). Thermal and kinematic evolution of the Eastern Cordillera fold and thrust belt, Colombia, in Swennen, R., Roure, F., and Granath, J.W., eds., Deformation, fluid flow, and reservoir appraisal in foreland fold and thrust belt Colombia. America Association of Petroleum Geologists Hedberg Series, 1: 79– 115.        [ Links ]

Uba, C., Heubeck. C., & Hulka, C. (2005). Facies analysis and basin architecture of the Neogene Subandean synorogenic wedge, southern Bolivia. Sedimentary Geology, 180 (3-4): 91-123.        [ Links ]

Vail, P.R., Mitchum, R.M.J., & Thompson, S.I. (1977). Seismic stratigraphy and global changes of sea level, part 3, relative changes of sea level from coastal onlap, in Payton, C.E., ed., Seismic stratigraphy - applications to hydrocarbon exploration. American Association of Petroleum Geologists Memoir, 26: 63-81.        [ Links ]

Van Wagoner, J., Mitchum, R.M.J., Campion, K., & Rahmanain, V. (1990). Siliciclastic sequence stratigraphy, in well logs, cores, and outcrops; concepts for high resolution correlation of time and facies. Methods in exploration series 7, Colombia. America Association of Petroleum Geologists, 55.        [ Links ]

Vergara, L., & Rodríguez, G. (1996). Consideraciones sobre la petrografía y diagénesis de los grupo Guadalupe (Cordillera Oriental) y Palmichal (Piedemonte Llanero). Geología Colombiana, 21: 41-63.        [ Links ]

Warren, E.A., & Pulham, A.J. (2001). Anomalous porosity and permeability preservation in deeply buried tertiary and mesozoic sandstones in the cusiana field, llanos foothills, colombia. Journal of Sedimentary Research, 71 (1): 2-14.         [ Links ]

(Received Aug. 23, 2007; Accepted Oct. 12, 2007)