A compact lightweight multipurpose ground-penetrating radar for glaciological applications

We describe a compact lightweight impulse radar for radio-echo sounding of subsurface structures designed specifically for glaciological applications. The radar operates at frequencies between 10 and 75 MHz. Its main advantages are that it has a high signal-to-noise ratio and a corresponding wide dynamic range of 132 dB due mainly to its ability to perform real-time stacking (up to 4096 traces) as well as to the high transmitted power (peak voltage 2800 V). The maximum recording time window, 40 ?s at 100 MHz sampling frequency, results in possible radar returns from as deep as 3300 m. It is a versatile radar, suitable for different geophysical measurements (common-offset profiling, common midpoint, transillumination, etc.) and for different profiling set-ups, such as a snowmobile and sledge convoy or carried in a backpack and operated by a single person. Its low power consumption (6.6 W for the transmitter and 7.5 W for the receiver) allows the system to operate under battery power for mayor que7 hours with a total weight of menor que9 kg for all equipment, antennas and batteries.

Observation and Detection of Underground Cellars in the Duero Basin By GNSS, Lidar and Ground Penetrating Radar Techniques

The underground cellars that appear in different parts of Spain are part of an agricultural landscape dispersed, sometimes damaged, others at risk of disappearing. This paper studies the measurement and display of a group of wineries located in Atauta (Soria), in the Duero River corridor. It is a unique architectural complex, facing rising, built on a smooth hillock as shown in Fig. 1. These constructions are excavated in the ground. The access to the cave or underground cellar has a shape of a narrow tube or down gallery. Immediately after, this space gets wider. There, wine is produced and stored [1]. Observation and detection of the underground cellar, both on the outside and underground, it is essential to make an inventory of the rural patrimony [2]. The geodetection is a noninvasive technique, adequate to accurately locate buried structures in the ground. Works undertaken include topographic work with the LIDAR techniques and integration with data obtained by GNSS and GPR.

Ice volume changes (1936-1990-2007) and ground-penetrating radar studies of Ariebreen, Hornsund, Spitsbergen

Ariebreen is a small (0.37 km2) valley glacier located in southern Spitsbergen. Our ground-penetrating radar surveys of the glacier show that it is less than 30 m thick on average, with a maximum thickness of 82 m, and it appears to be entirely cold. By analysing digital terrain models of the ice surface from different dates, we determine the area and volume changes during two periods, 1936-1990 and 1990-2007. The total ice volume of the glacier has decreased by 73% during the entire period 1936-2007, which is equivalent to a mean mass balance rate of -0.6190.17 m/yr w.eq. The glacier thinning rate has increased markedly between the first and second periods, from -0.5090.22 to -0.9590.17 m/yr w.eq.

Ice volume estimates from ground-penetrating radar surveys, western Nordenskiöld Land glaciers, Svalbard

As part of ongoing work within the SvalGlac project aimed to obtain a reliable estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea level rise, in this contribution we present volume calculations, with detailed error estimates, for ten glaciers on western Nordenskiöld Land, central Spitsbergen, Svalbard. The volume estimates are based upon a dense net of GPR-retrieved ice thickness data collected over several field campaigns spanning the period 1999-2012, all of them except one within 2010-2012. The total area and volume of the ensemble are 113.38±0.09 km2 and 10.439±0.185 km3, respectively, while the individual areas, volumes and average ice thickness lie within 2.5-49.1 km2, 0.08-5.48 km3 and 29-108 m, respectively. The maximum recorded ice thickness, 265±15 m, corresponds to Fridtjovbreen, which has also the largest average thickness (108±1m). Available empirical formulae for Svalbard glaciers overestimate the total volume of these glaciers by 24% with respect to our calculation. On the basis of the pattern of scattering in the radargrams, we also analyse the hydrothermal structure of these glaciers. Nine out of ten are polythermal, while only one is entirely cold.

Ice volume estimates from ground-penetrating radar surveys, Wedel Jarlsberg Land glaciers, Svalbard

