21 ground-penetrating radar lines of Area 1 - Beach ridge plain (low wave energy) from Potter Cove, King George Island, Antarctic Peninsula

During two field campaigns (Austral springs 2011 and 2012) the sedimentary architecture of a polar gravel-beach system at the western coast of Potter Peninsula (Area 1) was revealed using ground-penetrating radar (GPR, Geophysical Survey Systems, Inc. SIR-3000). 21 profiles were collected using a mono-static 200 MHz antenna operated in common offset mode. Trace increment was set to 0.05 m. A differential global-positioning system (dGPS, Leica GS09) was used to obtain topographical information along the GPR lines. GPR data are provided in RADAN-Format, dGPS coordinates are provided in ascii format; projection is UTM (WGS 84, zone 21S).

49 ground-penetrating radar lines of Area 2 - Gravel cuspate spit with gravel beach ridges from Potter Cove, King George Island, Antarctic Peninsula

During two field campaigns (Austral springs 2011 and 2012) the sedimentary architecture of a polar gravel-beach system at the southwestern coast of Potter Peninsula (Area 2) was revealed using ground-penetrating radar (GPR, Geophysical Survey Systems, Inc. SIR-3000). 49 profiles were collected using a mono-static 200 MHz antenna operated in common offset mode. Trace increment was set to 0.05 m. A differential global-positioning system (dGPS, Leica GS09) was used to obtain topographical information along the GPR lines. GPR data are provided in RADAN-Format, dGPS coordinates are provided in ascii format; projection is UTM (WGS 84, zone 21S).

Sensitivity of a distributed temperature-radiation index melt model based on AWS observations and surface energy balance fluxes, Hurd Peninsula glaciers, Livingston Island, Antarctica

We use an automatic weather station and surface mass balance dataset spanning four melt seasons collected on Hurd Peninsula Glaciers, South Shetland Islands, to investigate the point surface energy balance, to determine the absolute and relative contribution of the various energy fluxes acting on the glacier surface and to estimate the sensitivity of melt to ambient temperature changes. Long-wave incoming radiation is the main energy source for melt, while short-wave radiation is the most important flux controlling the variation of both seasonal and daily mean surface energy balance. Short-wave and long-wave radiation fluxes do, in general, balance each other, resulting in a high correspondence between daily mean net radiation flux and available melt energy flux. We calibrate a distributed melt model driven by air temperature and an expression for the incoming short-wave radiation. The model is calibrated with the data from one of the melt seasons and validated with the data of the three remaining seasons. The model results deviate at most 140 mm w.e. from the corresponding observations using the glaciological method. The model is very sensitive to changes in ambient temperature: a 0.5 ◦ C increase results in 56 % higher melt rates.

Caracterización dinámica del glaciar Hurd combinando observaciones de campo y simulaciones numéricas

