Modelling hydrodynamics in the Rio Parana, Argentina : an evaluation and inter-comparison of reduced-complexity and physics based models applied to a large sand-bed river

Depth-averaged velocities and unit discharges within a 30 km reach of one of the world's largest rivers, the Rio Parana, Argentina, were simulated using three hydrodynamic models with different process representations: a reduced complexity (RC) model that neglects most of the physics governing fluid flow, a two-dimensional model based on the shallow water equations, and a three-dimensional model based on the Reynolds-averaged Navier-Stokes equations. Row characteristics simulated using all three models were compared with data obtained by acoustic Doppler current profiler surveys at four cross sections within the study reach. This analysis demonstrates that, surprisingly, the performance of the RC model is generally equal to, and in some instances better than, that of the physics based models in terms of the statistical agreement between simulated and measured flow properties. In addition, in contrast to previous applications of RC models, the present study demonstrates that the RC model can successfully predict measured flow velocities. The strong performance of the RC model reflects, in part, the simplicity of the depth-averaged mean flow patterns within the study reach and the dominant role of channel-scale topographic features in controlling the flow dynamics. Moreover, the very low water surface slopes that typify large sand-bed rivers enable flow depths to be estimated reliably in the RC model using a simple fixed-lid planar water surface approximation. This approach overcomes a major problem encountered in the application of RC models in environments characterised by shallow flows and steep bed gradients. The RC model is four orders of magnitude faster than the physics based models when performing steady-state hydrodynamic calculations. However, the iterative nature of the RC model calculations implies a reduction in computational efficiency relative to some other RC models. A further implication of this is that, if used to simulate channel morphodynamics, the present RC model may offer only a marginal advantage in terms of computational efficiency over approaches based on the shallow water equations. These observations illustrate the trade off between model realism and efficiency that is a key consideration in RC modelling. Moreover, this outcome highlights a need to rethink the use of RC morphodynamic models in fluvial geomorphology and to move away from existing grid-based approaches, such as the popular cellular automata (CA) models, that remain essentially reductionist in nature. In the case of the world's largest sand-bed rivers, this might be achieved by implementing the RC model outlined here as one element within a hierarchical modelling framework that would enable computationally efficient simulation of the morphodynamics of large rivers over millennial time scales. (C) 2012 Elsevier B.V. All rights reserved.

Scales and causes of heterogeneity in bars in a large multi-channel river: Rio Parana, 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 Rio Parana, Argentina. The investigated bars are located between 30km upstream and 540km downstream of the Rio Parana - Rio 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 600m and thicknesses up to 12m, 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 Rio Parana 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 Rio Paraguay. Relative to other controls on downstream fining, the tributary input of fine-grained suspended material from the Rio 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 5m 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 Rio Parana are composed of dune sets, even directly downstream from the Rio 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.

Late Jurassic continental flood basalt doleritic dykes in northwestern Cuba: remnants of the Gulf of Mexico opening

Accreted terranes, comprising a wide variety of Late Jurassic and Early Cretaceous igneous and sedimentary rocks are an important feature of Cuban geology. Their characterization is helpful for understanding Caribbean paleogeography. The Guaniguanico terrane (western Cuba) is formed by upper Jurassic platform sediments intruded by microgranular dolerite dykes. The geochemical characteristics of the dolerite whole rock samples and their minerals (augitic clinopyroxene, labradorite and andesine) are consistent with a tholeiitic affinity. Major and trace element concentrations as well as Nd, Sr and Pb isotopes show that these rocks also have a continental affinity. Sample chemistry indicates that these lavas are similar to a low Ti-P2O5 (LTi) variety of continental flood basalts (CFB) similar to the dolerites of Ferrar (Tasmania). They derived from mixing of a lithospheric mantle Source and an asthenopheric component similar to E-MORB with minor markers of crustal contamination and sediment assimilation. However, the small quantity of Cuban magmatic rocks, similarly to Tasmania, Antarctica and Siberia differs from other volumetrically important CFB occurrences Such as Parana and Deccan. These dolerites are dated as 165-150 Ma and were emplaced during the separation of the Yucatan block from South America. They could in fact be part of the Yucatan-South America margin through which the intrusive system was emplaced and which was later accreted to the Cretaceous arc of central Cuba and to the Palaeogene arc of eastern Cuba. These samples could therefore reflect the pre-rift stage between North and South America and the opening of the gulf of Mexico.

Climate and sea-level variations along the northwestern Tethyan margin during the Valanginian C-isotope excursion: Mineralogical evidence from the Vocontian Basin (SE France)

A high resolution mineralogical study (bulk-rock and clay-fraction) was carried out upon the hemipelagic strata of the Angles section (Vocontian Basin, SE France) in which the Valanginian positive C-isotope excursion occurs. To investigate sea-level fluctuations and climate change respectively, a Detrital Index (DI: (phyllosilicates and quartz)/calcite) and a Weathering Index (WI: kaolinite/(illite + chlorite)) were established and compared to second-order sea-level fluctuations. In addition, the mineralogical data were compared with the High Nutrient Index (HNI, based on calcareous nannofossil taxa) data obtained by Duchamp-Alphonse et al. (2007), in order to assess the link between the hydrolysis conditions recorded on the surrounding continents and the trophic conditions inferred for the Vocontian Basin. It appears that the mineralogical distribution along the northwestern Tethyan margin is mainly influenced by sea-level changes during the Early Valanginian (Pertransiens to Stephanophorus ammonite Zones) and by climate variations from the late Early Valanginian to the base of the Hauterivian (top of the Stephanophorus to the Radiatus ammonite Zones). The sea-level fall observed in the Pertransiens ammonite Zone (Early Valanginian) is well expressed by an increase in detrital inputs (an increase in the DI) associated with a more proximal source and a shallower marine environment, whereas the sea-level rise recorded in the Stephanophorus ammonite Zone corresponds to a decrease in detrital influx (a decrease in the DI) as the source becomes more distal and the environment deeper. Interpretation of both DI and WI, indicates that the positive C-isotope excursion (top of the Stephanophorus to the Verrucosum ammonite Zones) is associated with an increase of detrital inputs under a stable, warm and humid climate, probably related to greenhouse conditions, the strongest hydrolysis conditions being reached at the maximum of the positive C-isotope excursion. From the Verrucosum ammonite Zone to the base of the Hauterivian (Radiatus ammonite Zone) climatic conditions evolved from weak hydrolysis conditions and, most likely, a cooler climate (resulting in a decrease in detrital inputs) to a seasonal climate in which more humid seasons alternated with more arid ones. The comparison of the WI to the HNI shows that the nutrification recorded al: the Angles section from the top of the Stephanophorus to the Radiatus ammonite Zones (including the positive C-isotope shift), is associated with climatic changes in the source areas. At that time, increased nutrient inputs were generally triggered by increased weathering processes in the source areas due to acceleration in the hydrological cycle under greenhouse conditions This scenario accords with the widely questioned palaeoenvironmental model proposed by Lini et al., (1992) and suggests that increasing greenhouse conditions are the main factor that drove the palaeoenvironmental changes observed in the hemipelagic realm of the Vocontian Basin, during the Valanginian positive C-isotope shift. This high-resolution mineralogical study highlights short-term climatic changes during the Valanginian, probably associated to rapid changes in the C-cycle. Coeval Massive Parana-Etendeka flood basalt eruptions may explain such rapid perturbations. (C) 2011 Elsevier B.V. All rights reserved.