Initiation and early evolution of a forearc basin: Georgia Basin,Canada

The lower Nanaimo Group was deposited in the (forearc) Georgia Basin, Canada and records the basin's initiation and early depositional evolution. Nanaimo Group strata are subdivided into 11 lithostratigraphic units, which are identified based on lithology, paleontology, texture and position relative to both the basal nonconformity and to each other. Significant topography on the basal nonconformity, however, has resulted in assignment of lithostratigraphic units that are not time correlative, and hence, cannot reliably be used to accurately reconstruct basin evolution. Herein, we present a sequence stratigraphic framework for lower Nanaimo Group strata in the Comox Sub-Basin (northern Georgia Basin) that integrates both facies analysis and maximum depositional ages (MDAs) derived from detrital zircon. This stratigraphic framework is used to define significant sub-basin-wide surfaces that bound depositional units and record the evolution of the basin during its early stages of development. Seven distinct depositional phases are identified in the lower 700 m of the lower Nanaimo Group. Depositional phases are separated by marine flooding surfaces, regressive surfaces, or disconformities. The overall stratigraphy reflects net transgression manifested as an upwards transition from braided fluvial conglomerates to marine mudstones. Transgression was interrupted by periods of shoreline progradation, and both facies analysis and MDAs reveal a disconformity in the lowermost part of the Nanaimo Group in the Comox Sub-Basin. Stratigraphic reconstruction of the Comox Sub-Basin reveals two dominant depocenters (along depositional strike) for coarse clastics (sandstones and conglomerates) during early development of the Georgia Basin. The development and position of these depocenters is attributed to subduction/tectonism driving both subsidence in the north-northwest and uplift in the central Comox Sub-Basin. Our work confirms that in its earliest stages of development, the Georgia Basin evolved from an underfilled, ridged forearc basin that experienced slow and stepwise drowning to a shoal-water ridged forearc basin that experienced rapid subsidence. We also propose that the Georgia Basin is a reasonable analogue for ridged forearc basins globally, as many ridged forearcs record similar depositional histories during their early evolution.     

Quantitative analysis of Miocene to Recent forearc basin evolution along the Colombian convergent margin

The Colombian accretionary complex forms the active convergent margin of the North Andes block of South America beneath which the east Panama Basin of the Nazca plate is subducted at a rate of 50–64 km Myr^?1. Multichannel seismic reflection data, collected as part of RRS Charles Darwin cruise CD40, image a series of well-developed forearc basins along the length of the margin, bounded on their oceanward side by an active accretionary complex and on their landward side by oceanward-dipping continental basement. Sedimentary sequences within the forearc basins indicate successive landward migration of the basin depocentre as the structural high bounding its oceanward edge is forced upward and landward by continued growth of the accretionary complex. Uplift beneath the oceanward side of the basins has resulted in progressive landward rotation of the older sedimentary sequences. Prominent seismic reflectors across the basins show a complex onlap–offlap relationship between successive sequences that reflects the interplay between tectonic uplift, sediment supply, differential sediment compaction and basement subsidence due to loading. A numerical model has been devised to investigate how Miocene to Recent forearc basin stratigraphy is controlled by progressive growth of the accretionary complex resulting in elevation of the outer-arc high and landward motion of the rear of the complex up the seaward-dipping backstop formed by the leading edge of the continental lithosphere. The effects of sediment accretion are modelled by treating the accretionary complex as a doubly vergent, noncohesive Coulomb wedge, where forces exerted by the proto- and retro-wedges must be balanced for the system to be in equilibrium. The model generates synthetic basin-fill architecture over a series of steps, each of which represents the deposition of individual sedimentary sequences and their subsequent deformation due to wedge growth. The model accounts for differential sediment compaction and the flexural response of the underlying lithosphere to changes in load distribution over time. Forearc basin evolution is simulated for various rates of sediment supply to the forearc and accretionary complex growth until the synthetic basin-fill geometry matches the observed geometry. The model enables either the rate of accretionary wedge growth or the rate of sediment supply to the forearc basin to be established. The technique is generally applicable to those convergent margins with forearc basins that have developed between an actively accreting wedge and a seaward-dipping backstop. Other examples include Peru, S. Chile, Sumatra and Barbados.     

Sediment‐budget modelling of multi‐sourced basin fills: application to recent deposits of the western Adriatic mud wedge (Italy)

