Debris-flow activity in abandoned channels of the Manival torrent reconstructed with LiDAR and tree-ring data

Hydrogeomorphic processes are a major threat in many parts of the Alps, where they periodically damage infrastructure, disrupt transportation corridors or even cause loss of life. Nonetheless, past torrential activity and the analysis of areas affected during particular events remain often imprecise. It was therefore the purpose of this study to reconstruct spatio-temporal patterns of past debris-flow activity in abandoned channels on the forested cone of the Manival torrent (Massif de la Chartreuse, French Prealps). A Light Detecting and Ranging (LiDAR) generated Digital Elevation Model (DEM) was used to identify five abandoned channels and related depositional forms (lobes, lateral levees) in the proximal alluvial fan of the torrent. A total of 156 Scots pine trees (Pinus sylvestris L.) with clear signs of debris flow events was analyzed and growth disturbances (GD) assessed, such as callus tissue, the onset of compression wood or abrupt growth suppression. In total, 375 GD were identified in the tree-ring samples, pointing to 13 debris-flow events for the period 1931–2008. While debris flows appear to be very common at Manival, they have only rarely propagated outside the main channel over the past 80 years. Furthermore, analysis of the spatial distribution of disturbed trees contributed to the identification of four patterns of debris-flow routing and led to the determination of three preferential breakout locations. Finally, the results of this study demonstrate that the temporal distribution of debris flows did not exhibit significant variations since the beginning of the 20th century.

A high-resolution soil erosion risk map of Switzerland as strategic policy support system

Soil erosion models and soil erosion risk maps are often used as indicators to assess potential soil erosion in order to assist policy decisions. This paper shows the scientific basis of the soil erosion risk map of Switzerland and its application in policy and practice. Linking a USLE/RUSLE-based model approach (AVErosion) founded on multiple flow algorithms and the unit contributing area concept with an extremely precise and high-resolution digital terrain model (2 m × 2 m grid) using GIS allows for a realistic assessment of the potential soil erosion risk, on single plots, i.e. uniform and comprehensive for the agricultural area of Switzerland (862,579 ha in the valley area and the lower mountain regions). The national or small-scale soil erosion prognosis has thus reached a level heretofore possible only in smaller catchment areas or single plots. Validation was carried out using soil loss data from soil erosion damage mappings in the field from long-term monitoring in different test areas. 45% of the evaluated agricultural area of Switzerland was classified as low potential erosion risk, 12% as moderate potential erosion risk, and 43% as high potential erosion risk. However, many of the areas classified as high potential erosion risk are located at the transition from valley to mountain zone, where many areas are used as permanent grassland, which drastically lowers their current erosion risk. The present soil erosion risk map serves on the one hand to identify and prioritise the high-erosion risk areas, and on the other hand to promote awareness amongst farmers and authorities. It was published on the internet and will be made available to the authorities in digital form. It is intended as a tool for simplifying and standardising enforcement of the legal framework for soil erosion prevention in Switzerland. The work therefore provides a successful example of cooperation between science, policy and practice.

Automatic extraction of the mid-facial plane for cranio-maxillofacial surgery planning

Recently developed computer applications provide tools for planning cranio-maxillofacial interventions based on 3-dimensional (3D) virtual models of the patient's skull obtained from computed-tomography (CT) scans. Precise knowledge of the location of the mid-facial plane is important for the assessment of deformities and for planning reconstructive procedures. In this work, a new method is presented to automatically compute the mid-facial plane on the basis of a surface model of the facial skeleton obtained from CT. The method matches homologous surface areas selected by the user on the left and right facial side using an iterative closest point optimization. The symmetry plane which best approximates this matching transformation is then computed. This new automatic method was evaluated in an experimental study. The study included experienced and inexperienced clinicians defining the symmetry plane by a selection of landmarks. This manual definition was systematically compared with the definition resulting from the new automatic method: Quality of the symmetry planes was evaluated by their ability to match homologous areas of the face. Results show that the new automatic method is reliable and leads to significantly higher accuracy than the manual method when performed by inexperienced clinicians. In addition, the method performs equally well in difficult trauma situations, where key landmarks are unreliable or absent.

