Rapid Volume Loss from Two East Greenland Outlet Glaciers Quantified Using Repeat Stereo Satellite Imagery

The coastal portions of Kangerdlugssuaq and Helheim glaciers in southeast Greenland lost at least 51 +/- 8 km(-3) yr(-1) of ice between 2001-2006 due to thinning and retreat, according to an analysis of sequential digital elevation models (DEMs) derived from stereo ASTER satellite imagery. The dominant contribution to this ice loss was dynamic thinning caused by the acceleration in flow of both glaciers. Peak rates of change, including thinning rates of similar to 90 m yr(-1), coincided with the rapid increases in flow speed. Extrapolation of the measured data to the ice divides yields an estimated combined catchment volume loss of similar to 122 +/- 30 km(-3) yr(-1), which accounts for half the total mass loss from the ice sheet reported in recent studies. These catchment-wide volume losses contributed similar to 0.31 +/- 0.07 mm yr(-1) to global sea level rise over the 5-year observation period with the coastal regions alone contributing at least 0.1 +/- 0.02 mm yr(-1).

Controls on the Basal Water Pressure in Subglacial Channels Near the Margin of the Greenland Ice Sheet

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66 degrees N and 66 degrees 30' N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Rothlisberger's equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.

Patterns of landscape form in the upper Rhône basin, Central Swiss Alps, predominantly show lithologic controls despite multiple glaciations and variations in rock uplift rates

The development of topography depends mainly on the interplay between uplift and erosion. These processes are controlled by various factors including climate, glaciers, lithology, seismic activity and short-term variables, such as anthropogenic impact. Many studies in orogens all over the world have shown how these controlling variables may affect the landscape's topography. In particular, it has been hypothesized that lithology exerts a dominant control on erosion rates and landscape morphology. However, clear demonstrations of this influence are rare and difficult to disentangle from the overprint of other signals such as climate or tectonics. In this study we focus on the upper Rhône Basin situated in the Central Swiss Alps in order to explore the relation between topography, possible controlling variables and lithology in particular. The Rhône Basin has been affected by spatially variable uplift, high orographically driven rainfalls and multiple glaciations. Furthermore, lithology and erodibility vary substantially within the basin. Thanks to high-resolution geological, climatic and topographic data, the Rhône Basin is a suitable laboratory to explore these complexities. Elevation, relief, slope and hypsometric data as well as river profile information from digital elevation models are used to characterize the landscape's topography of around 50 tributary basins. Additionally, uplift over different timescales, glacial inheritance, precipitation patterns and erodibility of the underlying bedrock are quantified for each basin. Results show that the chosen topographic and controlling variables vary remarkably between different tributary basins. We investigate the link between observed topographic differences and the possible controlling variables through statistical analyses. Variations of elevation, slope and relief seem to be linked to differences in long-term uplift rate, whereas elevation distributions (hypsometry) and river profile shapes may be related to glacial imprint. This confirms that the landscape of the Rhône Basin has been highly preconditioned by (past) uplift and glaciation. Linear discriminant analyses (LDAs), however, suggest a stronger link between observed topographic variations and differences in erodibility. We therefore conclude that despite evident glacial and tectonic conditioning, a lithologic control is still preserved and measurable in the landscape of the Rhône tributary basins.

Creating Useful Products From Connecticut's 2000 LIDAR Data Set JHR 08-314 Project 07-2

The State of Connecticut owns a LIght Detection and Ranging (LIDAR) data set that was collected in 2000 as part of the State’s periodic aerial reconnaissance missions. Although collected eight years ago, these data are just now becoming ready to be made available to the public. These data constitute a massive “point cloud”, being a long list of east-north-up triplets in the State Plane Coordinate System Zone 0600 (SPCS83 0600), orthometric heights (NAVD 88) in US Survey feet. Unfortunately, point clouds have no structure or organization, and consequently they are not as useful as Triangulated Irregular Networks (TINs), digital elevation models (DEMs), contour maps, slope and aspect layers, curvature layers, among others. The goal of this project was to provide the computational infrastructure to create a first cut of these products and to serve them to the public via the World Wide Web. The products are available at http://clear.uconn.edu/data/ct_lidar/index.htm.

Snow and ice surfaces on King George Island (Antarctic Peninsula) with high-resolution TerraSAR-X time series, with links to shapfiles and GeoTIFFs

