Scales and causes of heterogeneity in bars in a large multi‐channel river: Río Paraná, Argentina

To date, published studies of alluvial bar architecture in large rivers have been restricted mostly to case studies of individual bars and single locations. Relatively little is known about how the depositional processes and sedimentary architecture of kilometre‐scale bars vary within a multi‐kilometre reach or over several hundreds of kilometres downstream. This study presents Ground Penetrating Radar and core data from 11, kilometre‐scale bars from the Río Paraná, Argentina. The investigated bars are located between 30 km upstream and 540 km downstream of the Río Paraná – Río Paraguay confluence, where a significant volume of fine‐grained suspended sediment is introduced into the network. Bar‐scale cross‐stratified sets, with lengths and widths up to 600 m and thicknesses up to 12 m, enable the distinction of large river deposits from stacked deposits of smaller rivers, but are only present in half the surface area of the bars. Up to 90% of bar‐scale sets are found on top of finer‐grained ripple‐laminated bar‐trough deposits. Bar‐scale sets make up as much as 58% of the volume of the deposits in small, incipient mid‐channel bars, but this proportion decreases significantly with increasing age and size of the bars. Contrary to what might be expected, a significant proportion of the sedimentary structures found in the Río Paraná is similar in scale to those found in much smaller rivers. In other words, large river deposits are not always characterized by big structures that allow a simple interpretation of river scale. However, the large scale of the depositional units in big rivers causes small‐scale structures, such as ripple sets, to be grouped into thicker cosets, which indicate river scale even when no obvious large‐scale sets are present. The results also show that the composition of bars differs between the studied reaches upstream and downstream of the confluence with the Río Paraguay. Relative to other controls on downstream fining, the tributary input of fine‐grained suspended material from the Río Paraguay causes a marked change in the composition of the bar deposits. Compared to the upstream reaches, the sedimentary architecture of the downstream reaches in the top ca 5 m of mid‐channel bars shows: (i) an increase in the abundance and thickness (up to metre‐scale) of laterally extensive (hundreds of metres) fine‐grained layers; (ii) an increase in the percentage of deposits comprised of ripple sets (to >40% in the upper bar deposits); and (iii) an increase in bar‐trough deposits and a corresponding decrease in bar‐scale cross‐strata (<10%). The thalweg deposits of the Río Paraná are composed of dune sets, even directly downstream from the Río Paraguay where the upper channel deposits are dominantly fine‐grained. Thus, the change in sedimentary facies due to a tributary point‐source of fine‐grained sediment is primarily expressed in the composition of the upper bar deposits.     

Sources, distribution and storage of heavy metals in the Río Pilcomayo, Bolivia

The Río Pilcomayo rises on the Cerro Rico de Potosí precious metal-polymetallic tin deposits of Bolivia, and flows in a southeasterly direction for ca. 600 km to Bolivia's southern border with Argentina. Mining of the Potosí deposits has occurred continuously since 1545, generating large quantities of waste materials in the headwater of the basin. In addition, a tailings dam breach at the Porco mine in 1996 released an estimated 235 000 m³ of tailings and fluid into the upper reaches of the Río Pilaya, the largest tributary to the Pilcomayo.Concentrations of As, Sb, Cd, Cu, Pb, Hg, Ag, Tl and Zn in contemporary channel sediments upstream of the Pilaya confluence are significantly elevated above background values. Elevated levels appear to be associated with pyrite- and other sulphide mineral-bearing tailings materials transported more than 200 km downstream of the Potosí mines. Significant downstream declines in elemental concentrations occur within 15 km, and again between 150 and 200 km, from the mines. The initial decrease in concentrations is due to the rapid dilution of nearly pure tailings effluent released to the river from milling facilities near Potosí. The latter decrease results from a combination of geomorphic processes including the storage of sediment-borne metals within the channel bed and the influx of ‘clean’ sediment from several large tributaries. Downstream of the Pilaya confluence, concentrations of Cu, Pb, Hg and Zn are only slightly elevated above background values, and Ag, Cd, Sb and Tl cannot be distinguished from background levels. These data suggest that while the Porco tailings spill may have had a significant short-term impact on sediment and water quality along the lower reaches of the Río Pilcomayo, its longer-term impacts were limited. Metals stored and eroded from alluvial deposits of historical age in upstream reaches appear to be an important source of metals to the river today. An additional, and perhaps more significant source, is the release of tailings effluent to the river from modern milling operations. The transport of these contaminants downstream of Icla (203 km from Potosí) appears to be restricted by aggradational processes occurring in the vicinity of Puente Sucre. In addition, downstream of the confluence of the Río Pilaya, inputs of large amounts of ‘clean’ sediment have caused dilution of the metal contaminants. Data from other studies where similar geomorphic processes have occurred suggest that the metals in the upper Pilcomayo may eventually be moved downvalley as the aggradational processes are reversed and channel stabilisation occurs. Thus, the most significant impacts of metal contamination may not be realised in downstream areas for decades.     