El objetivo fundamental de esta tesis es la caracterización de la morfología y del estado de deformaciones y tensiones del Glaciar Hurd (Isla Livingston, Archipiélago de las Shetland del Sur, Antártida), mediante una combinación de observaciones de campo, registros de georradar y simulaciones numéricas. La morfología y el estado de deformaciones y tensiones actuales son la expresión de la evolución dinámica del glaciar desde tiempos pretéritos hasta recientes, y su análisis nos dará las pautas con las cuales ser capaces de predecir, con el apoyo de las simulaciones numéricas, su evolución futura. El primer aspecto que se aborda es el estudio de las estructuras que pueden observarse en la superficie del glaciar. Describimos las distintas técnicas utilizadas (medidas de campo, fotointerpretación de ortofotografías, análisis geoquímico de cenizas volcánicas, etc.) y presentamos el análisis e interpretación de los resultados morfo-estructurales, así como la correlación, mediante análisis geoquímicos (fluorescencia de rayos X), entre las cenizas volcánicas que extruyen en la superficie del Glaciar Hurd y las del volcán Decepción, origen de las cenizas. Esto nos permite realizar una datación de las mismas como Tefra 1, correspondiente a la erupción de 1970, Tefra 2, correspondiente a las erupciones pre-1829, y el conjunto Tefra 3, asociado a las erupciones más antiguas. En segundo lugar nos ocupamos de las estructuras presentes en el interior del glaciar, cuya herramienta de detección fundamental es el georradar. Identificadas estas estructuras internas, las vinculamos con las observadas en la superficie del glaciar. También hemos estudiado la estructura hidrotérmica del glaciar, obteniendo una serie de evidencias adicionales de su carácter politérmico. Entre éstas se contaban, hasta ahora, las basadas en el valor del parámetro de rigidez de la relación constitutiva del hielo determinada por ajuste de modelos dinámicos y observaciones realizados por Otero (2008) y las basadas en las velocidades de las ondas de radar en el hielo determinadas con el método de punto medio común por Navarro y otros (2009). Las evidencias adicionales que aportamos en esta tesis son: 1) la presencia de estructuras típicas de régimen compresivo en la zona terminal del glaciar y de cizalla en los márgenes del mismo, y 2) la presencia de un estrato superficial de hielo frío (por encima de otro templado) en la zona de ablación de los tres lóbulos del Glaciar Hurd –Sally Rocks, Argentina y Las Palmas–, que alcanzan espesores de 70, 50 y 40 m, respectivamente. Este estrato de hielo frío está probablemente congelado al lecho subglaciar en la zona terminal (Molina y otros, 2007; esta tesis). Por último, nos ocupamos de la simulación numérica de la dinámica glaciar. Presentamos el modelo físico-matemático utilizado, discutimos sus condiciones de contorno y cómo éstas se miden en los trabajos de campo, y describimos el procedimiento de resolución numérica del sistema de ecuaciones parciales del modelo. Presentamos los resultados para los campos de velocidades, deformaciones y tensiones, comparando estos resultados con las estructuras observadas. También incluimos el análisis de las elipses de deformación acumulativa, que proporcionan información sobre las estructuras a las que puede dar lugar la evolución del estado de deformaciones y tensiones a las que se ve sometido el hielo según avanza, lentamente, desde la cabecera hasta la zona terminal del glaciar, con tiempos de tránsito de hasta 1.250 años, recogiendo así la historia de deformaciones en el glaciar. Concluyendo, ponemos de manifiesto en esta tesis que las medidas de campo de las estructuras y niveles de cenizas, las medidas de georradar y las simulaciones numéricas de la dinámica glaciar, realizadas de forma combinada, permiten caracterizar el régimen actual de velocidades, deformaciones y tensiones del glaciar, entender su evolución en el pasado y predecir su evolución futura. ABSTRACT The main objective of this thesis is to characterize the morphology and the state of strains and stresses of Hurd Glacier (Livingston Island, South Shetland Islands archipelago, Antarctica) through a combination of field observations, ground-penetrating radar measurements and numerical simulations. The morphology and the current state of strain and stresses are the expression of the dynamic evolution of the glacier from the past to recent times, and their analysis gives us the guidelines to be able to predict, with the support of numerical simulations, its future evolution. The first subject addressed is the study of structures that can be observed on the glacier surface. We describe the different techniques used (field measurements, photointerpretation of orthophotos, geochemical analysis of volcanic ashes, etc.) and we present the analysis and interpretation of the morpho-structural results, as well as the correlation with geochemical analysis (XRF) between the volcanic ashes extruded to the surface of Hurd Glacier and those of Deception Island volcano, from which the ashes originate. This allows us dating the ashes as Tephra 1, corresponding to the 1970 eruption, Tephra 2, corresponding to the pre-1829 eruptions, and the Tephra 3 group, associated with older eruptions. Secondly we focus on the study of the structures present within the glacier, which are detected with the help of ground-penetrating radar. Once identified, we link these internal structures with those observed on the glacier surface. We also study the hydrothermal structure of the glacier, getting a series of additional evidences of its polythermal structure. Among the evidences available so far, we can mention those based on the value of the stiffness parameter of the constitutive relation of ice, determined by fitting dynamic models to observations, as done by Otero (2008), and those based on the velocity of propagation of the radar waves through the glacier ice, measured using the common midpoint method, as done by Navarro et al. (2009). The additional evidences that we provide in this thesis are: 1) the presence of structures typical of compressive regime in the terminal zone of the glacier, together with shear at its margins, and 2) the presence of a surface layer of cold ice (overlying a layer of temperate ice) in the ablation zone of the three lobes of Hurd Glacier –Sally Rocks, Argentina and Las Palmas–, reaching thicknesses of 70, 50 and 40 m, respectively. This cold layer is probably frozen to the subglacial bed in the terminal zone (Molina and others 2007; this thesis). Finally, we deal with the numerical simulation of glacier dynamics. We present the physical-mathematical model, discuss its boundary conditions and how they are measured in the field work, and describe the method of numerical solution of the model’s partial differential equations. We present the results for the velocity, strain and stress fields, comparing these results with the observed structures. We also include an analysis of the ellipses of cumulative deformation, which provide information about the structures that can result from the evolution of the strain and stress regime of the glacier ice as it moves slowly from the head to the snout of the glacier, with transit times of up to 1,250 years, so picking the history of deformation of the glacier. Summarizing, we show in this thesis that field measurements of structures and ash layers, ground-penetrating radar measurements and numerical simulations of glacier dynamics, performed in combination, allow us to characterize the current regime of velocities, strains and stresses of the glacier, to understand its past evolution and to predict its future evolution.

