Sink‐to‐source: Influence of offshore dynamics on upstream processes

When we model fluvial sedimentation and the resultant alluvial stratigraphy, we typically focus on the effects of local parameters (e.g., sediment flux, water discharge, grain size) and the effects of regional changes in boundary conditions applied in the source region (i.e., climate, tectonics) and at the shoreline (i.e., sea level). In recent years this viewpoint has been codified into the “source‐to‐sink” paradigm, wherein major shifts in sediment flux, grain‐size fining trends, channel‐stacking patterns, floodplain deposition and larger stratigraphic systems tracts are interpreted in terms of (1) tectonic and climatic signals originating in the hinterland that propagate downstream; and (2) eustatic fluctuation, which affects the position of the shoreline and dictates the generation of accommodation. Within this paradigm, eustasy represents the sole means by which downstream processes may affect terrestrial depositional systems. Here, we detail three experimental cases in which coastal rivers are strongly influenced by offshore and slope transport systems via the clinoform geometries typical of prograding sedimentary bodies. These examples illustrate an underdeveloped, but potentially important “sink‐to‐source” influence on the evolution of fluvial‐deltaic systems. The experiments illustrate the effects of (1) submarine hyperpycnal flows, (2) submarine delta front failure events, and (3) deformable substrates within prodelta and offshore settings. These submarine processes generate (1) erosional knickpoints in coastal rivers, (2) increased river channel occupancy times, (3) rapid rates of shoreline movement, and (4) localized zones of significant offshore sediment accumulation. Ramifications for coastal plain and deltaic stratigraphic patterns include changes in the hierarchy of scour surfaces, fluvial sand‐body geometries, reconstruction of sea‐level variability and large‐scale stratal geometries, all of which are linked to the identification and interpretation of sequences and systems tracts.     

Allogenic forcing of the late Quaternary Rhine–Meuse fluvial record: the interplay of sea-level change, climate change and crustal movements

The Rhine–Meuse system in the west‐central Netherlands is a continental‐scale fluvial system bordered by an extremely wide continental shelf. Consequently, late Quaternary eustatic sea‐level changes have resulted in dramatic shoreline displacements, by as much as 800 km. In addition, changes in climate have been severe, given the latitudinal and palaeogeographic setting of the Rhine–Meuse system. We investigated the relative importance of these allogenic controls on fluvial aggradation and incision during the last two glacial–interglacial cycles. We used optical dating of quartz from ～30 samples in a cross‐section perpendicular to the palaeoflow direction, allowing us to correlate periods of aggradation and incision with independent records of sea‐level change, climate change and glacio‐isostatic crustal movements. We found the long‐term aggradation rate to be ～8 cm kyr^?1, a value similar to previous estimates of tectonic subsidence rates in the study area. Several excursions from this long‐term aggradation trend could be identified for the last glacial–interglacial cycle. Dry climatic conditions with relatively high sediment supply induced aggradation during oxygen‐isotope stages (OIS) 4 and 3. Build‐up of a glacio‐isostatic forebulge during OIS 2 is a likely cause of incision around the Last Glacial Maximum, followed by an aggradation phase during forebulge collapse. Sea‐level highstands during OIS 5 have likely resulted in the aggradation of coastal prisms, but only minor, basal estuarine deposits have been preserved because these coastal prisms were prone to erosion during ensuing sea‐level falls. Overall, the sedimentary record is dominated by strata formed during time intervals when the study area was completely unaffected by sea‐level control, and our evidence shows that the falling‐stage systems tract has the highest preservation potential. Our study highlights the importance of considering the complex interplay of both upstream and downstream controls to obtain a comprehensive understanding of the evolution of basin‐margin successions.     

Macrogeomorphic evolution of the post-Triassic Appalachian mountains determined by deconvolution of the offshore basin sedimentary record

