The influence of water column variability and seafloor topography on multibeam bathymetric mapping

Senior thesis written for Oceanography 445

Batimetria fluvial estimada com dados orbitais: Estudo de caso no alto curso do rio Paraguai com o sensor ASTER

The Paraguay River is the main tributary of the Paraná River and has an extension of 1.693 km in Brazilian territory. The navigability conditions are very important for the regional economy because most of the central-west Brazilian agricultural and mineral production is transported by the Paraguay waterway. Increased sedimentation along the channel requires continuous dredging to waterway maintenance. Systematic bathymetric surveys are periodically carried out in order to check depth condition along the channel using echo-sounding devices. In this paper, digital image processing and geostatistical analysis methods were used to analyze the applicability of the ASTER sensor to estimate channel depths in a segment of the upper Paraguay River. The results were compared with field data in order to choose the band with better adjustment and to evaluate the standard deviation. Comparing the VNIR bands, the best fit was presented by the red wavelength (band 2; 0,63 - 0,69 μm), showing a good representation of the channel depths shallow than 1,7 m. Applying geostatistical methods, the model accuracy was enhanced from 43 cm to 36 cm and undesired components were slacked. It was concluded that the digital number of band 2, converted to bathymetry information allows a good estimation of river depths and channel morphology.

A snapshot of suspended sediment and fluid mud occurrence in a mixed-energy embayment, Tijucas Bay, Brazil

Along the southern Brazilian coast, Tijucas Bay is known for its unique muddy tidal flats associated with chenier plains. Previous field observations pointed to very high suspended sediment concentrations (SSCs) in the inner parts of the bay, and in the estuary of the Tijucas River, suggesting the presence of fluid mud. In this study, the occurrences of suspended sediments and fluid mud were examined during a larger-scale, high-resolution 2-day field campaign on 1-2 May 2007, encompassing survey lines spanning nearly 80 km, 75 water sampling stations for near-bottom density estimates, and ten sediment sampling stations. Wave refraction modeling provided qualitative wave energy estimates as a function of different incidence directions. The results show that SSC increases toward the inner bay near the water surface, but seaward near the bottom. This suggests that suspended sediment is supplied by the local rivers, in particular the Tijucas. Near-surface SSCs were of the order of 50 mg l(-1) close to the shore, but exceeded 100 mg l(-1) near the bottom in the deeper parts of the bay. Fluid mud thickness and location given by densimetry and echo-sounding agreed in some places, although being mostly discordant. The best agreement was observed where wave energy was high during the campaign. The discrepancy between the two methods may be an indication for the existence of fluid mud, which is recorded by one method but not the other. Agreement is considered to be an indication of fluidization, whereas disagreement indicates more consolidation. Wave modeling suggests that waves from the ENE and SE are the most effective in supplying energy to the inner bay, which may induce the liquefaction of mud deposits to form fluid mud. Nearshore mud resuspension and weak horizontal currents result in sediment-laden offshore flow, which explains the higher SSCs measured in the deeper parts of the bay, besides providing a mechanism for fine-sediment export to the inner shelf.

Evidence for a Frozen Bed, Byrd Glacier, Antarctica

Ice thickness, computed within the fjord region of Byrd Glacier on the assumptions that Byrd Glacier is in mass-balance equilibrium and that ice velocity is entirely due to basal sliding, are on average 400 m less than measured ice thicknesses along a radio-echo profile. We consider four explanations for these differences: (1) active glacier ice is separated from a zone of stagnant ice near the base of the glacier by a shear zone at depth; (2) basal melting rates are some 8 m/yr; (3) internal shear occurs with no basal sliding in much of the region above the grounding zone; or (4) internal creep and basal sliding contribute to the flow velocity in varying proportions above the grounding zone. Large gradients of surface strain rate seem to invalidate the first explanation. Computed values of basal shear stress (140 to 200 kPa) provide insufficient frictional heat to melt the ice demanded by the second explanation. Both the third and fourth explanations were examined by making simplifying assumptions that prevented a truly quantitative evaluation of their merit. Nevertheless, there is no escaping the qualitative conclusion that internal shear contributes strongly to surface velocities measured on Byrd Glacier, as is postulated in both these explanations.