Spatial Assessment of Erosion and its Impact on Soil Fertility in the Tajik Foothills

Efficient planning of soil conservation measures requires, first, to understand the impact of soil erosion on soil fertility with regard to local land cover classes; and second, to identify hot spots of soil erosion and bright spots of soil conservation in a spatially explicit manner. Soil organic carbon (SOC) is an important indicator of soil fertility. The aim of this study was to conduct a spatial assessment of erosion and its impact on SOC for specific land cover classes. Input data consisted of extensive ground truth, a digital elevation model and Landsat 7 imagery from two different seasons. Soil spectral reflectance readings were taken from soil samples in the laboratory and calibrated with results of SOC chemical analysis using regression tree modelling. The resulting model statistics for soil degradation assessments are promising (R2=0.71, RMSEV=0.32). Since the area includes rugged terrain and small agricultural plots, the decision tree models allowed mapping of land cover classes, soil erosion incidence and SOC content classes at an acceptable level of accuracy for preliminary studies. The various datasets were linked in the hot-bright spot matrix, which was developed to combine soil erosion incidence information and SOC content levels (for uniform land cover classes) in a scatter plot. The quarters of the plot show different stages of degradation, from well conserved land to hot spots of soil degradation. The approach helps to gain a better understanding of the impact of soil erosion on soil fertility and to identify hot and bright spots in a spatially explicit manner. The results show distinctly lower SOC content levels on large parts of the test areas, where annual crop cultivation was dominant in the 1990s and where cultivation has now been abandoned. On the other hand, there are strong indications that afforestations and fruit orchards established in the 1980s have been successful in conserving soil resources.

TanDEM-X High Resolution DEMs and Their Applications to Flow Modeling

Lava flow modeling can be a powerful tool in hazard assessments; however, the ability to produce accurate models is usually limited by a lack of high resolution, up-to-date Digital Elevation Models (DEMs). This is especially obvious in places such as Kilauea Volcano (Hawaii), where active lava flows frequently alter the terrain. In this study, we use a new technique to create high resolution DEMs on Kilauea using synthetic aperture radar (SAR) data from the TanDEM-X (TDX) satellite. We convert raw TDX SAR data into a geocoded DEM using GAMMA software [Werner et al., 2000]. This process can be completed in several hours and permits creation of updated DEMs as soon as new TDX data are available. To test the DEMs, we use the Harris and Rowland [2001] FLOWGO lava flow model combined with the Favalli et al. [2005] DOWNFLOW model to simulate the 3-15 August 2011 eruption on Kilauea's East Rift Zone. Results were compared with simulations using the older, lower resolution 2000 SRTM DEM of Hawaii. Effusion rates used in the model are derived from MODIS thermal infrared satellite imagery. FLOWGO simulations using the TDX DEM produced a single flow line that matched the August 2011 flow almost perfectly, but could not recreate the entire flow field due to the relatively high DEM noise level. The issues with short model flow lengths can be resolved by filtering noise from the DEM. Model simulations using the outdated SRTM DEM produced a flow field that followed a different trajectory to that observed. Numerous lava flows have been emplaced at Kilauea since the creation of the SRTM DEM, leading the model to project flow lines in areas that have since been covered by fresh lava flows. These results show that DEMs can quickly become outdated on active volcanoes, but our new technique offers the potential to produce accurate, updated DEMs for modeling lava flow hazards.

An integrated system for 3D hip joint reconstruction from 2D X-rays: a preliminary validation study