Changes of glaciers and snow cover in polar regions affect a wide range of physical and ecosystem processes on land and in the adjacent marine environment. In this study, we investigate the potential of 11-day repeat high-resolution satellite image time series from the TerraSAR-X mission to derive glaciological and hydrological parameters on King George Island, Antarctica during the period Oct/25/2010 to Apr/19/2011. The spatial pattern and temporal evolution of snow cover extent on ice-free areas can be monitored using multi-temporal coherence images. SAR coherence is used to map glacier extent of land terminating glaciers with an average accuracy of 25 m. Multi-temporal SAR color composites identify the position of the late summer snow line at about 220 m above sea level. Glacier surface velocities are obtained from intensity feature-tracking. Surface velocities near the calving front of Fourcade Glacier were up to 1.8 ± 0.01 m/d. Using an intercept theorem based on fundamental geometric principles together with differential GPS field measurements, the ice discharge of Fourcade Glacier was estimated to 20700 ± 5500 m**3/d (corresponding to ~19 ± 5 kt/d). The rapidly changing surface conditions on King George Island and the lack of high-resolution digital elevation models for the region remain restrictions for the applicability of SAR data and the precision of derived products.

The Austrian Glacier Inventories GI 1 (1969), GI 2 (1998), GI 3 (2006), and GI LIA in ArcGIS (shapefile) format

Glacier inventories provide the basis for further studies on mass balance and volume change, relevant for local hydrological issues as well as for global calculation of sea level rise. In this study, a new Austrian glacier inventory has been compiled, updating data from 1969 (GI 1) and 1998 (GI 2) based on high-resolution lidar digital elevation models (DEMs) and orthophotos dating from 2004 to 2012 (GI 3). To expand the time series of digital glacier inventories in the past, the glacier outlines of the Little Ice Age maximum state (LIA) have been digitalized based on the lidar DEM and orthophotos. The resulting glacier area for GI 3 of 415.11 ± 11.18 km**2 is 44% of the LIA area. The annual relative area losses are 0.3%/yr for the ~119-year period GI LIA to GI 1 with one period with major glacier advances in the 1920s. From GI 1 to GI 2 (29 years, one advance period of variable length in the 1980s) glacier area decreased by 0.6% yr?1 and from GI 2 to GI 3 (10 years, no advance period) by 1.2%/yr. Regional variability of the annual relative area loss is highest in the latest period, ranging from 0.3 to 6.19%/yr. The mean glacier size decreased from 0.69 km**2 (GI 1) to 0.46 km**2 (GI 3), with 47% of the glaciers being smaller than 0.1 km**2 in GI 3 (22%).

Modelos digitales de elevación generados por captura directa e indirecta de datos fuente : evaluación de su calidad

La información básica sobre el relieve de una cuenca hidrográfica, mediante metodologías analítico-descriptivas, permite a quienes evalúan proyectos relacionados con el uso de los recursos naturales, tales como el manejo integrado de cuencas, estudios sobre impacto ambiental, degradación de suelos, deforestación, conservación de los recursos hídricos, entre otros, contar para su análisis con los parámetros físicos necesarios. Estos procesos mencionados tienen un fuerte componente espacial y el empleo de Sistemas de Información Geográfica (SIG) son de suma utilidad, siendo los Modelos Digitales de Elevación (DEM) y sus derivados un componente relevante de esta base de datos. Los productos derivados de estos modelos, como pendiente, orientación o curvatura, resultarán tan precisos como el DEM usado para derivarlos. Por otra parte, es fundamental maximizar la habilidad del modelo para representar las variaciones del terreno; para ello se debe seleccionar una adecuada resolución (grilla) de acuerdo con los datos disponibles para su generación. En este trabajo se evalúa la calidad altimétrica de seis DEMs generados a partir de dos sistemas diferentes de captura de datos fuente y de distintas resoluciones de grilla. Para determinar la exactitud de los DEMs habitualmente se utiliza un grupo de puntos de control considerados como "verdad de campo" que se comparan con los generados por el modelo en la misma posición geográfica. El área seleccionada para realizar el estudio está ubicada en la localidad de Arrecifes, provincia de Buenos Aires (Argentina) y tiene una superficie de aproximadamente 120 ha. Los resultados obtenidos para los dos algoritmos y para los tres tamaños de grilla analizados presentaron los siguientes resultados: el algoritmo DEM from contourn, un RMSE (Root Mean Squared Error) de ± 0,11 m (para grilla de 1 m), ± 0,11 m (para grilla de 5 m) y de ± 0,15 m (para grilla de 10 m). Para el algoritmo DEM from vector/points, un RMSE de ± 0,09 m (para grilla de 1 m), ± 0,11 m (para grilla de 5 m) y de ± 0,11 m (para grilla de 10 m). Los resultados permiten concluir que el DEM generado a partir de puntos acotados del terreno como datos fuente y con el menor tamaño de grilla es el único que satisface los valores enumerados en la bibliografía, tanto nacional como internacional, lo que lo hace apto para proyectos relacionados con recursos naturales a nivel de ecotopo (predial). El resto de los DEMs generados presentan un RMSE que permite asegurar su aptitud para la evaluación de proyectos relacionados con el uso de los recursos naturales a nivel de unidad de paisaje (conjunto de ecotopos).