We analysed modern mass‐accumulation patterns on the western Adriatic mud wedge (Italy), an elongated belt of shelf mud formed by coalesced prodeltas of the Adige, Po, and Apennine rivers, as part of an integrated strategy aimed at producing a quantitative sediment budget model for muddy continental shelves sourced by multiple compositionally distinct fluvial systems. Sediment provenance and source‐specific accumulation rates of surface sediments were quantified by combining results of grain‐size analysis and geochemical analysis of specific size fractions with bulk mass accumulation rates. Statistical classification algorithms adapted to compositional data were used to partition the total (geochemical) variation of sediment properties into size‐related and provenance‐specific factors. We identified geochemically distinct fluvial end‐member sediment types in two different grain‐size fractions, which were grouped into sediments derived from the Apennine rivers, and sediments derived from the Po and Adige rivers. Compositional fingerprints (end‐member compositions) of each source area were estimated by taking into account relative rates of fluvial sediment supply from rivers as predicted by numerical modelling. The end members allow us to explain geochemical compositional variation of mud‐wedge surface sediments in terms of provenance and size‐selective dispersal, and map mass accumulation rates of sediments from individual source areas (grain size<63 μm), as well as bulk sand accumulation rates (grain size>63 μm) across the western Adriatic mud wedge. The source‐specific rates of fine‐grained sediment supply derived from geostatistical estimates of mass‐accumulation rates were used to calibrate the numerical model of sediment supply to present‐day conditions.     

Impact of structural and autocyclic basin-floor topography on the depositional evolution of the deep-water Valparaiso forearc basin, central Chile

The Valparaiso Basin constitutes a unique and prominent deep‐water forearc basin underlying a 40‐km by 60‐km mid‐slope terrace at 2.5‐km water depth on the central Chile margin. Seismic‐reflection data, collected as part of the CONDOR investigation, image a 3–3.5‐km thick sediment succession that fills a smoothly sagged, margin‐parallel, elongated trough at the base of the upper slope. In response to underthrusting of the Juan Fernández Ridge on the Nazca plate, the basin fill is increasingly deformed in the seaward direction above seaward‐vergent outer forearc compressional highs. Syn‐depositional growth of a large margin‐parallel monoclinal high in conjunction with sagging of the inner trough of the basin created stratal geometries similar to those observed in forearc basins bordered by large accretionary prisms. Margin‐parallel compressional ridges diverted turbidity currents along the basin axis and exerted a direct control on sediment depositional processes. As structural depressions became buried, transverse input from point sources on the adjacent upper slope formed complex fan systems with sediment waves characterising the overbank environment, common on many Pleistocene turbidite systems. Mass failure as a result of local topographic inversion formed a prominent mass‐flow deposit, and ultimately resulted in canyon formation and hence a new focused point source feeding the basin. The Valparaiso Basin is presently filled to the spill point of the outer forearc highs, causing headward erosion of incipient canyons into the basin fill and allowing bypass of sediment to the Chile Trench. Age estimates that are constrained by subduction‐related syn‐depositional deformation of the upper 700–800 m of the basin fill suggest that glacio‐eustatic sea‐level lowstands, in conjunction with accelerated denudation rates, within the past 350 ka may have contributed to the increase in simultaneously active point sources along the upper slope as well as an increased complexity of proximal depositional facies.     

Timing of Xunhua and Guide basin development and growth of the northeastern Tibetan Plateau,China

