Planform evolution of deltas with graded alluvial topsets: Insights from three‐dimensional tank experiments,geometric considerations and field applications

The profile of a river that conveys sediment without net deposition and net erosion is referred to as ‘graded’ with respect to vertical aggradation of the river segment. Three experimental series, designed in terms of the autostratigraphic view of alluvial grade, were conducted to clarify the diagnostic spatial behaviour of graded alluvial–deltaic rivers: an ‘R series’, which utilized a moving boundary setting with a stationary base level; an ‘F series’ in a fixed boundary setting with a stationary base level to produce ‘forced grade’; and an ‘M series’ in a moving boundary setting with constant base‐level fall to produce ‘autogenic grade.’ The results of the three experimental series, combined with geometrical modelling of the effects of basin water depth and other experimental data, suggest the following: (i) in a graded alluvial–deltaic system, lateral shifting and avulsing of active distributary channels are suppressed regardless of whether the downstream boundary of the deltaic system is fixed; (ii) in a delta with a downstream‐fixed boundary, the graded streams are stabilized within a valley that is incised in the axial part of the delta plain, whereby the alluvial plain outside the valley is abandoned and terraced; (iii) in moving boundary settings, the graded river simply extends basinward as a linearly elongated channel and lobe system without cutting a valley; and (iv) a modern forced‐graded alluvial river is most likely to be found in a valley incised into a fan delta in front of very deep water, and the stratigraphic signal of fossil autogenic‐graded rivers will be found in deltaic successions that accumulated in the outer to marginal areas of deltaic continental shelves during sea‐level falls. This renewed autostratigraphic view of alluvial grade suggests a thorough reconsideration of the conventional understanding that an alluvial river feeding a progradational delta is graded with a stationary base level.     

A three‐dimensional numerical model of sediment transport,erosion and deposition within a network of channel belts,floodplain and hill slope: extrinsic and intrinsic controls on floodplain dynamics and alluvial architecture

A three‐dimensional numerical model of sediment transport, erosion and deposition within a network of channel belts and associated floodplain is described. Sediment and water supply are defined at the upstream entry point, and base level is defined at the downstream edge of the model. Sediment and water are transported through a network of channels according to the diffusion equation, and each channel has a channel belt with a width that increases in time. The network of channels evolves as a result of channel bifurcation and abandonment (avulsion). The timing and location of channel bifurcation is controlled stochastically as a function of the cross‐valley slope of the floodplain adjacent to the channel belt relative to the down‐valley slope, and of annual flood discharge. A bifurcation develops into an avulsion when the discharge of one of the distributaries falls below a threshold value. The floodplain aggradation rate decreases with distance from the nearest active channel belt. Channel‐belt degradation results in floodplain incision. Extrinsic (extrabasinal, allogenic) and intrinsic (intrabasinal, autogenic) controls on floodplain dynamics and alluvial architecture were modelled, and sequence stratigraphy models were assessed. Input parameters were chosen based on data from the Rhine–Meuse delta. To examine how the model responds to extrinsic controls, the model was run under conditions of changing base level and increasing sediment supply. Rises and falls in base level and increases in sediment supply occurred over 10 000 years. Rising base level caused a wave of aggradation to move up‐valley, until aggradation occurred over the entire valley. Frequency of bifurcations and avulsions increased with rate of base‐level rise and aggradation rate. Channel‐belt width varied with water discharge and the lifespan of the channel belt. Wide, connected channel belts (and high channel‐deposit proportion) occurred around the upstream inflow point because of their high discharge and longevity. Less connected, smaller channel belts occurred further down‐valley. Such alluvial behaviour and architecture is also found in the Rhine–Meuse delta. During base‐level fall, valley erosion occurred, and the incised valley contained a single wide channel belt. During subsequent base‐level rise, a wave of aggradation moved up‐valley, filling the incised valley. Bifurcation and avulsion sites progressively moved upstream. Relatively thin, narrow channel belts bordered and cut into the valley fill. These results differ substantially from existing sequence stratigraphy models. The increase in sediment supply from upstream resulted in an alluvial fan. Most bifurcations and avulsions occurred at the fan apex (nodal avulsion), and channel belts were the widest and the thickest here (giving high channel‐deposit proportion) due to their high discharge and longevity. The width and thickness of channel belts decreased down‐valley due to decreased discharge, longevity and aggradation rate. This behaviour occurs in modern alluvial fans. Intrinsic controls also affect floodplain dynamics and alluvial architecture. Variation of aggradation rate, bifurcation frequency and number of coexisting channel belts occurred over periods of 500 to 2000 years, compared with 10 000 years for extrinsic controls. This variation is partly related to local aggradation and degradation of channel belts around bifurcation points. Channel belts were preferentially clustered near floodplain margins, because of low floodplain aggradation rate and topography there.     

Autogenic dynamics of alluvial fans in endorheic basins: Outcrop examples and stratigraphic significance