A perplexing macrogeomorphic problem is the persistence of topography in mountain ranges that were initially formed by orogenic events hundreds of millions of years old. In this paper, we deconvolve the post-Triassic macrogeomorphic history of a portion of one of these ranges, the central and northern Appalachians, using a well-documented offshore isopach sedimentary record of the US Atlantic margin. Topography is an important signature of tectonic, eustatic and/or geomorphic processes that produces changes in the erodible thickness of the crust (ETC). We define ETC as the total thickness of crust that would have to be consumed by erosion to reduce the mean elevation of a landscape to sea level. We use the term ‘source flux’, designated by ν˙, to describe the rate of change in ETC attributed to deep-seated geological processes such as crustal thickening, crustal extension, magmatic intrusions or dynamic flow in the mantle. In a mountain belt, the rate of change of mean elevation with respect to a base level, designated by ? ′, can be represented as ? ′ = c(ν˙ ? k _d z ′ ?; ? _c ) ?& hairsp;l˙ , where k _d is a proportionality constant relating the mean mechanical erosion rate to mean elevation, ? _c is the mean chemcial erosion rate, c  is a compensation ratio (held constant for Airy isostasy at 0.21) and l˙  is the rate of eustatic sea-level change. This equation describes the sum of constructive source terms, two destructive erosion terms and a eustatic sea-level term. We use this simple linear equation to develop a landscape evolution model based on the concept of a unit response function. The unit response function is analogous to a unit hydrograph which describes the relationship between input (rainfall) and output (discharge) in a hydrological system. In our case, we solve for the general relationship between the source flux into the mountain belt and the erosional flux out of the belt. Offshore sediment volumes are used to determine the erosional flux. Drainage basin area is treated as either a constant (pinned drainage divide) or as increasing through time (migrating drainage divide). We use a third-order polynomial fit to a global sea-level curve to account for long-term eustatically driven changes in ETC and in drainage basin area. Chemical erosion is treated as a constant fixed at 5 m Myr^?1. We consider two end-member models. The first is a ‘tectonic’ model in which the source flux is allowed to vary while k _d is assumed to be constant over geological time and equal to its mean Pleistocene value of about 0.07 Myr^?1. The second is an ‘erodibility’ model in which k _d is allowed to vary, reflecting changes in climatic, climatic variables and rock-type erodibility, while the source flux is held constant at zero. The ‘tectonic’ model reveals four important increases in the source flux, ranging from 200 to 2000 m Myr^?1 that occur over short (<10 Myr) time spans, followed by a protracted period (>25 Myr) where ν˙ drops below zero to values of ?1000 to ?6000 m Myr^?1. The ‘erodibility’ model produces a topography that decays in a step-like fashion from an initial mean elevation that ranges between ～1800 and 2300 m. These models cannot unequivocally distinguish the relative importance of tectonic vs. climatic processes in the macrogeomorphic evolution of the post-rift Appalachians, but they do provide some first-order quantitative prediction about these two end-member options. In light of existing stratigraphic, geological and thermochronological data, we favour the tectonic model because most of the events correlate well in time and form with known syn- and post-rift magmatic events. Nevertheless, the most recent episode of increased sediment flux to the offshore basins during the Miocene remains difficult to explain by either model. Limited evidence suggests that this event may reflect asthenospheric flow-driven uplift and the development of dynamically supported topography at a time when mechanical erosion rates were increasing in response to a cooling terrestrial climate.     

Fluvial response to sea-level changes: a quantitative analogue, experimental approach

ABSTRACT Quantitative evaluation of fluvial response to allogenic controls is crucial for further progress in understanding the stratigraphic record in terms of processes that control landscape evolution. For instance, without quantitative insight into time lags that are known to exist between sea‐level change and fluvial response, there is no way to relate fluvial stratigraphy to the sea‐level curve. It is difficult to put firm constraints on these time‐lag relationships on the basis of empirical studies. Therefore, we have started to quantify time‐averaged erosion and deposition in the fluvial and offshore realms in response to sea‐level change by means of analogue modelling in a 4 × 8‐m flume tank. The rate of sea‐level change was chosen as an independent variable, with other factors such as sediment supply, discharge and initial geometry kept constant over the course of 18 experiments. Our experimental results support the common view that neither fall nor rise in sea level affects the upstream fluvial system instantaneously. An important cause for the delayed fluvial response is that a certain amount of time is required to connect initial incisions on the newly emergent shelf (canyons) with the fluvial valley. Lowering of the fluvial longitudinal profile starts only after the connection of an active shelf canyon with the fluvial valley; until that moment the profile remains steady. We quantified the process of connection and introduced the quantity ‘connection rate’. It controlled, in conjunction with the rate of sea‐level fall: (1) the amount of fluvial degradation during sea‐level fall; (2) the total sediment volume that bypasses the shelf edge; (3) the percentage of fluvial relative to shelf sediment in the lowstand delta; (4) the volume of the transgressive systems tract and (5) the amount of diachroneity along the sequence boundary. Our experiments demonstrate also that the sequence‐stratigraphic concept is difficult to apply to continental successions, even when these successions have been deposited within the influence of sea level.     