Surface velocity and mass balance of Livingston Island ice cap, Antarctica

The mass budget of the ice caps surrounding the Antarctica Peninsula and, in particular, the partitioning of its main components are poorly known. Here we approximate frontal ablation (i.e. the sum of mass losses by calving and submarine melt) and surface mass balance of the ice cap of Livingston Island, the second largest island in the South Shetland Islands archipelago, and analyse variations in surface velocity for the period 2007–2011. Velocities are obtained from feature tracking using 25 PALSAR-1 images, and used in conjunction with estimates of glacier ice thicknesses inferred from principles of glacier dynamics and ground-penetrating radar observations to estimate frontal ablation rates by a flux-gate approach. Glacier-wide surface mass-balance rates are approximated from in situ observations on two glaciers of the ice cap. Within the limitations of the large uncertainties mostly due to unknown ice thicknesses at the flux gates, we find that frontal ablation (−509 ± 263 Mt yr−1, equivalent to −0.73 ± 0.38 m w.e. yr−1 over the ice cap area of 697 km2) and surface ablation (−0.73 ± 0.10 m w.e. yr−1) contribute similar shares to total ablation (−1.46 ± 0.39 m w.e. yr−1). Total mass change (δM = −0.67 ± 0.40 m w.e. yr−1) is negative despite a slightly positive surface mass balance (0.06 ± 0.14 m w.e. yr−1). We find large interannual and, for some basins, pronounced seasonal variations in surface velocities at the flux gates, with higher velocities in summer than in winter. Associated variations in frontal ablation (of ~237 Mt yr−1; −0.34 m w.e. yr−1) highlight the importance of taking into account the seasonality in ice velocities when computing frontal ablation with a flux-gate approach.

Tab. 1: Zircon isotopic data and rock ages from the Antarctic Peninsula

Morphological and U-Pb isotope studies on sedimentary zircons reflect the orogenic evolution of their former host rocks. The orogenic history of detrital zircons from the Trinity Peninsula Formation (TPF) defines the former geological surrounding of the sedimentation basin of the TPF. Same few weil rounded, polycyclic zircons of Precambrian age and Cambrian overprint give hints for an old cratonic source rock. Because of their very low frequency compared with euhedral types, the contribution of an cratonic shield area to the bulk of the sedimentary debris is neglectable low. Euhedral zircons of granitoid origin and Carboniferous age indicate a derivation from an area of widespread Carboniferous intrusions. Except for southern South America and unsurveyed regions in the Antarctic Peninsula itself, no region could deliver zircons with a Carboniferous age record. The only acceptable explanation for the origin of these zircons is a position of the Antarctic Peninsula during the sedimentation of the TPF approximately southwest of southern South America.

