Reactivation of basement faults: interplay of ice‐sheet advance,glacial lake formation and sediment loading

The Emme Delta is a small glacilacustrine delta, which developed on the southern flank of the Wesergebirge Mountains in NW Germany. Shallow shear‐wave seismic surveys allow a detailed assessment of the structural style of the delta body. Two different fault systems are developed within the delta, both showing syn‐sedimentary activity. The faults have planar to slightly listric geometries and show vertical offsets in a range of 2–15 m. They form small graben and half‐graben systems, which locally show roll‐over structures. The fill of the half‐grabens has a wedge‐shaped geometry, with the greatest sediment thickness close to the fault. The fault system in the upper portion of the Emme Delta is restricted to the delta body and probably gravity induced. In the lower portion of the delta, normal faults occur that originate in the underlying Jurassic basement rocks and penetrate into the delta deposits. The grid of seismic lines shows that the normal faults are trending E–W. This fits to a late Triassic–early Jurassic deformation phase in the Central European Basin System. We hypothese that these faults were reactivated during the Pleistocene by the advancing ice‐sheet, water and sediment loading. Based on the seismic data set, an overall model for the reactivation of the basement fault was developed. The advancing ice‐sheet caused far field extension, which might have reactivated pre‐existing normal faults. Later, the fault activity was enhanced due to sediment and water loading. In addition, high pore pressure due to lake formation might have supported the slip processes along the faults. After glacial unloading and lake drainage, the fault activity stopped.     

Middle Miocene–Pliocene siliciclastic influx across a carbonate shelf and influence of deltaic sedimentation on shelf construction,Northern Carnarvon Basin,Northwest Shelf of Australia

Middle Miocene to Pliocene siliciclastics of the Bare Formation represent a long‐lived (ca. 11 Myr) break in the otherwise carbonate‐dominated shelf of the Northern Carnarvon Basin, Northwest Shelf of Australia. The quartz‐sandstone interval is correlated with the appearance of spectacular clinoform sets mapped on 3D and dense 2D seismic data. Twenty‐seven clinoform sets are interpreted as delta lobes primarily based on their plan‐view morphology (strike‐elongate to lobate features) and their 40–100‐m‐high clinoform amplitudes. The delta lobes were deposited on outer‐shelf to shelf‐edge positions, and the older deltas show evidence of a higher degree wave reworking than the younger deltas. Measurements of the along‐strike (migration) and down‐dip (progradation) movement of these deltas are compared with relative sea‐level behaviour inferred from shelf‐edge trajectory analysis. Delta lobes exhibit greater lateral shifting during relative sea‐level rise, whereas delta lobes are more restricted to dip‐oriented fairways during sea‐level fall, although no major incised valleys have been identified. Long‐term (cumulative) progradation of this delta system and subsequent backstepping correlates with long‐term sea‐level fall and rise during the late middle and late Miocene. In addition, a long‐term northeastward migration trend for these delta lobes was likely a result of localized uplift of an inversion anticline in the Rosemary–Legendre Trend; the growth of this anticline probably steered the fluvial source for the delta system towards the northeast. The Bare Formation siliciclastic influx correlates with other middle Miocene increases in siliciclastic sediment supply worldwide. Global cooling and a shift to more arid conditions, negatively influencing vegetation cover, may have combined with more seasonally variable rainfall to generate the high sediment supply that built the deltas. Retreat of the siliciclastics could correlate with ice‐sheet growth in the Northern Hemisphere and/or increase in the Indonesian Throughflow and Leeuwin Current (ca. 1.6 Ma), which might have modified climate regionally.     

Late Ordovician,deep‐water gravity‐flow deposits,palaeogeography and tectonic setting,Tarim Basin,Northwest China

The Upper Ordovician in the Tarim Basin contains 5000–7000 m of siliciclastic and calciclastic deep‐water, gravity‐flow deposits. Their depositional architecture and palaeogeographical setting are documented in this investigation based on an integrated analysis of seismic, borehole and outcrop data. Six gravity‐flow depositional–palaeogeomorphological elements have been identified as follows: submarine canyon or deeply incised channels, broad and shallow erosional channels, erosional–depositional channel and levee–overbank complexes, frontal splays‐lobes and nonchannelized sheets, calciclastic lower slope fans and channel lobes or sheets, and debris‐flow complexes. Gravity‐flow deposits of the Sangtamu and Tierekeawati formations comprise a regional transgressive‐regressive megacycle, which can be further classified into six sequences bounded by unconformities and their correlative conformities. A series of incised valleys or canyons and erosional–depositional channels are identifiable along the major sequence boundaries which might have been formed as the result of global sea‐level falls. The depositional architecture of sequences varies from the upper slope to abyssal basin plain. Palaeogeographical patterns and distribution of the gravity‐flow deposits in the basin can be related to the change in tectonic setting from a passive continental margin in the Cambrian and Early to Middle Ordovician to a retroarc foreland setting in the Late Ordovician. More than 3000 m of siliciclastic submarine‐fan deposits accumulated in south‐eastern Tangguzibasi and north‐eastern Manjiaer depressions. Sedimentary units thin onto intrabasinal palaeotopographical highs of forebulge origin and thicken into backbulge depocentres. Sediments were sourced predominantly from arc terranes in the south‐east and the north‐east. Slide and mass‐transport complexes and a series of debris‐flow and turbidite deposits developed along the toes of unstable slopes on the margins of the deep‐water basins. Turbidite sandstones of channel‐fill and frontal‐splay origin and turbidite lobes comprise potential stratigraphic hydrocarbon reservoirs in the basin.     