The Xunhua, Guide and Tongren intermontane basin system in the NE Tibetan Plateau, situated near the Xining basin to the N and the Linxia basin to the E, is bounded by thrust fault‐controlled ranges. These include to the N, the Riyue Shan, Laji Shan and Jishi Shan ranges, and to the S the northern West Qinling Shan (NWQ). An integrated study of the structural geology, sedimentology and provenance of the Cenozoic Xunhua and Guide basins provides a detailed record of the growth of the NE Tibetan Plateau since the early Eocene. The Xining Group (ca. 52–21 Ma) is interpreted as consisting of unified foreland basin deposits which were controlled by the bounding thrust belt of the NWQ. The Xunhua, Guide and Xining subbasins were interconnected prior to later uplift and damming by the Laji Shan and Jishi Shan ranges. Their sediment source, the NWQ, is constrained by strong unidirectional paleocurrent trends towards the N, a northward fining lithology, distinct and recognizable clast types and detrital zircon ages. Collectively, formation of this mountain–basin system indicates that the Tibetan Plateau expanded into the NWQ at a time roughly coinciding with Eocene to earliest Miocene continental collision between India and Eurasia. The Guide Group (ca. 21–1.8 Ma) is inferred to have been deposited in the separate Xunhua, Guide and Tongren broken foreland basins. Each basin was filled by locally sourced alluvial fans, braided streams and deltaic‐lacustrine systems. Structural, paleogeographic, paleocurrent and provenance data indicate that thrust faulting in the NWQ stepped northward to the Laji Shan from ca. 21 to 16 Ma. This northward shift was accompanied by E–W shortening related to nearly N–S‐striking thrust faulting in Jishi Shan after 11–13 Ma. A lower Pleistocene conglomerate (1.8–1.7 Ma) was deposited by a through‐flowing river system in the overfilled and connected Guide and Xunhua basins following the termination of thrust activity. All of the basin–mountain zones developed along the Tibetan Plateau's NE margin since Indian–Tibetan continental collision may have been driven by collision‐induced basal drag of old slab remnants in the manner of N‐dipping and flat‐slab subduction, and their subsequent sinking into the deep mantle.     

Sedimentology,provenance and geochronology of the upper Cretaceous–lower Eocene western Xigaze forearc basin,southern Tibet

Located on the southern margin of the Lhasa terrane in southern Tibet, the Xigaze forearc basin records Cretaceous to lower Eocene sedimentation along the southern margin of Asia, prior to and during the initial stages of continental collision with the Tethyan Himalaya in the Early Eocene. We present new measured stratigraphic sections, totalling 4.5 km stratigraphic thickness, from a 60 km E–W segment of the western portion of the Xigaze forearc basin, northeast of the Lopu Kangri Range (29.8007° N, 84.91827° E). In addition, we apply U–Pb detrital zircon geochronology to constrain the provenance and maximum depositional ages of investigated strata. Stratigraphic ages range between ca. 88 and ca. 54 Ma and sedimentary facies indicate a shoaling‐upward trend from deep‐marine turbidites to fluvial deposits. Depositional environments of coeval Cretaceous strata along strike include deep‐marine distal turbidites, slope‐apron debris‐flow deposits and marginal marine carbonates. This along‐strike variability in facies suggests an irregular paleogeography of the Asian margin prior to collision. Paleocene–Eocene strata are composed of shallow marine carbonates with abundant foraminifera such as Nummulites‐Discocyclina and Miscellanea‐Daviesina and transition into fluvial deposits dated at ca. 54 Ma. Sandstone modal analyses, conglomerate clast compositions and detrital zircon U–Pb geochronology indicate that forearc detritus in this region was derived solely from the Gangdese magmatic arc to the north. In addition, U–Pb detrital zircon age spectra within the upper Xigaze forearc stratigraphy are similar to those from Eocene foreland basin strata south of the Indus‐Yarlung suture near Sangdanlin, suggesting that the Xigaze forearc was a possible source of Sangdanlin detritus by ca. 55 Ma. We propose a model in which the Xigaze forearc prograded south over the accretionary prism and onto the advancing Tethyan Himalayan passive margin between 58 and 54 Ma, during late stage evolution of the forearc basin and the beginning of collision with the Tethyan Himalaya. The lack of documented forearc strata younger than ca. 51 Ma suggests that sedimentation in the forearc basin ceased at this time owing to uplift resulting from continued continental collision.     

The Astrakhan Arch of the Pricaspian basin: Geothermal analysis and modelling

