An investigation of geological and geochemical characteristics of late-Quaternay sediments in the Georgian Bay Region, Southern Ontario

Core samples of postglacial sediments and sediment surface samples from Shepherd Lake on the Bruce Peninsula, Harts Lake on the Canadian Shield, and two cores from Georgian Bay (core P-l in the western deep part and core P-7 in the eastern shallow part) have been analyzed for pH, grain size distribution, water content, bulk density, loss on ignition at 4500C and 11000 C, major oxides (Si02 ,A1203,!FeO,MgO,CaO, Na20,K20,Ti02 ,MnO and P205) and trace elements (Ba,Zr,Sr,y,S, Zn,Cu,Ni,Ce and Rb). The sediment in Georgian Bay are generally fine grained (fine silt to very fine silty clay) and the grain size decreases from the Canadian Shield (core p-7) towards the Bruce Peninsula (core P-l) along the assumed direction of sediment transport. This trend coincides with a decrease in sorting coefficient and an increase in roundness. Other physical characteristics, such as water content, bulk density and loss on ignition are positively correlated with the composition of sediments and their compaction, as well as with the energy of the depositional environment. Analyses of sediment surface samples from Shepherd Lake and Harts Lake indicate the influence of bedrock and surficial deposits in the watershed on pH condition that is also influenced by the organic matter content and probably I ! I man's activities. Organic matter content increases significantly in the surface sediment in these small lakes as a result of either natural eutrophication or anthropogenic organic loading. The extremely high organic matter content in Shepherd Lake sediment indicates rapid natural eutrophication in this closed basin and high biological productivity during postglacial time, probably due to high nutrient levels and shallow depth. The chemical composition of the Canadian Shield bedrock is positively correlated with the chemical characteristics of predominantly inorganic lake sediments that were derived from the Shield rocks by glacial abrasion and by postglacial weathering and erosion of both bedrock and surficial deposits. High correlation coefficients were found between organic matter in lake sediments and major oxides (Si02,AI203,.~FeO, MgO,CaO,K20 and MnO) , as well as some trace elements (Ba,Y, S,Zn,Cu,Ni and Rb). The chemical composition of sediments in Harts Lake and core P-7 in Georgian Bay on the Canadian Shield differs from the chemistry of sediments in Shepherd Lake and core P-l in Georgian Bay on the Bruce Peninsula. The difference between cores P-l and P-7 is indicated by values of Si02 , AI203 ,:LFeo,Mgo,CaO,Ba,Zr,Sr,y and S, and also by the organic matter content. This study indicates that the processes of sediment transport, depositional environment, weathering of the rocks and surficial deposits in the watershed, as well as chemical composition of source rocks all affect the chemical characteristics of lake sediments. The stratigraphic changes and variations in lake sediment chemistry with regard to major oxides, trace elements, and organic matter content are probably related to the history of glacial and postglacial lake stages of the Georgian Bay Region and, therefore, the geochemical data can make a useful contribution to a better understanding of the Late-Quaternary history of the Great Lakes.

Numerical modelling of the impact of climate change on the morphology of Saint-Lawrence tributaries