Quantitative progradation dynamics and stratigraphic architecture of ancient shallow‐marine clinoform sets: a new method and its application to the Upper Jurassic Sognefjord Formation,Troll Field,offshore Norway

This article presents a new numerical inversion method to estimate progradation rates in ancient shallow‐marine clinoform sets, which is then used to refine the tectono‐stratigraphic and depositional model for the Upper Jurassic Sognefjord Formation reservoir in the super‐giant Troll Field, offshore Norway. The Sognefjord Formation is a 10–200‐m thick, coarse‐grained clastic wedge, that was deposited in ca. 6 Myr by a fully marine, westward‐prograding, subaqueous delta system sourced from the Norwegian mainland. The formation comprises four, 10–60‐m thick, westerly dipping, regressive clinoform sets, which are mapped for several tens of kilometres along strike. Near‐horizontal trajectories are observed in each clinoform set, and the sets are stacked vertically. Clinoform age and progradation rates are constrained by: (i) regionally correlatable bioevents, tied to seismically mapped clinoforms and clinoform set boundaries that intersect wells, (ii) exponential age–depth interpolations between bioevent‐dated surfaces and a distinctive foreset‐to‐bottomset facies transition within each well, and (iii) distances between wells along seismic transects that are oriented perpendicular to the clinoform strike and tied to well‐based stratigraphic correlations. Our results indicate a fall in progradation rate (from 170–500 to 10–65 km Myr^?1) and net sediment flux (from 6–14 to ≤1 km² Myr^?1) westwards towards the basin, which is synchronous with an overall rise in sediment accumulation rate (from 7–16 to 26–102 m Myr^?1). These variations are attributed to progradation of the subaqueous delta into progressively deeper waters, and a concomitant increase in the strength of alongshore currents that transported sediment out of the study area. Local spatial and temporal deviations from these overall trends are interpreted to reflect a subtle structural control on sedimentation. This method provides a tool to improve the predictive potential of sequence stratigraphic and clinoform trajectory analyses and offers a greater chronostratigraphic resolution than traditional approaches.     

Low‐accommodation and backwater effects on sequence stratigraphic surfaces and depositional architecture of fluvio‐deltaic settings (Cretaceous Mesa Rica Sandstone,Dakota Group,USA)

The adequate documentation and interpretation of regional‐scale stratigraphic surfaces is paramount to establish correlations between continental and shallow marine strata. However, this is often challenged by the amalgamated nature of low‐accommodation settings and control of backwater hydraulics on fluvio‐deltaic stratigraphy. Exhumed examples of full‐transect depositional profiles across river‐to‐delta systems are key to improve our understanding about interacting controlling factors and resultant stratigraphy. This study utilizes the ~400 km transect of the Cenomanian Mesa Rica Sandstone (Dakota Group, USA), which allows mapping of down‐dip changes in facies, thickness distribution, fluvial architecture and spatial extent of stratigraphic surfaces. The two sandstone units of the Mesa Rica Sandstone represent contemporaneous fluvio‐deltaic deposition in the Tucumcari sub‐basin (Western Interior Basin) during two regressive phases. Multivalley deposits pass down‐dip into single‐story channel sandstones and eventually into contemporaneous distributary channels and delta‐front strata. Down‐dip changes reflect accommodation decrease towards the paleoshoreline at the Tucumcari basin rim, and subsequent expansion into the basin. Additionally, multi‐storey channel deposits bound by erosional composite scours incise into underlying deltaic deposits. These represent incised‐valley fill deposits, based on their regional occurrence, estimated channel tops below the surrounding topographic surface and coeval downstepping delta‐front geometries. This opposes criteria offered to differentiate incised valleys from flood‐induced backwater scours. As the incised valleys evidence relative sea‐level fall and flood‐induced backwater scours do not, the interpretation of incised valleys impacts sequence stratigraphic interpretations. The erosional composite surface below fluvial strata in the continental realm represents a sequence boundary/regional composite scour (RCS). The RCS’ diachronous nature demonstrates that its down‐dip equivalent disperses into several surfaces in the marine part of the depositional system, which challenges the idea of a single, correlatable surface. Formation of a regional composite scour in the fluvial realm throughout a relative sea‐level cycle highlights that erosion and deposition occur virtually contemporaneously at any point along the depositional profile. This contradicts stratigraphic models that interpret low‐accommodation settings to dominantly promote bypass, especially during forced regressions. Source‐to‐sink analyses should account for this in order to adequately resolve timing and volume of sediment storage in the system throughout a complete relative sea‐level cycle.     