One of the aims of the SvalGlac project is to obtain an improved estimate, with reliable error estimates, of the volume of Svalbard glaciers and their potential contribution to sea level rise. As part of this work, we present volume calculations, with detailed error estimates, for eight glaciers on Wedel Jarlsberg Land, southern Spitsbergen, Svalbard. The volume estimates are based upon a dense net of GPR-retrieved ice thickness data collected over several field campaigns spanning the period 2004-2011. The total area and volume of the ensemble are 502.9±18.6 km2 and 80.72±2.85 km3, respectively. Excluding Ariebreen (a tiny glacier, menor que 0.4 km2 in area), the individual areas, volumes and average ice thickness lie within 4.7-141.0 km2, 0.30-25.85 km3 and 64-183 m, respectively. The maximum recorded ice thickness, ca. 619±13 m, is found in Austre Torellbreen. To estimate the ice volume of small non-echo-sounded tributary glaciers, we used a function providing the best fit to the ice thickness along the centre line of a collection of such tributaries where echo-soundings were available, and assuming parabolic cross-sections. We did some tests on the effect on the measured ice volumes of the distinct radio-wave velocity (RWV) of firn as compared to ice, and cold versus temperate ice, concluding that the changes in volume implied by such corrections were within the error bounds of our volume estimate using a constant RWV for the entire glacier inferred from common mid-point measurements on the upper ablation area.

Ground-penetrating radar surveys and ice volume estimates of Wedel Jarlsberg Land glaciers, Svalbard

We present volume calculations, with detailed error estimates, for eight glaciers on Wedel Jarlsberg Land, southern Spitsbergen, Svalbard, and compare them to those obtained from area-volume scaling relationships. The volume estimates are based upon a dense net of GPR-retrieved ice thickness data collected over several field campaigns spanning the period 2004-2011.

Ice volume estimates from ground-penetrating radar surveys, Wedel Jarlsberg Land glaciers, Svalbard

We present ground-penetrating radar (GPR)—based volume calculations, with associated error estimates, for eight glaciers on Wedel Jarlsberg Land, southwestern Spitsbergen, Svalbard, and compare them with those obtained from volume-area scaling relationships. The volume estimates are based upon GPR ice-thickness data collected during the period 2004–2013. The total area and volume of the ensemble are 502.91 ± 18.60 km2 and 91.91 ± 3.12 km3, respectively. The individual areas, volumes, and average ice thickness lie within 0.37–140.99 km2, 0.01–31.98 km3, and 28–227 m, respectively, with a maximum recorded ice thickness of 619 ± 13 m on Austre Torellbreen. To estimate the ice volume of unsurveyed tributary glaciers, we combine polynomial cross-sections with a function providing the best fit to the measured ice thickness along the center line of a collection of 22 surveyed tributaries. For the time-to-depth conversion of GPR data, we test the use of a glacierwide constant radio-wave velocity chosen on the basis of local or regional common midpoint measurements, versus the use of distinct velocities for the firn, cold ice, and temperate ice layers, concluding that the corresponding volume calculations agree with each other within their error bounds.

Evaluación de la profundidad de la fisuración superficial descendente en pavimentos asfálticos mediante técnicas de ultrasonidos. Validación teórico-práctica y modelos