Cette thèse examine les impacts sur la morphologie des tributaires du fleuve Saint-Laurent des changements dans leur débit et leur niveau de base engendrés par les changements climatiques prévus pour la période 2010–2099. Les tributaires sélectionnés (rivières Batiscan, Richelieu, Saint-Maurice, Saint-François et Yamachiche) ont été choisis en raison de leurs différences de taille, de débit et de contexte morphologique. Non seulement ces tributaires subissent-ils un régime hydrologique modifié en raison des changements climatiques, mais leur niveau de base (niveau d’eau du fleuve Saint-Laurent) sera aussi affecté. Le modèle morphodynamique en une dimension (1D) SEDROUT, à l’origine développé pour des rivières graveleuses en mode d’aggradation, a été adapté pour le contexte spécifique des tributaires des basses-terres du Saint-Laurent afin de simuler des rivières sablonneuses avec un débit quotidien variable et des fluctuations du niveau d’eau à l’aval. Un module pour simuler le partage des sédiments autour d’îles a aussi été ajouté au modèle. Le modèle ainsi amélioré (SEDROUT4-M), qui a été testé à l’aide de simulations à petite échelle et avec les conditions actuelles d’écoulement et de transport de sédiments dans quatre tributaires du fleuve Saint-Laurent, peut maintenant simuler une gamme de problèmes morphodynamiques de rivières. Les changements d’élévation du lit et d’apport en sédiments au fleuve Saint-Laurent pour la période 2010–2099 ont été simulés avec SEDROUT4-M pour les rivières Batiscan, Richelieu et Saint-François pour toutes les combinaisons de sept régimes hydrologiques (conditions actuelles et celles prédites par trois modèles de climat globaux (MCG) et deux scénarios de gaz à effet de serre) et de trois scénarios de changements du niveau de base du fleuve Saint-Laurent (aucun changement, baisse graduelle, baisse abrupte). Les impacts sur l’apport de sédiments et l’élévation du lit diffèrent entre les MCG et semblent reliés au statut des cours d’eau (selon qu’ils soient en état d’aggradation, de dégradation ou d’équilibre), ce qui illustre l’importance d’examiner plusieurs rivières avec différents modèles climatiques afin d’établir des tendances dans les effets des changements climatiques. Malgré le fait que le débit journalier moyen et le débit annuel moyen demeurent près de leur valeur actuelle dans les trois scénarios de MCG, des changements importants dans les taux de transport de sédiments simulés pour chaque tributaire sont observés. Ceci est dû à l’impact important de fortes crues plus fréquentes dans un climat futur de même qu’à l’arrivée plus hâtive de la crue printanière, ce qui résulte en une variabilité accrue dans les taux de transport en charge de fond. Certaines complications avec l’approche de modélisation en 1D pour représenter la géométrie complexe des rivières Saint-Maurice et Saint-François suggèrent qu’une approche bi-dimensionnelle (2D) devrait être sérieusement considérée afin de simuler de façon plus exacte la répartition des débits aux bifurcations autour des îles. La rivière Saint-François est utilisée comme étude de cas pour le modèle 2D H2D2, qui performe bien d’un point de vue hydraulique, mais qui requiert des ajustements pour être en mesure de pleinement simuler les ajustements morphologiques des cours d’eau.

Évaluation de l’effet des vagues de bateau sur les conditions hydrauliques près des berges en milieu fluvial

Les vagues de bateau ajoutent une pression supplémentaire sur les berges de rivières et doivent être considérées dans les modèles de prédiction des taux de recul des berges. L’objectif de cette étude est d’examiner le rôle des vagues de bateau sur l’écoulement et le transport en suspension le long des berges en milieu fluvial. Pour atteindre cet objectif, nous utilisons un transect perpendiculaire à la berge de quatre courantomètres électromagnétiques (ECMs) mesurant deux dimensions de l’écoulement et deux turbidimètres (OBSs) placés dos à dos, orientés vers la berge et le large pour mesurer les conditions moyennes et turbulentes de l’écoulement longitudinal et vertical ainsi que les flux de sédiments en suspension provoqués par les vagues. Une chaloupe à moteur de 16 pieds, équipée d’un moteur 40 hp, a été utilisée afin de générer des vagues. Nous avons mesuré l’effet de trois distances à partir de la berge (5, 10, 15 m) et trois vitesses de bateau (5, 15 et 25 km/h) et cinq répliques de chaque combinaison de distance et de vitesse ont été réalisées, totalisant 45 passages. Nous avons caractérisé la variabilité des conditions d’écoulement, de vagues et de transport de sédiments et nous avons réalisé des analyses spectrales afin de séparer les portions oscillatoire et turbulente de l’écoulement généré par les vagues de bateau. L’effet de la distance et de la vitesse du bateau sur le transport de sédiments est non-linéaire et la réponse sédimentaire induite par les passages de bateau montre une variabilité importante entre les répliques et les deux sondes OBS, ce qui suggère un changement morphologique induit par les vagues de bateau. Les corrélations entre les variables d’écoulement et de transport montrent l’importance des relations entre le cisaillement et la puissance de la portion turbulente de l’écoulement avec le transport de sédiments. Cette étude a permis de quantifier les relations entre la dynamique des vagues et les flux de concentrations de sédiments en suspension, ce qui représente une contribution importante au développement de mesures de mitigation dans les environnements fluviaux où les berges sont fragilisées par le trafic plaisancier.

