Intraformational unconformities as a record of late Miocene eustatic falls of sea level in the Pisco Formation (southern Peru)

ABSTRACT Field mapping and sedimentological study of outcrop sections exposed along the Ica River valley permitted the establishment of a regional allostratigraphic framework for the upper Miocene portion of the Pisco Formation. The stratigraphy of the studied interval is illustrated using a new 1:20,000-scale geological map which reveals that this formation is a cyclical sedimentary succession composed of three fining-upward allomembers. The bounding surfaces defining each allomember are transgressively modified subaerial unconformities. They converge and merge landward into a single composite surface representing the time-transgressive lower boundary of the Pisco Formation. Accordingly, the extent of the stratigraphic gap associated with the basal unconformity varies significantly throughout the basin and increases toward the basin margins. The timing of allomember-bounding surfaces coincides with that of major oxygen-isotope maxima in the deep-sea oxygen isotopic record and matches the ages of eustatic sequence boundaries identified elsewhere, indicating glacio-eustatic falls due to the growth of Antarctica ice sheets as a viable mechanism for their development.

In the last few years, the stratigraphic architecture of the Pisco Formation has been documented in local studies Di Celma, Malinverno, Cantalamessa, et al., 2016;Di Celma, Malinverno, Gariboldi, et al., 2016) and constrained with biostratigraphic and geochronologic studies (Bosio et al., in press Gariboldi, 2016;Gariboldi et al., 2017). This paper forms part of this ongoing attempt to create a refined chronostratigraphic model for the vertebrate fossil localities of the Pisco Formation  and illustrates a high-resolution outcrop-based allostratigraphic framework for upper Miocene sediments of the Pisco Formation exposed along the western bank of the Ica River. Such a detailed stratigraphic analysis provides both the physical and temporal framework within which a number of broader issues concerning patterns of marine vertebrate evolution, speciation, and extinction events can be addressed.
Within strata of the Pisco Formation, some of the most important shifts in composition and diversity of the marine vertebrate community occur across a series of basin-wide intraformational unconformities . At present, the causes that drove these turnover events are uncertain and a causal link between them and the development of the unconformities has yet to be established. Here we investigate the nature of the unconformities and explore a possible relationship between them and eustatic sea-level oscillations related to global cooling. During the middle Miocene, indeed, Earth's climate transitioned from a prolonged phase of global warmth (the middle Miocene Climatic Optimum) into a Middle to late Miocene stepwise climate cooling phase (the middle Miocene Climate Transition) (Flower & Kennett, 1994;Holbourn, Kuhnt, Clemens, Prell, & Andersen, 2013). This prominent cooling phase represents one of the major steps in Cenozoic climatic evolution (e.g. Zachos, Shackleton, Revenaugh, Pälike, & Flower, 2001) and was punctuated by several short-lived episodes of glaciations, the Mi-Events of Miller, Wright, and Fairbanks (1991), and concomitant eustatic lowerings reflecting increased continental ice volume on Antarctica.
The objectives of the present contribution are: (i) to identify Miocene unconformities in the Pisco Formation; (ii) to trace them, and the sedimentary packages in between, over the entire studied area; and (iii) to assess the nature of the basin-wide discontinuities that have been used to outline the Pisco allomembers by comparing their ages with the timing of eustatic sea-level lowstands inferred from Middle and late Miocene oxygen isotope records.