The Astrakhan Arch (ASAR) region contains one of the largest sub‐salt carbonate structures of the Pricaspian salt basin (located to the northwest of the Caspian Sea), where prospects for hydrocarbon generation and accumulation in the Devonian to Carboniferous deposits are considered to be high. We evaluate the regional vertical temperature gradient within stratigraphic units based on the analysis of 34 boreholes drilled in the region. To show that the thermal gradient is altered in the vicinity of salt diapirs, we study measured temperatures in six deep boreholes. We develop a three‐dimensional geothermal model of the ASAR region constrained by temperature measurements, seismic stratigraphic and lithological data. The temperatures of the sub‐salt sediments predicted by the geothermal model range from about 100 °C to 200 °C and are consistent with the temperatures obtained from the analysis of vitrinite reflectivity and from previous two‐dimensional geothermal models. Temperature anomalies are positive in the uppermost portions of salt diapirs as well as within the salt‐withdrawal basins at the depth of 3.5 km depth and are negative beneath the diapirs. Two areas of positive temperature anomalies in the sub‐salt sediments are likely to be associated with the deep withdrawal basins above and with the general uplift of salt/sub‐salt interface in the southern part of the study region. This implies an enhancement of thermal maturity of any organically rich source rocks within these areas. The surface heat flux in the model varies laterally from about 40 to 55 mW m^?2. These variations in the heat flux are likely to be associated with structural heterogeneities of the sedimentary rocks and with the presence of salt diapirs. The results of our modelling support the hypothesis of oil and gas condensate generation in the Upper Carboniferous to Middle Devonian sediments of the ASAR region.     

汾河流域地貌发育对构造运动和气候变化的响应

Tectonic movements and climate changes are two main controllers on the development of landfrorm.In order to reconstruct the history of the evolution of the landform in the Fenhe drainage basin during middle-late Quaternary comprehensively,this paper has provided a variety of geomorphological and geologic evidences to discuss how tetonic movements and climate changes worked together to influence the landform processes,According to the features of the lacustrine and alluvial terraces in this drainage basin,it is deduced that it was the three tectonic uplifts that resulted in the three great lake-regressions with an extent of about 40-60 m and the formation of the three lacustrine terraces.The times when the tectonic uplifts took place are 0.76 MaBP,0.55 MaBP and 0.13 MaBP respectively,synchronous with the formation of paleosol units S8,S5 and S1 respectively.During the intervals between two tectonic uplifts when tectonic movement was very weak ,climate changes played a major role in the evolution of the paleolakes and caused frequent fluctuations of lake levels.The changes of the features of lacustrine sediment in the grabens show the extent of such fluctuations of lake level is about 2-3m.     

Responses of landform development to tectonic movements and climate change during Quaternary in Fenhe drainage basin

Tectonic movements and climate changes are two main controllers on the development of landform. In order to reconstruct the history of the evolution of the landform in the Fenhe drainage basin during middle-late Quaternary comprehensively, this paper has provided a variety of geomorphological and geologic evidences to discuss how tectonic movements and climate changes worked together to influence the landform processes. According to the features of the lacustrine and alluvial terraces in this drainage basin, it is deduced that it was the three tectonic uplifts that resulted in the three great lake-regressions with an extent of about 40-60 m and the formation of the three lacustrine terraces. The times when the tectonic uplifts took place are 0.76 MaBP, 0.55 MaBP and 0.13 MaBP respectively, synchronous with the formation of paleosol units S8, S5 and S1 respectively. During the intervals between two tectonic uplifts when tectonic movement was very weak, climate changes played a major role in the evolution of the paleolakes and caused frequent fluctuations of lake levels. The changes of the features of lacustrine sediment in the grabens show the extent of such fluctuations of lake level is about 2-3 m.     

Neogene supradetachment basin development on Crete (Greece) during exhumation of the South Aegean core complex

Tertiary extension in the Aegean region has led to extensional detachment faulting, along which metamorphic core complexes were exhumed, among which is the Early to Middle Miocene South Aegean core complex. This paper focuses on the supradetachment basin developed during the final stages of exhumation of the South Aegean core complex along the Cretan detachment, plus the Late Miocene to Pliocene basin development and palaeogeography associated with the southward motion of Crete during the opening of the Aegean arc. For the latter purpose, the sedimentary and palaeobathymetric evolutions of a large number of Middle Miocene to Late Pliocene sequences exposed on Crete, Gavdos and Koufonisi were studied. The supradetachment basin development of Crete is characterised by a break‐up of the hanging wall of the Cretan detachment into extensional klippen and subsequent migration of laterally coexisting sedimentary systems, and finally the deformation of the exhumed core complex by processes related to the opening of the Aegean arc. Hence, three main tectonic phases are recognised: (1) Early to Middle Miocene N–S extension formed during the Cretan detachment, exhumed in the South Aegean core complex. The Cretan detachment remained active until 11–10 Ma, based on the oldest sediments that unconformably overlie the metamorphic rocks. Successions older than 11–10 Ma unconformably overlie only the hanging wall of the Cretan detachment, and do not contain fragments of the footwall rocks; they therefore predate the oldest exposure of the metamorphic rocks of the footwall. The hanging wall rocks and Middle Miocene sediments form isolated blocks on top of the exhumed metamorphic rocks, which are interpreted as extensional klippen. (2) From approximately 10 Ma onward, southward migration of the area that presently covers Crete was accompanied by E–W extension, and the opening of the Sea of Crete to the north. Contemporaneously, large folds with WNW–ESE striking, NNE dipping axial planes developed, possibly in response to sinistral transpression. (3) During the Pliocene, Crete emerged and tilted to the NNW, probably as a result of left‐lateral transpression in the Hellenic fore‐arc, induced by the collision with the African promontory.     