Determinación de la precisión de DEMs generados con puntos digitalizados a partir de cartografía existente : relación con el algoritmo utilizado

La información básica sobre el relieve de cuencas hidrográficas, mediante metodologías analítico-descriptivas permite contar con los parámetros físicos necesarios para el estudio de los recursos naturales tales como manejo integrado de cuencas, impacto ambiental, degradación de suelos, deforestación, conservación de los recursos hídricos, entre otros. Todos estos procesos ligados a una fuerte componente espacial permiten el uso de sistemas de información geográfica. Los Modelos Digitales de Elevación (DEM) y sus derivados son un componente relevante de estas fuentes de datos. En este trabajo se planteó evaluar los DEMs, generados por captura indirecta de datos fuente, (puntos digitalizados sobre un plano, en formato papel, con curvas de nivel) para determinar si su calidad altimétrica varía con la resolución espacial del mismo (tamaño de grilla). Asimismo se evaluó la posible existencia de alguna interacción entre el software utilizado y los tamaños de grilla, que podría influir sobre la calidad altimétrica del DEM. Se concluyó que la resolución espacial afecta la calidad altimétrica del DEM. Se constató, asimismo, que es posible encontrar un tamaño de grilla ideal en cuanto a calidad altimétrica y peso informático del modelo digital de elevaciones generado. Respecto de la existencia de una interacción entre el tamaño de grilla y el software, que pudieran afectar la calidad de los modelos obtenidos en función de la metodología de adquisición de datos fuente, se demostró que existe esta posibilidad cuando la base de datos se genera a partir de puntos digitalizados.

Ice thickness distribution and glacier bed of Hallstätter Gletscher

Hallstätter Glacier is the northernmost glacier of Austria. Appendant to the northern Limestone Alps, the glacier is located at 47°28'50'' N, 13°36'50'' E in the Dachstein-region. At the same time with its advance linked to the Little Ice Age (LIA), research on changes in size and mass of Hallstätter glacier was started in 1842 by Friedrich Simony. He observed and documented the glacier retreat related to its last maximum extension in 1856. In addition, Hallstätter Glacier is a subject to scientific research to date. In this thesis methods and results of ongoing mass balance measurements are presented and compared to long term volume changes and meteorological observations. The current mass balance monitoring programm using the direct glaciological method was started 2006. In this context, 2009 the ice thickness was measured with ground penetrating radar. The result are used with digital elevation models reconstucted from historical maps and recent digital elevation models to calculate changes in shape and volume of Hallstätter Glacier. Based on current meteorological measurements near the glacier and longtime homogenized climate data provided by HISTALP, time series of precipitation and temperature beginning at the LIA are produced. These monthly precipitation and monthly mean temperature data are used to compare results of a simple degree day model with the volume change calculated from the difference of the digital elevation models. The two years of direct mass balance measurements are used to calibrate the degree day model. A number of possible future scenarios are produced to indicate prospective changes. Within the 150-year-period between 1856 and 2007 the Hallstätter Glacier lost 1940 meters of its length and 2.23 km**2 in area. 37% of the initial volume of 1856 remained. This retreat came along with a change in climate. The application of a running avarage of 30 years shows an increase in precipitation of 18.5% and a warming of 1.3°C near the glacier between 1866 and 1993. The mass loss was continued in the hydrological years 2006/2007 and 2007/2008 showing mean specific mass balance of -376 mm and -700 mm, respectively. Applying a temperature correction for the different minimum elevations of the glacier, the degree day approach based on the two measured mass balances can reproduce sign and order of magnitude of the volume change of Hallstätter Glacier since 1856. Nevertheless, the relative deviation is significant. Future scenarios show, that 30% of the entire glacier volume remains after subtracting the elevation changes between the digital elevation models of 2002 and 2007 ten times from the surface of 2007. The past and present mass changes of Hallstätter Glacier are showing a retreating glacier as a consequence of rising temperatures. Due to high precepitation, increased with previous warming, the Hallstätter Glacier can and will exist in lower elevation compared to inner alpine glaciers.