Geological and stratigraphic setting
The Peruvian margin is generally recognised as a longlived convergent margin, with the oceanic Farallon/ Nazca Plate subducting below the western margin of South America at least since the beginning of the Jurassic (Coira, Davidson, Mpodozis, & Ramos, 1982). The subduction of normal oceanic lithosphere and collision of the aseismic Nazca Ridge beneath this part of the forearc zone have been associated with a wide range of effects on the leading edge of the overriding South America continent, including long-term subsidence driven by basal tectonic erosion, rapid crustal uplift and exhumation, and the development of a composite transform-convergent margin dominated by extensional and strike-slip tectonic activity (Bishop et al., 2017;Clift, Pecher, Kukowski, & Hampel, 2003;Dunbar, Marty, & Baker, 1990;Hampel, Kukowski, Bialas, Huebscher, & Heinbockel, 2004;Hsu, 1992;Kulm, Resig, Thornburg, & Schrader, 1982;León, Aleman, Torres, Rosell, & De La Cruz, 2008;Macharé & Ortlieb, 1992). According to Thornburg and Kulm (1981), two trench-parallel structural highs were formed on the continental shelf and upper slope of the Peruvian margin during late Cretaceous-early Paleogene time, namely the Outer Shelf High and the Upper Slope Ridge. These two structural highs subdivided the Peruvian offshore into a chain of fault-bounded sedimentary basins, which formed in response to prolonged extension and subsidence since the middle-late Eocene (Viveen & Schlunegger, 2018) and may be separated into an inner set of shelf basins and a seaward set of slope basins (Figure 1). In southern Peru, the onshore (East) Pisco Basin is separated from the adjacent West Pisco Basin, which is situated on the slope, by the Outer Shelf High.
The basin fill is variously folded and faulted and deformation is generally less intense as one moves upward through the sedimentary succession. In the Pisco Formation faults display dip angles mostly comprised between 75°and 85°. Kinematic indicators suggest that all of the faults are either pure normal faults or have a minor strike-slip component (Rustichelli, Di Celma, Tondi, & Bianucci, 2016). These faults display two main ranges of orientation, with the large majority of them having a NW-SE orientation and the remainder generally striking NE. Along these faults the displacement is generally minimal (metres to tens of metres) and also varies along strike, occasionally tipping out into a fracture with no displacement. Most of the fault zones are marked by an abundance of gypsum filled veins (Rustichelli, Di Celma, Tondi, Baud, & Vinciguerra, 2016).

Methods
The stratigraphic architecture of the study interval is illustrated using a new geological map (Main Map), which encompasses an area of about 200 km 2 , and two northwest-southeast and south-north trending geological cross-sections. Field mapping has been carried out at 1:10,000 scale, however the final map has been compiled at 1:20,000 scale. The allostratigraphic approach employed for the purpose of the present study is based on recommendations of the North American Commission on Stratigraphic Nomenclature (NACSN, 2005). However, the allomembers defined herein are considered to be informal stratigraphic units, designated by letters and numbers rather than formally defined names. By emphasizing the importance of chronostratigraphically significant discontinuities, the allostratigraphic approach has facilitated correlation of discontinuity-bounded stratal units and permitted an appreciation of the genetic relationships between coeval depositional settings. The bounding discontinuities that have been used to define the Pisco allomembers are represented by composite subaerial unconformity/marine-generated transgressive surfaces that, in general, are characterised by an abrupt transition from finer-grained offshore diatomites to coarser-grained shoreface sandstones. Age control for Pisco strata and bounding unconformities relies on the integration of recently published biostratigraphic data and 40 Ar/ 39 Ar age estimates .
The allostratigraphic scheme presented in this paper builds upon the original allostratigraphy established by Di Celma, Malinverno, Cantalamessa, et al. (2016) and Di Celma, Malinverno, Gariboldi, et al. (2016) for the Pisco Formation at Cerro Colorado and Cerro Los Quesos. In order to identify the depositional settings, to define stratigraphic surfaces, and to delineate genetic units, four detailed sedimentary sections (totalling some 710 m of measured stratigraphy) were logged at the decimetre scale through the exposed succession. For every single bed, lithology, colour, composition, texture, sedimentary structures, amount and type of bioturbation, palaeocurrent directions, and macrofossil content have been described. These four sections were successfully correlated by walking out a series of distinctive physical stratigraphic markers and bounding erosion surfaces in the field, by visually tracing them along outcrop photopanels, and by using high-resolution Google Earth imagery.