La fisuración iniciada en la superficie de los pavimentos asfálticos constituye uno de los más frecuentes e importantes modos de deterioro que tienen lugar en los firmes bituminosos, como han demostrado los estudios teóricos y experimentales llevados a cabo en la última década. Sin embargo, este mecanismo de fallo no ha sido considerado por los métodos tradicionales de diseño de estos firmes. El concepto de firmes de larga duración se fundamenta en un adecuado seguimiento del proceso de avance en profundidad de estos deterioros y la intervención en el momento más apropiado para conseguir mantenerlos confinados como fisuras de profundidad parcial en la capa superficial más fácilmente accesible y reparable, de manera que pueda prolongarse la durabilidad y funcionalidad del firme y reducir los costes generalizados de su ciclo de vida. Por lo tanto, para la selección de la estrategia óptima de conservación de los firmes resulta esencial disponer de metodologías que posibiliten la identificación precisa in situ de la fisuración descendente, su seguimiento y control, y que además permitan una determinación fiable y con alto rendimiento de su profundidad y extensión. En esta Tesis Doctoral se presentan los resultados obtenidos mediante la investigación sistemática de laboratorio e in situ llevada a cabo para la obtención de datos sobre fisuración descendente en firmes asfálticos y para el estudio de procedimientos de evaluación de la profundidad de este tipo de fisuras empleando técnicas de ultrasonidos. Dichos resultados han permitido comprobar que la metodología no destructiva propuesta, de rápida ejecución, bajo coste y sencilla implementación (principalmente empleada hasta el momento en estructuras metálicas y de hormigón, debido a las dificultades que introduce la naturaleza viscoelástica de los materiales bituminosos) puede ser aplicada con suficiente fiabilidad y repetibilidad sobre firmes asfálticos. Las medidas resultan asimismo independientes del espesor total del firme. Además, permite resolver algunos de los inconvenientes frecuentes que presentan otros métodos de diagnóstico de las fisuras de pavimentos, tales como la extracción de testigos (sistema destructivo, de alto coste y prolongados tiempos de interrupción del tráfico) o algunas otras técnicas no destructivas como las basadas en medidas de deflexiones o el georradar, las cuales no resultan suficientemente precisas para la investigación de fisuras superficiales. Para ello se han realizado varias campañas de ensayos sobre probetas de laboratorio en las que se han estudiado diferentes condiciones empíricas como, por ejemplo, distintos tipos de mezclas bituminosas en caliente (AC, SMA y PA), espesores de firme y adherencias entre capas, temperaturas, texturas superficiales, materiales de relleno y agua en el interior de las grietas, posición de los sensores y un amplio rango de posibles profundidades de fisura. Los métodos empleados se basan en la realización de varias medidas de velocidad o de tiempo de transmisión del pulso ultrasónico sobre una única cara o superficie accesible del material, de manera que resulte posible obtener un coeficiente de transmisión de la señal (mediciones relativas o autocompensadas). Las mediciones se han realizado a bajas frecuencias de excitación mediante dos equipos de ultrasonidos diferentes dotados, en un caso, de transductores de contacto puntual seco (DPC) y siendo en el otro instrumento de contacto plano a través de un material especialmente seleccionado para el acoplamiento (CPC). Ello ha permitido superar algunos de los tradicionales inconvenientes que presenta el uso de los transductores convencionales y no precisar preparación previa de las superficies. La técnica de autocalibración empleada elimina los errores sistemáticos y la necesidad de una calibración local previa, demostrando el potencial de esta tecnología. Los resultados experimentales han sido comparados con modelos teóricos simplificados que simulan la propagación de las ondas ultrasónicas en estos materiales bituminosos fisurados, los cuales han sido deducidos previamente mediante un planteamiento analítico y han permitido la correcta interpretación de dichos datos empíricos. Posteriormente, estos modelos se han calibrado mediante los resultados de laboratorio, proporcionándose sus expresiones matemáticas generalizadas y gráficas para su uso rutinario en las aplicaciones prácticas. Mediante los ensayos con ultrasonidos efectuados en campañas llevadas a cabo in situ, acompañados de la extracción de testigos del firme, se han podido evaluar los modelos propuestos. El máximo error relativo promedio en la estimación de la profundidad de las fisuras al aplicar dichos modelos no ha superado el 13%, con un nivel de confianza del 95%, en el conjunto de todos los ensayos realizados. La comprobación in situ de los modelos ha permitido establecer los criterios y las necesarias recomendaciones para su utilización sobre firmes en servicio. La experiencia obtenida posibilita la integración de esta metodología entre las técnicas de auscultación para la gestión de su conservación. Abstract Surface-initiated cracking of asphalt pavements constitutes one of the most frequent and important types of distress that occur in flexible bituminous pavements, as clearly has been demonstrated in the technical and experimental studies done over the past decade. However, this failure mechanism has not been taken into consideration for traditional methods of flexible pavement design. The concept of long-lasting pavements is based on adequate monitoring of the depth and extent of these deteriorations and on intervention at the most appropriate moment so as to contain them in the surface layer in the form of easily-accessible and repairable partial-depth topdown cracks, thereby prolonging the durability and serviceability of the pavement and reducing the overall cost of its life cycle. Therefore, to select the optimal maintenance strategy for perpetual pavements, it becomes essential to have access to methodologies that enable precise on-site identification, monitoring and control of top-down propagated cracks and that also permit a reliable, high-performance determination of the extent and depth of cracking. This PhD Thesis presents the results of systematic laboratory and in situ research carried out to obtain information about top-down cracking in asphalt pavements and to study methods of depth evaluation of this type of cracking using ultrasonic techniques. These results have demonstrated that the proposed non-destructive methodology –cost-effective, fast and easy-to-implement– (mainly used to date for concrete and metal structures, due to the difficulties caused by the viscoelastic nature of bituminous materials) can be applied with sufficient reliability and repeatability to asphalt pavements. Measurements are also independent of the asphalt thickness. Furthermore, it resolves some of the common inconveniences presented by other methods used to evaluate pavement cracking, such as core extraction (a destructive and expensive procedure that requires prolonged traffic interruptions) and other non-destructive techniques, such as those based on deflection measurements or ground-penetrating radar, which are not sufficiently precise to measure surface cracks. To obtain these results, extensive tests were performed on laboratory specimens. Different empirical conditions were studied, such as various types of hot bituminous mixtures (AC, SMA and PA), differing thicknesses of asphalt and adhesions between layers, varied temperatures, surface textures, filling materials and water within the crack, different sensor positions, as well as an ample range of possible crack depths. The methods employed in the study are based on a series of measurements of ultrasonic pulse velocities or transmission times over a single accessible side or surface of the material that make it possible to obtain a signal transmission coefficient (relative or auto-calibrated readings). Measurements were taken at low frequencies by two short-pulse ultrasonic devices: one equipped with dry point contact transducers (DPC) and the other with flat contact transducers that require a specially-selected coupling material (CPC). In this way, some of the traditional inconveniences presented by the use of conventional transducers were overcome and a prior preparation of the surfaces was not required. The auto-compensating technique eliminated systematic errors and the need for previous local calibration, demonstrating the potential for this technology. The experimental results have been compared with simplified theoretical models that simulate ultrasonic wave propagation in cracked bituminous materials, which had been previously deduced using an analytical approach and have permitted the correct interpretation of the aforementioned empirical results. These models were subsequently calibrated using the laboratory results, providing generalized mathematical expressions and graphics for routine use in practical applications. Through a series of on-site ultrasound test campaigns, accompanied by asphalt core extraction, it was possible to evaluate the proposed models, with differences between predicted crack depths and those measured in situ lower than 13% (with a confidence level of 95%). Thereby, the criteria and the necessary recommendations for their implementation on in-service asphalt pavements have been established. The experience obtained through this study makes it possible to integrate this methodology into the evaluation techniques for pavement management systems.