Alluvial fans are relatively simple depositional systems, due to the direct coupling of sediment sources and adjacent accumulation areas. Nonetheless, general models of alluvial‐fan evolution and stratigraphy remain elusive, due to the great sensitivity of such systems to allogenic controls and their strongly case‐specific responses. Autogenic processes intrinsic to alluvial‐fan dynamics can complicate stratigraphic architectures, with effects not easily distinguishable from those of allogenic forcing. A distinction is made here between lateral autogenic dynamics, tied to spatial sediment distribution over fan surfaces, and vertical autogenic dynamics, related to independent incision‐aggradation cycles. Autogenic mechanisms have been highlighted recently by modelling studies, but remain poorly constrained in field‐based studies. Examples are presented here from the margins of the Cenozoic Teruel and Ebro basins (Spain), where alluvial fans accumulated thick successions during phases of basin topographic closure and endorheic drainage which promoted forced aggradation. Fan successions consist of conformable architectures of stacked clastic sheets, laterally continuous and with no evidence of internal unconformities, inset architectures, fan segmentation or preserved incised channels. Continuous aggradation in these closed basins strongly inhibited ‘vertical’ autogenic dynamics in the form of fan head and through fan incision, due to the forced rise in geomorphic base level and the creation of positive accommodation. Furthermore, the lack of incised channels favoured widespread sediment transport and aggradation over broad fan sectors in relatively short time spans, in contrast to the typical occurrence of active lobes and abandoned fan surfaces caused by ‘lateral’ autogenic dynamics. Stratigraphic records of alluvial fans developed in endorheic basins are essentially complete and largely unaffected by autogenic processes. The latter characteristic implies that they can be more unambiguously interpreted in terms of allogenic forcing, because stratigraphic signatures are not complicated by the effects of complex fan autodynamics.     

A late Ordovician (Hirnantian) karstic surface in a submarine channel,recording glacio‐eustatic sea‐level changes: Meifod,central Wales

The growth and decay of the end‐Ordovician Gondwanan glaciation is globally reflected by facies changes in sedimentary sequences, which record a major eustatic fall and subsequent rise in the Hirnantian Stage at the end of the Ordovician. However, there are different reported estimates of the magnitude and pattern of sea‐level change. Particularly good evidence for end‐Ordovician sea‐level change comes from a sequence at Meifod in central Wales, which has a karstified limestone unit within a channel incised into marine shelf sediments. Pre‐glacial (Rawtheyan) mudstones have a diverse fauna suggesting a mid‐to‐deep‐shelf water depth of c. 60 m. The channel, 20 m deep, was incised into these mudstones and partially filled with a mixture of fine sand and detrital carbonate. The taphonomy of bioclasts and intraclasts indicates that many had a long residence time on the sea floor or suffered diagenesis after shallow burial before being resedimented into the channel. The presence of carbonates on the Welsh shelf is atypical and they are interpreted as having accumulated as patches during a minor regression prior to the main glacio‐eustatic fall. Comparison of the carbon stable‐isotopic values of the bioclast material with the global isotopic record confirms that most of the material is of Rawtheyan age, but that some is Hirnantian. The resedimented carbonates lithified rapidly and formed a limestone, several metres thick, in the deepest parts of the channel. As sea‐level fell, this limestone was exposed and eroded into karstic domes and pillars with a relief of over 2 m. The overall, glacio‐eustatic, sea‐level fall is estimated to be in excess of 80 m. A succeeding sea‐level rise estimated to be 40–50 m is recorded in the laminated crust that mantles the karstic domes and pillars. The crust is formed of encrusting bryozoans, associated cystoids, crinoid holdfasts and clusters of the brachiopod Paromalomena, which is normally associated with mid‐shelf environments. Fine sands buried the karst topography and accumulated to fill the channel. In the sandstones at the base of the channel there is a Hirnantia fauna, while in the sandstones high in the channel‐sequence there is cross‐stratification characteristic of mid‐shoreface environments. This would indicate a fall of sea‐level of c. 30 m. The subsequent major transgression marking the end of the glaciation is not recorded at the Meifod locality, but nearby exposures of mudstones suggest a return to mid‐to‐deep‐shelf environments, similar to those that prevailed before the Hirnantian regression. The Meifod sequence provides strong evidence for the magnitude of the Hirnantian sea‐level changes and by implication confirm larger estimates for the size of the ice sheets. Smaller oscillations in relative sea‐level seen at Meifod may be local phenomena or may reflect eustatic changes that have not been widely reported elsewhere. Copyright © 2005 John Wiley & Sons, Ltd.     

Tectonically driven deposition and landscape evolution within upland incised valleys: Ambra Valley fill,Pliocene–Pleistocene,Tuscany, Italy