Role of autoretreat and A/S changes in the understanding of deltaic shoreline trajectory: a semi-quantitative approach

ABSTRACT There is continued interest in how the rate of relative sea‐level rise [A ( > 0)] and the rate of sediment supply [S] function during the growth and evolution of deltaic shorelines. The theory of shoreline autoretreat, recently corroborated in flume experiments, claims that (1) A( > 0) and S can never be in equilibrium, and (2) shoreline or shelf‐edge progradation inevitably turns to retrogradation, when relative sea level is rising even modestly and even if A/S = const (> 0). Autoretreat arises because the area of the clinoform surface of the delta (or shelf edge) per kilometer of shoreline must increase as the relative sea level rises, and the delta (or shelf edge) progrades into deeper water. A finite sediment supply rate is thus liable to become inadequate to sustain progradation. The problem increases further as a rising sea level also greatly increases the delta‐plain volume that needs to be filled, further limiting the progradation of the system. The fundamental trajectory of shoreline migration is thus one characterized by a concave‐landward shape, even under the steady forcing of the basin. The magnitudes of A (> 0) and S, or A/S do not determine whether the landward turnaround of the shoreline is realized or not, but affect merely the length and height of the fundamental trajectory curve. Thus, any attempt to detect and interpret temporal changes in A and S from the observed stratigraphic record of shoreline trajectory needs first to take full account of the inbuilt autoretreat mechanism. We develop here a simple, semi‐quantitative method of reconstructing the basin conditions (A and S) from the stratigraphic record of prograding deltaic shorelines (or prograding shelf‐margin clinoforms) on the basis of the theory of shoreline autoretreat. The deterministic nature of the autoretreat theory is advantageous in managing this latter issue, because any expected or unexpected change emerges as some discrepancy from a trajectory that was predicted for the initial conditions. The autoretreat theory also provides a convenient graphical method of dealing with the uncertainty of the field data, and with evaluating the accuracy of any reconstruction. Our methodology has been developed to deal with the behaviour of deltaic shorelines, but is basically applicable to any clinoform system, the development of which is affected by relative sea level. The suggested method is applied to an Early Eocene (Ypresian) regressive shoreline succession in the Central Tertiary Basin on Spitsbergen. The studied regressive wedge developed as a delta‐driven, progradational shelf‐margin system under a regime of overall (i.e. long‐term) rise of relative sea level, but also suffered short‐term sea‐level falls associated with valley incisions on the coastal plain and shelf. On the assumption that S was constant or was steadily decreasing, the analysis of field data obtained from three sites within the basin suggests that the initial water depth in the basin was around 0.45 km, and that the overall relative sea‐level rise (c. 0.80 km) happened largely during an early time period and was followed by a longer period of much lower rate of rise. This pattern of relative sea‐level rise is consistent with the Palaeogene tectonic subsidence trend of the basin which was determined independently through a geohistory analysis. The uncertainty of the field data does not negate our reconstruction. The combined effects of autoretreat and A/S changes on a deltaic shoreline trajectory are confirmed through the development of an autoretreat‐based methodology. Conventional sequence stratigraphic models that assume a possible equilibrium condition between A and S are both conceptually misleading and insufficient to analyse basin conditions quantitatively. Sequence stratigraphic analyses of shorelines need to incorporate the autoretreat concept.     

Stream-terrace genesis: implications for soil development