Soil fragmentation study applying different tillage systems

Runoff generation depends on rainfall, infiltration, interception, and surface depressional storage. Surface depressional storage depends on surface microtopography, usually quantified trough soil surface roughness (SSR). SSR is subject to spatial and temporal changes that create a high variability. In an agricultural environment, tillage operations produce abrupt changes in roughness. Subsequent rainfall gradually decreases roughness. Beside it, local variation in soil properties and hydrology cause its SSR to vary spatially at different scales. The methods commonly used to measure it involve collecting point elevations in regular grids using laser profilers or scanners, digital close range stereo-photogrammetry and terrestrial laser scanning or LIDAR systems. In this case, a laser-scanning instrument was used to obtain representative digital elevation models (DEMs) at a grid resolution of 7.2x7.2mm that cover an area of 0.9x0.9m. The DEMs were obtained from two study sites with different soils. The first study site was an experimental field on which five conventional tillage methods were applied. The second study site was a large olive orchard with trees planted at 7.5x5.0m and bare soils between rows. Here, three tillage treatments were applied. In this work we have evaluated the spatial variability of SSR at several scales studying differences in height calculated from points separated by incremental distances h were raised to power values q (from 0 to 4 in steps of 0.1). The q = 2 data were studied as a semivariogram model. The logarithm of average differences plotted vs. log h were characterized by their slope, ?(q). Structure functions [?(q) vs. q] were fitted showing that data had nonlinear structure functions typical of multiscale phenomena. Comparisson of the two types of soil in their respective structure functions are shown.

Modelización de la zona de inundación del río Forcón-Barranca a su paso por San Miguel de Reinante

El presente proyecto muestra el estudio realizado para evaluar la capacidad de drenaje de la canalización construida en el río Forcón-Barranca, en la población de San Miguel de Reinante (Lugo), a consecuencia de las severas inundaciones acaecidas los días 9 y 10 de junio de 2010. Para evaluar y analizar la peligrosidad se ha seguido un método hidrometeorológico agregado, debido a la falta de aforos, pluviómetros y datos, y se han usado Modelos Digitales de Elevaciones como punto de partida. Todo ello ha llevado a la conclusión de la ineficacia del canal para un período de retorno de 500 años. ABSTRACT After severe flooding in San Miguel de Reinante (Lugo) the 9th and 10th June 2010, a concrete channel was built. This project assesses the drainage capacity of such channel, performing a flood hazard analysis. Because no flow measures or rain gauge data were available a lumped model was followed in order to evaluate flood hazard, using Digital Elevation Models as a starting point. This has led to the conclusion that the channel cannot effectively drain a flow after a rainfall with a return period of 500 years.

Nesting In The Clouds: Evaluating And Predicting Sea Turtle Nesting Beach Parameters From Lidar Data

Humans' desire for knowledge regarding animal species and their interactions with the natural world have spurred centuries of studies. The relatively new development of remote sensing systems using satellite or aircraft-borne sensors has opened up a wide field of research, which unfortunately largely remains dependent on coarse-scale image spatial resolution, particularly for habitat modeling. For habitat-specialized species, such data may not be sufficient to successfully capture the nuances of their preferred areas. Of particular concern are those species for which topographic feature attributes are a main limiting factor for habitat use. Coarse spatial resolution data can smooth over details that may be essential for habitat characterization. Three studies focusing on sea turtle nesting beaches were completed to serve as an example of how topography can be a main deciding factor for certain species. Light Detection and Ranging (LiDAR) data were used to illustrate that fine spatial scale data can provide information not readily captured by either field work or coarser spatial scale sources. The variables extracted from the LiDAR data could successfully model nesting density for loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) sea turtle species using morphological beach characteristics, highlight beach changes over time and their correlations with nesting success, and provide comparisons for nesting density models across large geographic areas. Comparisons between the LiDAR dataset and other digital elevation models (DEMs) confirmed that fine spatial scale data sources provide more similar habitat information than those with coarser spatial scales. Although these studies focused solely on sea turtles, the underlying principles are applicable for many other wildlife species whose range and behavior may be influenced by topographic features.