Description of map units
The study area is located along the eastern, internal margin of the basin, where the infilling sedimentary succession becomes stratigraphically incomplete and only the Miocene Chilcatay and Pisco formations occur at the surface, resting directly on the eroded surface of the pre-Cenozoic basement. The Eocene-Oligocene strata of the Otuma Formation are widely exposed outside of the southern and south-southeastern boundaries of the map area, but their extension in the subsurface of the mapped area is virtually unknown and remains tentative. As a consequence, the description provided below will focus only on the exposed stratigraphic units.

Basement
All the pre-Cenozoic units are grouped together under this title. Within the study area, the basement rocks consist of a complex assemblage of lower Palaeozoic gabbroic to granitoid rocks forming the San Nicolás batholiths (Mukasa & Henry, 1990), intruding a Precambrian metamorphic complex known as the Arequipa massif and covered by Jurassic volcanosedimentary rocks of the Guaneros Formation. The Guaneros Formation consists of volcaniclastic rocks with andesitic to rhyolitic composition deriving from a Jurassic calcalkaline volcanic arc, interbedded with calcareous sediments (León et al., 2008). In the study area, exposures of the basement rocks are limited and only occur between Cerro la Bruja and Cerro Blanco.

Chilcatay Formation
Sediments of the Chilcatay Formation rest either with angular unconformity on the Otuma Formation or directly on the undulatory surface of the basement rocks. Over much of the study area, Chilcatay strata dip gently to the northeast and are subdivided by an intraformational unconformity into two distinctive sediment wedges (Figure 2; Di Celma et al., 2018). The lower unit displays a stratigraphic thickness in excess of 65 m and comprises a package of massive or weakly bedded, heavily bioturbated, sandy siltstones with disperse comminuted skeletal debris. These sandy siltstones underlie and landward interfinger (from southwest to northeast) with a package of seaward-dipping, 20 m-thick clinoformed grainstones. The clinostratified beds display decimetric thickness and, generally, their dip directions indicate a distinct mode southwestward. Reworked barnacle plates, pectinids and oysters account for a large portion of the highly degraded skeletal fraction, emphasising the non-tropical character of these mixed siliciclastic-carbonate rocks (compositional mixing sensu Chiarella, Longhitano, & Tropeano, 2017) and suggesting their development within a shallow-marine carbonate factory. Similar heterozoan carbonate wedges with seaward dipping, steep clinoforms have been documented by Pomar and Tropeano (2001) and Massari and D'Alessandro (2012). These authors interpreted the clinoforms as the product of offshore-directed storm waves and currents that shed off the skeletal materials from a shoreface carbonate factory and accumulated them along a steep clinoformed ramp, within an offshore-transition setting. In this scenario, the underlying and laterally adjacent sub-horizontal sandy siltstones are interpreted to reflect deposition of the suspended load within a proximal offshore setting.
The clinostratified body of mixed siliciclastic-carbonate rocks is unconformably overlain by a finingupward sedimentary succession that is about 10 m thick at Cerro las Tres Piramides and about 25 m thick at Cerro Yesera de Amara. This sedimentary unit comprises a 3.5 m-thick package of massive, medium-grained sandstones overlain by a package of massive siltstones that represent shoreface and outer shelf depositional settings, respectively. A volcanic ash layer sampled in the upper portion of massive siltstones, just 1 m below the erosional contact with the overlying Pisco Formation, provided an 40 Ar/ 39 Ar age of 18.02 ± 0.07 Ma, whereas biostratigraphic data from the same stratigraphic interval seem to indicate a slightly younger age .