The acquisition of conventional X-ray radiographs remains the standard imaging procedure for the diagnosis of hip-related problems. However, recent studies demonstrated the benefit of using three-dimensional (3D) surface models in the clinical routine. 3D surface models of the hip joint are useful for assessing the dynamic range of motion in order to identify possible pathologies such as femoroacetabular impingement. In this paper, we present an integrated system which consists of X-ray radiograph calibration and subsequent 2D/3D hip joint reconstruction for diagnosis and planning of hip-related problems. A mobile phantom with two different sizes of fiducials was developed for X-ray radiograph calibration, which can be robustly detected within the images. On the basis of the calibrated X-ray images, a 3D reconstruction method of the acetabulum was developed and applied together with existing techniques to reconstruct a 3D surface model of the hip joint. X-ray radiographs of dry cadaveric hip bones and one cadaveric specimen with soft tissue were used to prove the robustness of the developed fiducial detection algorithm. Computed tomography scans of the cadaveric bones were used to validate the accuracy of the integrated system. The fiducial detection sensitivity was in the same range for both sizes of fiducials. While the detection sensitivity was 97.96% for the large fiducials, it was 97.62% for the small fiducials. The acetabulum and the proximal femur were reconstructed with a mean surface distance error of 1.06 and 1.01 mm, respectively. The results for fiducial detection sensitivity and 3D surface reconstruction demonstrated the capability of the integrated system for 3D hip joint reconstruction from 2D calibrated X-ray radiographs.

3D volumetric intensity reconsturction from 2D x-ray images using partial least squares regression

Reconstruction of shape and intensity from 2D x-ray images has drawn more and more attentions. Previously introduced work suffers from the long computing time due to its iterative optimization characteristics and the requirement of generating digitally reconstructed radiographs within each iteration. In this paper, we propose a novel method which uses a patient-specific 3D surface model reconstructed from 2D x-ray images as a surrogate to get a patient-specific volumetric intensity reconstruction via partial least squares regression. No DRR generation is needed. The method was validated on 20 cadaveric proximal femurs by performing a leave-one-out study. Qualitative and quantitative results demonstrated the efficacy of the present method. Compared to the existing work, the present method has the advantage of much shorter computing time and can be applied to both DXA images as well as conventional x-ray images, which may hold the potentials to be applied to clinical routine task such as total hip arthroplasty (THA).

A New Velocity Map for Byrd Glacier, East Antarctica, from Sequential Aster Satellite Imagery

New ice-velocity measurements are obtained for the main trunk of Byrd Glacier, East Antarctica, using recently acquired Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. The velocities are derived from the application of a cross-correlation technique to sequential images acquired in 2000 and 2001. Images were co-registered and ortho-rectified with the aid of a digital elevation model (DEM) generated from ASTER stereo imagery. This paper outlines the process of DEM generation, image co-registration and correction, and the application of the cross-correlation technique to obtain ice velocities. Comparison of the new velocity map with earlier measurements of velocity from 1978 indicates that the glacier has undergone a substantial deceleration between observations. Portions of the glacier flowing at speeds of similar to 850 m a(-1) in 1978/79 were flowing at similar to 650 m a(-1) in 2000/01. The cause of this change in ice dynamics is not known, but the observation shows that East Antarctic outlet glaciers can undergo substantial changes on relatively short timescales.

Fast Algorithms Using Minimal Data Structures for Common Topological Relationships in Large, Irregularly-spaced Topographic Data Sets

Digital terrain models (DTM) typically contain large numbers of postings, from hundreds of thousands to billions. Many algorithms that run on DTMs require topological knowledge of the postings, such as finding nearest neighbors, finding the posting closest to a chosen location, etc. If the postings are arranged irregu- larly, topological information is costly to compute and to store. This paper offers a practical approach to organizing and searching irregularly-space data sets by presenting a collection of efficient algorithms (O(N),O(lgN)) that compute important topological relationships with only a simple supporting data structure. These relationships include finding the postings within a window, locating the posting nearest a point of interest, finding the neighborhood of postings nearest a point of interest, and ordering the neighborhood counter-clockwise. These algorithms depend only on two sorted arrays of two-element tuples, holding a planimetric coordinate and an integer identification number indicating which posting the coordinate belongs to. There is one array for each planimetric coordinate (eastings and northings). These two arrays cost minimal overhead to create and store but permit the data to remain arranged irregularly.

The glacier inventory of Salzburg, Austria, 2007/2009