A data driven study of relationships between relief and farmland abandonment in a Mediterranean region

The environmental, cultural and socio-economic causes and consequences of farmland abandonment are issues of increasing concern for researchers and policy makers. In previous studies, we proposed a new methodology for selecting the driving factors in farmland abandonment processes. Using Data Mining and GIS, it is possible to select those variables which are more significantly related to abandonment. The aim of this study is to investigate the application of the above mentioned methodology for finding relationships between relief and farmland abandonment in a Mediterranean region (SE Spain).We have taken into account up to 28 different variables in a single analysis, some of them commonly considered in land use change studies (slope, altitude, TWI, etc), but also other novel variables have been evaluated (sky view factor, terrain view factor, etc). The variable selection process provides results in line with the previous knowledge of the study area, describing some processes that are region specific (e.g. abandonment versus intensification of the agricultural activities). The European INSPIRE Directive (2007/2/EC) establishes that the digital elevation models for land surfaces should be available in all member countries, this means that the research described in this work can be extrapolated to any European country to determine whether these variables (slope, altitude, etc) are important in the process of abandonment.

Tectonic geomorphology assessment of the Ken River catchment near Panna, Madhya Pradesh, India

This thesis focuses on tectonic geomorphology and the response of the Ken River catchment to postulated tectonic forcing along a NE-striking monocline fold in the Panna region, Madhya Pradesh, India. Peninsular India is underlain by three northeast-trending paleotopographic ridges of Precambrian Indian basement, bounded by crustal-scale faults. Of particular interest is the Pokhara lineament, a crustal scale fault that defines the eastern edge of the Faizabad ridge, a paleotopographic high cored by the Archean Bundelkhand craton. The Pokhara lineament coincides with the monocline structure developed in the Proterozoic Vindhyan Supergroup rocks along the Bundelkhand cratonic margin. A peculiar, deeply incised meander-like feature, preserved along the Ken River where it flows through the monocline, may be intimately related to the tectonic regime of this system. This thesis examines 41 longitudinal stream profiles across the length of the monocline structure to identify any tectonic signals generated from recent surface uplift above the Pokhara lineament. It also investigates the evolution of the Ken River catchment in response to the generation of the monocline fold. Digital Elevation Models (DEM) from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) were used to delineate a series of tributary watersheds and extract individual stream profiles which were imported into MATLAB for analysis. Regression limits were chosen to define distinct channel segments, and knickpoints were defined at breaks between channel segments where there was a discrete change in the steepness of the channel profile. The longitudinal channel profiles exhibit the characteristics of a fluvial system in transient state. There is a significant downstream increase in normalized steepness index in the channel profiles, as well as a general increase in concavity downstream, with some channels exhibiting convex, over-steepened segments. Normalized steepness indices and uppermost knickpoint elevations are on average much higher in streams along the southwest segment of the monocline compared to streams along the northeast segment. Most channel profiles have two to three knickpoints, predominantly exhibiting slope-break morphology. These data have important implications for recent surface uplift above the Pokhara lineament. Furthermore, geomorphic features preserved along the Ken River suggest that it is an antecedent river. The incised meander-like feature appears to be the abandoned river valley of a former Ken River course that was captured during the evolution of the landscape by what is the present day Ken River.

Calculating the Uncertainty of a Structure from Motion (SfM) Model, Cadman Quarry, Monroe, Washington

Structure from Motion (SfM) is a new form of photogrammetry that automates the rendering of georeferenced 3D models of objects using digital photographs and independently surveyed Ground Control Points (GCPs). This project seeks to quantify the error found in Digital Elevation Models (DEMs) produced using SfM. I modeled a rockslide found at the Cadman Quarry (Monroe, Washington) because the surface is vegetation-free, which is ideal for SfM and Terrestrial LiDAR Scanner (TLS) surveys. By using SfM, TLS, and GPS positioning at the same time, I attempted to find the deviation in the SfM model from the TLS model and GPS points. Using the deviation, I found the Root-Mean-Square Error (RMSE) between the SfM DEM and GPS positions. The RMSE of the SfM model when compared to surveyed GPS points is 17cm. I propagated the uncertainty of the GPS points with the RMSE of the SfM model to find the uncertainty of the SfM model compared to the NAD 1984 datum. The uncertainty of the SfM model compared to the NAD 1984 is 27cm. This study did not produce a model from the TLS that had sufficient resolution on horizontal surfaces to compare to surveyed GPS points.