Bathymetry of Lake Lisan controls late Pleistocene and Holocene stream incision in response to base level fall

This paper examines the millennial-scale evolution of the longitude profile of Nahal (Wadi) Zin in the Dead Sea basin in the northern Arava valley, Israel. Nahal Zin has incised ~ 50 m into relatively soft late Pleistocene Lake Lisan sediments. Incision was forced by the regressive (> 10 km) lake level fall of a total of > 200 m of Lake Lisan from its highest stand at ~ 25 ka and exposure of the lake-floor sediments to fluvial and coastal processes. Alluvial cut terraces of the incising channel are well preserved along the 17.5 km of the lowermost reach of Nahal Zin. At its outlet into the Dead Sea basin, Nahal Zin deposited a Holocene alluvial fan at the base of a 10–80 m high escarpment in unconsolidated sediments. The escarpment is associated with the Amazyahu fault, which forms the southern structural boundary of the present Dead Sea basin. Geomorphic mapping, optically stimulated luminescence (OSL) ages, and soil stratigraphy allowed correlation of terrace remnants and reconstruction of several past longitudinal profiles of Nahal Zin and its incision history. Together with the published lake level chronology, these data provide an opportunity to examine stream incision related to base level lowering at a millennial scale. OSL ages of the terraces fit relatively well with the established lake level chronology and follow its regression and fall. For a few thousands of years the longitudinal profile response to the lake level fall was downstream lengthening onto the exposed former lake bed. Most of the incision (~ 40 m) occurred later, when the lake level reached the top of the Amazyahu fault escarpment and continued to drop. The incision was a relatively short episode at about 17 ka and cut through this escarpment almost to its base. The fast incision, its timing, and the profiles of the incising channels indicate that the escarpment was an underwater feature and was not formed after the lake retreated.This fairly simple scenario of regressive lake level fall and knickpoint exposure and incision is modeled here using a one-dimensional numerical incision model based on a linear diffusion equation. The calculated diffusion coefficient fits earlier results and data obtained from other streams in the area and confirms the upscaling of this simple model to the millennial scale.     

The role of climate variation in delta architecture: lessons from analogue modelling

Sequence‐stratigraphic models for fourth to sixth order, glacio‐eustatic sequences based only on relative sea‐level variations result in simplified and potentially false interpretations. Glacio‐eustatic sea‐level variations form only one aspect of cyclic climate variation; other aspects, such as variations in fluvial water discharge, vegetation cover, weathering and sediment supply can lead to variable sediment yield, thus adding complexity to sequence‐stratigraphic patterns normally attributed to sea‐level variations. Analogue flume models show a significant impact of water discharge on the timing and character of sequence boundaries, and on changes in the relative importance of systems tracts, as expressed in sediment volumes. Four deltas, generated under the influence of an identical sea‐level curve, and affected by different water‐discharge cycles were generated in the Eurotank facility: (1) constant discharge; (2) high‐frequency discharge variations (HFD); (3) discharge leading sea level by a quarter phase; (4) discharge lagging sea level by a quarter phase. HFD shift the parasequence stacking pattern consistently but do not alter large‐scale delta architecture. Water‐discharge changes that lead sea‐level changes result in high sediment yield during sea‐level rise and in the poor development of maximum flooding surfaces. Delta‐front erosion during sea‐level fall is expressed by multiple, small channels related to upstream avulsions, and does not result in an incised valley that efficiently routs sediment to the shelf edge. When water‐discharge changes lag sea‐level changes, sediment yield is high during falling sea level and results in rapid progradation during forced regression. Erosion from incised valleys is strong on the proximal delta top and dissipates towards the delta front. The combination of high discharge and sea‐level fall provides the most efficient mode of valley incision and sediment transport to the shelf edge. During sea‐level rise, low water discharge results in sediment starvation and well‐developed maximum flooding surfaces. Water‐discharge variations thus alter sequence‐stratigraphic patterns and provide an alternative explanation to the amplitude of sea‐level fall for generating either type 1 or 2 erosional unconformities.     

Sequence stratigraphy of a Paleogene coal bearing rim syncline: interplay of salt dynamics and sea‐level changes,Schöningen,Germany