Late Cretaceous basement foundering of the Rosario embayment, Peninsular Ranges forearc basin: backstripping of the El Gallo Formation, Baja California Norte, Mexico

Excellent exposure, well-controlled palaeobathymetry, and tightly-spaced, high-precision radiometric age control in the El Gallo Fm. permit rigorous quantitative analysis. Backstripping of these proximal nonmarine, forearc basin deposits reveals that, during the Late Cretaceous, the Rosario embayment of the Peninsular Ranges forearc was undergoing an episode of rapid tectonic subsidence. This subsidence had several marked effects on the sedimentology of the Rosario embayment: formation of a broad alluvial plain consisting of coarse-grained clastics; rapid (∼ 600 m Myr^-1) aggradation of sediments; and a retrogradational succession of facies, capped by a marine transgression, as deposition failed to keep pace with eustatic rise and subsidence.
Long-term sedimentation is driven by some combination of two allocyclic mechanisms: tectonic subsidence and eustatic sea-level rise. In order to evaluate which force predominated during deposition of the El Gallo Fm., the processes of sedimentation, compaction, and isostasy are evaluated through the interval in question. A sensitivity analysis is performed, in which the maximum tectonic and maximum eustatic contributions are estimated, along with the best-fit model. These results are qualitatively the same: tectonic subsidence was the major driving force of sedimentation in the Rosario embayment in late Campanian time. Regional sedimentological similarities suggest that this tectonic subsidence may have characterized the Peninsular Ranges forearc margin at this time, reflecting an episode of active down-faulting during the Late Cretaceous.     

Formation and fragmentation of a late Miocene supradetachment basin in central Crete: implications for exhumation mechanisms of high‐pressure rocks in the Aegean forearc

We present a new lithostratigraphy and chronology for the Miocene on central Crete, in the Aegean forearc. Continuous sedimentation started at ～10.8 Ma in the E–W trending fluvio‐lacustrine Viannos Basin, formed on the hangingwall of the Cretan detachment, which separates high‐pressure (HP) metamorphic rocks from very low‐grade rocks in its hangingwall. Olistostromes including olistoliths deposited shortly before the Viannos Basin submerged into the marine Skinias Basin between 10.4 and 10.3 Ma testifies to significant nearby uplift. Uplift of the Skinias Basin between 9.7 and 9.6 Ma, followed by fragmentation along N–S and E–W striking normal faults, marks the onset of E–W arc‐parallel stretching superimposed on N–S regional Aegean extension. This process continued between 9.6 and 7.36 Ma, as manifested by tilting and subsidence of fault blocks with subsidence events centred at 9.6, 8.8, and 8.2 Ma. Wholesale subsidence of Crete occurred from 7.36 Ma until ～5 Ma, followed by Pliocene uplift and emergence. Subsidence of the Viannos Basin from 10.8 to 10.4 Ma was governed by motion along the Cretan detachment. Regional uplift at ～10.4 Ma, followed by the first reworking of HP rocks (10.4–10.3 Ma) is related to the opening and subsequent isostatic uplift of extensional windows exposing HP rocks. Activity of the Cretan detachment ceased sometime between formation of extensional windows around 10.4 Ma, and high‐angle normal faulting cross‐cutting the detachment at 9.6 Ma. The bulk of exhumation of the Cretan HP‐LT metamorphic rocks occurred between 24 and 12 Ma, before basin subsidence, and was associated with extreme thinning of the hangingwall (by factor ～10), in line with earlier inferences that the Cretan detachment can only explain a minor part of total exhumation. Previously proposed models of buyoant rise of the Cretan HP rocks along the subducting African slab provide an explanation for extension without basin subsidence.     

