Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Little ice age advance and retreat sediment budgets for an outlet glacier in western Norway

Bødalsbreen is an outlet glacier of the Jostedalsbreen Ice Field in western Norway. Nine moraine ridges formed during and after the maximum extent of the Little Ice Age (LIA). The stratigraphy of proglacial sediments in the Bødalen basin inside the LIA moraines is examined, and corresponding sediment volumes are calculated based on georadar surveys and seismic profiling. The total erosion rates (etot) by the glacier are determined for the periods AD 1650?1930 and AD 1930?2005 as 0.8 ± 0.4 mm/yr and 0.7 ± 0.3 mm/yr, respectively. These rates are based on the total amount of sediment delivered to the glacier margin. The values are almost one order of magnitude higher than total erosion rates previously calculated for Norwegian glaciers. This is explained by the large amount of pre-existing sediment that was recycled by Bødalsbreen. Thus, the total erosion rate must be considered as a composite of eroded bedrock and of removed pre-existing sediments. The total erosion rate is likely to vary with time owing to a decreasing volume of easily erodible, unconsolidated sediment and till under the glacier. A slight increase in the subglacial bedrock erosion is expected owing to the gradually increasing bedrock surface area exposed to subglacial erosion.

Biogeochemical weathering in sedimentary chronosequences of the Rhône and Oberaar Glaciers (Swiss Alps): rates and mechanisms of biotite weathering