Salt rim synclines contain important hydrocarbon and coal resources in central Europe. The Schöningen salt rim syncline is filled with >300 m of Early to Middle Eocene unconsolidated clastics with interbedded lignitic coal seams that are mined at the surface. In this study, 357 lithologic logs are integrated with measured outcrop sections and paleo‐botanical data to interpret the depositional environments and sequence stratigraphic framework of the rim syncline fill. As salt withdrew, it generated an elongate mini‐basin that mimicked an incised valley. The sustained accommodation and slow broadening of the syncline affected the stratigraphic architecture and contributed to the preservation of coal units. The clastic units in the syncline filled in seven depositional stages: (1) tidally influenced fluvial estuarine channels; (2) mixed tide‐ and wave‐ dominated estuaries; (3) prograding wave dominate deltas; (4) transgressive shoreline deposits; (5) braided fluvial channels; (6) estuaries; and (7) prograding tide‐dominated channels. The succession defines four 3rd order sequences and several higher order sequences that are possibly related to Milankovitch cycles. The higher order sequences are dominantly characterized by stacked transgressive cycles of thick, lowstand coals overlain by estuarine sands. The nearly continuous warm and wet Eocene climate was conducive to continuous peat production with a climatic overprint recorded in the mire type: ombrotrophic mires developed in wetter times and rheotrophic mires developed in relatively drier conditions pointing to the presence of orbitally controlled seasonality. Both mire types were impacted by the interplay of subsidence and base‐level. The continuous dropping of the mires below base‐level via subsidence protected the mires against erosion and may account for the absence of coals outside of the rim synclines in the region.     

Former glacial lakes of the Bunger Hills,Antarctica

The Bunger Hills in East Antarctica occupy a land area of approximately 400 km². They have been exposed by Holocene retreat of the Antarctic ice sheet and its outlet glaciers. The accompanying sea level rise flooded the marine inlets that now separate the northern islands and peninsulas from the major part of the hills. During deglaciation the continental ice sheet margin retreated south‐eastwards with several temporary halts, during which ice‐dammed lakes were formed in some valleys. These lakes were maintained long enough to permit formation of beaches of sand and gravel, and for the erosion of shore platforms and low cliffs in bedrock. Around the western end of Fish Tail Bay impressive shoreline features 20 m above sea level define a former ice‐dammed lake that was 5.5 km long. A similar 7 km long former ice‐dammed lake was formed at Lake Dolgoe. The more extensive and deeper glacial lake is revealed by well‐developed and preserved shoreline features cut at 29 m which is 16 m above present lake level. In addition, several small ice‐dammed lakes existed temporarily near Lake Shchel and in the valley to the west. Lake Fish Tail existed more than 6,900 ¹⁴C years ago and Lake Shchel probably more than 6,680 ¹⁴C years ago. It is inferred that the shore platforms and beaches were formed by lake ice and wave action over considerable periods when the lakes were impounded by steep cold ice margins. There appears to have been a balance between meltwater input and evaporative loss from the lakes in the cold dry continental climate. There is no evidence for rapid lake level fluctuations, and there was very little input of clastic sediment. This resulted in poor development of deltaic and rhythmically laminated lake floor deposits. It is suggested that such deposits are more characteristic of ice‐dammed lakes formed in association with wet‐based temperate ice than those associated with dry‐based polar ice.     

Assessment of the Evolution in Velocity of Two Debris‐Covered Valley Glaciers in Nepal and New Zealand

Feature tracking of orthorectified pairs of Advanced Spaceborne Thermal Emission and Reflection Radiometer satellite images is used to calculate velocities for the Tasman Glacier, New Zealand (2002–2014) and the Khumbu Glacier, Nepal (2001–2008). Velocities in the middle and upper ablation zones of both glaciers show a long‐term decrease of ～10–20%, while the terminus of Khumbu Glacier has remained near stagnation throughout the study period. In contrast, there has been a recent acceleration of the lower terminus of Tasman Glacier, from ～5 m a^–1 in 2002 to 40 m a^–1 in 2014. Both of these glaciers have an extensive supraglacial debris cover across their lower ablation regions, with the Khumbu Glacier terminating on land and the Tasman Glacier terminating in a proglacial lake. The rapid recent increase in velocity of the terminus of Tasman Glacier is closely correlated with the increase in size of its proglacial lake. These results indicate the complex dynamic changes that mountain valley glaciers may undergo in response to long‐term negative mass balance.     

100 Myr record of sequences, sedimentary facies and sea level change from Ocean Drilling Program onshore coreholes, US Mid-Atlantic coastal plain