Genesis of three distinct types of stream terraces can be understood through application of the concepts of tectonically induced downcutting, base level of erosion, complex response, threshold of critical power, diachronous and synchronous response times, and static and dynamic equilibrium. Climatic and tectonic stream terraces are major terraces below which flights of minor complex-response degradation terraces can form.These three types of terraces can be summarized by describing a downcutting-aggradation-renewed downcutting sequence for streams with gravell bedload. By tectonically induced downcutting, streams degrade to achieve and maintain a dynamic equilibrium longitudinal profile at the base level of erosion. Lateral erosion bevels bedrock beneath active channels to create major straths that are the fundamental tectonic stream-terrace landform. Aggradation events record brief reversals of long-term tectonically induced downcutting because they raise active channels. They may be considered as major (the result of climatic perturbations) or minor (the result of complex-response model types of perturbations). Climatically controlled aggradation followed by degradation leaves an aggradation surface; this type of fill-terrace tread is the fundamental climatic stream-terrace landform. Aggradation surfaces may be buried by subsequent episodes of deposition unless intervening tectonically induced downcutting is sufficient for younger aggradation surfaces to form below older surfaces. Raising of the active channel by either tectonic uplift or by climatically induced aggradation provides the vertical space for degradation terraces to form; first in alluvial fill and then in underlying bedrock along tectonically active streams. These are complex-response terraces because they result from interactions of dependent variables within a given fluvial system. Pauses in degradation to a new base level of erosion, and/or minor episodes of backfilling, lead to formation of complex-response fill-cut and strath, or of fill terraces. Fill-cut terraces are formed in alluvium; they are complex-response terraces because they are higher than the base level of erosion. Good exposures and dating are needed to distinguish static equilibrium complex-response minor strath terraces from dynamic equilibrium tectonic (major) straths. Strath terraces may be regarded as complex-response terraces where degradation rates between times terrace-tread formation exceed the long-term uplift rate for the reach based on ages and positions of tectonic terraces.Late Quaternary global climatic changes control aggradation events and even the times of cutting of major (tectonic) straths, because the base level of erosion can not be attained during times of climatically driven aggradation-degradation events.Most terrace soils form on treads of climatic and complex-response terraces. Aggradation surfaces may provide an ideal flight of terraces on which to study a soils chronosequence. Each aggradation event is recorded by a single relict soil where tectonically induced downcutting is sufficient to provide clear altitudinal separation of the terrace treads. Multiple paleosols are typical of tectonically stable regions where younger aggradation events spread alluvium over treads of older climatic terraces. Pedons on a climatic terrace in a small fluvial system commonly are roughly synchronous - variations of soil properties that can be attributed to temporal differences will be minor compared to altitudinally controlled climatic factors. Climatic terraces of adjacent watersheds also should be roughly synchronous (correlatable) - variations of soil properties that can be attributed to temporal differences will be minor compared to lithologic and climatic factors between different watersheds. Such generalizations may not apply to basins with sufficient relief that geomorphic responses to climatic changes occur at different and overlapping times, and to large rivers whose widely separated reaches are characterized by different response times to climatic perturbations. Soils on climatic terraces of distant watershedswill not be synchronous if their respective aggradation events occur during full-glacial times and interglacial times. Soils on some complex-response terraces may be diachronous within a given fluvial system, and typically are diachronous between watersheds.     

Stratigraphic signatures of translation of thrust-sheet top basins over low-angle detachment faults

Abstract Low‐angle detachment faults and thrust‐sheet top basins are common features in foreland basins. However, in stratigraphic analysis their influence on sequence architecture is commonly neglected. Usually, only eustatic sea level and changing flexural subsidence are accounted for, and when deformation is considered, the emphasis is on the generation of local thrust‐flank unconformities. This study analyses the effects of detachment angle and repetitive detachment activation on stratigraphic stacking patterns in a large thrust‐sheet top basin by applying a three‐dimensional numerical model. Model experiments show that displacement over low‐angle faults (2–6°) at moderate rates (～5.0 m kyr^?1) results in a vertical uplift component sufficient to counteract the background flexural subsidence rate. Consequently, the basin‐wide accommodation space is reduced, fluvio‐deltaic systems carried by the thrust‐sheet prograde and part of the sediment supply is spilled over towards adjacent basins. The intensity of the forced regression and the interconnectedness of fluvial sheet sandstones increases with the dip angle of the detachment fault or rate of displacement. In addition, the delta plain is susceptible to the formation of incised valleys during eustatic falls because these events are less compensated by regional flexural subsidence, than they would be in the absence of fault displacement.     

Drainage integration and sediment dispersal in active continental rifts: A numerical modelling study of the central Italian Apennines

Progressive integration of drainage networks during active crustal extension is observed in continental areas around the globe. This phenomenon is often explained in terms of headward erosion, controlled by the distance to an external base‐level (e.g. the coast). However, conclusive field evidence for the mechanism(s) driving integration is commonly absent as drainage integration events are generally followed by strong erosion. Based on a numerical modelling study of the actively extending central Italian Apennines, we show that overspill mechanisms (basin overfilling and lake overspill) are more likely mechanisms for driving drainage integration in extensional settings and that the balance between sediment supply vs. accommodation creation in fault‐bounded basins is of key importance. In this area drainage integration is evidenced by lake disappearance since the early Pleistocene and the transition from internal (endorheic) to external drainage, i.e. connected to the coast. Using field observations from the central Apennines, we constrain normal faulting and regional surface uplift within the surface process model CASCADE (Braun & Sambridge, 1997, Basin Research, 9, 27) and demonstrate the phenomenon of drainage integration, showing how it leads to the gradual disappearance of lakes and the transition to an interconnected fluvial transport system over time. Our model results show that, in the central Apennines, the relief generated through both regional uplift and fault‐block uplift produces sufficient sediment to fill the extensional basins, enabling overspill and individual basins to eventually become fluvially connected. We discuss field observations that support our findings and throw new light upon previously published interpretations of landscape evolution in this area. We also evaluate the implications of drainage integration for topographic development, regional sediment dispersal and offshore sediment supply. Finally, we discuss the applicability of our results to other continental rifts (including those where regional uplift is absent) and the importance of drainage integration for transient landscape evolution.     