Pisco Formation
The Pisco Formation is a regionally extensive stratal unit recording the final widespread marine incursion into the Pisco Basin. Within the study area, its strata form a northeastward dipping monocline and consists mainly of marine diatomites with subordinate siltstones and sandstones. The base of the Pisco Formation has a composite nature and includes at least three regionally mappable unconformities ( Figure  3a), termed PE0.0 through PE0.2 in ascending order, which converge and merge landward (i.e. to the northeast) to form a single diachronous surface (PE0). The stratal configuration below and above this composite surface varies with position within the basin and is most commonly that of an angular unconformity where it overlies older sediments of the Chilcatay Formation, or it is a nonconformity where it rests directly on crystalline basement rocks (i.e. along the basin margins). A high-resolution allostratigraphic approach, based on mapping of three basin-wide unconformities, revealed that the investigated interval is composed of three allomembers (named alphanumerically P0, P1, P2 from older to younger). The three allomembers are thought to be equivalent to depositional sequences separated by transgressively-modified subaerial unconformities, associated with shifts in base level .
Internally, each allomember preserves a record of water-depth change during deposition and displays an overall fining-upward trend of facies associations marking a gradual vertical shift from shoreface to offshore deposits. None of the allomember-bounding discontinuities shows evidence of subaerial exposure, although each marks an abrupt shallowing and seaward shift in facies. These surfaces are penetrated by a Thalassinoides-dominated Glossifungites ichnofacies (Figure 3b), indicating the erosional exposure of sedimentary firmgrounds that were subsequently colonised by burrowing animals (e.g. MacEachern, Raychaudhuri, & Pemberton, 1992). Directly overlying the Glossifungites-demarcated surfaces is a relatively thin (up to 0.4 m thick), coarse-grained layer including a variable mixture of pebble-size phosphatic nodules, well-rounded pebble-to boulder-size basement clasts, cobble-size clasts of dolomite-cemented mudstone showing bivalve borings, internal moulds of gastropods and articulated bivalves, shark teeth, and polished bone fragments (Figure 3a and b). This coarse-grained layer is regarded as a condensed transgressive lag produced in the shoreface by low net sedimentation rates during shoreline retreat (e.g. Boessenecker, Perry, Schmitt, & Farke, 2014;Di Celma, Ragaini, Cantalamessa, & Curzio, 2002;Föllmi, 2016;Grimm, 2000). A sand-prone facies association is invariably found immediately above the aforementioned transgressive lag. It comprises yellow-weathering, intensely bioturbated, fineto very fine-grained sandstones containing elongated, U-shaped, sandfilled gutter casts (Figure 3c). These structures characterise sedimentation during storms and are of particular sedimentological significance for the interpretation of nearshore depositional environments (i.e. Chiocci & Clifton, 1991;Leckie & Walker, 1982;Myrow, 1992). The basal sandstones grade up into finely laminated gray-white diatomaceous mudstones recording deposition by suspension settling of diatom frustules in a low-energy offshore setting (Figure 3d). Subordinate lithologies include volcanic ash layers, more resistant dolomite-cemented mudstones, fineto medium-grained sandstones, and phosphatic nodule beds forming laterally continuous ledges in steep slopes. Traced northeastward to outcrops between Cerro la Bruja and Cerro Blanco, the diatomaceous mudstones gradually give way to laterally adjacent storm-influenced shoreface sandstones, displaying swaley cross-stratification, large-scale unidirectional cross-bedding, and small-scale, gravelfilled gutter casts oriented perpendicular to the palaeoshoreline (Figure 3e-g).