Sedimentation in the upstream reaches of incised valleys is predominantly of alluvial origin and, in most cases, independent from relative sea‐level or lake‐level oscillations. Preserved facies distributions record the depositional response to a combination of allogenic factors, including tectonics, climate and landscape evolution. Tectonics drive fluvial aggradation and degradation through local changes in gradient, both longitudinal and transverse to the valley slope. This article deals with a Pliocene–Pleistocene fluvial valley fill developed in the north‐eastern shoulder of the Siena Basin (Northern Apennines, Italy). Evolution of the valley was not influenced by sea‐level or lake‐level changes and morphological and depositional evolution of valley resulted from extensional tectonics that gave rise to normal and oblique‐slip faults orthogonal and parallel to the valley axis. Data from both field observations and geophysical study are interpreted to develop a comprehensive tectono‐sedimentary model of coeval longitudinal and lateral tilting of the developing alluvial plain. Longitudinal tilting was generated by a transverse, upstream‐dipping normal fault that controlled the aggradation of fining‐upward strata sets. Upstream of the fault zone, valley back‐filling generated an architecture similar to that of classic, sea‐level‐controlled, coastal incised valleys. Downstream of the fault zone, valley down‐filling was related to an overwhelming sediment supply sourced and routed from the active fault zone itself. Lateral tilting was promoted by the activity of a fault oriented parallel to the valley axis, as well as by different offsets along near orthogonal faults. As a result, the valley trunk system experienced complex lateral shifts, which were governed by interacting fault‐generated subsidence and by the topographic confinement of progradational, flank‐sourced alluvial fans.     

Upper Ordovician–Lower Silurian transgressive–regressive cycles of the Central Iberian Zone (NE Jaén,Spain)

A thick Upper Ordovician shelf sequence was developed in the northern Gondwana margin (southernmost exposures of the Central Iberian Zone). Integrated sedimentologic and stratigraphic studies allow distinction between pedogenetic processes (Facies association C), shoreline deposits (Facies association S), proximal to distal shelf (Facies association L, H₁, H₂, H₃) and outer shelf zone or open marine environments (Facies association M, Mo). The vertical distribution of facies is characterized by the presence of regressive high frequency sequences (partial shelf progradational sequences), affected by the presence of catastrophic phenomena (storms). These sequences, in turn, can be classified into higher‐order transgressive (T)–regressive (R) cycles. Two second‐order T‐R megacycles (MC. Ord‐2 and MC. Sil‐1) limited by a major sequence boundary are identified. Traces of emersion (palaeokarsts and palaeosols) are detected along the sequence boundary, and these are related to the eustatic sea‐level fall that occurred during the Ashgillian. The MC. Ord‐2 and MC. Sil‐1 megacycles extend respectively from the Middle Arenig to the Ashgillian and from Late Ashgillian to the Late Llandovery. Major transgressive peaks occurred at the Llanvirn and at the Middle Llandovery (Aeronian). The vertical distribution of the facies delineates successive genetically related units in relation to relative sea‐level changes. Within the upper part of the first megacycle (MC. Ord‐2) six third‐order cycles are proposed (Lla‐1, Car‐1, Car‐2, Car‐3, Car‐4, Ash‐1), in which a transgressive and a regressive interval can be distinguished. Within the lower part of the second megacycle (MC. Sil‐1) two transgressive–regressive third‐order cycles are proposed (Lly‐1, Lly‐2). Copyright © 2005 John Wiley & Sons, Ltd.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Aggradation and carbonate accumulation of Holocene Norwegian cold‐water coral reefs

Cold‐water coral ecosystems present common carbonate factories along the Atlantic continental margins, where they can form large reef structures. There is increasing knowledge on their ecology, molecular genetics, environmental controls and threats available. However, information on their carbo‐nate production and accumulation is still very limited, even though this information is essential for their evaluation as carbonate sinks. The aim of this study is to provide high‐resolution reef aggradation and carbonate accumulation rates for Norwegian cold‐water coral reefs from various settings (sunds, inner shelf and shelf margin). Furthermore, it introduces a new approach for the evaluation of the cold‐water coral preservation within cold‐water coral deposits by computed tomography analysis. This approach allows the differentiation of various kinds of cold‐water coral deposits by their macrofossil clast size and orientation signature. The obtained results suggest that preservation of cold‐water coral frameworks in living position is favoured by high reef aggradation rates, while preservation of coral rubble prevails by moderate aggradation rates. A high degree of macrofossil fragmentation indicates condensed intervals or unconformities. The observed aggradation rates with up to 1500 cm kyr^?1 exhibit the highest rates from cold‐water coral reefs so far. Reef aggradation within the studied cores was restricted to the Early and Late Holocene. Available datings of Norwegian cold‐water corals support this age pattern for other fjords while, on the shelf, cold‐water coral ages are reported additionally from the early Middle Holocene. The obtained mean carbonate accumulation rates of up to 103 g cm^?2 kyr^?1 exceed previous estimates of cold‐water coral reefs by a factor of two to three and by almost one order of magnitude to adjacent sedimentary environments (shelf, slope and deep sea). Only fjord basins locally exhibit carbonate accumulation rates in the range of the cold‐water coral reefs. Furthermore, cold‐water coral reef carbonate accumulation rates are in the range of tropical reef carbonate accumulation rates. These results clearly suggest the importance of cold‐water coral reefs as local, maybe regional to global, carbonate sinks.     