Development of transient carbonate ramps in an evolving foreland basin

This study of Eocene carbonate succession in the Dinaric Foreland Basin of northern Dalmatia, Croatia, integrates palaeontological and sedimentological data to document a range of carbonate ramps formed intermittently during the basin tectonic development. The end‐Cretaceous basal erosional unconformity records the coupling of Adria and Eurasia crustal plates, with an antiformal uplift along their suture zone. The overlying late Ypresian carbonate ramp, spanning biozones SBZ 11–12, developed on the forebulge flank of a shallow‐marine early synclinal basin. Basal grainstone/packstone facies, dominated by encrusting foraminifers with alveolinids and miliolids, pass upwards into packstones dominated by miliolids and rotaliids with bryozoan and echinoid fragments, indicating an increased bathymetry of the retreating forebulge flank. Deposition of grainstone facies preceded an end‐Ypresian (SBZ 12/13 transition) subaerial exposure due to post‐subductional isostatic uplift. The younger, middle to late Eocene carbonate ramps (SBZ 13–19) formed episodically as perched isolated features on blind‐thrust anticlines in a bathymetrically diversified wedge‐top basin, where phases of clastic and skeletal biogenic sedimentation alternated due to disharmonic thrusting and relative sea‐level changes. Clastic sedimentation reflects anticline crest erosion and a forced‐regressive progradation of gravelly foreshore and sandy shoreface facies over heterolithic offshore‐transition and muddy offshore facies on the anticline flank. Biogenic sedimentation represents inner‐ to middle‐ramp environments, with the latter terminating bluntly in muddy offshore environment. An outer‐ramp environment, known from classic ramp models, was lacking due to bathymetric threshold. Analysis of larger benthic foraminifers (LBF), as biostratigraphic age indicators and palaeobathymetric proxies, helped distinguish systems tracts and determine their time span. A comparison of local and global sea‐level changes allowed the interplay of tectonic and eustatic forcing to be deciphered for the study area.     

CO2‐warming,rising sea‐level and retreating coasts: review and critique

A CO₂‐warming atmospheric scenario, whereby increased concentrations of ‘greenhouse’ gases result in higher temperatures that either melt near‐polar ice or cause thermal expansion of ocean waters, thus leading to increased sea‐levels and exacerbated coastal erosion, assumes fundamental but unproven cause‐and‐effect relationships. General circulation models have reinforced claims of an accelerated warming and indirectly given support to the complete scenario, but ignore the point that global climate and not just air temperatures have changed over the past century. Indeed, it is difficult to prove that air temperatures have warmed naturally outside of urban centres over this period. To attribute recent temperature increases to anthropogenic factors and to extrapolate these trends to the future also ignores the historic variability of climate. What is more, an eustatic rise in sea‐level cannot be discerned from the background noise of technically or climatically induced changes. Even if sea‐level was rising, coastal erosion may be accounted for better by a suite of inter‐related climatic factors including changes in rainfall regimes, hemispheric circulation and storminess.     

Origin of the in situ stress field in south-eastern Australia

The in situ stress field of south‐eastern Australia inferred from earthquake focal mechanisms and bore‐hole breakouts is unusual in that it is characterised by large obliquity between the maximum horizontal compressive stress orientation (S_Hmax) and the absolute plate motion azimuth. The evolution of the neotectonic strain field deduced from historical seismicity and both onshore and offshore faulting records is used to address the origin of this unusual stress field. Strain rates derived from estimates of the seismic moment release rate (up to ～10^?16 s^?1) are compatible with Quaternary fault–slip rates. The record of more or less continuous tectonic activity extends back to the terminal Miocene or early Pliocene (10–5 Ma). Terminal Miocene tectonic activity was characterised by regional‐scale tilting and local uplift and erosion, now best preserved by unconformities in offshore basins. Plate‐scale stress modelling suggests the in situ stress field reflects increased coupling of the Australian and Pacific Plate boundary in the late Miocene, associated with the formation of the Southern Alps in New Zealand.     