We analysed the composition of phyllosilicate minerals in sediments deposited by the Rhone and Oberaar glaciers (Swiss Alps), in order to identify processes and rates of biogeochemical weathering in relation to glacial erosion. The investigated sediments are part of chronosequences consisting of (A) suspended, "fresh" sediment in melt water; (B) terminal moraines from the Little Ice Age (LIA; approximately 1560-1850); and (C) tilts of the Younger Dryas interval (YD; approximately 11'500y BP). Secondary weathering products associated with the suspended sediment have not been observed: we therefore exclude intermittent subglacial storage and weathering of this material and assume that the suspended sediment is directly derived from mechanically abraded bedrock. This implies that biogeochemical weathering processes started once the glacially-derived sediment was deposited in the proglacial area. The combination of a developing vegetation cover, the generally high permeability allowing the percolation of precipitation, and the chemical reactivity related to the dominance of fine-grained material (<63 pm) drives the weathering process and the initial Umbrepts present in LIA profiles undergo podzolisation and lead to the formation of Humods observed in YD profiles. Systematic XRD analyses of these chronosequences show a progressive decrease in biotite contents and a concomitant increase in pedogenically formed vermiculite with increasing sediment age. Biotite contents decrease by 25-50% in the upper 30 cm of the moraines after 145-275 yr in the proglacial environment. Biotite weathering rates are calculated using the difference in the biotite content between unweathered and weathered glacial sediments within the investigated profiles. The reactive mineral surface area is estimated geometrically, both with regards to the total relative surface (WRT) as well as to the relative edge surface (WRE). WRT Biotite weathering rates are estimated as 10(-13)-10-(15) mol(biotite) m(biotite)(-2) s(-1). WRE Biotite weathering rates are on the order of 10(-13)-10(-14) mol(biotite) m(biotite)(-2) s(-1). Biotite weathering rates obtained by this study are in the order of one magnitude higher in comparison to other published field-based weathering rates. Using biotite as an indicator, we therefore suggest that glacially-derived material in the area of the Oberaar and Rhone glaciers is generally subjected to enhanced biogeochemical weathering, starting immediately after deposition in the proglacial zone and subsequently continuing for thousands of years after glacier retreat.

Weathering and the mobility of phosphorus in the catchments and forefields of the Rhone and Oberaar glaciers, central Switzerland: Implications for the global phosphorus cycle on glacial-interglacial timescales

In this study we evaluate the dynamics of the biophile element phosphorus (P) in the catchment and proglacial areas of the Rhone and Oberaar glaciers (central Switzerland). We analysed erosion and dissolution rates of P-containing minerals in the subglacial environment by sampling water and suspended sediment in glacier outlets during three ablation and two accumulation seasons. We also quantified biogeochemical weathering rates of detrital P in proglacial sedimentary deposits using two chronosequences of samples of fresh, suspended, material obtained from the Oberaar and Rhone water outlets, Little-Ice-Age (LIA) moraines and Younger Dryas (YD) tills in each catchment. Subglacial P weathering is mainly a physical process and detrital P represents more than 99%, of the precipitation-corrected total P denudation flux (234 and 540 kg km(-2) yr(-1) for the Rhone and Oberaar catchments, respectively). The calculated detrital P flux rates are three to almost five times higher than the world average flux. The precipitation-corrected soluble reactive P (SRP) flux corresponds to 1.88-1.99 kg km(-2) yr(-1) (Rhone) and 2.12-2.44 kg km(-2) yr(-1) (Oberaar), respectively. These fluxes are comparable to those of tropical rivers draining transport-limited, tectonically inactive weathering areas. In order to evaluate the efficiency of detrital P weathering in the Rhone and Oberaar proglacial areas, we systematically graded apatite grains extracted from the chronosequence in each catchment relative to weathering-induced changes in their surface morphologies (grades 1-4). Fresh apatite grains are heavily indented and dissolution rounded (grade 1). LIA grains from two 0-10 cm deep moraine samples show extensive dissolution etching, similar to surface grains from the YD profile (mean grades 2.7, 3.5 and 3.5, respectively). In these proglacial deposits, the weathering front deepens progressively as a function of time due to biocorrosion in the evolving acidic pedosphere, with mechanical indentations on grains acting as sites of preferential dissolution. We also measured iron-bound, organic and detrital P concentrations in the chronosequence and show that organic and iron-bound P has almost completely replaced detrital P in the top layers of the YD profiles. Detrital P weathering rates are calculated as 3 10 and 280 kg km(-2) yr(-1) for LIA moraines and 10 kg km(-2) yr(-1) for YD tills. During the first 300 years of glacial sediment exposure P dissolution rates are shown to be approximately 70 times higher than the mean global dissolved P flux from ice-free continents. After 11.6 kyr the flux is 2.5 times the global mean. These data strengthen the argument for substantial changes in the global dissolved P flux on glacial-interglacial timescales. A crude extrapolation from the data described here suggests that the global dissolved P flux may increase by 40-45% during the first few hundred years of a deglaciation phase

Topographic forcing of flow partition and flow structures at river bifurcations

This paper presents the predicted flow dynamics from the application of a Reynolds-averaged NavierStokes model to a series of bifurcation geometries with morphologies measured during previous flume experiments. The topography of the bifurcations consists of either plane or bedform-dominated beds which may or may not possess discordance between the two bifurcation distributaries. Numerical predictions are compared with experimental results to assess the ability of the numerical model to reproduce the division of flow into the bifurcation distributaries. The hydrodynamic model predicts: (1) diverting fluxes in the upstream channel which direct water into the distributaries; (2) super-elevation of the free surface induced at the bifurcation edge by pressure differences; and (3) counter-rotating secondary circulation cells which develop upstream of the apex of the bifurcation and move into the downstream channels, with water converging at the surface and diverging at the bed. When bedforms are not present, weak transversal fluxes characterize the upstream channel for almost its entire length, associated with clearly distinguishable secondary circulation cells, although these may be under-estimated by the turbulence model used in the solution. In the bedform dominated case, the same hydrodynamic conditions were not observed, with the bifurcation influence restricted and depth scale secondary circulation cells not forming. The results also demonstrate the dominant effect bed discordance has upon flow division between the two distributaries. Finally, results indicate that in bedform dominated rivers. Consequently, we suggest that sand-bed river bifurcations are more likely to have an influence that extends much further upstream and have a greater impact upon water distribution. This may contribute to observed morphological differences between sand-bedded and gravel-bedded braided river networks. Copyright (C) 2012 John Wiley & Sons, Ltd.

Étude des sédiments quaternaires, de la molasse et sa tectonique, dans le Grand Lac (Léman) à partir de données sismiques 2D et 3D

RESUME L'Institut de Géophysique de l'Université de Lausanne a développé au cours de ces dernières années un système d'acquisition de sismique réflexion multitrace à haute résolution 2D et 3D. L'objectif de cette thèse était de poursuivre ce développement tout améliorant les connaissances de la géologie sous le lac Léman, en étudiant en particulier la configuration des grands accidents sous-lacustres dans la Molasse (Tertiaire) qui forme l'essentiel du substratum des formations quaternaires. En configuration 2D, notre système permet d'acquérir des profils sismiques avec une distance inter-CDP de 1,25 m. La couverture varie entre 6 et 18 selon le nombre de traces et la distance inter-tir. Le canon à air (15/15 eu. in.), offre une résolution verticale de 1,25 ni et une pénétration maximale de 300 m sous le fond de l'eau. Nous avons acquis au total plus de 400 km de sections 2D dans le Grand Lac et le Haut Lac entre octobre 2000 et juillet 2004. Une campagne de sismique 3D a fourni des données au large d'Evian sur une surface de 442,5 m sur 1450 m, soit 0,64 km2. La navigation ainsi que le positionnement des hydrophones et de la source ont été réalisés avec des GPS différentiels. Nous avons utilisé un traitement sismique conventionnel, sans appliquer d'AGC et en utilisant une migration post-stack. L'interprétation du substratum antéquaternaire est basée sur l'identification des sismofaciès, sur leurs relations avec les unités géologiques adjacentes au lac, ainsi que sur quelques données de forages. Nous obtenons ainsi une carte des unités géologiques dans le Grand Lac. Nous précisons la position du chevauchement subalpin entre la ville de Lausanne, sur la rive nord, et le bassin de Sciez, sur la rive sud. Dans la Molasse de Plateau, nous avons identifié les décrochements de Pontarlier et de St. Cergue ainsi que plusieurs failles non reconnues jusqu'ici. Nous avons cartographié les accidents qui affectent la Molasse subalpine ainsi que le plan de chevauchement du flysch sur la Molasse près de la rive sud du lac. Une nouvelle carte tectonique de la région lémanique a ainsi pu être dressée. L'analyse du substratum ne montre pas de failles suggérant une origine tectonique de la cuvette lémanique. Par contre, nous suggérons que la forme du creusement glaciaire, donc de la forme du lac Léman, a été influencée par la présence de failles dans le substratum antéquaternaire. L'analyse des sédiments quaternaires nous a permis de tracer des cartes des différentes interfaces ou unités qui les composent. La carte du toit du substratum antéquaternaire montre la présence de chenaux d'origine glaciaire dont la profondeur maximale atteint la cote -200 ni. Leur pente est dirigée vers le nord-est, à l'inverse du sens d'écoulement actuel des eaux. Nous expliquons cette observation par l'existence de circulations sous-glaciaires d'eau artésienne. Les sédiments glaciaires dont l'épaisseur maximale atteint 150 ni au centre du lac ont enregistré les différentes récurrences glaciaires. Dans la zone d'Evian, nous mettons en évidence la présence de lentilles de sédiments glaciolacustres perchées sur le flanc de la cuvette lémanique. Nous avons corrélé ces unités avec des données de forage et concluons qu'il s'agit du complexe inférieur de la pile sédimentaire d'Evian. Celui-ci, âgé de plus de 30 000 ans, serait un dépôt de Kame associé à un lac périglaciaire. La sismique réflexion 3D permet de préciser l'orientation de l'alimentation en matériel détritique de l'unité. La finesse des images obtenues nous permet également d'établir quels types d'érosion ont affecté certaines unités. Les sédiments lacustres, dont l'épaisseur maximale imagée atteint plus de 225 m et sans doute 400 ni sous le delta du Rhône, indiquent plusieurs mécanismes de dépôts. A la base, une mégaturbidite, épaisse d'une trentaine de mètres en moyenne, s'étend entre l'embouchure de la Dranse et le delta du Rhône. Au-dessus, la décantation des particules en suspension d'origine biologique et détritique fournit l'essentiel des sédiments. Dans la partie orientale du lac, les apports détritiques du Rhône forment un delta qui prograde vers l'ouest en s'imbriquant avec les sédiments déposés par décantation. La structure superficielle du delta a brutalement évolué, probablement à la suite de l'évènement catastrophique du Tauredunum (563 A.D.). Sa trace probable se marque par la présence d'une surface érosive que nous avons cartographiée. Le delta a ensuite changé de géométrie, avec notamment un déplacement des chenaux sous-lacustres. Sur l'ensemble de nos sections sismiques, nous n'observons aucune faille dans les sédiments quaternaires qui attesterait d'une tectonique postglaciaire du substratum. ABSTRACT During the last few years the institute of Geophysics of the University of Lausanne cleveloped a 2D and 3D high-resolution multichannel seismic reflection acquisition system. The objective of the present work was to carry on this development white improving our knowledge of the geology under Lake Geneva, in particular by studying the configuration of the large accidents affecting the Tertiary Molasse that makes up the basement of most Quaternary deposits. In its 2D configuration, our system makes it possible to acquire seismic profiles with a CDP interval of 1.25 m. The fold varies from 6 to 18 depending on the number of traces and the shooting interval. Our air gun (15/15 cu. in.) provides a vertical resolution of 1.25 m and a maximum penetration depth of approximately 300 m under water bottom. We acquired more than 400 km of 2D sections in the Grand Lac and the Haut Lac between October 2000 and July 2004. A 3D seismic survey off the city of Evian provided data on a surface of 442.5 m x 1450 m (0.64 km2). Ship's navigation as well as hydrophone- and source positioning were carried out with differential GPS. The seismic data were processed following a conventional sequence without .applying AGC and using post-stack migration. The interpretation of the pre-Quaternary substratum is based on sismofacies, on their relationships with terrestrial geological units and on some borehole data. We thus obtained a map of the geological units in the Grand Lac. We defined the location of the subalpine thrust from Lausanne, on the north shore, to the Sciez Basin, on the south shore. Within the Molasse de Plateau, we identified the already know Pontarlier and St Cergue transforms Fault as well as faults. We mapped faults that affect subalpine Molasse as well as the thrust fault plane between alpine flysch and Molasse near the lake's south shore. A new tectonic map of the Lake Geneva region could thus be drawn up. The substratum does not show faults indicating a tectonic origin for the Lake Geneva Basin. However, we suggest that the orientation of glacial erosion, and thus the shape of Lake Geneva, vas influenced by the presence of faults in the pre-Quaternary basement. The analysis of Quaternary sediments enabled us to draw up maps of various discontinuities or internal units. The top pre-Quaternary basement map shows channels of glacial origin, the deepest of them reaching an altitude of 200 m a.s.l. The channel's slopes are directed to the North-East, in opposite direction of the present water flow. We explain this observation by the presence of artesian subglacial water circulation. Glacial sediments, the maximum thickness of which reaches 150 m in the central part of the lake, record several glacial recurrences. In the Evian area, we found lenses of glacio-lacustrine sediments set high up on the flank of the Lake Geneva Bassin. We correlated these units with on-land borehole data and concluded that they represent the lower complex of the Evian sedimentary pile. The lower complex is aider than 30 000 years, and it could be a Kame deposit associated with a periglacial lake. Our 3D seismic reflexion survey enables us to specify the supply direction of detrital material in this unit. With detailed seismic images we established how some units were affected by different erosion types. The lacustrine sediments we imaged in Lake Geneva are thicker than 225 m and 400 m or more Linder the Rhone Delta. They indicate several depositional mechanisms. Their base is a major turbidite, thirty meters thick on average, that spreads between the Dranse mouth and the Rhone delta. Above this unit, settling of suspended biological and detrital particles provides most of the sediments. In the eastern part of the lake, detrital contribution from the Rhone builds a delta that progrades to the west and imbricates with the settling sediments. The shallow structure of the Rhone delta abruptly evolved, probably after the catastrophic Tauredunum event (563 A.D.). It probably coincides with an erosive surface that we mapped. As a result, the delta geometry changed, in particular associated with a displacement of water bottom channels. In all our seismic sections, we do not observe fault in the Quaternary sediments that would attest postglacial tectonic activity in the basement.

Lidar measurement of surface melt for a temperate Alpine glacier at the seasonal and hourly scales

This study shows how a new generation of terrestrial laser scanners can be used to investigate glacier surface ablation and other elements of glacial hydrodynamics at exceptionally high spatial and temporal resolution. The study area is an Alpine valley glacier, Haut Glacier d'Arolla, Switzerland. Here we use an ultra-long-range lidar RIEGL VZ-6000 scanner, having a laser specifically designed for measurement of snow- and ice-cover surfaces. We focus on two timescales: seasonal and daily. Our results show that a near-infrared scanning laser system can provide high-precision elevation change and ablation data from long ranges, and over relatively large sections of the glacier surface. We use it to quantify spatial variations in the patterns of surface melt at the seasonal scale, as controlled by both aspect and differential debris cover. At the daily scale, we quantify the effects of ogive-related differences in ice surface debris content on spatial patterns of ablation. Daily scale measurements point to possible hydraulic jacking of the glacier associated with short-term water pressure rises. This latter demonstration shows that this type of lidar may be used to address subglacial hydrologic questions, in addition to motion and ablation measurements.