Climate control on the evolution of Late Pleistocene alluvial‐fan and aeolian sand‐sheet systems in NW Germany

The Late Pleistocene was characterized by rapid climate oscillations with alternation of warm and cold periods that lasted up to several thousand years. Although much work has been carried out on the palaeoclimate reconstruction, a direct correlation of ice‐core, marine and terrestrial records is still difficult. Here we present new data from late Middle Pleniglacial to Lateglacial alluvial‐fan and aeolian sand‐sheet deposits in northwestern Germany. Records of Late Pleniglacial alluvial fans in central Europe are very rare, and OSL dating is used to determine the timing of fan aggradation. In contrast to fluvial systems that commonly show a delay between climate change and incision/aggradation, the small alluvial‐fan systems of the Senne area responded rapidly to climatic changes and therefore act as important terrestrial climate archives for this time span. The onset of alluvial‐fan deposition correlates with the climate change from warm to cold at the end of MIS 3 (29.3±3.2 ka). Strong fan progradation started at 24.4±2.8 ka and may be related to a period of higher humidity. The vertical stacking pattern of sedimentary facies and channel styles indicate a subsequrent overall decrease in water and sediment supply, with less sustained discharges and more sporadic runoffs from the catchment area, corresponding to an increasing aridity in central Europe during the Late Pleniglacial. Major phases of channel incision and fan aggradation may have been controlled by millennial‐scale Dansgaard–Oeschger cycles. The incision of channel systems is attributed to unstable climate phases at cold–warm (dry–wet) or warm–cold (wet–dry) transitions. The alluvial‐fan deposits are bounded by an erosion surface and are overlain by aeolian sand‐sheets that were periodically affected by flash‐floods. This unconformity might be correlated with the Beuningen Gravel Bed, which is an important marker horizon in deposits of the Late Pleniglacial resulting from deflation under polar desert conditions. The deposition of aeolian sand‐sheet systems (19.6±2.1 to 13.1±1.5 ka) indicates a rapid increase in aridity at the end of the Late Pleniglacial. Intercalated flash‐floods deposits and palaeosols (Finow type) point to temporarily wet conditions during the Lateglacial. The formation of an ephemeral channel network probably marks the warm‐cold transition from the Allerød to the Younger Dryas.     

Correlating Alpine glaciation with Adriatic sea‐level changes through lake and alluvial stratigraphy

We compare lake and alluvial stratigraphy along a frame connecting the southern Alpine foothills and the Adriatic Sea, with the aim of matching the effects of Alpine glaciation and sea‐level changes on sedimentation during the last glacial cycle. The palynostratigraphy of Lake Fimon provided proxies for regional vegetation and climate change and was coupled with sediment petrography, loss on ignition and magnetic susceptibility, disentangling alluvial phases from fluvioglacial activity related to culminations of the southeastern Alpine glaciers. The Fimon area was not reached by alluvial fans during the penultimate glacial maximum, nor by the sea transgression during the last interglacial, but a closed lake soon developed at the Eemian onset due to enhanced rainfall. Sea‐level fall at glacial inception triggered the entrenchment of the drainage network in the plain reaching the outer Fimon Basin. Slow aggradation, but no sign of fluvioglacial activity, lasted to 38.2 ± 1.45 cal. ka BP, when a major forest withdrawal took place, coeval to the spread of alluvial fans. By 27.5 ± 0.5 cal. ka BP the Fimon Basin was dammed by the Brenta outwash system. The main step of of forest recovery commenced at around (15.8) cal. ka BP, when apex trenching of the outwash fans was triggered by the glacier's decay. Copyright © 2011 John Wiley & Sons, Ltd.     

Fluvial response to rapid high‐amplitude lake‐level changes during the Late Weichselian and early Holocene,Ain River valley,Jura, France