Silica diagenesis in Cenozoic mudstones of the North Viking Graben: physical properties and basin modelling

Silica diagenesis can significantly change physical properties of the host strata and release large volumes of water. Predicting these changes and their timing is essential to understanding compaction, fluid flow and rock deformation in sedimentary basins. In this paper, the influence of silica diagenesis (opal‐A/CT transformation) on physical properties is determined, the sediment volume affected by these changes is mapped, and a new technique to model silica diagenesis is introduced. A petrophysical analysis of 16 exploration wells shows that the opal‐A/CT transformation leads to a porosity reduction of c.20% (from 49 to 29%) in Cenozoic mudstones of the North Viking Graben. Using three‐dimensional seismic reflection data, it is shown that the c.50 m thick opal‐A/CT transformation zone covers an area of >1500 km², equating to a minimum volume of 75 km³. The spatial and temporal evolution of opal‐A/CT transformation is simulated using an innovative basin modelling approach, the results of which indicate that the transformation started around Middle‐to‐Late Eocene times and then migrated upwards until it gradually fossilised between the Miocene and present. These findings are important, as they help understanding how these sediments compact and when fluids are released by diagenesis.     

On the importance of rock thermal conductivity and heat flow density in basin and petroleum system modelling

Basin and petroleum systems are routinely modelled to provide qualitative and quantitative assessments of a hydrocarbon play. The importance of the rock thermal properties and heat flow density in thermal modelling the history of a basin are well-known, but little attention is paid to assumptions of the thermal conductivity, present-day heat flow density and thermal history of basins. Assumed values are often far from measured values when data are available to check parameters, and effective thermal conductivity models prescribed in many basin simulators require improvement. The reconstructed thermal history is often justified by a successful calibration to present-day temperature and vitrinite reflectance data. However, a successful calibration does not guarantee that the reconstruction history is correct. In this paper, we describe the pitfalls in setting the thermal conductivity and heat flow density in basin models and the typical uncertainties in these parameters, and we estimate the consequences by means of a one-dimensional model of the super-deep Tyumen SG-6 well area that benefits from large amounts of reliable input and calibration data. The results show that the entire approach to present-day heat flow evaluations needs to be reassessed. Unreliable heat flow density data along with a lack of measurements of rock thermal properties of cores can undermine the quality of basin and petroleum system modelling.     

Mid-Palaeocene palaeogeography of the eastern North Sea basin: integrating geological evidence and 3D geodynamic modelling

The Mid‐Palaeocene palaeogeography of Denmark and the surrounding areas have been reconstructed on the basis of published geological data integrated with 3D geodynamic modelling. The use of numerical modelling enables quantitative testing of scenarios based on geological input alone and thus helps constrain likely palaeo‐water depths in areas where the geological data are inconclusive or incomplete. The interpretation of large‐scale erosional valleys and small‐scale circular depressions at the Mid‐Palaeocene Top Chalk surface in the Norwegian–Danish basin as either submarine or subaerial features is enigmatic and has strong implications for palaeogeographical reconstructions of the eastern North Sea basin. A 3D thermo‐mechanical model is employed in order to constrain the likely palaeo‐water depths of the eastern North Sea basin during the Palaeocene. The model treats the lithosphere as an elasto‐visco‐plastic continuum and models the lithospheric response to the regional stress field and thermal structure. The model includes the effects of sea‐level change, sedimentation and erosion, from the Mid Cretaceous to the present. Modelling results reproduce to first order geological data such as present sediment isopachs and palaeo‐water depths. It is concluded that the Mid Palaeocene water depths in the Norwegian–Danish basin were about 250 m. The erosional valleys and circular depressions at the top of the Upper Cretaceous‐Danian Chalk Group are thus interpreted to have formed in relatively deep water rather than due to subaerial exposure. Likely interpretations of the structures are therefore submarine valleys and pockmarks.     

Three-dimensional analogue modelling of an alluvial basin margin affected by hydrological cycles: processes and resulting depositional sequences