Estudio de la variación del contenido de humedad en el Pinus pinaster Ait. por medio de la técnica no destructiva del georradar / Assessment of the variation of the moisture content in the Pinus pinaster Ait using the non destructive gpr technique

La variación en el contenido de humedad (CH) tiene una influencia significativa tanto en las propiedades físico- químicas de la madera, como en sus propiedades electromagnéticas y por tanto afecta a las características de la propagación de las ondas. En concreto, en este trabajo se estudia la capacidad del georradar (GR) empleando una antena de 1.6GHz de frecuencia central para registrar las variaciones que se producen en la velocidad y en la amplitud de las ondas electromagnéticas cuando se propagan en unas viguetas de madera de Pinus pinaster Ait de uso estructural cuyo CH va disminuyendo. Se ha comprobado como cuando el CH descendía la velocidad de propagación y las amplitudes, tanto de la onda directa como de la reflejada aumentaba. Los altos factores de correlación encontrados demuestran que el GR es una técnica capaz de evaluar, de forma no destructiva, el CH de la madera de uso estructural. The moisture content variations in wood have a significant influence in wood?s physicochemical properties, as well as in its electromagnetic properties and to specific effects upon waves? characteristics. In particular, this paper focuses on the analysis of the Ground-penetrating Radar?s (GPR) using an antenna of 1.6 GHz central frequency capacity to register the velocity and the amplitude of the electromagnetic waves? variation during the drying process of Pinus pinaster Ait timber joists. The results showed that when timber MC descended, the propagation velocity and amplitude of both the direct and the reflected wave increased. The high correlation found between the variables studied demonstrates GPR efficiency and the innovative application of this technique as a non-destructive evaluation tool for timber structures, particularly when studying its moisture content.