Investigating the surface process response to fault interaction and linkage using a numerical modelling approach

In order to better understand the evolution of rift‐related topography and sedimentation, we present the results of a numerical modelling study in which elevation changes generated by extensional fault propagation, interaction and linkage are used to drive a landscape evolution model. Drainage network development, landsliding and sediment accumulation in response to faulting are calculated using CASCADE, a numerical model developed by Braun and Sambridge, and the results are compared with field examples. We first show theoretically how the ‘fluvial length scale’, L_f, in the fluvial incision algorithm can be related to the erodibility of the substrate and can be varied to mimic a range of river behaviour between detachment‐limited (DL) and transport‐limited (TL) end‐member models for river incision. We also present new hydraulic geometry data from an extensional setting which show that channel width does not scale with drainage area where a channel incises through an area of active footwall uplift. We include this information in the coupled model, initially for a single value of L_f, and use it to demonstrate how fault interaction controls the location of the main drainage divide and thus the size of the footwall catchments that develop along an evolving basin‐bounding normal fault. We show how erosion by landsliding and fluvial incision varies as the footwall area grows and quantify the volume, source area, and timing of sediment input to the hanging‐wall basin through time. We also demonstrate how fault growth imposes a geometrical control on the scaling of river discharge with downstream distance within the footwall catchments, thus influencing the incision rate of rivers that drain into the hanging‐wall basin. Whether these rivers continue to flow into the basin after the basin‐bounding fault becomes fully linked strongly depends on the value of L_f. We show that such rivers are more likely to maintain their course if they are close to the TL end member (small L_f); as a river becomes progressively more under supplied, i.e. the DL end member (large L_f), it is more likely to be deflected or dammed by the growing fault. These model results are compared quantitatively with real drainage networks from mainland Greece, the Italian Apennines and eastern California. Finally, we infer the calibre of sediments entering the hanging‐wall basin by integrating measurements of erosion rate across the growing footwall with the variation in surface processes in space and time. Combining this information with the observed structural control of sediment entry points into individual hanging‐wall depocentres we develop a greater understanding of facies changes associated with the rift‐initiation to rift‐climax transition previously recognised in syn‐rift stratigraphy.     

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.     

Controls on landscape and drainage evolution in regions of distributed normal faulting: Perachora Peninsula,Corinth Rift,Central Greece

Zones of distributed faulting with narrow (2–3 km) across‐strike spacing form a common structural style within rifts, especially in accommodation zones, and contrast with crustal‐scale half‐grabens, where strain is localised on normal faults spaced 10–30 km apart. These contrasting styles are likely to have a significant impact on geomorphic development, sediment routing and the stratigraphic record. Perachora Peninsula, in the eastern part of the active Corinth Rift, Greece, is one such zone of distributed faulting. We analyse the topography and drainage networks developed around these closely spaced normal faults, and compare our results with published studies from crustal‐scale half‐grabens. We subdivide the Perachora Peninsula into a series of drainage domains and examine the tectono‐geomorphic evolution of three domains that best represent the range of topographic characteristics, base levels and drainage network styles. We interpret that the perched, endorheic nature of the Asprokampos domain developed due to uplift and backtilt on offshore faults. The Pisia West domain, which drains the valley between the Skinos and Pisia Faults and responds to a perched base level, is interpreted to have experienced a complex base‐level history with episodic connections to sea level. The Skinos Relay domain drains to sea level, lying on the relay ramp between the closely spaced Kamarissa and Skinos Faults. Here, interaction between the displacement fields associated with each of the closely spaced faults controls the rate and style of landscape evolution. In contrast to crustal‐scale half‐grabens, observations from Perachora Peninsula suggest that zones of distributed faulting may be characterised by: (i) perched, internal sediment sinks at different elevations, responding to multiple base levels; (ii) minimal fault‐transverse sediment transport; (iii) interaction of uplift and subsidence fields associated with closely spaced faults, which modulate the rate and style of landscape response; and (iv) complex erosion and sedimentation histories, the evidence for which may have low preservation potential in the stratigraphic record.     

Shoreline trajectories on a glacially influenced stable margin – insight from the Barents Sea Shelf, NW Russia