Architecture and development of the alluvial sediments of the Upper Jurassic Tordillo Formation in the Cañada Ancha Valley, northern Neuquén Basin, Argentina

The Upper Jurassic Tordillo Formation at Cañada Ancha area, northern Neuquén Basin, Argentina, comprises a multi-stage suit of predominantly alluvial sediments that is heterolithic in nature. In that suit, several lithofacies, architectural elements, and bounding surfaces of different order have been identified and their lateral and vertical distribution characterized. This analysis allowed the differentiation of 3 main units (lower, middle and upper), 20 subunits (C-1 to C-20), and the characterization of their alluvial styles.The lower unit (which comprises subunits C-1 to C-4) is mainly formed by fine- to medium-grained sandstones, which become medium- to coarse-grained towards the top. These sandstones characterize settings ranging from floodplains with isolated, unconfined flows, to more complex, vertically stacked, multi-storey sheet sandstones of braided fluvial systems. The middle unit (C-5 to C-10) is dominated by pale brown-grey fine-to coarse-grained sands and medium size subangular to angular conglomerates, which reflect amalgamated complexes of sandstone sheets and downstream accretion macroforms. Remarkably, this alluvial sedimentation was episodically punctuated by volcaniclastic flows. The upper unit (C-11 to C-20) consists of finer sediments, mainly pink to white fine-to medium grained sandstones and red to green siltstones. Towards the top, bioturbation becomes important, and also the presence of volcanosedimentary flows is noticeable. Fluvial settings include braided sheet sandstones with waning flood deposits evolving to isolated high-sinuosity fluvial systems, with flash flood deposits. At the top of this unit, facies may suggest marine influence.Vertical changes in the fluvial style result from both climatic and tectonic controls. A semiarid to arid climate and the active tectonism linked to the eastward migration of the Andean volcanic arc determined major bounding surfaces, fluvial style evolution and the presence of the volcano-sedimentary deposits. Different stages of high and low subsidence rates has been deduced from the vertical stacking of sediments.     

Fine‐grained meandering systems of the Lower Permian Clear Fork Formation of north‐central Texas,USA: Lateral and oblique accretion on an arid plain

Facies models that adequately represent the diverse range of fine‐grained fluvial systems are currently lacking from the literature. In this paper, the spectrum of these systems on the arid plains of western equatorial Pangea is explored, as well as the source and nature of the fine‐grained sediments. Eight fluvial elements in the Early Permian Clear Fork Formation of north‐central Texas represent channel systems up to 7 m deep with coarse basal deposits, three types of lateral‐accretion deposits and sandstone sheets, with laminated, disrupted and massive mudstones laid down in abandoned channels and on floodplains. The three fine‐grained fluvial styles represent a continuum between two end‐members: sustained lateral accretion of bedload composed of quartzose sediments and mud aggregates on point bars, and oblique accretion of suspended sediment on steep accretionary benches and banks with limited lateral migration. This spectrum is controlled, in part, by grain size and the proportion of suspended to bedload sediments. The presence of rarely documented swept ripples on exhumed accretion surfaces is attributed to rapid decline in water levels and downstream re‐entry of overbank floodwaters into the channel. Rill casts, roots and disrupted mudstones low down in channel bodies indicate periods of near‐dryness. Laterally extensive sheet sandstones were formed by episodic flows in broad, sandbed channels. The fluvial sediments were primarily intrabasinally sourced with extrabasinal sediments brought in during major floods from upland source areas or reworked from local storage in the basin, representing a supply limited system. The upward change in cement composition from mainly calcite and ankerite to dolomite and gypsum with minor celestine implies increasingly saline groundwater and progressive aridification, supporting Late Palaeozoic palaeoclimatic models. By integrating petrographic data with sedimentology, a plethora of information about ancient landscapes and climate is provided, allowing a fuller comparison between the Clear Fork Formation and modern dryland alluvial plains.     