Chemistry of volcanic ash layers from the southern Atlantic

Petrographic and geochemical investigations were carried out on 21 ash layers from four sites of ODP Legs 113 and 114 in the southern Atlantic Ocean. With the help of geochemical data and petrographic characterization three rock series can be distinguished for stratigraphically different ash layers from Site 701 (Leg 114) located east of the South Sandwich Island Arc, whereas the Leg 113 tephras from the southern slope of the South Orkney Microcontinent belong to another magmatic series. Geochemical correlation of the Leg 113 tephras with possible source areas indicates that they were probably erupted from the Antarctic Peninsula. The Miocene ashes from Site 701 are probably derived from the now-extinct Discovery Arc, the precursor of the South Sandwich Islands. The Pliocene ashes from the site show some affinity with the South Shetland Islands, although the available data do not permit a clear correlation. The Quaternary ashes from Site 701 display a chemistry typical of island-arc tholeiites and are therefore most probably derived from eruptions on the South Sandwich Islands. Because of their distant position the southern Andes seem to be rather improbable as a potential source region for the tephra layers investigated.

(Table 1) Occurrences of Epimeria georgiana species complex (Amphipoda) in the Weddel and Scotia Sea

DNA barcoding revealed four well-supported clades among amphipod specimens that keyed out to Epimeria georgiana Schellenberg, 1931, three clades with specimens from the southern Scotia Arc and one clade with specimens from the Weddell Sea. Detailed morphological investigations of sequenced specimens were conducted, through light and scanning electron microscopy. High magnification (500-2,000 fold) revealed features such as comb-scales on the first antenna and trich bearing pits on the fourth coxal plate to be similar for all specimens in the four clades. Consistent microstructure character differences in the Weddell Sea specimens combined with high genetic distances (COI divergence>20%) allowed the description of Epimeria angelikae, a species new to science. Specimens of E. georgiana in the other three COI clades from the Scotia Arc were morphologically indistinguishable. Representative specimens of clade A are also illustrated in detail. Our results on the high genetic divergences in epimeriid amphipods support the theory of the southern Scotia Arc being a centre of Antarctic diversification.

Relative sea-level history of Marguerite Bay, Antarctic Peninsula derived from optically stimulated luminescence-dated beach cobbles

Calmette Bay within Marguerite Bay along the western side of the Antarctic Peninsula contains one of the most continuous flights of raised beaches described to date in Antarctica. Raised beaches extend to 40.8 m above sea level (masl) and are thought to reflect glacial isostatic adjustment due to the retreat of the Antarctic Peninsula Ice Sheet. Using optically stimulated luminescence (OSL), we dated quartz extracts from cobble surfaces buried in raised beaches at Calmette Bay. The beaches are separated into upper and lower beaches based on OSL ages, geomorphology, and sedimentary fabric. The two sets of beaches are separated by a prominent scarp. One of our OSL ages from the upper beaches dates to 9.3 thousand years ago (ka; as of 1950) consistent with previous extrapolation of sea-level data and the time of ice retreat from inner Marguerite Bay. However, four of the seven ages from the upper beaches date to the timing of glaciation. We interpret these ages to represent reworking of beaches deposited prior to the Last Glacial Maximum (LGM) by advancing and retreating LGM ice. Ages from the lower beaches record relative sea-level fall due to Holocene glacial-isostatic adjustment. We suggest a Holocene marine limit of 21.7 masl with an age of 5.5-7.3 ka based on OSL ages from Calmette Bay and other sea-level constraints in the area. A marine limit at 21.7 masl implies half as much relative sea-level change in Marguerite Bay during the Holocene as suggested by previous sea-level reconstructions. No evidence for a relative sea-level signature of neoglacial events, such as a decrease followed by an increase in RSL fall due to ice advance and retreat associated with the Little Ice Age, is found within Marguerite Bay indicating either: (1) no significant neoglacial advances occurred within Marguerite Bay; (2) rheological heterogeneity allows part of the Antarctic Peninsula (i.e. the South Shetland Islands) to respond to rapid ice mass changes while other regions are incapable of responding to short-lived ice advances; or (3) the magnitude of neoglacial events within Marguerite Bay is too small to resolve through relative sea-level reconstructions. Although the application of reconstructing sea-level histories using OSL-dated raised beach deposits provides a better understanding of the timing and nature of relative sea-level change in Marguerite Bay, we highlight possible problems associated with using raised beaches as sea-level indices due to post-depositional reworking by storm waves.