Application of ocean colour data to study the oceanographic and atmospheric features off the southwest coast of India with special reference to upwelling

Satellite remote sensing is being effectively used in monitoring the ocean surface and its overlying atmosphere. Technical growth in the field of satellite sensors has made satellite measurement an inevitable part of oceanographic and atmospheric research. Among the ocean observing sensors, ocean colour sensors make use of visible band of electromagnetic spectrum (shorter wavelength). The use of shorter wavelength ensures fine spatial resolution of these parameters to depict oceanographic and atmospheric characteristics of any region having significant spaio-temporal variability. Off the southwest coast of India is such an area showing very significant spatio-temporal oceanographic and atmospheric variability due to the seasonally reversing surface winds and currents. Consequently, the region is enriched with features like upwelling, sinking, eddies, fronts, etc. Among them, upwelling brings nutrient-rich waters from subsurface layers to surface layers. During this process primary production enhances, which is measured in ocean colour sensors as high values of Chl a. Vertical attenuation depth of incident solar radiation (Kd) and Aerosol Optical Depth (AOD) are another two parameters provided by ocean colour sensors. Kd is also susceptible to undergo significant seasonal variability due to the changes in the content of Chl a in the water column. Moreover, Kd is affected by sediment transport in the upper layers as the region experiences land drainage resulting from copious rainfall. The wide range of variability of wind speed and direction may also influence the aerosol source / transport and consequently AOD. The present doctoral thesis concentrates on the utility of Chl a, Kd and AODprovided by satellite ocean colour sensors to understand oceanographic and atmospheric variability off the southwest coast of India. The thesis is divided into six Chapters with further subdivisions

Dynamics of hydrothermal plumes in Lake Banyoles

Aquesta tesi permet adquirir una millor comprensió de la dinàmica dels plomalls hidrotèrmics en l'estany de Banyoles i correlacionar la meteorologia amb els processos de fluïdització que s'hi produeixen. S'han trobat els patrons atmosfèrics que generen els esdeveniments de les fluïditzacions així com la seva freqüència. Això és crucial per a determinar la qualitat de l'aigua de l'estany. S'han trobat inhomogeneïtats espacials tant en el plomall crònic que es desenvolupa en la cubeta B1, com en el plomall episòdic de la cubeta B2. S'ha caracteritzat l'estructura del plomall generat a la cubeta B2 i s'ha comparat amb el desenvolupat a la cubeta B1. Finalment, s'han realitzat simulacions numèriques mitjançant un model numèric: MIT General Circulation Model, en el qual s'ha hagut de fer canvis en les condicions de contorn per a situar la font convectiva al fons de la columna d'aigua. Les simulacions s'han comparat amb els resultats experimentals trobats mitjançant campanyes.

Abrasion control on dune colour: Muleshoe Dunes, SW USA

The Muleshoe Dunes, an east-west trending dunefield on the border separating Texas and New Mexico, consist of two distinct components: a white (carbonate rich) component and an overlying pink (quartz rich) component. The pink component exhibits significant spatial variation in redness. The reddest sands, in the western part of the dunefield, decrease in redness towards the east. This gradient is thought to result from abrasion of all iron-rich, red clay coating as the sediments were transported eastward by Late Quaternary aeolian processes. The effects of aeolian abrasion on the spectral signature and surface texture of the sediments were examined using laboratory abrasion experiments. Changes in spectral reflectance of abrasion samples from the laboratory were compared to field samples that were abraded naturally because of sediment transport. The changes resulting from increased time of abrasion are similar to those observed with increased distance downwind in the dunefield. These results suggest that downwind abrasion can explain the pattern of dune colour in the Muleshoe Dunes, although this does not preclude other possible causes. (C) 2008 Elsevier B.V. All rights reserved.

Abrasion control on dune colour: Muleshoe Dunes, SW USA

The Muleshoe Dunes, an east-west trending dunefield on the border separating Texas and New Mexico, consist of two distinct components: a white (carbonate rich) component and an overlying pink (quartz rich) component. The pink component exhibits significant spatial variation in redness. The reddest sands, in the western part of the dunefield, decrease in redness towards the east. This gradient is thought to result from abrasion of all iron-rich, red clay coating as the sediments were transported eastward by Late Quaternary aeolian processes. The effects of aeolian abrasion on the spectral signature and surface texture of the sediments were examined using laboratory abrasion experiments. Changes in spectral reflectance of abrasion samples from the laboratory were compared to field samples that were abraded naturally because of sediment transport. The changes resulting from increased time of abrasion are similar to those observed with increased distance downwind in the dunefield. These results suggest that downwind abrasion can explain the pattern of dune colour in the Muleshoe Dunes, although this does not preclude other possible causes. (C) 2008 Elsevier B.V. All rights reserved.

