ABSTRACT. Talos Dome is an ice dome on the edge of the East Antarctic plateau; because accumulation is higher here than in other domes of East Antarctica, the ice preserves a good geochemical and palaeoclimatic record. A new map of the Talos Dome area locates the dome summit using the global positioning system (GPS) (72˚47' 14''S, 159˚04' 2'' E; 2318.5 m elevation (WGS84)). A surface strain network of nine stakes was measured using GPS. Data indicate that the stake closest to the summit moves south-southeast at a few cm a Airborne radar measurements indicate that the bedrock at the Talos Dome summit is about 400 m in elevation, and that it is covered by about 1900 m of ice. Snow radar and GPS surveys show that internal layering is continuous and horizontal in the summit area (15 km radius). The depth distribution analysis of snow radar layers reveals that accumulation decreases downwind of the dome (north-northeast) and increases upwind (south-southwest).The palaeomorphology of the dome has changed during the past 500 years, probably due to variation in spatial distribution of snow accumulation, driven by wind sublimation. In order to calculate a preliminary age vs depth profile for Talos Dome, a simple onedimensional steady-state model was formulated. This model predicts that the ice 100 m above the bedrock may cover one glacial-interglacial period.

SUMMARYAn accurate, efficient and robust numerical method for the solution of the section-averaged De St. Venant equations of open channel flow is presented and discussed. The method consists in a semi-implicit, finitevolume discretization of the continuity equation capable to deal with arbitrary cross-section geometry and in a semi-implicit, finite-difference discretization of the momentum equation. By using a proper semiLagrangian discretization of the momentum equation, a highly efficient scheme that is particularly suitable for subcritical regimes is derived. Accurate solutions are obtained in all regimes, except in presence of strong unsteady shocks as in dam-break cases. By using a suitable upwind, Eulerian discretization of the same equation, instead, a scheme capable of describing accurately also unsteady shocks can be obtained, although this scheme requires to comply with a more restrictive stability condition. The formulation of the two approaches allows a unified implementation and an easy switch between the two. The code is verified in a wide range of idealized test cases, highlighting its accuracy and efficiency characteristics, especially for long time range simulations of subcritical river flow. Finally, a model validation on field data is presented, concerning simulations of a flooding event of the Adige river.

The problem of describing ice dynamics has been faced by many researchers; in this paper a fully three-dimensional model for ice dynamics is presented and tested. Using an approach followed by other researchers, ice is considered a non-linear incompressible viscous fluid so that a fluid-dynamic approach can be used. The model is based on the full three-dimensional Stokes equations for the description of pressure and velocity fields, on the Saint-Venant equation for the description of the freesurface time evolution and on a constitutive law derived from Glen's law for the description of ice viscosity. The model computes the complete pressure field by considering both the hydrostatic and hydrodynamic pressure components; it is time-evolutive and uses high-order numerical approximation for equations and boundary conditions. Moreover it can deal with both constant and variable viscosity. Three theoretical tests and two applications to Priestley Glacier, Antarctica, are presented in order to evaluate the performance of the model and to investigate important phenomena of ice dynamics such as the influence of viscosity on pressure and velocity fields, basal sliding and flow over perturbed bedrocks. All these applications demonstrate the importance of treating the complete pressure and stress fields.

SUMMARYA semi-implicit numerical method for the three-dimensional incompressible Navier-Stokes equations is presented. The method describes the velocity and the non-hydrostatic pressure ÿeld and the free surface evolution in time. The governing equations are discretized by means of the ÿnite volume method on a structured non-uniform grid; for this reason both local and global mass conservation are guaranteed. Convective and di usive uxes on the control volume faces, as well as boundary conditions, are approximated by high-order formulae.