Evaluación de la madera estructural para la rehabilitación de las Naves Cros (Valencia, España) por medio de diferentes técnicas no destructivas

El proyecto de rehabilitación de una de las naves del complejo fabril de la industria química ?CROS? en Valencia se llevó a cabo con el criterio de mantener, en la medida de lo posible, los elementos estructurales presentes en la nave. Con este objetivo se realizaron una serie de ensayos no destructivos (END) in situ. Estos ensayos permitieron evaluar la calidad de la madera, determinar qué elementos estructurales debían ser sustituidos y comprobar la aptitud de los que iban a ser reutilizados. Los END empleados en este estudio fueron los siguientes: (1) Identificación de la especie por técnicas anatómicas, (2) Clasificación resistente por método visual, (3) Estimación de humedad por la técnica de resistencia eléctrica; (4) Obtención de velocidades de propagación ultrasónicas (5) Resistógrafía y (6) Alteración de la propagación de ondas electromagnéticas por medio de Georradar. Para la calibración de estos END se tomó una muestra de piezas y se hicieron ensayos destructivos bajo condiciones controladas en laboratorio. En el trabajo que aquí se presenta se muestra la metodología empleada durante el proceso de toma de datos, de análisis de resultados y de cruce de la información obtenida a partir de cada uno de los ensayos hasta llegar a un diagnóstico para los elementos analizados. The assessment of structural timber was requested in the rehabilitation project of the Naves of the chemical industry "CROS". The criterion was to maintain as much as possible timber of the structure and to make only partial replacements. In order not to damage the existing structure and to assess the quality of the existing timber, a series of non-destructive testing (NDT) in the entire structure were performed: (1) Identification of the species by anatomical techniques, (2) Strength grading by visual method, (3) Estimation of moisture content by the technique of electrical resistance, (4) Acquisition of ultrasonic propagation velocities (5) Resistography and (6) Record of the propagation of electromagnetic waves by means of Ground-penetrating radar. Following, a sample group was chose to carry out destructive testing in the lab and compare the NDT results with those obtained with the standard UNE-EN408 (modules of strength, stiffness and density). In the present work, the results provided by each of the NDT techniques are detailed and above all, what is more important, the validity of these after they have been contrasted with the destructive standard tests.

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.

Integración de modelos numéricos de glaciares y procesado de datos de georradar en un sistema de información geográfica