In the present paper the effects of rapid, high‐amplitude base‐level changes during the last glacial‐interglacial transition were studied for the Ain River in eastern France. During the Würm glacial maximum (MIS 2) rapid aggradation by deep‐water Gilbert‐type deltas and shallow‐water fan deltas occurred at the margins of a 20 to 50 m deep proglacial lake. A temporal high‐amplitude lake‐level fall of 60 m resulted in gravel deposition by forced‐regressive deltas, followed by rapid lake‐level rise and fine‐grained glaciolacustrine deposition. During the final deglaciation, a rapid base‐level fall of 40 m resulted in a complex fluvial response. Knickpoint formation and headward incision of the highstand deltas and concomitant deposition of gravel sheets by forced‐regressive deltas and braided systems occurred in several depocentres on the former glacial lake floor. Preservation of highstand and falling‐stage deposits and terrace formation in the incised valley depended on vertical incision and lateral channel migration. Terraces are well developed in the former lake‐floor depressions, whereas vertical incision was dominant in the higher lake‐floor areas. The Ain terrace staircase was likely formed by autogenic processes during a single allogenic base‐level fall. This case study possibly offers an analogue for the preservation of interglacial highstand coastal deltas during sea‐level fall at warm‐to‐cold climate transitions, although the rates of base‐level fall are different.     

Stratigraphic organization and predictability of mixed coarse‐grained and fine‐grained successions in an upper slope Pleistocene turbidite system of the Peri‐Adriatic basin

The early Pleistocene clastic succession of the Peri‐Adriatic basin, eastern central Italy, records the filling of a series of piggyback sub‐basins that formed in response to the development of the eastward‐verging Apennine fold‐thrust belt. During the Gelasian (2·588 to 1·806 Ma), large volumes of Apennine‐derived sediments were routed to these basins through a number of slope turbidite systems. Using a comprehensive outcrop‐based dataset, the current study documents the depositional processes, stratigraphic organization, foraminiferal age and palaeodepth, and stratigraphic evolution of one of these systems exposed in the surroundings of the Castignano village. Analysis of foraminiferal assemblages consistently indicates Gelasian deposition in upper bathyal water depths. Sediments exposed in the study area can be broken into seven main lithofacies, reflecting specific gravity‐induced depositional elements and slope background deposition: (i) clast‐supported conglomerates (conglomerate channel‐fill); (ii) amalgamated sandstones (late stage sandstone channel‐fill); (iii) medium to thick‐bedded tabular sandstones (frontal splay sandstones); (iv) thin to thick‐bedded channelized sandstones (sandy channel‐fill); (v) medium to very thin‐bedded sandstones and mudstones (levée‐overbank deposits); (vi) pebbly mudstones and chaotic beds (mudstone‐rich mass‐transport deposits); and (vii) massive mudstones (hemipelagic deposits). Individual lithofacies combine vertically and laterally to form decametre‐scale, disconformably bounded, fining‐upward lithofacies successions that, in turn, stack to form slope valley fills bounded by deeply incised erosion surfaces. A hierarchical approach to the physical stratigraphy of the slope system indicates that it has evolved through multiple cycles of waxing then waning flow energy at multiple scales and that its packaging can be described in terms of a six‐fold hierarchy of architectural elements and bounding surfaces. In this scheme, the whole system (sixth‐order element) is comprised of three distinct fifth‐order stratigraphic cycles (valley fills), which define sixth‐order initiation, growth and retreat phases of slope deposition, respectively; they are separated by discrete periods of entrenchment that generated erosional valleys interpreted to record fifth‐order initiation phases. Backfilling of individual valleys progressed through deposition of two vertically stacked lithofacies successions (fourth‐order elements), which record fifth‐order growth and retreat phases. Fourth‐order initiation phases are represented by erosional surfaces bounding lithofacies successions. The component lithofacies (third‐order element) record fourth‐order growth and retreat phases. Map trends of erosional valleys and palaeocurrent indicators converge to indicate that the sea floor bathymetric expression of a developing thrust‐related anticline markedly influenced the downslope transport direction of gravity currents and was sufficient to cause a major diversion of the turbidite system around the growing structure. This field‐based study permits the development of a sedimentological model that predicts the evolutionary style of mixed coarse‐grained and fine‐grained turbidite slope systems, the internal distribution of reservoir and non‐reservoir lithofacies within them, and has the potential to serve as an analogue for seismic or outcrop‐based studies of slope valley fills developed in actively deforming structural settings and under severe icehouse regimes.     

Studying the effect of non‐spherical micro‐particles on Hoek–Brown strength parameter mi using numerical true triaxial compressive tests

The strength parameter m_i in the Hoek–Brown strength criterion is empirical and was developed by trial and error. To better understand the fundamental relationship between m_i and the physical characteristics of intact rock, this paper presents a systematic study of m_i by representing intact rock as a densely packed cemented particle material and simulating its mechanical behavior using particle flow modeling. Specifically, the three‐dimensional particle flow code (PFC3D) was used to conduct numerical true triaxial compression tests on intact rock and to investigate the effect of non‐spherical micro‐particle parameters on m_i. To generate numerical intact rock specimens containing non‐spherical micro‐particles, a new genesis process was proposed, and a specific loop algorithm was used based on the efficiency of the process and the acceptability of generated specimens. Four main parameters—number, aspect ratio, size, and shape—of non‐spherical micro‐particles were studied, and the results indicated that they all have great effect on m_i. The strength parameter m_i increases when the number, aspect ratio, or size is larger or the shape becomes more irregular, mainly as a result of the higher level of interlocking between particles. This confirms the observations from engineering experience and laboratory experiments. To simulate the right strength parameter m_i, it is important to use appropriate non‐spherical micro‐particles by controlling these four parameters. This is further demonstrated by the simulation of two widely studied rocks, Lac du Bonnet granite and Carrara marble. Copyright © 2014 John Wiley & Sons, Ltd.     