We analyzed the latest Early Cretaceous to Miocene sections (～110–7 Ma) in 11 New Jersey and Delaware onshore coreholes (Ocean Drilling Program Legs 150X and 174AX). Fifteen to seventeen Late Cretaceous and 39–40 Cenozoic sequence boundaries were identified on the basis of physical and temporal breaks. Within‐sequence changes follow predictable patterns with thin transgressive and thick regressive highstand systems tracts. The few lowstands encountered provide critical constraints on the range of sea‐level fall. We estimated paleowater depths by integrating lithofacies and biofacies analyses and determined ages using integrated biostratigraphy and strontium isotopic stratigraphy. These datasets were backstripped to provide a sea‐level estimate for the past ～100 Myr. Large river systems affected New Jersey during the Cretaceous and latest Oligocene–Miocene. Facies evolved through eight depositional phases controlled by changes in accommodation, long‐term sea level, and sediment supply: (1) the Barremian–earliest Cenomanian consisted of anastomosing riverine environments associated with warm climates, high sediment supply, and high accommodation; (2) the Cenomanian–early Turonian was dominated by marine sediments with minor deltaic influence associated with long‐term (10⁷ year) sea‐level rise; (3) the late Turonian through Coniacian was dominated by alluvial and delta plain systems associated with long‐term sea‐level fall; (4) the Santonian–Campanian consisted of marine deposition under the influence of a wave‐dominated delta associated with a long‐term sea‐level rise and increased sediment supply; (5) Maastrichtian–Eocene deposition consisted primarily of starved siliciclastic, carbonate ramp shelf environments associated with very high long‐term sea level and low sediment supply; (6) the late Eocene–Oligocene was a starved siliciclastic shelf associated with moderately high sea‐level and low sediment supply; (7) late early–middle Miocene consisted of a prograding shelf under a strong wave‐dominated deltaic influence associated with major increase in sediment supply and accommodation due to local sediment loading; and (8) over the past 10 Myr, low accommodation and eroded coastal systems were associated with low long‐term sea level and low rates of sediment supply due to bypassing.     

Post glacial lacustrine event sedimentation in an ancient mountain setting: Carboniferous Lake Malanzán (Western Argentina)

Lacustrine deposits of the Malanzán Formation record sedimentation in a small and narrow mountain paleovalley. Lake Malanzán was one of several water bodies formed in the Paganzo Basin during the Late Carboniferous deglaciation. Five sedimentary facies have been recognized. Facies A (Dropstones-bearing laminated mudstones) records deposition from suspension fall-out and probably underflow currents coupled with ice-rafting processes in a basin lake setting. Facies B (Ripple cross-laminated sandstones and siltstones) was deposited from low density turbidity currents in a lobe fringe environment. Facies C (Massive or graded sandstones) is thought to represent sedimentation from high and low density turbidity currents in sand lobes. Facies D (Folded sandstones and siltstones) was formed from slumping in proximal lobe environments. Facies E (Wave-rippled sandstones) records wave reworking of sands supplied by turbidity currents above wave base level.The Lake Malanzán succession is formed by stacked turbidite sand lobe deposits. These lobes were probably formed in proximal lacustrine settings, most likely relatively high gradient slopes. Paleocurrents indicate a dominant direction from cratonic areas to the WSW. Although the overall sequence shows a regressive trend from basin fine-grained deposits to deltaic and braided fluvial facies, individual lobe packages lack of definite vertical trends in bed thickness and grain size. This fact suggests aggradation from multiple-point sources, rather than progradation from single-point sources. Sedimentologic and paleoecologic evidence indicate high depositional rate and sediment supply. Deposition within the lake was largely dominated by event sedimentation. Low diversity trace fossil assemblages of opportunistic invertebrates indicate recolonization of event beds under stressed conditions.Three stages of lake evolutionary history have been distinguished. The vertical replacement of braided fluvial deposits by basinal facies indicates high subsidence and a lacustrine transgressive episode. This flooding event was probably linked to a notable base level rise during postglacial times. The second evolutionary stage was typified by the formation of sand turbidite lobes from downslope mass-movements. Lake history culminates with the progradation of deltaic and braided fluvial systems     

Quantifying faulting and base level controls on syn‐rift sedimentation using stratigraphic architectures of coeval,adjacent Early‐Middle Pleistocene fan deltas in Lake Corinth,Greece

Quantification of allogenic controls in rift basin‐fills requires analysis of multiple depositional systems because of marked along‐strike changes in depositional architecture. Here, we compare two coeval Early‐Middle Pleistocene syn‐rift fan deltas that sit 6 km apart in the hangingwall of the Pirgaki‐Mamoussia Fault, along the southern margin of the Gulf of Corinth, Greece. The Selinous fan delta is located near the fault tip and the Kerinitis fan delta towards the fault centre. Selinous and Kerinitis have comparable overall aggradational stacking patterns. Selinous comprises 15 cyclic stratal units (ca. 25 m thick), whereas at Kerinitis 11 (ca. 60 m thick) are present. Eight facies associations are identified. Fluvial and shallow water facies dominate the major stratal units in the topset region, with shelfal fine‐grained facies constituting ca. 2 m thick intervals between major topset units and thick conglomeratic foresets building down‐dip. It is possible to quantify delta build times (Selinous: 615 kyr; Kerinitis: >450 kyr) and average subsidence and equivalent sedimentation rates (Selinous: 0.65 m/kyr; Kerinitis: >1.77 m/kyr). The presence of sequence boundaries at Selinous, but their absence at Kerinitis, enables sensitivity analysis of the most uncertain variables using a numerical model, ‘Syn‐Strat’, supported by an independent unit thickness extrapolation method. Our study has three broad outcomes: (a) the first estimate of lake level change amplitude in Lake Corinth for the Early‐Middle Pleistocene (10–15 m), which can aid regional palaeoclimate studies and inform broader climate‐system models; (b) demonstration of two complementary methods to quantify faulting and base level signals in the stratigraphic record—forward modelling with Syn‐Strat and a unit thickness extrapolation—which can be applied to other rift basin‐fills; and (c) a quantitative approach to the analysis of stacking patterns and key surfaces that could be applied to stratigraphic pinch‐out assessment and cross‐hole correlations in reservoir analysis.     