Los modelos de termomecánica glaciar están definidos mediante sistemas de ecuaciones en derivadas parciales que establecen los principios básicos de conservación de masa, momento lineal y energía, acompañados por una ley constitutiva que define la relación entre las tensiones a las que está sometido el hielo glaciar y las deformaciones resultantes de las mismas. La resolución de estas ecuaciones requiere la definición precisa del dominio (la geometría del glaciar, obtenido a partir de medidas topográficas y de georradar), así como contar con un conjunto de condiciones de contorno, que se obtienen a partir de medidas de campo de las variables implicadas y que constituyen un conjunto de datos geoespaciales. El objetivo fundamental de esta tesis es desarrollar una serie de herramientas que nos permitan definir con precisión la geometría del glaciar y disponer de un conjunto adecuado de valores de las variables a utilizar como condiciones de contorno del problema. Para ello, en esta tesis se aborda la recopilación, la integración y el estudio de los datos geoespaciales existentes para la Península Hurd, en la Isla Livingston (Antártida), generados desde el año 1957 hasta la actualidad, en un sistema de información geográfica. Del correcto tratamiento y procesamiento de estos datos se obtienen otra serie de elementos que nos permiten realizar la simulación numérica del régimen termomecánico presente de los glaciares de Península Hurd, así como su evolución futura. Con este objetivo se desarrolla en primer lugar un inventario completo de datos geoespaciales y se realiza un procesado de los datos capturados en campo, para establecer un sistema de referencia común a todos ellos. Se unifican además todos los datos bajo un mismo formato estándar de almacenamiento e intercambio de información, generándose los metadatos correspondientes. Se desarrollan asimismo técnicas para la mejora de los procedimientos de captura y procesado de los datos, de forma que se minimicen los errores y se disponga de estimaciones fiables de los mismos. El hecho de que toda la información se integre en un sistema de información geográfica (una vez producida la normalización e inventariado de la misma) permite su consulta rápida y ágil por terceros. Además, hace posible efectuar sobre ella una serie de operaciones conducentes a la obtención de nuevas capas de información. El análisis de estos nuevos datos permite explicar el comportamiento pasado de los glaciares objeto de estudio y proporciona elementos esenciales para la simulación de su comportamiento futuro. ABSTRACT Glacier thermo-mechanical models are defined by systems of partial differential equations stating the basic principles of conservation of mass, momentum and energy, accompanied by a constitutive principle that defines the relationship between the stresses acting on the ice and the resulting deformations. The solution of these equations requires an accurate definition of the model domain (the geometry of the glacier, obtained from topographical and ground penetrating radar measurements), as well as a set of boundary conditions, which are obtained from measurements of the variables involved and define a set of geospatial data. The main objective of this thesis is to develop tools able to provide an accurate definition of the glacier geometry and getting a proper set of values for the variables to be used as boundary conditions of our problem. With the above aim, this thesis focuses on the collection, compilation and study of the geospatial data existing for the Hurd Peninsula on Livingston Island, Antarctica, generated since 1957 to present, into a geographic information system. The correct handling and processing of these data results on a new collection of elements that allow us to numerically model the present state and the future evolution of Hurd Peninsula glaciers. First, a complete inventory of geospatial data is developed and the captured data are processed, with the aim of establishing a reference system common to all collections of data. All data are stored under a common standard format, and the corresponding metadata are generated to facilitate the information exchange. We also develop techniques for the improvement of the procedures used for capturing and processing the data, such that the errors are minimized and better estimated. All information is integrated into a geographic information system (once produced the standardization and inventory of it). This allows easy and fast viewing and consulting of the data by third parties. Also, it is possible to carry out a series of operations leading to the production of new layers of information. The analysis of these new data allows to explain past glacier behavior, and provides essential elements for explaining its future evolution.

Characterization of Cavities Using the GPR, LIDAR and GNSS Techniques

The study of the many types of natural and manmade cavities in different parts of the world is important to the fields of geology, geophysics, engineering, architectures, agriculture, heritages and landscape. Ground-penetrating radar (GPR) is a noninvasive geodetection and geolocation technique suitable for accurately determining buried structures. This technique requires knowing the propagation velocity of electromagnetic waves (EM velocity) in the medium. We propose a method for calibrating the EM velocity using the integration of laser imaging detection and ranging (LIDAR) and GPR techniques using the Global Navigation Satellite System (GNSS) as support for geolocation. Once the EM velocity is known and the GPR profiles have been properly processed and migrated, they will also show the hidden cavities and the old hidden structures from the cellar. In this article, we present a complete study of the joint use of the GPR, LIDAR and GNSS techniques in the characterization of cavities. We apply this methodology to study underground cavities in a group of wine cellars located in Atauta (Soria, Spain). The results serve to identify construction elements that form the cavity and group of cavities or cellars. The described methodology could be applied to other shallow underground structures with surface connection, where LIDAR and GPR profiles could be joined, as, for example, in archaeological cavities, sewerage systems, drainpipes, etc.

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.

Estimate of the total volume of Svalbard glaciers, and their potential contribution to sea-level rise, using new regionally based scaling relationships

We present a set of new volume scaling relationships specific to Svalbard glaciers, derived from a sample of 60 volume–area pairs. Glacier volumes are computed from ground-penetrating radar (GPR)-retrieved ice thickness measurements, which have been compiled from different sources for this study. The most precise scaling models, in terms of lowest cross-validation errors, are obtained using a multivariate approach where, in addition to glacier area, glacier length and elevation range are also used as predictors. Using this multivariate scaling approach, together with the Randolph Glacier Inventory V3.2 for Svalbard and Jan Mayen, we obtain a regional volume estimate of 6700 ± 835 km3, or 17 ± 2 mm of sea-level equivalent (SLE). This result lies in the mid- to low range of recently published estimates, which show values as varied as 13 and 24 mm SLE. We assess the sensitivity of the scaling exponents to glacier characteristics such as size, aspect ratio and average slope, and find that the volume of steep-slope and cirque-type glaciers is not very sensitive to changes in glacier area.