Deep drilling of glaciers: Russian projects in the Arctic (1975-1995)

The data collection "Deep Drilling of Glaciers: Soviet-Russian projects in Arctic, 1975-1995" was collected by the following basic considerations: - compilation of deep (>100 m) drilling projects on Arctic glaciers, using data of (a) publications; (b) archives of IGRAN; (c) personal communication of project participants; - documentation of parameters, references. Accuracy of data and techniques applied to determine different parameters are not evaluated. The accuracy of some geochemical parameters (up to 1984 and heavy metalls) is uncertain. Most reconstructions of ice core age and of annual layer thickness are discussed; - digitizing of published diagrams (in case, when original numerical data were lost) and subsequent data conversion to equal range series and adjustment to the common units. Therefore, the equal-range series were calculated from original data or converted from digitized chart values as indicated in the metadata. For the methodological purpose, the equal-range series obtained from original and reconstructed data were compared repeatedly; the systematic difference was less then 5-7%. Special attention should be given to the fact, that the data for individual ice core parameters varies, because some parameters were originally measured or registered. Parameters were converted in equal-range series using 2 m steps; - two or more parameter values were determined, then the mean-weighted (i.e. accounting the sample length) value is assigned to the entire interval; - one parameter value was determined, measured or registered independently from the parameter values in depth intervals which over- and underlie it, then the value is assigned to the entire interval; - one parameter value was determined, measured or registered for two adjoining depth intervals, then the specific value is assigned to the depth interval, which represents >75% of sample length ; if each of adjoining depth intervals represents <75% of sample length, then the correspondent parameter value is assigned to both intervals of depth. This collection of ice core data (version 2000) was made available through the EU funded QUEEN project by S.M. Arkhipov, Moscow.

Penetration velocity and deceleration of penetrometer stations N1 to N97 during cruise Senckenberg_11_2008

In-situ geotechnical measurements of surface sediments were carried out along large subaqueous dunes in the Knudedyb tidal inlet channel in the Danish Wadden Sea using a small free-falling penetrometer. Vertical profiles showed a typical stratification pattern with a resolution of ~1 cm depicting a thin surface layer of low sediment strength and a stiffer substratum below (quasi-static bearing capacity equivalent: 1-3 kPa in the top layer, 20-140 kPa in the underlying sediment; thickness of the top layer ca. 5-8 cm). Observed variations in the thickness and strength of the surface layer during a tidal cycle were compared to mean current velocities (measured using an acoustic Doppler current profiler, ADCP), high-resolution bathymetry (based on multibeam echo sounding, MBES) and qualitative estimates of suspended sediment distributions in the water column (estimated from ADCP backscatter intensity). The results revealed an ebb dominance in sediment remobilization, and a general accretion of the bed towards low water. A loose top layer occurred throughout the tidal cycle, likely influenced by bedload transport and small events of suspended sediment resettlement (thickness: 6 +-2 cm). Furthermore, this layer showed a significant increase in thickness (e.g. from 8 cm to 16 cm) related to periods of overall deposition. These findings imply that dynamic penetrometers can conveniently serve to (1) quantify potentially mobile sediments by determining the thickness of a loose sediment surface layer, (2) unravel sediment strength development in potentially mobile sediments and (3) identify sediment accumulation. Such data are an important complement and add a new geotechnical perspective during investigations of sediment remobilization processes in highly dynamic coastal environments.

Active methane seepage in the Southern Ocean, off the sub-Antarctic island of South Georgia

An extensive submarine cold-seep area was discovered on the northern shelf of South Georgia during R/V Polarstern cruise ANT-XXIX/4 in spring 2013. Hydroacoustic surveys documented the presence of 133 gas bubble emissions, which were restricted to glacially-formed fjords and troughs. Video-based sea floor observations confirmed the sea floor origin of the gas emissions and spatially related microbial mats. Effective methane transport from these emissions into the hydrosphere was proven by relative enrichments of dissolved methane in near-bottom waters. Stable carbon isotopic signatures pointed to a predominant microbial methane formation, presumably based on high organic matter sedimentation in this region. Although known from many continental margins in the world's oceans, this is the first report of an active area of methane seepage in the Southern Ocean. Our finding of substantial methane emission related to a trough and fjord system, a topographical setting that exists commonly in glacially-affected areas, opens up the possibility that methane seepage is a more widespread phenomenon in polar and sub-polar regions than previously thought.