Fission‐track dating of a tephra layer related to Poti‐Malal and Seguro drifts in the Río Grande basin,Mendoza, Argentina

Two old drift units called Poti‐Malal and Seguro have been differentiated in the Río Grande basin based on relative‐age criteria, stratigraphical relationships, morphology and fission‐track dating. A tephra dated at 0.226 ± 0.025 Ma was deposited on the Poti‐Malal drift and underlies the Seguro outwash, which is inferred to equate with marine oxygen isotope stage 6. The stratigraphical position and age suggest that the tephra post‐dates the Poti‐Malal glaciation and that it is older than the Seguro drift. The Poti‐Malal glaciation must be at least as old as Early–Middle Pleistocene, and the Seguro glaciation is assigned to the penultimate glaciation. The tephra unit may have been deposited during marine oxygen isotope stage 7. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

On the relationship between flow and suspended sediment transport over the crest of a sand dune, Río Paraná, Argentina

The links between large‐scale turbulence and the suspension of sediment over alluvial bedforms have generated considerable interest in the last few decades, with past studies illustrating the origin of such turbulence and its influence on flow resistance, sediment transport and bedform morphology. In this study of turbulence and sediment suspension over large sand dunes in the Río Paraná, Argentina, time series of three‐dimensional velocity, and at‐a‐point suspended sediment concentration and particle‐size, were measured with an acoustic Doppler current profiler and laser in situ scattering transmissometer, respectively. These time series were decomposed using wavelet analysis to investigate the scales of covariation of flow velocity and suspended sediment. The analysis reveals an inverse relationship between streamwise and vertical velocities over the dune crest, where streamwise flow deceleration is linked to the vertical flux of fluid towards the water surface in the form of large turbulent fluid ejections. Regions of high suspended sediment concentration are found to correlate well with such events. The frequencies of these turbulent events have been assessed from wavelet analysis and found to concentrate in two zones that closely match predictions from empirical equations. Such a finding suggests that a combination and interaction of vortex shedding and wake flapping/changing length of the lee‐side separation zone are the principal contributors to the turbulent flow field associated with such large alluvial sand dunes. Wavelet analysis provides insight upon the temporal and spatial evolution of these coherent flow structures, including information on the topology of dune‐related turbulent flow structures. At the flow stage investigated, the turbulent flow events, and their associated high suspended sediment concentrations, are seen to grow with height above the bed until a threshold height (ca 0·45 flow depth) is reached, above which they begin to decay and dissipate.     

Changing lake dynamics and sequence stratigraphy of synrift lacustrine strata in a half-graben: an example from the Triassic Ischigualasto–Villa Unión Basin, Argentina

Well‐exposed Triassic rift strata from the Ischigualasto–Villa Unión Basin (NW Argentina) include a 80 to ca 515 m thick lacustrine‐dominated package that can be correlated across a half‐graben using key stratigraphic surfaces (sequence boundaries, lacustrine flooding surfaces and forced regressive surfaces). The characteristics of the synrift lacustrine fill in different parts of the half‐graben have been examined and the mechanisms controlling sedimentation inferred. A variety of sedimentary environments are recognized including; volcaniclastic floodplain, mildly saline lake and playa lake, offshore lacustrine, delta front to fluvial‐dominated and wave‐dominated deltas, distributary and fluvial channel, and interdistributary bay. The succession can be divided into four stratigraphic sequences (SS1 to SS4), the oldest of which (SS1) contains volcaniclastic, fluvial and saline lake deposits; it is thickest close to the western border fault zone, reflecting more rapid subsidence here. Accommodation exceeded sediment and water input during SS1. The second and third sequences (SS2 and SS3) mark the onset of widespread lacustrine sedimentation, reflecting a balance between accommodation creation and water and sediment fluxes. Sequences SS2 and SS3 are represented by offshore meromictic lacustrine and deltaic deposits, the latter mostly sourced from the flexural and southern axial margins of the half‐graben. The presence of stacked parasequences bound by lacustrine flooding surfaces is related to climatically induced lake‐level fluctuations superimposed on variable rates of subsidence on the controlling rift border fault zone. The youngest sequence (SS4) is represented by the deposits of littoral lacustrine and shallow shelf deltas distinguished by a change in lithofacies, palaeocurrents and sandstone composition, suggesting a switch in sediment supply to the footwall margin to the NW. The change in the sediment source is related to reduced footwall uplift, the possible presence of a relay ramp and/or supply from a captured antecedent drainage network. During SS4, the rate of creation of accommodation was exceeded by the sediment and water discharge. The stratigraphic evolution of lacustrine strata in the half‐graben was mainly controlled by tectonic processes, including subsidence rate and the growth and evolution of the border fault zone, but changing climate (inducing changes in water balance and lake level) and autocyclic processes (delta lobe switching) were also important.     