The dynamics between sediment erosion and accumulation at an alluvial basin margin affected by changes in the surface hydrology are explored using scaled analogue models produced in a flume. The presented results differ from previous counterparts in that accumulation or erosion has not been forced at a spreading outlet, but occurred at a slope change produced by previously accumulated sediment. Cyclical upstream incision produced by increased stream discharge generated incised valleys, and these were subsequently filled by sediment carried by less efficient streams generated during the low discharge period. High resolution mapping using 2.5 mm contour maps allowed the study of sediment accumulation and terrain modelling. The results of three selected experiments are analysed. The only variable explored was discharge. The basin margin was simulated by a ramp inserted in a low sloping flume, consisting of two segments of different slopes selected to emulate high and low efficiency flume fans produced elsewhere. Water and fine‐medium sand entered the ramp along a narrow (0.1 m) channel and flow expanded but without occupying the complete 1.2 m flume width. Flows were highly concentrated and noncohesive. Fan‐like accumulation (slope: 0.11) began during low discharge (LD) periods at the ramp slope break, and proceeded upstream, onlapping quickly at first, but shifting to mostly progradation at the end of the period. High discharges (HD) usually generated two or three incised channels at the beginning of the period, but one of them prevailed and rapidly eroded parts of the LD fan and moved the sediment to a more distal low‐sloping fan (slope: 0.045). Both LD and HD fans passed downstream into a system of small parallel channels resembling a braided alluvial plain ending in sediment lobes. The mapping of the accumulated sediment during the various periods allowed calculation of sediment budgets for the entire flume. The stratal architecture of the deposits was investigated along five parallel trenches cut after experiment termination. The regression analysis of depositional profiles at fan‐like features (expanding flow) and at braided plains (parallel flow) indicated that these fan‐like systems are linear and dependent on applied discharge, while the latter showed an exponential decrease of slope downstream, with a starting value set up by the fan slope. Two main types of stratigraphic units were generated, the LDST and HDST (system tracts). The LDST has a nonerosive base over ‘bedrock’ and the previous HDST, filling proximal erosional topography and prograding as well, generating an onlap–downlap array. Its geometry is highly variable and dependent on pre‐existing topography. The HDST base is an important erosive surface comparable to sequence boundaries. However, there are places without erosion due to a marginal position with respect to the main stream. Indeed, the results suggest that the three‐dimensional variability of erosion and depositional processes might produce very different architectures along the same basin margin.     

Numerical modelling of Quaternary terrace staircase formation in the Ebro foreland basin,southern Pyrenees,NE Iberia

The southern foreland basin of the Pyrenees (Ebro basin) is an exorheic drainage basin since Late Miocene times. Remnants of an early exorheic Ebro drainage system are not preserved, but morphology provides evidence for the Pliocene–Quaternary drainage development. The incision history of the Ebro system is denoted by (i) extensive, low gradient pedimentation surfaces which are associated with the denudation of the southern Pyrenean piedmont around the Pliocene–Quaternary transition and (ii) deeply entrenched Quaternary river valleys. Presumably since the Middle Pleistocene fluvial incision intensified involving the formation of extensive terrace staircase in the Ebro basin. Terrace exposure dating in major Ebro tributary rivers indicates climate‐triggered terrace formation in response to glacial–interglacial climate and glacier fluctuations in the Pyrenean headwaters. The overall (semi)parallel longitudinal terrace profiles argue for progressive base level lowering for the whole Ebro drainage network. The landscape evolution model, TISC, is used to evaluate climatic, tectonic and base level scenarios for terrace staircase formation in the Ebro drainage system. Model simulations are compared with morpho‐climatic, tectonic and chronologic data. Results show that climatic fluctuations cause terrace formation, but the incision magnitudes and convergent terrace profiles predicted by this climate model scenario are not consistent with the (semi)parallel terraces in the Ebro basin. A model including previous (late Pliocene) uplift of the lower Ebro basin results in rapid base‐level lowering and erosion along the drainage network, small late stage incision magnitudes and terrace convergence, which are not in agreement with observations. Instead, continuous Quaternary uplift of both the Pyrenees and the Ebro foreland basin triggers (semi)parallel terrace staircase formation in southern Pyrenean tributary rivers in consistency with the observed longitudinal terrace profiles and Middle–Late Pleistocene incision magnitudes. Forward model simulations indicate that the present Ebro drainage system is actively incising, providing further evidence for uplift.