Petrography and geochemical composition of the Atlantis II Fracture Zone

SeaBeam echo sounding, seismic reflection, magnetics, and gravity profiles were run along closely spaced tracks (5 km) parallel to the Atlantis II Fracture Zone on the Southwest Indian Ridge, giving 80% bathymetric coverage of a 30- * 170-nmi strip centered over the fracture zone. The southern and northern rift valleys of the ridge were clearly defined and offset north-south by 199 km. The rift valleys are typical of those found elsewhere on the Southwest Indian Ridge, with relief of more than 2200 m and widths from 22 to 38 km. The ridge-transform intersections are marked by deep nodal basins lying on the transform side of the neovolcanic zone that defines the present-day spreading axis. The walls of the transform generally are steep (25°-40°), although locally, they can be more subdued. The deepest point in the transform is 6480 m in the southern nodal basin, and the shallowest is an uplifted wave-cut terrace that exposes plutonic rocks from the deepest layer of the ocean crust at 700 m. The transform valley is bisected by a 1.5-km-high median tectonic ridge that extends from the northern ridge-transform intersection to the midpoint of the active transform. The seismic survey showed that the floor of the transform contains up to 0.5 km of sediment. Piston-coring at two locations on the transform floor recovered more than 1 m of sand and gravel, which appears to be turbidites shed from the walls of the fracture zone. Extensive dredging showed that more than two-thirds of the crust exposed in the transform valley and its walls were plutonic rocks, principally gabbros and residual mantle peridotites. In contrast, based on dredging and seafloor morphology, only relatively undisrupted pillow basalt flows have been exposed on crust of the same age spreading away from the transform. Magnetic anomalies are well defined out to 11 m.y. over the flanking transverse ridges and transform valley, even where layer 2 appears to be absent. The total opening rate is 1.6 cm/yr, but the arrangement of the anomalies indicates that the spreading for each ridge is asymmetric, with the ridge flanks facing the transform spreading at a rate of 1.0 cm/yr. Such an asymmetric spreading pattern requires that both the northern and southern ridges migrate away from each other at 0.2 cm/yr, thus lengthening the transform at 0.4 cm/yr for the last 11 m.y. To the north, the fracture zone valley is oriented differently from the present-day transform, indicating a paleospreading direction change at 17 m.y. from N10°E to due north-south. This change placed the transform into extension for the 11-m.y. period required for simple orthogonal ridge-transform geometry to be reestablished and produced a large transtensional basin within the transform valley. This basin was split by continued transform slip after 11 m.y., with the larger half moving to the north with the African Plate.

Gridded bathymetry from multibeam echosounder EM122 data of the cruise M84/2, off Eregli (2011)

We investigated gas bubble emissions at the Don-Kuban paleo-fan in the northeastern Black Sea regarding their geological setting, quantities as well as spatial and temporal variabilities during three ship expeditions between 2007 and 2011. About 600 bubble-induced hydroacoustic anomalies in the water column (flares) originating from the seafloor above the gas hydrate stability zone (GHSZ) at ~700 m water depth were found. At about 890 m water depth a hydrocarbon seep area named "Kerch seep area" was newly discovered within the GHSZ. We propose locally domed sediments ('mounds') discovered during ultra-high resolution bathymetric mapping with an autonomous underwater vehicle (AUV) to result from gas hydrate accumulation at shallow depths. In situ measurements indicated spatially limited temperature elevations in the shallow sediment likely induced by upward fluid flow which may confine the local GHSZ to a few meters below the seafloor. As a result, gas bubbles are suspected to migrate into near-surface sediments and to escape the seafloor through small-scale faults. Hydroacoustic surveys revealed that several flares originated from a seafloor area of about 1 km**2 in size. The highest flare disappeared in about 350 m water depth, suggesting that the released methane remains in the water column. A methane flux estimate, combining data from visual quantifications during dives with a remotely operated vehicle (ROV) with results from ship-based hydroacoustic surveys and gas analysis revealed that between 2 and 87 x 10**6 mol CH4 yr-1 escaped into the water column above the Kerch seep area. Our results show that the finding of the Kerch seep area represents a so far underestimated type of hydrocarbon seep, which has to be considered in methane budget calculations.

Estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea-level rise = Estimación del volumen de hielo de Svalbard y su contribución potencial a la subida del nivel del mar

El objetivo final de las investigaciones recogidas en esta tesis doctoral es la estimación del volumen de hielo total de los ms de 1600 glaciares de Svalbard, en el Ártico, y, con ello, su contribución potencial a la subida del nivel medio del mar en un escenario de calentamiento global. Los cálculos más exactos del volumen de un glaciar se efectúan a partir de medidas del espesor de hielo obtenidas con georradar. Sin embargo, estas medidas no son viables para conjuntos grandes de glaciares, debido al coste, dificultades logísticas y tiempo requerido por ellas, especialmente en las regiones polares o de montaña. Frente a ello, la determinación de áreas de glaciares a partir de imágenes de satélite sí es viable a escalas global y regional, por lo que las relaciones de escala volumen-área constituyen el mecanismo más adecuado para las estimaciones de volúmenes globales y regionales, como las realizadas para Svalbard en esta tesis. Como parte del trabajo de tesis, hemos elaborado un inventario de los glaciares de Svalbard en los que se han efectuado radioecosondeos, y hemos realizado los cálculos del volumen de hielo de más de 80 cuencas glaciares de Svalbard a partir de datos de georradar. Estos volúmenes han sido utilizados para calibrar las relaciones volumen-área desarrolladas en la tesis. Los datos de georradar han sido obtenidos en diversas campañas llevadas a cabo por grupos de investigación internacionales, gran parte de ellas lideradas por el Grupo de Simulación Numérica en Ciencias e Ingeniería de la Universidad Politécnica de Madrid, del que forman parte la doctoranda y los directores de tesis. Además, se ha desarrollado una metodología para la estimación del error en el cálculo de volumen, que aporta una novedosa técnica de cálculo del error de interpolación para conjuntos de datos del tipo de los obtenidos con perfiles de georradar, que presentan distribuciones espaciales con unos patrones muy característicos pero con una densidad de datos muy irregular. Hemos obtenido en este trabajo de tesis relaciones de escala específicas para los glaciares de Svalbard, explorando la sensibilidad de los parámetros a diferentes morfologías glaciares, e incorporando nuevas variables. En particular, hemos efectuado experimentos orientados a verificar si las relaciones de escala obtenidas caracterizando los glaciares individuales por su tamaño, pendiente o forma implican diferencias significativas en el volumen total estimado para los glaciares de Svalbard, y si esta partición implica algún patrón significativo en los parámetros de las relaciones de escala. Nuestros resultados indican que, para un valor constante del factor multiplicativo de la relacin de escala, el exponente que afecta al área en la relación volumen-área decrece según aumentan la pendiente y el factor de forma, mientras que las clasificaciones basadas en tamaño no muestran un patrón significativo. Esto significa que los glaciares con mayores pendientes y de tipo circo son menos sensibles a los cambios de área. Además, los volúmenes de la población total de los glaciares de Svalbard calculados con fraccionamiento en grupos por tamaño y pendiente son un 1-4% menores que los obtenidas usando la totalidad de glaciares sin fraccionamiento en grupos, mientras que los volúmenes calculados fraccionando por forma son un 3-5% mayores. También realizamos experimentos multivariable para obtener estimaciones óptimas del volumen total mediante una combinación de distintos predictores. Nuestros resultados muestran que un modelo potencial simple volumen-área explica el 98.6% de la varianza. Sólo el predictor longitud del glaciar proporciona significación estadística cuando se usa además del área del glaciar, aunque el coeficiente de determinación disminuye en comparación con el modelo más simple V-A. El predictor intervalo de altitud no proporciona información adicional cuando se usa además del área del glaciar. Nuestras estimaciones del volumen de la totalidad de glaciares de Svalbard usando las diferentes relaciones de escala obtenidas en esta tesis oscilan entre 6890 y 8106 km3, con errores relativos del orden de 6.6-8.1%. El valor medio de nuestras estimaciones, que puede ser considerado como nuestra mejor estimación del volumen, es de 7.504 km3. En términos de equivalente en nivel del mar (SLE), nuestras estimaciones corresponden a una subida potencial del nivel del mar de 17-20 mm SLE, promediando 19_2 mm SLE, donde el error corresponde al error en volumen antes indicado. En comparación, las estimaciones usando las relaciones V-A de otros autores son de 13-26 mm SLE, promediando 20 _ 2 mm SLE, donde el error representa la desviación estándar de las distintas estimaciones. ABSTRACT The final aim of the research involved in this doctoral thesis is the estimation of the total ice volume of the more than 1600 glaciers of Svalbard, in the Arctic region, and thus their potential contribution to sea-level rise under a global warming scenario. The most accurate calculations of glacier volumes are those based on ice-thicknesses measured by groundpenetrating radar (GPR). However, such measurements are not viable for very large sets of glaciers, due to their cost, logistic difficulties and time requirements, especially in polar or mountain regions. On the contrary, the calculation of glacier areas from satellite images is perfectly viable at global and regional scales, so the volume-area scaling relationships are the most useful tool to determine glacier volumes at global and regional scales, as done for Svalbard in this PhD thesis. As part of the PhD work, we have compiled an inventory of the radio-echo sounded glaciers in Svalbard, and we have performed the volume calculations for more than 80 glacier basins in Svalbard from GPR data. These volumes have been used to calibrate the volume-area relationships derived in this dissertation. Such GPR data have been obtained during fieldwork campaigns carried out by international teams, often lead by the Group of Numerical Simulation in Science and Engineering of the Technical University of Madrid, to which the PhD candidate and her supervisors belong. Furthermore, we have developed a methodology to estimate the error in the volume calculation, which includes a novel technique to calculate the interpolation error for data sets of the type produced by GPR profiling, which show very characteristic data distribution patterns but with very irregular data density. We have derived in this dissertation scaling relationships specific for Svalbard glaciers, exploring the sensitivity of the scaling parameters to different glacier morphologies and adding new variables. In particular, we did experiments aimed to verify whether scaling relationships obtained through characterization of individual glacier shape, slope and size imply significant differences in the estimated volume of the total population of Svalbard glaciers, and whether this partitioning implies any noticeable pattern in the scaling relationship parameters. Our results indicate that, for a fixed value of the factor in the scaling relationship, the exponent of the area in the volume-area relationship decreases as slope and shape increase, whereas size-based classifications do not reveal any clear trend. This means that steep slopes and cirque-type glaciers are less sensitive to changes in glacier area. Moreover, the volumes of the total population of Svalbard glaciers calculated according to partitioning in subgroups by size and slope are smaller (by 1-4%) than that obtained considering all glaciers without partitioning into subgroups, whereas the volumes calculated according to partitioning in subgroups by shape are 3-5% larger. We also did multivariate experiments attempting to optimally predict the volume of Svalbard glaciers from a combination of different predictors. Our results show that a simple power-type V-A model explains 98.6% of the variance. Only the predictor glacier length provides statistical significance when used in addition to the predictor glacier area, though the coefficient of determination decreases as compared with the simpler V-A model. The predictor elevation range did not provide any additional information when used in addition to glacier area. Our estimates of the volume of the entire population of Svalbard glaciers using the different scaling relationships that we have derived along this thesis range within 6890-8106 km3, with estimated relative errors in total volume of the order of 6.6-8.1% The average value of all of our estimates, which could be used as a best estimate for the volume, is 7,504 km3. In terms of sea-level equivalent (SLE), our volume estimates correspond to a potential contribution to sea-level rise within 17-20 mm SLE, averaging 19 _ 2 mm SLE, where the quoted error corresponds to our estimated relative error in volume. For comparison, the estimates using the V-A scaling relations found in the literature range within 13-26 mm SLE, averaging 20 _ 2 mm SLE, where the quoted error represents the standard deviation of the different estimates.