P0 allomember
The P0 allomember is bounded at the base by PE0.0 and at the top by PE0.1 (Figure 4a and b). This stratal unit is confined to the southern part of the study area, where it attains its greatest thickness (about 40 m) at Cerro las Tres Piramides. Based on physical stratigraphic relationships, it pinches out toward the north and the west accompanied by the merging of the PE0.0 and PE0.1 bounding surfaces.
At present, due to the lack of direct biostratigraphic and radiometric data, the age of P0 is constrained by the 40 Ar/ 39 Ar radiometric ages for volcanic ash beds in the youngest sediments of the underlying Chilcatay Formation at Cerro Yesera de Amara (18.02 ± 0.07 Ma, Di Celma et al., 2018) and for the onset of deposition of the overlying P1 allomember strata at Cerro Colorado (9.10 ± 0.04 Ma, Gariboldi et al., 2017) and at Cerros Cadenas de los Zanjones (9.00 ± 0.02 Ma, Di Celma et al., 2017). These three unconformities converge and merge landward into a single surface (informally referred to as PE0) representing the composite lower boundary of the Pisco Formation. The resulting allomembers (Amb) have been named for their lower bounding surface. Thus, the P0 allomember refers to strata lying between the PE0.0 unconformity and the next higher unconformity (PE0.1). The vertical scale is only indicative of thickness of sediment packages between unconformity surfaces. Figure 3. (a) Detail of the unconformity between the P0 and P1 allomembers at Cerro las Tres Piramides (black dotted line). This surface places coarse-grained sediments composed of tightly packed phosphate pebbles, basement clasts, shark teeth, and polished bones surrounded by a sandy matrix (above) in direct and erosive contact with diatomite mudstones (below); (b) below the surface, frequent Thalassinoides burrows are filled with small phosphate nodules; (c) close up of sand-filled gutter casts in the nearshore facies at the base of the P2 allomember. They occur parallel to each other as discrete scours that are usually less than 1.50 m in width and 0.5 m deep. In cross section, these erosional structures display a gently (concave up) curved base and outlines that flare upward. Internally, they are filled with slightly sagging sandy laminae that either drape the margins concordantly or pinch out discordantly against them (Cadenas de los Zanjones); (d) finely laminated diatomites (P2, Cerro Hueco la Zorra); (e) close up view of multiple gravel-filled gutter casts at the base of a gravel bed in nearshore facies. Typically, these linear features are 0.20-0.25 m wide and 0.15 m deep and laterally connected. Lateral spacing is regular, at distances of about 0.3 m. Their crosssectional shape is generally symmetrical and some have walls with stepped outlines. Gutters show a strongly preferred NE orientation, nearly perpendicular to the inferred shoreline trend. The infill consists of pebble-size phosphatic nodules set in a well-sorted, medium-to coarse-grained sand (base of P2, Cerro Blanco); (f) swaley cross-stratification in fine-grained sandstone produced by storm-induced oscillatory and combined flows (base of P2, Cerro Blanco); (g) oblique view of trough cross-stratification produced in medium-to coarse-grained bioclastic sandstones by migration of lunate bedforms (base of P1, Cerro la Bruja).

P1 allomember
The P1 allomember is bounded at the base by PE0.1 and at the top by PE0.2 (Figure 4b-d). The geological mapping suggests that this unit onlaps onto the underlying P0 strata and basement rocks and progressively oversteps them from southwest to northeast. The P1 allomember is exposed throughout the entire study area and is about 100 m thick at Cerros Cadenas de los Zanjones, where a continuous stratigraphic section can be measured. To help stratigraphically subdividing this unit, two marker beds, namely P1-1 and P1-2, have been correlated across part of the study area (Figure 4d).

P2 allomember
This allomember, which is more than 210 m thick at Cerro la Bruja, has 10 ledge-forming marker beds 2) at Cadenas de los Zanjones; (d) crystalline basement nonconformably overlain by P1 strata and ledge-forming marker beds within P1 and P2 (Cerro la Bruja); (e) crystalline basement nonconformably overlain by P2 strata and ledge-forming marker beds serving as formidable correlation tools over wide areas within P2 (Cerro Hueco la Zorra).
(designated P2-1 through P2-10 in ascending order, Figure 4d and e) being traceable across most of the central and northern portions of the study area and, at least in part, westward into the sedimentary succession exposed at Cerro los Quesos . Careful mapping of these marker beds between Cerro la Bruja and Cerro Blanco shows that the shallowwater strata of the P2 sequence extend beyond the northern limit of the present study area and onlap an unconformity surface cut across the P1 sequence, which they overstep north of Cerro el Brujito to onlap onto basement crystalline rocks. Lateral tracing of the same marker beds in the southern portion of the study area, between Cerros Cadenas de los Zanjones and Cerro la Mama y la Hija, is hampered by the distance of these localities and the presence of extensive sediment cover.