Remote sensing of intertidal morphological change in Morecambe Bay, U.K., between 1991 and 2007

Tidal Flats are important examples of extensive areas of natural environment that remain relatively unaffected by man. Monitoring of tidal flats is required for a variety of purposes. Remote sensing has become an established technique for the measurement of topography over tidal flats. A further requirement is to measure topographic changes in order to measure sediment budgets. To date there have been few attempts to make quantitative estimates of morphological change over tidal flat areas. This paper illustrates the use of remote sensing to measure quantitative and qualitative changes in the tidal flats of Morecambe Bay during the relatively long period 1991–2007. An understanding of the patterns of sediment transport within the Bay is of considerable interest for coastal management and defence purposes. Tidal asymmetry is considered to be the dominant cause of morphological change in the Bay, with the higher currents associated with the flood tide being the main agency moulding the channel system. Quantitative changes were measured by comparing a Digital Elevation Model (DEM) of the intertidal zone formed using the waterline technique applied to satellite Synthetic Aperture Radar (SAR) images from 1991–1994, to a second DEM constructed from airborne laser altimetry data acquired in 2005. Qualitative changes were studied using additional SAR images acquired since 2003. A significant movement of sediment from below Mean Sea Level (MSL) to above MSL was detected by comparing the two Digital Elevation Models, though the proportion of this change that could be ascribed to seasonal effects was not clear. Between 1991 and 2004 there was a migration of the Ulverston channel of the river Leven north-east by about 5 km, followed by the development of a straighter channel to the west, leaving the previous channel decoupled from the river. This is thought to be due to independent tidal and fluvial forcing mechanisms acting on the channel. The results demonstrate the effectiveness of remote sensing for measurement of long-term morphological change in tidal flat areas. An alternative use of waterlines as partial bathymetry for assimilation into a morphodynamic model of the coastal zone is also discussed.

Extraction of tidal channel networks from aerial photographs alone and combined with laser altimetry

Tidal channel networks play an important role in the intertidal zone, exerting substantial control over the hydrodynamics and sediment transport of the region and hence over the evolution of the salt marshes and tidal flats. The study of the morphodynamics of tidal channels is currently an active area of research, and a number of theories have been proposed which require for their validation measurement of channels over extensive areas. Remotely sensed data provide a suitable means for such channel mapping. The paper describes a technique that may be adapted to extract tidal channels from either aerial photographs or LiDAR data separately, or from both types of data used together in a fusion approach. Application of the technique to channel extraction from LiDAR data has been described previously. However, aerial photographs of intertidal zones are much more commonly available than LiDAR data, and most LiDAR flights now involve acquisition of multispectral images to complement the LiDAR data. In view of this, the paper investigates the use of multispectral data for semiautomatic identification of tidal channels, firstly from only aerial photographs or linescanner data, and secondly from fused linescanner and LiDAR data sets. A multi-level, knowledge-based approach is employed. The algorithm based on aerial photography can achieve a useful channel extraction, though may fail to detect some of the smaller channels, partly because the spectral response of parts of the non-channel areas may be similar to that of the channels. The algorithm for channel extraction from fused LiDAR and spectral data gives an increased accuracy, though only slightly higher than that obtained using LiDAR data alone. The results illustrate the difficulty of developing a fully automated method, and justify the semi-automatic approach adopted.

A hybrid data assimilation scheme for model parameter estimation: application to morphodynamic modelling

We present a novel algorithm for joint state-parameter estimation using sequential three dimensional variational data assimilation (3D Var) and demonstrate its application in the context of morphodynamic modelling using an idealised two parameter 1D sediment transport model. The new scheme combines a static representation of the state background error covariances with a flow dependent approximation of the state-parameter cross-covariances. For the case presented here, this involves calculating a local finite difference approximation of the gradient of the model with respect to the parameters. The new method is easy to implement and computationally inexpensive to run. Experimental results are positive with the scheme able to recover the model parameters to a high level of accuracy. We expect that there is potential for successful application of this new methodology to larger, more realistic models with more complex parameterisations.