(Table 2) Radiocarbon dates from lakes on Fildes Peninsula, King George Island

Precise relative sea level (RSL) data are important for inferring regional ice sheet histories, as well as helping to validate numerical models of ice sheet evolution and glacial isostatic adjustment. Here we develop a new RSL curve for Fildes Peninsula, South Shetland Islands (SSIs), a sub-Antarctic archipelago peripheral to the northern Antarctic Peninsula ice sheet, by integrating sedimentary evidence from isolation basins with geomorphological evidence from raised beaches. This combined approach yields not only a Holocene RSL curve, but also the spatial pattern of how RSL change varied across the archipelago. The curve shows a mid-Holocene RSL highstand on Fildes Peninsula at 15.5 m above mean sea level between 8000 and 7000 cal a BP. Subsequently RSL gradually fell as a consequence of isostatic uplift in response to regional deglaciation. We propose that isostatic uplift occurred at a non-steady rate, with a temporary pause in ice retreat ca. 7200 cal a BP, leading to a short-lived RSL rise of ~1 m and forming a second peak to the mid-Holocene highstand. Two independent approaches were taken to constrain the long-term tectonic uplift rate of the SSIs at 0.22-0.48 m/ka, placing the tectonic contribution to the reconstructed RSL highstand between 1.4 and 2.9 m. Finally, we make comparisons to predictions from three global sea level models.

Quarternary tephra layers in sediment cores from the Bellingshausen Sea and Scotia Sea

The Southern Ocean is a region of the world's ocean which is fundamental to the generation of cold deep ocean water which drives the global therrno-haline circulation. Previous investigations of deep-sea sediments south of the Polar Front have been significantly constrained by the lack of a suitable correlation and dating technique. In this study, deep-sea sediment cores from the Bellingshausen, Scotia and Weddell seas have been investigated for the presence of tephra layers. The major oxide and trace element composition of glass shards have been used to correlate tephra isochrons over distances in excess of 600 km. The source volcanoes for individual tephra layers have been identified. Atmospheric transport distances greater than 1500 km for >32 pm shards are reported. One megascopic tephra is identified and correlated across 7 sediment drifts on the continental rise in the Bellingshausen Sea. Its occurrence in a sedimentary unit that has been biostratigraphically dated to delta 18O substage 5e identifies it as a key regional marker horizon for that stage. An unusual bimodal megascopic ash layer erupted from Deception Island, South Shetland Islands, has been correlated between 6 sediment cores which form a 600 km NE-SW transect from the central Scotia Sea to Jane Basin. This megascopic ash layer has been 14C dated at c. 10,670 years BP. It represents the last significant input of tephra into the Scotia Sea or Jane Basin from that volcano and forms an important early Holocene marker horizon for the region. Five disseminated tephras can be correlated to varying extents across the central Scotia Sea cores. Together with the megascopic tephra they form a tephrostratigraphic framework that will greatly aid palaeoclimatic, palaeoenvironrnental and palaeoceanographic investigations in the region.