Origin and timing of erosional unconformities and correlation to global events
Diatom biostratigraphy combined with available 40 Ar/ 39 Ar radiometric ages indicate that strata of the P1 allomember in the study area were deposited between about 9.5 and 8.9 Ma, whereas those of the P2 allomember are younger than 8.5 Ma . A comparison between the estimated ages of allomember-bounding unconformities in the Pisco Formation and a glacio-eustatic proxy afforded by deep-sea oxygen isotopic records indicates that, within the resolution of our age model, PE0.2 approximates the drop in sea level recognised as the Mi7 event of Miller et al. (1991), an isotopic maximum dated at 8.7 Ma by Miller et al. (1998) and Westerhold, Bickert, and Röhl (2005). The PE0.1 unconformity correlates reasonably well with the Mi6 oxygen isotope maximum dated at 10.3 Ma by Miller et al. (1998) and astronomically dated at 10.4 Ma by Turco, Hilgen, Lourens, Shackleton, and Zachariasse (2001) and Westerhold et al. (2005). Assuming that this correlation between the studied PE0.1 and PE0.2 unconformities and the eustatic curve derived from deep-sea oxygen isotope records is correct and, therefore, that these breaks in deposition have an eustatic (global) rather than tectonic (local) origin, we speculate that the PE0.0 unconformity reflects the Mi5 event, an isotopic maximum dated at 11.7 Ma by Miller et al. (1998) and astronomically at 11.4 and 11.7 Ma by Turco et al. (2001) and Westerhold et al. (2005), respectively ( Figure 5). In this frame, intermittent periods of subaerial exposure during deposition of the Pisco Formation were likely produced by recurring third-order (0.5-3 m.y.) eustatic oscillations with the rate of eustatic fall exceeding that of long-term tectonic subsidence. As a main result, the several hundred metres thick sediments of the Pisco Formation reflect accommodation space produced by a combination of moderate subsidence rates and eustatic sea-level fluctuations.
An additional support to the eustatic origin for these allomember-bounding surfaces would be the similar timing of the unconformities in different ocean basins. Comparison of ages of depositional breaks in the Pisco Formation with ages of seismic sequence boundaries identified and dated in the Bahamas (Eberli, Anselmetti, Kroon, Sato, & Wright, 2002;ODP Leg 166) indicates that the timing of the PE0.1 and PE0.2 unconformities is remarkably similar with the timing of eustatic sequence boundaries lettered H (9.4 Ma) and G (8.7 Ma) along the Bahamas Transect, respectively.
Considering the above assumption and available chronological constraints, the time represented by the erosional unconformity interposed between the youngest sediments of the Chilcatay Formation (18.02 ± 0.07 Ma) and the oldest sediments of the Pisco Formation (11.7 Ma) is about 6.3 m.y. This long hiatus in the stratigraphic record indicates a prolonged phase of subaerial exposure and erosion of this marginal basin and possibly reflects a widespread pulse of uplift and compression related to the Quechua 1 tectonic event (Viveen & Schlunegger, 2018) exacerbating the effects of the long-term global lowering of sea level recorded during the the middle Miocene Climate Transition. The most recent estimates for the amplitude of the eustatic changes during this period indicate that the sea-level fell about 53-69 m between 16.5 and 13.9 Ma (John et al., 2011) and 50 ± 5 m between 13.6 and 11.4 Ma (John, Karner, & Mutti, 2004). As such, the Pisco Formation records a subsequent, widespread marine transgression and re-flooding of the forearc basin.

Allomember stacking pattern
The evidence of an eustatic origin for the allomemberbounding unconformities in the Pisco Formation does, obviously, not imply that the regional tectonic activity did not affect the overall stratigraphic architecture of this unit. Within the study area, the PE0.1 and PE0.2 surfaces onlap onto the older PE0.0 surface in a broadly northeasterly direction, merging into a composite PE0 unconformity. Accordingly, the intervening P0, P1, and P2 allomembers taper progressively to the northeast and overstep each other to lap out against the pre-Cenozoic crystalline basement, forming a set of retrogradationally stacked stratal wedges. This architectural pattern, suggesting a steady increase in the accommodation space, is indicative of deposition during a long-term period of relative sea-level rise punctuated by repeated sea-level falls and is interpreted as the product of third-order eustatic sea-level changes superimposed on prolonged extension and relatively uniform subsidence (Viveen & Schlunegger, 2018).