Data assimilation for state and parameter estimation: application to morphodynamic modelling

Data assimilation is predominantly used for state estimation; combining observational data with model predictions to produce an updated model state that most accurately approximates the true system state whilst keeping the model parameters fixed. This updated model state is then used to initiate the next model forecast. Even with perfect initial data, inaccurate model parameters will lead to the growth of prediction errors. To generate reliable forecasts we need good estimates of both the current system state and the model parameters. This paper presents research into data assimilation methods for morphodynamic model state and parameter estimation. First, we focus on state estimation and describe implementation of a three dimensional variational(3D-Var) data assimilation scheme in a simple 2D morphodynamic model of Morecambe Bay, UK. The assimilation of observations of bathymetry derived from SAR satellite imagery and a ship-borne survey is shown to significantly improve the predictive capability of the model over a 2 year run. Here, the model parameters are set by manual calibration; this is laborious and is found to produce different parameter values depending on the type and coverage of the validation dataset. The second part of this paper considers the problem of model parameter estimation in more detail. We explain how, by employing the technique of state augmentation, it is possible to use data assimilation to estimate uncertain model parameters concurrently with the model state. This approach removes inefficiencies associated with manual calibration and enables more effective use of observational data. We outline the development of a novel hybrid sequential 3D-Var data assimilation algorithm for joint state-parameter estimation and demonstrate its efficacy using an idealised 1D sediment transport model. The results of this study are extremely positive and suggest that there is great potential for the use of data assimilation-based state-parameter estimation in coastal morphodynamic modelling.

Integration of a 3D Variational data assimilation scheme with a coastal area morphodynamic model of Morecambe Bay

This paper describes the implementation of a 3D variational (3D-Var) data assimilation scheme for a morphodynamic model applied to Morecambe Bay, UK. A simple decoupled hydrodynamic and sediment transport model is combined with a data assimilation scheme to investigate the ability of such methods to improve the accuracy of the predicted bathymetry. The inverse forecast error covariance matrix is modelled using a Laplacian approximation which is calibrated for the length scale parameter required. Calibration is also performed for the Soulsby-van Rijn sediment transport equations. The data used for assimilation purposes comprises waterlines derived from SAR imagery covering the entire period of the model run, and swath bathymetry data collected by a ship-borne survey for one date towards the end of the model run. A LiDAR survey of the entire bay carried out in November 2005 is used for validation purposes. The comparison of the predictive ability of the model alone with the model-forecast-assimilation system demonstrates that using data assimilation significantly improves the forecast skill. An investigation of the assimilation of the swath bathymetry as well as the waterlines demonstrates that the overall improvement is initially large, but decreases over time as the bathymetry evolves away from that observed by the survey. The result of combining the calibration runs into a pseudo-ensemble provides a higher skill score than for a single optimized model run. A brief comparison of the Optimal Interpolation assimilation method with the 3D-Var method shows that the two schemes give similar results.

Late-Holocene changes in character and behaviour of land-terminating glaciers on James Ross Island, Antarctica

Virtually no information is available on the response of land-terminating Antarctic Peninsula glaciers to climate change on a centennial timescale. This paper analyses the topography, geomorphology and sedimentology of prominent moraines on James Ross Island, Antarctica, to determine geometric changes and to interpret glacier behaviour. The moraines are very likely due to a late-Holocene phase of advance and featured (1) shearing and thrusting within the snout, (2) shearing and deformation of basal sediment, (3) more supraglacial debris than at present and (4) short distances of sediment transport. Retreat of ∼100 m and thinning of 15–20 m has produced a loss of 0.1 km3 of ice. The pattern of surface lowering is asymmetric. These geometrical changes are suggested most simply to be due to a net negative mass balance caused by a drier climate. Comparisons of the moraines with the current glaciological surface structure of the glaciers permits speculation of a transition from a polythermal to a cold-based thermal regime. Small land-terminating glaciers in the northern Antarctic Peninsula region could be cooling despite a warming climate.

Controls of heavy minerals and grain size in a holocene regressive barrier (Ilha Comprida, southeastern Brazil)