From a river valley to a ria: Evolution of an incised valley (Ría de Ferrol,north‐west Spain) since the Last Glacial Maximum

The evolution of incised valleys is an important area of research due to the invaluable data it provides on sea‐level variations and depositional environments. In this article the sedimentary evolution of the Ría de Ferrol (north‐west Spain) from the Last Glacial Maximum to the present is reconstructed using a multidisciplinary approach, combining seismic and sedimentary facies, and supported by radiocarbon data and geochemical proxies to distinguish the elements of sedimentary architecture within the ria infill. The main objectives are: (i) to analyse the ria environment as a type of incised valley to evaluate the response of the system to the different drivers; (ii) to investigate the major controlling factors; and (iii) to explore the differentiation between rias and estuaries. As a consequence of the sea‐level rise subsequent to the Last Glacial Maximum (ca 20 kyr bp ), an extensive basin, drained by a braided palaeoriver, evolved into a tide‐dominated estuary and finally into a ria environment. Late Pleistocene and Holocene high‐frequency sea‐level variations were major factors that modulated the type of depositional environments and their evolution. Another major modulating factor was the antecedent morphology of the ria, with a rock‐incised narrow channel in the middle of the basin (the Ferrol Strait), which influenced the evolution of the ria as it became flooded during Holocene transgression. The strait acted as a rock‐bounded ‘tidal inlet’ enhancing the tidal erosion and deposition at both ends, i.e. with an ebb‐tidal delta in the outer sector and tidal sandbanks in the inner sector. The final step in the evolution of the incised valley into the modern‐defined ria system was driven by the last relative sea‐level rise (after 4 kyr bp ) when the river mouths retreated landward and a single palaeoriver was converted into minor rivers and streams with scattered mouths in an extensive coastal area.     

The role of tidal,wave and river currents in the evolution of mixed‐energy deltas: Example from the Lajas Formation (Argentina)

Many modern deltas show complex morphologies and architectures related to the interplay of river, wave and tidal currents. However, methods for extracting the signature of the individual processes from the stratigraphic architecture are poorly developed. Through an analysis of facies, palaeocurrents and stratigraphic stacking patterns in the Jurassic Lajas Formation, this paper: (i) separates the signals of wave, tide and river currents; (ii) illustrates the result of strong tidal reworking in the distal reaches of deltaic systems; and (iii) discusses the implications of this reworking for the evolution of mixed‐energy systems and their reservoir heterogeneities. The Lajas Formation, a sand‐rich, shallow‐marine, mixed‐energy deltaic system in the Neuquén Basin of Argentina, previously defined as a tide‐dominated system, presents an exceptional example of process variability at different scales. Tidal signals are predominantly located in the delta front, the subaqueous platform and the distributary channel deposits. Tidal currents vigorously reworked the delta front during transgressions, producing intensely cross‐stratified, sheet‐like, sandstone units. In the subaqueous platform, described for the first time in an ancient outcrop example, the tidal reworking was confined within subtidal channels. The intensive tidal reworking in the distal reaches of the regressive delta front could not have been predicted from knowledge of the coeval proximal reaches of the regressive delta front. The wave signals occur mainly in the shelf or shoreface deposits. The fluvial signals increase in abundance proximally but are always mixed with the other processes. The Lajas system is an unusual clean‐water (i.e. very little mud is present in the system), sand‐rich deltaic system, very different from the majority of mud‐rich, modern tide‐influenced examples. The sand‐rich character is a combination of source proximity, syndepositional tectonic activity and strong tidal‐current reworking, which produced amalgamated sandstone bodies in the delta‐front area, and a final stratigraphic record very different from the simple coarsening‐upward trends of river‐dominated and wave‐dominated delta fronts.