A third glacier inventory (GI3) is presented for the province of Salzburg where 173 glaciers are located in the seven mountain ranges: Ankogel (47°4'N, 13°14'E), Glockner, Granatspitz, Sonnblick (Goldberg), Hochkönig, Venediger and Zillertal (47°8'N, 12°7'E). The basis for the new GI3 are orthophotos of 2007 and 2009 and the digital elevation model (DEM) of the southern part of Salzburg. On the basis of former inventories, area- and volume changes have been calculated. The biggest relative loss of glacier area per mountain range was found in the Ankogel range and on Hochkönig as a result of the disrupted structure of their small and thin glaciers. In terms of absolute values, the largest changes took place in the Glockner- and Venediger range with an area loss of -10.1 km**2 and -9.7 km**2 during the period between GI1 (1969) and GI3 (2007/2009), respectively. Volume changes have been calculated for nearly half of the glacier area in Salzburg, where DEMs were available. The Glockner, Granatspitz and Sonnblick mountain ranges showed a volume loss of -0.481 km**3 which corresponds to a mean thickness change of -10.5 m. An extrapolation of these changes to all of the 173 glaciers in Salzburg results in a loss of about 1.04 km**3 between GI1 and GI3 and 0.44 km**3 between GI2 and GI3. Overall annual changes in the province of Salzburg between GI2 and GI3 were higher than between GI1 and GI2 and show likewise changes such as those of Tyrol.

Sediment volume variations inferred from DGPS measurements within Bossons glacier proglacial area (Mont-Blanc massif, France)

This dataset presents Differential Global Positioning System data (DGPS) acquired within the Bossons glacier proglacial area. Bossons glacier is a rapidly retreating glacier and its proglacial area is deglaciated for ~30 years. Bossons stream is one of the outlets of the subglacial drainage system. It starts as a 800 m steep cascade reach, then flows through an area with gentler slope : the Plan des Eaux (PdE). PdE is a 300 m long, 50 m wide proglacial alluvial plain with an increasing channel mobility in the downstream direction but decreasing slope gradient and incision. As it may act a sediment trap, studying periglacial and proglacial erosion processes in the Bossons catchment requires to quantify PdE sediment volume evolution. A several meter-sized block located within Bossons proglacial area was set up as GPS base : its location was measured by one antenna (Topcon Hyper Pro) by performing 600 consecutive measurements throughout one day. A second antenna (Topcon Hyper Pro) was then used to measure XYZ location of points in the proglacial area with a ~2 m grid. Radio communication between the two antennas allowed differential calculations to be automatically carried out on field using the Topcon FC-250 hand controller. This methodology yields 3 cm XY and 1.5 cm Z uncertainties. DGPS data have been acquired through 10 campaigns from 2004 to 2014; campaigns from 2004 to 2008 cover a smaller area than those from 2010 to 2014. Digital Elevation Model (DEM) have been interpolated from DGPS data and difference between two DEMs yields deposited and eroded volume within PdE. Maps of PdE volume variation between two campaigns show that incision mainly occurs in the upper and lower sections where as deposition dominates in the middle section. Deposition, denudation and net rate (deposition rate - denudation rate) are calculated by normalizing volumes by DEM areas. Deposition dominates results with a mean net rate of 29 mm/yr. However, strong inter-annual variability exists and some years are dominated by denudation : -36 mm/yr and -100 mm/yr for 2006 and 2011, respectively. Nonetheless, oldest campaigns (2004 to 2008) were carried out on the lower part part of the alluvial plain and ruling them out to keep only complete DEM (2010 to 2014) yields a mean net rate of ~15 mm/yr. This results is coherent with field observations of both strong deposition (e.g. flood deposits) and strong erosion (e.g. 30 cm incision) evidences. Bossons glacier proglacial area is thus dynamic with year-to-year geormorphological changes but may leans toward increasing its mean elevation through a deposition dominated system.

Probability inundation map for a scenario combining storm surge and sea-level rise