Ilha Comprida is a regressive barrier island located in southeastern Brazil that was formed essentially by Quaternary sandy sediments. Ilha Comprida sediments were analyzed to assess heavy mineral indices and grain size variables. The spatial variation of heavy minerals and grain size was interpreted in terms of the present barrier dynamics and the barrier`s evolution since the Middle Holocene. These analyses allowed for the identification of the main factors and processes that control the variation of heavy minerals and grain size on the barrier. Rutile and zircon (RZi) and tourmaline and hornblende (THi) are significantly sensitive to provenance and exhibit the contributions of the Ribeira de Iguape River sediments, which reach the coast next to the northeastern end of Ilha Comprida. In addition to the influence of provenance, TZi responds mainly to hydraulic sorting processes. This agrees with a sediment transport pattern characterized by a divergence of two resultant net alongshore drifts southwest of the barrier. The sediments from the Ribeira de Iguape River reach the barrier directly through the river mouth and indirectly after temporary storage in the inner shelf. The combination of grain size and heavy mineral analyses is a reliable method for determining sediment transport patterns and provenance. (C) 2010 Elsevier Ltd. All rights reserved.

Geração de sulcos em meio coesivo e caracterização hidrossedimentométrica de sua evolução

Na pesquisa em erosão, nas últimas décadas, está se formando um consenso de que é importante entender os processos básicos que regem o fenômeno. Uma alternativas para tentar compreender melhor as etapas do processo erosivo é separá-lo na fase de sulco (fluxo concentrado) e de entressulco. Dentro desse enfoque foi construído no Laboratório de Processos Erosivos e Deposicionais (LaPED) do IPH/UFRGS um canal de declividade para estudar o processo de incisão e o desenvolvimento dos sulcos de erosão. A estrutura experimental projetada e construída permite que seja controlada a vazão através de um medidor eletromagnético e que seja alterada a declividade do canal através de um sistema hidráulico associado a um nível digital. O solo colocado no canal foi um Latossolo Vermelho distrófico típico, as declividades de trabalho foram 3,0; 6,0 e 9,0% e a seqüência de vazões aplicadas foi 10,0; 18,5; 25,5; 38,5 e 51,0L.min-1. A estrutura experimental montada se mostrou de fácil operação e eficiente para permitir o avanço no entendimento dos processos de desagregação e de transporte de partículas sólidas pela ação do escoamento superficial, além de possibilitar a geração de sulco(s) de erosão na superfície do solo. O escoamento passou da condição de difuso para concentrado a partir do momento em que a velocidade superficial do fluxo alcançou 0,26m.s-1, a altura de lâmina atingiu 0,0102m, a velocidade de cisalhamento superou os 0,059m.s-1, a tensão de cisalhamento chegou a 3,50Pa e que a potência do escoamento atingiu pelo menos 0,22N.s-1. O processo de incisão iniciou-se com o canal experimental colocado em baixa declividade e em regime de escoamento sub-crítico e de transição. A velocidade de cisalhamento, no momento da incisão, foi, praticamente, o dobro daquela encontrada na literatura para solos siltosos e arenosos. Entretanto, para as três declividades a fase de sulco definido ocorreu somente em regime de escoamento turbulento. A tensão de cisalhamento foi o parâmetro que melhor descreveu a evolução da perda de solo. A potência do escoamento foi o parâmetro hidráulico que mostrou maior eficiência para separar as fases evolutivas dos sulcos. O desenvolvimento do(s) sulco(s) teve o seu início em uma condição de escoamento difuso (ausência de sulcos) e com a potência do escoamento oscilando entre 0,057 e 0,198N.s-1. O avanço do(s) sulco(s) começou com uma zona de transição (fase de incisão e de aprofundamento) onde a potência do escoamento varia entre 0,220 e 0,278N.s-1 e, logo em seguida, teve início a fase de sulco definido, com a potência do escoamento entre 0,314 e 0,544N.s-1. Na fase de escoamento concentrado foi preponderante o papel do processo de erosão regressiva para aumentar tanto o tamanho como o peso das partículas sólidas em transporte pelo escoamento superficial e assim fazer com que predominasse o transporte via fundo sobre o transporte via suspensão. As cargas de sedimento geradas nos solos de diferentes classes texturais foram separadas em grupos distintos em função da potência unitária do escoamento.