This study describes shoreline migration paths for late Quaternary sediments on the inner Barents Sea shelf between Kola and the Pechora Sea. The depositional geometries provide an example of stratigraphical architecture in a glacially influenced basin prone to isostatic movements as well as rapid and high-amplitude changes in eustatic sea level. The depositional geometries reflect asymmetrical relative sea level changes characterised by marine inundation upon deglaciation and prolonged forced regressions. Thus, all deposition occurs during the falling stage and lowstand systems tracts. The transgressive and highstand systems tracts are lacking and the maximum landward position of the shoreline is coinciding with the basal surface of forced regression. Shoreline migration is dominated by downward and seaward trajectories, but aggradation occurs on the falling limb of the relative sea level curve due to superimposed eustatic cycles of lower hierarchical order. Fluvial aggradation behind the shoreline takes place during the lowstand systems tract, but is also linked to high sediment supply and may also respond to superimposed lower order sea level fluctuations. Lateral variations in isostatic load due to asynchronous ice advances lead to regional variations in shoreline trajectories. Significant differences in sea level history exist across former ice margins leading to time-transgressive and laterally discontinuous stratigraphical surfaces. Sequence boundaries are not only diachronous along the depositional profile, but also laterally, and basal surfaces of forced regression are strongly diachronous across former ice margins. Absolute age control allows for estimates of the time differences along significant stratigraphical surfaces.     

Cenozoic uplift of Nuussuaq and Disko, West Greenland—elevated erosion surfaces as uplift markers of a passive margin

Remnants of a high plateau have been identified on Nuussuaq and Disko, central West Greenland. We interpret the plateau as an erosion surface (the summit erosion surface) formed mainly by a fluvial system and graded close to its former base level and subsequently uplifted to its present elevation. It extends over 150 km east–west, being of low relative relief, broken along faults, tilted westwards in the west and eastwards in the east, and having a maximum elevation of ca. 2 km in central Nuussuaq and Disko. The summit erosion surface cuts across Precambrian basement rocks and Paleocene–Eocene lavas, constraining its age to being substantially younger than the last rift event in the Nuussuaq Basin, which took place during the late Maastrichtian and Danian. The geological record shows that the Nuussuaq Basin was subjected to subsidence of several kilometres during Paleocene–Eocene volcanism and was transgressed by the sea later during the Eocene. By comparing with results from apatite fission track analysis and vitrinite reflectance maturity data, it is suggested that formation of the erosion surface was probably triggered by an uplift and erosion event starting between 40 and 30 Ma. Surface formation was completed prior to an uplift event that started between 11 and 10 Ma and caused valley incision. This generation of valleys graded to the new base level and formed a lower erosion surface, at most 1 km below the summit erosion surface, thus indicating the magnitude of its uplift. Formation of this generation of valleys was interrupted by a third uplift event also with a magnitude of 1 km that lifted the landscape to near its present position. Correlation with the fission-track record suggests that this uplift event started between 7 and 2 Ma. Uplift must have been caused initially by tectonism. Isostatic compensation due to erosion and loading and unloading of ice sheets has added to the magnitude of uplift but have not significantly altered the configuration of the surface. It is concluded that the elevations of palaeosurfaces (surfaces not in accordance with present climate or tectonic conditions) on West Greenland's passive margin can be used to define the magnitude and lateral variations of Neogene uplift events. The striking similarity between the landforms in West Greenland and those on many other passive margins is also noted.     

河流阶地形成过程及其驱动机制再研究

河流阶地的形成是在内因(河流内部动力变化)和外因(低频和高频气候变化、构造运动、基准面变化)共同作用下的结果。受单一气候变化制约的河流阶地发育模式可以解释由于沉积物通量和径流量变化引起的河流堆积-侵蚀过程,但它难以解释形成多级阶地的逐步(或间歇性)下切过程。多级阶地的形成可能同时受到构造抬升和周期性气候变化的制约。由于下切过程的滞后效应,侵蚀和冰川均衡抬升、河谷的侧向侵蚀过程等影响,山地的构造抬升与河谷的下切之间并非一种简单的线性关系,应当慎用河谷的下切速率来代表山地的抬升速率。     

Using fluvial terraces to determine Holocene coastal erosion and Late Pleistocene uplift rates: An example from northwestern Hawke Bay, New Zealand