Palaeoenvironmental reconstruction of a middle Miocene alluvial fan to cyclic shallow lacustrine depositional system in the Calatayud Basin (NE Spain)

ABSTRACT The middle Miocene sedimentary fill of the Calatayud Basin in north‐eastern Spain consists of proximal to distal alluvial fan‐floodplain and shallow lacustrine deposits. Four main facies groups characteristic of different sedimentary environments are recognized: (1) proximal and medial alluvial fan facies that comprise clast‐supported gravel and subordinate sandstone and mudstone, the latter exhibiting incipient pedogenic features; (2) distal alluvial fan facies, formed mainly of massive mudstone, carbonate‐rich palaeosols and local carbonate pond deposits; (3) lake margin facies, which show two distinct lithofacies associations depending on their distribution relative to the alluvial fan system, i.e. front (lithofacies A), comprising massive siliciclastic mudstone and tabular carbonates, or lateral (lithofacies B) showing laminated and/or massive siliciclastic mudstone alternating with tabular and/or laminated carbonate beds; and (4) mudflat–shallow lake facies showing a remarkable cyclical alternation of green‐grey and/or red siliciclastic mudstone units and white dolomitic carbonate beds. The cyclic mudflat–shallow lake succession, as exposed in the Orera composite section (OCS), is dominantly composed of small‐scale mudstone–carbonate/dolomite cycles. The mudstone intervals of the sedimentary cycles are interpreted as a result of sedimentation from suspension by distal sheet floods, the deposits evolving either under subaerial exposure or water‐saturated conditions, depending on their location on the lacustrine mudflat and on climate. The dolomite intervals accumulated during lake‐level highstands with Mg‐rich waters becoming increasingly concentrated. Lowstand to highstand lake‐level changes indicated by the mudstone/dolomite units of the small‐scale cycles reflect a climate control (from dry to wet conditions) on the sedimentation in the area. The spatial distribution of the different lithofacies implies that deposition of the small‐scale cycles took place in a low‐gradient, shallow lake basin located in an interfan zone. The development of the basin was constrained by gradual alluvial fan aggradation. Additional support for the palaeoenvironmental interpretation is derived from the isotopic compositions of carbonates from the various lithofacies that show a wide range of δ¹⁸O and δ¹³C values varying from ?7·9 to 3·0‰ PDB and from ?9·2 to ?1·7‰ PDB respectively. More negative δ¹⁸O and δ¹³C values are from carbonate‐rich palaeosols and lake‐margin carbonates, which extended in front of the alluvial fan systems, whereas more positive values correspond to dolomite beds deposited in the shallow lacustrine environment. The results show a clear trend of δ¹⁸O enrichment in the carbonates from lake margin to the centre of the shallow lake basin, thereby also demonstrating that the lake evolved under hydrologically closed conditions.     

SOME APPLICATIONS OF THEORY AND EXPERIMENT TO THE STUDY OF BEDDING GENESIS

The origin of many sedimentation units deposited in granular cohesionless materials can be rationally explained by utilizing the concept of the profile of equilibrium. This basic concept can be formally expressed in terms of the variables, or groups of variables, that characterize a sediment transport system. A change in one or more variables or a shift in local base level will generally cause a shift in the spatial position of the profile, resulting in either aggradation or degradation. As shown by flume experiments, the rate of shift of the profile is a critical factor in bedding genesis. A relatively large and rapid upward shift of the profile results in the deposition of a tabular or wedge-shaped unit of cross-bedding, i.e., a laboratory delta. On the other hand, a gradual upward shift of the profile results in the deposition of a sequence of horizontal bedding that is commonly associated with intercalations of ripple or dune cross-bedding. For the intermediate case of a moderately rapid shift of the profile, the depositional sequence includes trough units and poorly-defined tabular units of cross-bedding with numerous intercalations of horizontal bedding. The concept of the profile of equilibrium therefore provides a rationale for considering the depositional framework for sedimentary structures produced by current flow.     

Fluvial response to Holocene climate change in low‐order streams of central Mexico