Shallow-marine sequences as the building blocks of stratigraphy: insights from numerical modelling

Two types of depositional sequences can be defined within the sequence stratigraphic framework: the parasequence and the high‐frequency sequence. Both sequences consist of stacked regressive and transgressive deposits. However, a parasequence forms under conditions of overall sea‐level rise, whereas a high‐frequency sequence forms as the sea level oscillates which results in typical forced regressive deposits during sea‐level fall. Both depositional sequences may develop over comparable temporal (10–100 kyr) and spatial (1–20 km wide and 1–40 m thick) scales. Numerical modelling is used to compare the architecture, preservation potential, internal volumes, bounding surfaces, condensed and expanded sections and facies assemblages of parasequences and high‐frequency sequences. Deposits originating from transgression are less pronounced than their regressive counterparts and consist of either preserved backbarrier deposits or shelf deposits. Shoreface deposits are not preserved during transgression. The second half of the paper evaluates in detail the preservation potential of backbarrier deposits and proposes a mechanism that explains the occurrence of both continuous and discontinuous barrier retreat in terms of varying rates of sea‐level rise and sediment supply. The key to this mechanism is the maximum washover capacity, which plays a part in both barrier shoreline retreat and backbarrier‐lagoonal shoreline retreat. If these two shorelines are not balanced, then the retreat of the coastal system as a whole is discontinuous and in time barrier overstep may take place.     

Response of the St. Joseph River to lake level changes during the last 12,000 years in the Lake Michigan basin

The water level of the Lake Michigan basin is currently 177 m above sea level. Around 9,800 ¹⁴C years B.P., the lake level in the Lake Michigan basin had dropped to its lowest level in prehistory, about 70 m above sea level. This low level (Lake Chippewa) had profound effects on the rivers flowing directly into the basin. Recent studies of the St. Joseph River indicate that the extreme low lake level rejuvenated the river, causing massive incision of up to 43 m in a valley no more than 1.6 km wide. The incision is seen 25 km upstream of the present shoreline. As lake level rose from the Chippewa low, the St. Joseph River lost competence and its estuary migrated back upstream. Floodplain and channel sediments partially refilled the recently excavated valley leaving a distinctly non-classical morphology of steep sides with a broad, flat bottom. The valley walls of the lower St. Joseph River are 12–18 m tall and borings reveal up to 30 m of infill sediment below the modern floodplain. About 3 × 10⁸ m³ of sediment was removed from the St. Joseph River valley during the Chippewa phase lowstand, a massive volume, some of which likely resides in a lowstand delta approximately 30 km off-shore in Lake Michigan. The active floodplain below Niles, Michigan, is inset into an upper terrace and delta graded to the Calumet level (189 m) of Lake Chicago. In the lower portion of the terrace stratigraphy a 1.5–2.0 m thick section of clast-supported gravel marks the entry of the main St. Joseph River drainage above South Bend, Indiana, into the Lake Michigan basin. This gravel layer represents the consolidation of drainage that probably occurred during final melting out of ice-marginal kettle chains allowing stream piracy to proceed between Niles and South Bend. It is unlikely that the St. Joseph River is palimpsest upon a bedrock valley. The landform it cuts across is a glaciofluvial-deltaic feature rather than a classic unsorted moraine that would drape over pre-glacial topography. This is the fifth in a series of ten papers published in this special issue of Journal of Paleolimnology. These papers were presented at the 47th Annual Meeting of the International Association for Great Lakes Research (2004), held at the University of Waterloo, Waterloo, Ontario, Canada. P.F. Karrow and C.F.M. Lewis were guest editors of this special issue.     

PUNCTUATED DRAINAGE OF AN ICE-DAMMED QUATERNARY LAKE IN SOUTHERN SOUTH AMERICA