In order to make robust predictions of future coastal processes and hazards, historical rates of coastal processes such as coastal erosion need to be put into a long-term (Holocene) context. In this study a methodology is proposed that uses fluvial terraces to construct longitudinal profiles which can be projected offshore to infer paleo-coastline positions. From these positions, an average Holocene coastal erosion rate can be calculated. This study also shows how constraints can be placed on sea level changes and Late Pleistocene uplift rates using fluvial terraces, and by assuming the latter has been constant since  55–37 ka, these constraints feedback into the coastal erosion rate calculations. For the northwestern Hawke Bay (North Island, New Zealand) coastline, Late Pleistocene uplift rates of 0.6 ± 0.2, 0.6 ± 0.2, and − 0.1 ± 0.1 (i.e., stable or subsiding) mm/yr have been determined for the Waikari, Mohaka, and Waihua River mouths, respectively. These rates are consistent with previous interpretations of subsidence to the northeast and uplift being the result of regional, subduction-related processes. A Holocene coastal erosion rate of 0.5 ± 0.1 m/yr was determined for the Waikari River mouth, which is at the higher end of the calculated historical ( 1880–1980) rates (0.02–0.5 m/yr). If this difference is significant, then two possible reasons for this difference are: (i) the historical rate is affected by events such as the 1931 Napier earthquake, and (ii) the Holocene rate is the average of a steadily declining rate over the last 7.3 ka.     

Decoding temporal and spatial patterns of fault uplift using transient river long profiles

We present detailed observations of rivers crossing active normal faults in the Central Apennines, Italy, where excellent constraints exist on the temporal and spatial history of fault movement. We demonstrate that rivers with drainage areas > 10 km² and crossing faults that have undergone an increase in throw rate within the last 1 My, have significant long-profile convexities. In contrast, channels that cross faults that have had a constant-slip rate for 3 My have concave-up profiles and have similar concavities and steepness indices to rivers that do not cross any active fault structures. This trend is consistent across the Central Apennines and cannot be explained by appeal to lithology or regional base level change. The data challenge the belief that active faulting must always be reflected in river profiles; instead, the long-profile convexities are best explained as a transient response of the river system to a change in tectonic uplift rate. Moreover, for these rivers we demonstrate that the height of the profile convexity, as measured from the fault, scales with the magnitude of the uplift rate increase on the fault; and we establish that this relationship holds for throw rate variation along strike for the same fault segment, as well as between faults. These findings are shown to be consistent with predictions of channel response to changing uplift rate rates using a detachment-limited fluvial erosion model, and they illustrate that analysis of the magnitude of profile convexities has considerable predictive potential for extracting tectonic information. We also demonstrate that the migration rate of the profile convexities varies from 1.5–10 mm/y, and is a function of the slip rate increase as well as the drainage area. This is consistent with n > 1 for the slope exponent in a classical detachment-limited stream-power erosion law, but could potentially be explained by incorporating an erosion threshold or an explicit role for sediment in enhancing erosion rates. Finally, we show that for rivers in extensional settings, where the response times to tectonic perturbation are long (in this case > 1 My), attempts to extract tectonic uplift rates from normalised steepness indices are likely to be flawed because topographic steady state has not yet been achieved.     

Spatial relationship between topography and rock uplift patterns in asymmetric mountain ranges based on a stream erosion model

The spatial relationship between topography and rock uplift patterns in asymmetric mountain ranges was investigated using a stream erosion model in which the asymmetric rock uplift was given and erosion rates were proportional to the m-th power of the drainage area and the n-th power of the channel gradient. The model conditions were simple, and thus the effects of horizontal rock movement, diffusional processes, and erosion thresholds were neglected, and spatially uniform precipitation, lithology, and vegetation were assumed. In asymmetric mountain ranges, under realistic exponent conditions (m < n) and the above assumptions, the surface erosion rate is faster on the steeper side and slower on the gentler side. The topographic axis migrates away from the rock uplift axis toward the center of the mountain range owing to the contrast in erosion rates. This migration continues until the erosion is balanced with rock uplift. In a dynamic steady state, the topographic pattern is independent of the rock uplift rate as indicated by an analytical solution, and is prescribed by the rock uplift pattern and the exponents m and n. As the asymmetry of the rock uplift pattern increases, the topographic axis migrates a greater distance. The location of the topographic axis is related to the location of the rock uplift axis by a simple logarithmic function, for a wide range of m and n. The fit of the numerical results and the logarithmic function is particularly good when m = 0.5 and n = 1.0. If the rock uplift pattern in asymmetric mountain ranges is known, the value of n − 5m/4 can be constrained based on the logarithmic relation, assuming a dynamic steady state. On the other hand, if the value of n − 5m/4 is known in an asymmetric mountain range, the rock uplift pattern can be estimated directly from the topography. This relation was applied to the Suzuka Range in central Japan, and the value of n − 5m/4 was estimated for an assumed reverse fault motion.