Coastal managers require reliable spatial data on the extent and timing of potential coastal inundation, particularly in a changing climate. Most sea level rise (SLR) vulnerability assessments are undertaken using the easily implemented bathtub approach, where areas adjacent to the sea and below a given elevation are mapped using a deterministic line dividing potentially inundated from dry areas. This method only requires elevation data usually in the form of a digital elevation model (DEM). However, inherent errors in the DEM and spatial analysis of the bathtub model propagate into the inundation mapping. The aim of this study was to assess the impacts of spatially variable and spatially correlated elevation errors in high-spatial resolution DEMs for mapping coastal inundation. Elevation errors were best modelled using regression-kriging. This geostatistical model takes the spatial correlation in elevation errors into account, which has a significant impact on analyses that include spatial interactions, such as inundation modelling. The spatial variability of elevation errors was partially explained by land cover and terrain variables. Elevation errors were simulated using sequential Gaussian simulation, a Monte Carlo probabilistic approach. 1,000 error simulations were added to the original DEM and reclassified using a hydrologically correct bathtub method. The probability of inundation to a scenario combining a 1 in 100 year storm event over a 1 m SLR was calculated by counting the proportion of times from the 1,000 simulations that a location was inundated. This probabilistic approach can be used in a risk-aversive decision making process by planning for scenarios with different probabilities of occurrence. For example, results showed that when considering a 1% probability exceedance, the inundated area was approximately 11% larger than mapped using the deterministic bathtub approach. The probabilistic approach provides visually intuitive maps that convey uncertainties inherent to spatial data and analysis.

The International Bathymetric Chart of the Southern Ocean (IBCSO) Version 1.0

The International Bathymetric Chart of the Southern Ocean (IBCSO) Version 1.0 is a new digital bathymetric model (DBM) portraying the seafloor of the circum-Antarctic waters south of 60° S. IBCSO is a regional mapping project of the General Bathymetric Chart of the Oceans (GEBCO). IBCSO Version 1.0 DBM has been compiled from all available bathymetric data collectively gathered by more than 30 institutions from 15 countries. These data include multibeam and single beam echo soundings, digitized depths from nautical charts, regional bathymetric gridded compilations, and predicted bathymetry. Specific gridding techniques were applied to compile the DBM from the bathymetric data of different origin, spatial distribution, resolution, and quality. The IBCSO Version 1.0 DBM has a resolution of 500 x 500 m, based on a polar stereographic projection, and is publicly available together with a digital chart for printing from the project website (http://www.ibcso.org) and from the two data sets shown at the bottom of this page.

Plant Functional Type (PFT) map over Siberia derived from GlobCover 2005 dataset, link to dataset in NetCDF format

High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover data sets are required to both document these changes as well as to provide land-surface information for benchmarking and initializing Earth system models. Earth system models also require specific land cover classification systems based on plant functional types (PFTs), rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover data sets against one another and with auxiliary data to identify key uncertainties that contribute to variability in PFT classifications that would introduce errors in Earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, vegetation continuous fields (MODIS VCFs) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover data set, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFT maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT data set, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to represent the water and carbon cycles in northern latitudes better. Updated land cover data sets are critical for improving and maintaining the relevance of Earth system models for assessing climate and human impacts on biogeochemistry and biophysics.

Satellite image map of King George Island, Antarctica

A new topographic database for King George Island, one of the most visited areas in Antarctica, is presented. Data from differential GPS surveys, gained during the summers 1997/98 and 1999/2000, were combined with up to date coastlines from a SPOT satellite image mosaic, and topographic information from maps as well as from the Antarctic Digital Database. A digital terrain model (DTM) was generated using ARC/INFO GIS. From contour lines derived from the DTM and the satellite image mosaic a satellite image map was assembled. Extensive information on data accuracy, the database as well as on the criteria applied to select place names is given in the multilingual map. A lack of accurate topographic information in the eastern part of the island was identified. It was concluded that additional topographic surveying or radar interferometry should be conducted to improve the data quality in this area. In three case studies, the potential applications of the improved topographic database are demonstrated. The first two examples comprise the verification of glacier velocities and the study of glacier retreat from the various input data-sets as well as the use of the DTM for climatological modelling. The last case study focuses on the use of the new digital database as a basic GIS (Geographic Information System) layer for environmental monitoring and management on King George Island.