Conclusions
A systematic analysis of the sedimentological and stratigraphic features of the outcrop belt exposed along the western bank of the Ica River valley permitted the establishment of a regional allostratigraphic framework for the upper Miocene portion of the Pisco Formation. Coplanar subaerial unconformity/marine-generated transgressive surfaces, designated PE0.0, PE0.1 and PE0.2, can be traced throughout the entire study area and provide a basis for defining the presence of three informal allomembers, designated P0 through P2, in the order of a few tens to hundreds of metres thick. These allomembers display repetitive vertical facies successions representing a spectrum of storm- . Geochronologic data from these sections were supplemented by additional 40 Ar/ 39 Ar radiometric ages from Cerro los Quesos (CLQ, vertical dashed line). The Pisco Formation exhibits pronounced thinning to the northeast with component allomembers onlapping onto the basal composite surface PE0 and arranged in a retrogradational pattern, progressively offset to the northeast. The marine erosion surfaces that bound Pisco sequences, namely PE0.0, PE0.1, and PE0.2, were formed as a result of wave erosion in a transgressing shoreline following sea-level lowstands and, therefore, will be close to planar and inclined slightly seaward. Red lines indicate chronostratigraphic markers; (b) correlation of Pisco allomembers with late Miocene sea-level fluctuations by using integrated biostratigraphic and geochronologic data from Gariboldi et al. (2017). Diatom biostratigraphy revised following most recent publications on Equatorial Pacific diatom biostratigraphy (mostly Barron, 2003). Calculated sea-level curve from Westerhold et al. (2005).
influenced shoreface to offshore environments laid down during the relative rise of sea-level cycles.
In the study area, the Pisco Formation forms a northeastward-thinning wedge and its composite lower boundary steps up-section from the base of the P0 allomember in the southwest, through that of the overlying P1, to coincide with the base of P2 to the northeast (paleolandward). Accordingly, the extent of the stratigraphic gap associated with the basal unconformity varies significantly throughout the basin and increases toward the basin margins.
At this stage, the underlying control on these fluctuations of sea level is in part speculative. However, combined biostratigraphic and geochronologic age constraints indicate that at least two bounding unconformities developed near discrete glacial intervals inferred from deep-sea oxygen isotope records, suggesting a causal link between the observed depositional cyclicity and third-order late Miocene changes in global sea level caused by the waxing and waning of the Antarctic ice sheet. On the basis of available age control, unconformities PE0.1 and PE0.2 appear to correlate with isotopic maxima observed at about 10.4 and 8.7 Ma, respectively. We speculate that the PE0.0 unconformity at the base of the P0 allomember might correlate with the isotopic maxima observed at 11.7 Ma.
As a consequence of third-order glacio-eustatic fluctuations in sea-level acting in concert with relatively uniform subsidence rates, the study area was alternatively exposed and submerged, leading to development of bounding subaerial unconformities during glacial lowerings of sea level and deposition of allomembers during interglacial sea-level rises. Over the long-term, this combination of eustatic fluctuations and basin subsidence resulted into a punctuated relative sea-level rise, with P0, P1 and P2 allomembers stacked in a retrogradational, onlapping configuration onto the underlying Chilcatay Formation and the pre-Cenozoic basement.

Software
The geological map and associated geological sections were compiled by scanning hand drafts as black and white TIF files, and then digitising the linework using the Corel Draw X3 graphics package. By using the GIS Data processing application Global Mapper 12, contour lines for the 1:10,000 scale topographic base map were generated from digital elevation models (DEMs) based on the Shuttle Radar Topography Mission 26 (SRTM) as released by the United States Geological Survey (SRTM3 USGS version 2.1).