Alluvial sequences constitute a recognised source of information on past environmental change, but one that has scarcely been tapped in central Mexico. This paper reviews what is currently known about the Holocene alluvial stratigraphy of the region, focusing on the interplay between climate and the pace and style of sedimentation in the incised headwater reaches of stream networks. The records obtained in five different drainage basins – four in the state of Tlaxcala and one in Guanajuato – are presented and compared to published reconstructions of climate change. A near‐synchronous incision of all stream networks occurred close to 10 200 ¹⁴C a BP in response to an increase in precipitation and stream discharge. A spell of very humid but markedly seasonal conditions ensued, resulting in the formation of wet meadows along streams and the accumulation of thick fine‐textured valley fills dominated by cumulic soil A horizons. After 9100 ¹⁴C a BP a transition to a warmer and more arid climate provoked the thinning of vegetation cover on slopes, accelerated runoff and increased sediment delivery to streams. The aggradation of coarser‐textured valley fills poor in organic matter set in. It ceased or slowed down significantly after a few millennia as the studied stream reaches achieved a near‐graded condition adjusted to the relatively stable climate. Arid mid Holocene conditions are also reflected in the abundant precipitation of secondary carbonates in Guanajuato. At 3100 ¹⁴C a BP higher precipitation caused more frequent flooding and a resumption of aggradation. Shortly after that date sedentary farmers colonised Tlaxcala. Agriculture altered runoff and sediment delivery to streams and accelerated cut‐and‐fill cycles on a scale that masked the impact of any climatic fluctuations. Guanajuato was colonised later and its alluvial record suggests the persistence of a humid climate at least until 1000 ¹⁴C a BP. Copyright © 2009 John Wiley & Sons, Ltd.     

Late Quaternary marine deposition in New South Wales and southern Queensland — An evolutionary model

From new data on coastal and continental shelf morphology, sediments, stratigraphy and chronology, it is possible to formulate a general model of late Quaternary marine sedimentation, for New South Wales and southern Queensland. This model integrates various factors influencing deposition in coastal and shelf environments, in relation to glacio‐eustatic sea level oscillations.

The model involves several components, including (i) very slow to negligible continental margin subsidence during the Quaternary, (ii) an inherited geomorphic framework; (iii) oscillations of sea level of c 100 m amplitude every 100 000 years, with interglacial high sea levels being close to present and only the Last Interglacial being significantly higher; and (iv) a wave climate that induces a potential south to north littoral sand transport at all sea level positions.

Terrigenous sediment that is moved from the hinterland through embayments to the shelf is either stored as barrier, estuarine or inner shelf deposits, or lost to depositional sinks on the continental slope or into coastal dune fields. Over many glacial‐interglacial cycles, sand has been progressively moved northward and has accumulated in vast aeolian sand deposits in southern Queensland. Littoral sand transport was especially effective during sea levels lower than present. The relatively shallow and lower gradient shelf north of Newcastle (33°S) has encouraged preservation at the coast of a wide range of depositional morphologies, including Pleistocene barriers, whereas the steeper southern shelf has induced net sediment loss seawards and shoreline erosion, excpt in the Holocene. To account for Holocene barrier development in the southern region, the model invokes reworking of sand deposits stranded high on the inner shelf at the end of the Pleistocene Epoch. These were in disequilibrium with Postglacial marine processes that operated at a lower level of the sea than did those during the Last Interglacial maximum.     

Late Quaternary multiple incised valley systems: An unusually well‐preserved stratigraphic record of two interglacial valley‐fill successions from the Arno Plain (northern Tuscany,Italy)

Un‐fragmented stratigraphic records of late Quaternary multiple incised valley systems are rarely preserved in the subsurface of alluvial‐delta plains due to older valley reoccupation. The identification of a well‐preserved incised valley fill succession beneath the southern interfluve of the Last Glacial Maximum Arno palaeovalley (northern Italy) represents an exceptional opportunity to examine in detail evolutionary trends of a Mediterranean system over multiple glacial–interglacial cycles. Through sedimentological and quantitative meiofauna (benthic foraminifera and ostracods) analyses of two reference cores (80 m and 100 m long) and stratigraphic correlations, a mid‐Pleistocene palaeovalley, 5 km wide and 50 m deep, was reconstructed. Whereas valley filling is chronologically constrained to the penultimate interglacial (Marine Isotope Stage 7) by four electron spin resonance ages on bivalve shells (Cerastoderma glaucum), its incision is tentatively correlated with the Marine Isotope Stage 8 sea‐level fall. Above basal fluvial‐channel gravels, the incised valley fill is formed by a mud‐prone succession, up to 44 m thick, formed by a lower floodplain unit and an upper unit with brackish meiofauna that reflects the development of a wave‐dominated estuary. Subtle meiofauna changes towards less confined conditions record two marine flooding episodes, chronologically linked to the internal Marine Isotope Stage 7 climate‐eustatic variability. After the maximum transgressive phase, recorded by coastal sands, the interfluves were flooded around 200 ka (latest Marine Isotope Stage 7). The subsequent shift in river incision patterns, possibly driven by neotectonic activity, prevented valley reoccupation guiding the northward formation of the Last Glacial Maximum palaeovalley. The applied multivariate approach allowed the sedimentological characterization of the Marine Isotope Stage 7 and Marine Isotope Stage 1 palaeovalley fills, including shape, size and facies architecture, which revealed a consistent river‐coastal system response over two non‐consecutive glacial–interglacial cycles (Marine Isotope Stages 8 to 7 and Marine Isotope Stages 2 to 1). The recurring stacking pattern of facies documents a predominant control exerted on stratigraphy by Milankovitch and sub‐Milankovitch glacio‐eustatic oscillations across the late Quaternary period.     