A chain of east‐west elongated lakes fringes the eastern flank of the Andes Mountains between 39 and 52°S in southern Chile and Argentina. These deep lakes fill valleys left by the retreat of outlet glaciers of the Patagonian Icefield about 16000 years ago. Two of the lakes, Lago General Carrera/Buenos Aires and Lago Cochrane, represent the shrunken remnants of a much larger lake (called here the predecessor lake) that developed in a moat alongside the icefield. A series of seven stepped lacustrine braid deltas on the southeastern shore of Lago General Carrera/Buenos Aires mark the progressive, punctuated drainage of the predecessor lake. The deltas were formed by the build‐out of fans of sand and gravel into the lake. They have the form of Gilbert deltas, with virtually flat subaerial topset beds and steeply inclined subaqueous foreset beds. The exposed delta fronts are marked by a series of small terraces produced by wave erosion during falling lake levels. On either side of the fans are embayments with beaches formed by deposition in the lee of the pro‐grading fans. These embayments are filled with ridges of well‐sorted gravel deposited by wave‐driven long‐shore drift. About 13 000 years ago the regional drainage underwent a remarkable change. Rivers flowing eastwards into the Atlantic reversed their direction to flow westwards through the Andes Mountains and into the Pacific Ocean. The watershed moved some 200 km to the east. The cause of this change, and the falling lake levels, was melting of the Patagonian Icefield. As the ice retreated the icefield split into north and south components. This broke the ice dam that impounded the western end of the predecessor lake. The stepped, progressive fall of the lake level suggests that as the ice melted a series of lower overflow sills of hard rock were exposed. Each of the major episodes of sudden drop in lake level was followed by a long period of up to thousands of years when lake level remained constant. Terraces on the delta fronts indicate that each major episode of fall in the lake level was punctuated by up to ten minor episodes possibly representing major storm events or annual cycles of freezing and melting, blocking the overflow sills.     

Generic 3D interpolation of Holocene base‐level rise and provision of accommodation space,developed for the Netherlands coastal plain and infilled palaeovalleys

We present an interpolation model that describes Holocene groundwater level rise and the creation of accommodation space in 3D in the Rhine‐Meuse delta – the Netherlands. The model area (ca. 12 400 km²) covers two palaeovalleys of Late Pleistocene age (each 30 km wide) and the overlying Holocene deposits of the Rhine‐Meuse delta, the Holland coastal plain, and the Zuiderzee former lagoon. Water table rise is modelled from 10 800 to 1000 cal. BP, making use of age‐depth relations based on 384 basal peat index points, and producing output in the form of stacked palaeo groundwater surfaces, groundwater age‐depth curves, and voxel sets. These products allow to resolve (i) regional change and variations of inland water table slopes, (ii) spatial differences in the timing and pacing of transgression, and (iii) analysis of interplay of coastal, fluvial and subsidence controls on the provision of accommodation space. The interpolation model is a multi‐parameter trend function, to which a 3D‐kriging procedure of the residuals is added. This split design deploys a generic approach for modelling provision of accommodation space in deltas and coastal lowlands, aiming to work both in areas of intermediate data availability and in the most data‐rich environments. Major provision of accommodation space occurred from 8500 cal BP onwards, but a different evolution occurred in each of the two palaeovalleys. In the northern valley, creation of accommodation space began to stall at 7500 cal BP, while in the southern valley provision of new accommodation space in considerable quantities continued longer. The latter is due to the floodplain gradient that was maintained by the Rhine, which distinguishes the fluvial deltaic environment from the rest of the back‐barrier coastal plain. The interpolation results allow advanced mapping and investigation of apparent spatial differences in Holocene aggradation in larger coastal sedimentary systems. Furthermore, they provide a means to generate first‐order age information with centennial precision for 3D geological subsurface models of Holocene deltas and valley fills. As such, the interpolation is of use in studies into past and present land subsidence and into low land sedimentation.     

Microbial‐dominated carbonate platforms during the Ladinian rifting: sequence stratigraphy and evolution of accommodation in a fault‐controlled setting (Catalan Coastal Ranges,NE Spain)

The Upper Muschelkalk sedimentary record constitutes a major transgressive pulse of north‐eastern Iberia during the Ladinian. This record is arranged in two transgressive–regressive (T–R) sequences formed by two stepped microbial‐dominated carbonate ramp systems where accommodation was mainly controlled by extensional faults. This study seeks to gain new insights into how the evolution of syn‐rift subsidence controls the creation of accommodation space, the depositional styles and, especially, the palaeogeographical domains where specific microbialites developed (thrombolites and stromatolites). Thrombolite bodies (at least 40 m thick) display two types of architecture, biostromal and mud‐mounded and stromatolite bodies (at least 7 m thick) consist of tabular and domed, head‐shaped morphologies. Domed and mounded forms are usually developed during stages of increasing accommodation rates, low‐to flat‐nelief forms tend to grow in association with periods of low accommodation rates. A sea‐level fall of at least 50 m occurred at the end of the Early Ladinian leaving the platform subaerially exposed. As a result, a prominent karst with significant erosional incisions and profuse collapse breccia fillings was formed in the inner and middle ramp settings. The resultant subaerial unconformity bounds T–R sequences 1 and 2. Subsidence curves display two stages of rapid/decelerated total subsidence, constituting two discrete rift/post‐rift pulses in the large Triassic rifting period: (i) Buntsandstein – Middle Muschelkalk, and (ii) Late Muschelkalk – Imon Formation (Rhaetian). The second pulse is characterized by a rapid syn‐rift subsidence during the Late Muschelkalk, and a decelerated post‐rift subsidence throughout the deposition of Keuper facies and Imon Formation. The Late Muschelkalk rapid syn‐rift pulse of total subsidence produces gains in accommodation, which controls the development of the stromatolites and thrombolites (biostromes and mud‐mounds).     