Upper Cretaceous marine‐continental transition (Leonese Area,NW Spain) defined from integrated outcrop and seismic stratigraphy

The Upper Cretaceous succession of the Leonese Area (NW Spain) comprises mixed clastic and carbonate sediments. This succession is divided into two lithostratigraphic units, the Voznuevo Member and the Boñar Formation, which represent fluvial, shoreface, intertidal, subtidal and open‐shelf sedimentary environments. Regional seismic interpretation and sequence stratigraphic analysis have allowed the study of lateral and vertical changes in the sedimentary record and the definition of third‐order levels of stratigraphic cyclicity. On the basis of these data, the succession can be divided into two second‐order depositional sequences (DS‐1 and DS‐2), incorporating three system tracts in a lowstand to transgressive to highstand system tract succession (LST–TST–HST). These sequences are composed of fluvial systems at the base with palaeocurrents that flowed westward and south‐westward. The upper part of DS‐1 (Late Albian–Middle Turonian) shows evidence of intertidal to subtidal and offshore deposits. DS‐2 (Late Turonian–Campanian) comprises intertidal to subtidal, tidal flat, shallow marine and lacustrine deposits and interbedded fluvial deposits. Two regressive–transgressive cycles occurred in the area related to eustatic controls. The evolution of the basin can be explained by base‐level changes and associated shifts in depositional trends of successive retrogradational episodes. By using isobath and isopach maps, the main palaeogeographic features of DS‐1 and DS‐2 were constrained, namely coastline positions, the existence and orientation of corridors through which fluvial networks were channelled and the location of the main depocentres of the basin. Sedimentation on the Upper Cretaceous marine platform was mainly controlled by (i) oscillations of sea level and (ii) the orientation of Mesozoic faults, which induced sedimentation along depocentres. Copyright © 2013 John Wiley & Sons, Ltd.     

Redbed‐associated sabkhas and tidal flats at Shark Bay,Western Australia: Their significance for genetic models of stratiform Cu‐(Pb‐Zn) deposits

Interaction of metalliferous continental brines with biogenic sulphide is the basis of some syngenetic and early diagenetic models for the formation of Cu‐(Pb‐Zn) sulphides during a depositional cycle of carbonates in restricted marine environments. A variation of these models (an ‘evaporative concentration‐lateral groundwater flow’ model) is proposed, using hydrological, geochemical and biological data from low metal, but otherwise pertinent redbed‐associated, sabkha, tidal flat and subtidal environments at Nilemah Embayment, in Hamelin Pool (Shark Bay, Western Australia).

The model is constrained by: (i) the short time available for ore accumulation during a single depositional cycle; (ii) limitation of adequate rates of bacterial sulphate reduction for the formation of an ore deposit to near‐surface sediments; (iii) restriction of the most favourable ore‐forming sites to the intertidal zone and the littoral shelf; (iv) coincidence in these sites of laterally‐flowing marine/meteoric groundwater brine, and mosaics of in situ cyanobacterial mats and shallow erosional depressions containing detrital organic matter eroded from the mats. Under these conditions the metalliferous fluid would have to contain about 1000 ppm Cu and flow for 1000 years at a rate of 5 m/a through the intertidal/littoral shelf environment to produce an ore deposit.

Critical features of a model that could generate this combination of very high metal concentrations and flow rates are: (i) a highly permeable unconfined aquifer system comprising alluvial fans at the base of basaltic mountain ranges and continental redbeds beneath a broad coastal plain; (ii) mobilization, concentration and transport of the metals in this aquifer to intertidal/littoral shelf sites of ore deposition; (iii) effective concentration processes in the aquifer, involving evaporation and reflux of brines in groundwater discharge areas on the coastal plain and evaporation in marine‐continental and marine sabkhas bordering the sites of deposition; (iv) rapid lateral groundwater flow of the concentrated metalliferous brines under a strong seawards‐directed hydraulic gradient; and (v) discharge of the metalliferous brines into or through topographic depressions generated by erosion and shoaling in the peritidal and littoral shelf environments.

The model hydrodynamic processes and their magnitude are within the range observed in modern environments but they are most likely to be effective in coarse‐grained, topographically irregular carbonate sabkhas and tidal flats, which usually form under high‐energy conditions. Even under these conditions, the individual ore‐forming processes must combine in an optimum manner before the highly demanding metal concentrations and flow rates required for ore formation in a single marine depositional cycle can be met.