Unroofing the core of the central Andean fold‐thrust belt during focused late Miocene exhumation: evidence from the Tipuani‐Mapiri wedge‐top basin,Bolivia

As the highest part of the central Andean fold‐thrust belt, the Eastern Cordillera defines an orographic barrier dividing the Altiplano hinterland from the South American foreland. Although the Eastern Cordillera influences the climatic and geomorphic evolution of the central Andes, the interplay among tectonics, climate and erosion remains unclear. We investigate these relationships through analyses of the depositional systems, sediment provenance and ⁴⁰Ar/³⁹Ar geochronology of the upper Miocene Cangalli Formation exposed in the Tipuani‐Mapiri basin (15–16°S) along the boundary of the Eastern Cordillera and Interandean Zone in Bolivia. Results indicate that coarse‐grained nonmarine sediments accumulated in a wedge‐top basin upon a palaeotopographic surface deeply incised into deformed Palaeozoic rocks. Seven lithofacies and three lithofacies associations reflect deposition by high‐energy braided river systems, with stratigraphic relationships revealing significant (～500 m) palaeorelief. Palaeocurrents and compositional provenance data link sediment accumulation to pronounced late Miocene erosion of the deepest levels of the Eastern Cordillera. ⁴⁰Ar/³⁹Ar ages of interbedded tuffs suggest that sedimentation along the Eastern Cordillera–Interandean Zone boundary was ongoing by 9.2 Ma and continued until at least ～7.4 Ma. Limited deformation of subhorizontal basin fill, in comparison with folded and faulted rocks of the unconformably underlying Palaeozoic section, implies that the thrust front had advanced into the Subandean Zone by the 11–9 Ma onset of basin filling. Documented rapid exhumation of the Eastern Cordillera from ～11 Ma onward was decoupled from upper‐crustal shortening and coeval with sedimentation in the Tipuani‐Mapiri basin, suggesting climate change (enhanced precipitation) or lower crustal and mantle processes (stacking of basement thrust sheets or removal of mantle lithosphere) as possible controls on late Cenozoic erosion and wedge‐top accumulation. Regardless of the precise trigger, we propose that an abruptly increased supply of wedge‐top sediment produced an additional sedimentary load that helped promote late Miocene advance of the central Andean thrust front in the Subandean Zone.     

Tectono‐sedimentary evolution of the Plio‐Pleistocene Corinth rift,Greece

The onshore central Corinth rift contains a syn‐rift succession >3 km thick deposited in 5–15 km‐wide tilt blocks, all now inactive, uplifted and deeply incised. This part of the rift records upward deepening from fluviatile to lake‐margin conditions and finally to sub‐lacustrine turbidite channel and lobe complexes, and deep‐water lacustrine conditions (Lake Corinth) were established over most of the rift by 3.6 Ma. This succession represents the first of two phases of rift development – Rift 1 from 5.0–3.6 to 2.2–1.8 Ma and Rift 2 from 2.2–1.8 Ma to present. Rift 1 developed as a 30 km‐wide zone of distributed normal faulting. The lake was fed by four major N‐ to NE‐flowing antecedent drainages along the southern rift flank. These sourced an axial fluvial system, Gilbert fan deltas and deep lacustrine turbidite channel and lobe complexes. The onset of Rift 2 and abandonment of Rift 1 involved a 30 km northward shift in the locus of rifting. In the west, giant Gilbert deltas built into a deepening lake depocentre in the hanging wall of the newly developing southern border fault system. Footwall and regional uplift progressively destroyed Lake Corinth in the central and eastern parts of the rift, producing a staircase of deltaic and, following drainage reversal, shallow marine terraces descending from >1000 m to present‐day sea level. The growth, linkage and death of normal faults during the two phases of rifting are interpreted to reflect self‐organization and strain localization along co‐linear border faults. In the west, interaction with the Patras rift occurred along the major Patras dextral strike‐slip fault. This led to enhanced migration of fault activity, uplift and incision of some early Rift 2 fan deltas, and opening of the Rion Straits at ca. 400–600 ka. The landscape and stratigraphic evolution of the rift was strongly influenced by regional palaeotopographic variations and local antecedent drainage, both inherited from the Hellenide fold and thrust belt.