Science.gov

Sample records for elevation models lahar

  1. Lahar Hazard Modeling at Tungurahua Volcano, Ecuador

    NASA Astrophysics Data System (ADS)

    Sorensen, O. E.; Rose, W. I.; Jaya, D.

    2003-04-01

    lahar-hazard-zones using a digital elevation model (DEM), was used to construct a hazard map for the volcano. The 10 meter resolution DEM was constructed for Tungurahua Volcano using scanned topographic lines obtained from the GIS Department at the Escuela Politécnica Nacional, Quito, Ecuador. The steep topographic gradients and rapid downcutting of most rivers draining the edifice prevents the deposition of lahars on the lower flanks of Tungurahua. Modeling confirms the high degree of flow channelization in the deep Tungurahua canyons. Inundation zones observed and shown by LAHARZ at Baños yield identification of safe zones within the city which would provide safety from even the largest magnitude lahar expected.

  2. Hydraulic modeling for lahar hazards at cascades volcanoes

    USGS Publications Warehouse

    Costa, J.E.

    1997-01-01

    The National Weather Service flood routing model DAMBRK is able to closely replicate field-documented stages of historic and prehistoric lahars from Mt. Rainier, Washington, and Mt. Hood, Oregon. Modeled time-of-travel of flow waves are generally consistent with documented lahar travel-times from other volcanoes around the world. The model adequately replicates a range of lahars and debris flows, including the 230 million km3 Electron lahar from Mt. Rainier, as well as a 10 m3 debris flow generated in a large outdoor experimental flume. The model is used to simulate a hypothetical lahar with a volume of 50 million m3 down the East Fork Hood River from Mt. Hood, Oregon. Although a flow such as this is thought to be possible in the Hood River valley, no field evidence exists on which to base a hazards assessment. DAMBRK seems likely to be usable in many volcanic settings to estimate discharge, velocity, and inundation areas of lahars when input hydrographs and energy-loss coefficients can be reasonably estimated.

  3. Modeling of the Guagua Pichincha volcano (Ecuador) lahars

    NASA Astrophysics Data System (ADS)

    Canuti, Paolo; Casagli, Nicola; Catani, Filippo; Falorni, Giacomo

    Lahars, here defined as debris flows of volcanic origin, are rapid mass movements that pose a serious threat to cities located in the vicinity of many volcanoes. Quito, capital city of Ecuador and placed at the foot of the Pichincha volcano complex, is exposed to serious inundation hazard as part of the city is built on numerous deposits of large lahars that have occurred in the last 10,000 years. The objective of this paper is to model the potential lahars of the Pichincha volcano to predict inundation areas within the city of Quito. For this purpose two models that apply different approaches were utilized and their results were compared. The programs used were LAHARZ, a semi-empirical model conceived by the United States Geological Survey (USGS), and FLO-2D, a hydraulic model distributed by FLO Software Inc. LAHARZ is designed as a rapid, objective and reproducible automated method for mapping areas of potential lahar inundation (Proc. First Int. Conf. on Debris Flow Hazards Mitigation, San Francisco, USA, ASCE, 1998, p. 176). FLO-2D is a two-dimensional flood routing model for simulating overland flow on complex surfaces such as floodplains, alluvial fans or urbanized areas (FLO-2D Users manual, version 99.2). Both models run within geographical information systems (GIS). Fieldwork was focused on collecting all available information involved in lahar modeling. A total of 49 channel cross-sections were measured along the two main streams and stratigraphic investigations were carried out on the fans to estimate the volume of previous events. A global positioning system was utilized to determine the coordinates of each cross-section. Further data collection concerned topography, rainfall characteristics and ashfall thicknesses. All fieldwork was carried out in cooperation with the Instituto Geofisico of the Escuela Politecnica Nacional. Modeling in a GIS environment greatly aided the exportation of results for the creation of thematic maps and facilitated model

  4. Risk modeling, assessment, and management of lahar flow threat.

    PubMed

    Leung, M F; Santos, J R; Haimes, Y Y

    2003-12-01

    The 1991 eruption of Mount Pinatubo in the Philippines is considered one of the most violent and destructive volcanic activities in the 20th century. Lahar is the Indonesian term for volcanic ash, and lahar flows resulting from the massive amount of volcanic materials deposited on the mountain's slope posed continued post-eruption threats to the surrounding areas, destroying lives, homes, agricultural products, and infrastructures. Risks of lahar flows were identified immediately after the eruption, with scientific data provided by the Philippine Institute of Volcanology, the U.S. Geological Survey, and other research institutions. However, competing political, economic, and social agendas subordinated the importance of scientific information to policy making. Using systemic risk analysis and management, this article addresses the issues of multiple objectives and the effective integration of scientific techniques into the decision-making process. It provides a modeling framework for identifying, prioritizing, and evaluating policies for managing risk. The major considerations are: (1) applying a holistic approach to risk analysis through hierarchical holographic modeling, (2) applying statistical methods to gain insight into the problem of uncertainty in risk assessment, (3) using multiobjective trade-off analysis to address the issue of multiple decisionmakers and stakeholders in the decision-making process, (4) using the conditional expected value of extreme events to complement and supplement the expected value in quantifying risk, and (5) assessing the impacts of multistage decisions. Numerical examples based on ex post data are formulated to illustrate applications to various problems. The resulting framework from this study can serve as a general baseline model for assessing and managing risks of natural disasters, which the Philippines' lead agency-the National Disaster Coordinating Council (NDCC)-and other related organizations can use for their decision

  5. Challenges of modeling current very large lahars at Nevado del Huila Volcano, Colombia

    NASA Astrophysics Data System (ADS)

    Worni, Raphael; Huggel, Christian; Stoffel, M.; Pulgarín, B.

    2012-03-01

    Nevado del Huila, a glacier-covered volcano in the South of Colombia's Cordillera Central, had not experienced any historical eruptions before 2007. In 2007 and 2008, the volcano erupted with phreatic and phreatomagmatic events which produced lahars with flow volumes of up to about 300 million m3 causing severe damage to infrastructure and loss of lives. The magnitude of these lahars and the prevailing potential for similar or even larger events, poses significant hazards to local people and makes appropriate modeling a real challenge. In this study, we analyze the recent lahars to better understand the main processes and then model possible scenarios for future events. We used lahar inundation depths, travel duration, and flow deposits to constrain the dimensions of the 2007 event and applied LAHARZ and FLO-2D for lahar modeling. Measured hydrographs, geophone seismic sensor data and calculated peak discharges served as input data for the reconstruction of flow hydrographs and for calibration of the models. For model validation, results were compared with field data collected along the Páez and Simbola Rivers. Based on the results of the 2007 lahar simulation, we modeled lahar scenarios with volumes between 300 million and 1 billion m3. The approach presented here represents a feasible solution for modeling high-magnitude flows like lahars and allows an assessment of potential future events and related consequences for population centers downstream of Nevado del Huila.

  6. The use of FLO2D numerical code in lahar hazard evaluation at Popocatépetl volcano: a 2001 lahar scenario

    NASA Astrophysics Data System (ADS)

    Caballero, L.; Capra, L.

    2014-12-01

    Lahar modeling represents an excellent tool for designing hazard maps. It allows the definition of potential inundation zones for different lahar magnitude scenarios and sediment concentrations. Here, we present the results obtained for the 2001 syneruptive lahar at Popocatépetl volcano, based on simulations performed with FLO2D software. An accurate delineation of this event is needed, since it is one of the possible scenarios considered if magmatic activity increases its magnitude. One of the main issues for lahar simulation using FLO2D is the calibration of the input hydrograph and rheological flow properties. Here, we verified that geophone data can be properly calibrated by means of peak discharge calculations obtained by the superelevation method. Digital elevation model resolution also resulted as an important factor in defining the reliability of the simulated flows. Simulation results clearly show the influence of sediment concentrations and rheological properties on lahar depth and distribution. Modifying rheological properties during lahar simulation strongly affects lahar distribution. More viscous lahars have a more restricted aerial distribution and thicker depths, and resulting velocities are noticeably smaller. FLO2D proved to be a very successful tool for delimitating lahar inundation zones as well as generating different lahar scenarios not only related to lahar volume or magnitude, but also taking into account different sediment concentrations and rheologies widely documented as influencing lahar-prone areas.

  7. Modelling 2001 lahars at Popocatépetl volcano using FLO2D numerical code

    NASA Astrophysics Data System (ADS)

    Caballero, L.; Capra, L.

    2013-12-01

    Popocatépetl volcano is located on the central part of the Transmexican Volcanic Belt. It is one of the most active volcanoes in Mexico and endanger more than 25 million people that lives in its surroundings. In the last months, the renewal of its volcanic activity put into alert scientific community. One of the possible scenarios is the 2001 explosive activity, which was characterized by a 8 km eruptive column and the subsequent formation of pumice flows up to 4 km from the crater. Lahars were generated few hours after, remobilizing the new deposits towards NE flank of the volcano, along Huiloac Gorge, almost reaching Santiago Xalitzintla town (Capra et al., 2004). The occurrence of a similar scenario makes very important to reproduce this event to delimitate accurately lahar hazard zones. In this work, 2001 lahar deposit is modeled using FLO2D numerical code. Geophone data is used to reconstruct initial hydrograph and sediment concentration. Sensitivity study of most important parameters used by this code like Manning, and α and β coefficients was conducted in order to achieve a good simulation. Results obtained were compared with field data and demonstrated a good agreement in thickness and flow distribution. A comparison with previously published data with laharZ program (Muñoz-Salinas, 2009) is also made. Additionally, lahars with fluctuating sediment concentrations but with similar volume are simulated to observe the influence of the rheological behavior on lahar distribution.

  8. Soil development along elevational transects on granite, andesitic lahar and basalt in the western Sierra Nevada, California

    NASA Astrophysics Data System (ADS)

    Dahlgren, R. A.; Rasmussen, C.; Southard, R. J.

    2009-04-01

    Soil development along three elevational transects, consisting of granite, andesitic lahar and basalt, were investigated on the western slopes of the Sierra Nevada, California to assess the effects of climate on soil properties and processes. The transects, each consisting of four to seven soils, spanned elevations between 150 to 2900 m with mean annual temperatures (3-17 C) decreasing and precipitation (33-150 cm) increasing with increasing elevation. All sites were characterized by a Mediterranean climate with warm to hot, dry summers and cool to cold, wet winters. Vegetation progressed from oak woodland/annual grasslands at low elevations to mixed conifer forest at mid elevations and subalpine mixed conifer forest at high elevations. Soil pH and base saturation decreased with increasing elevation with the largest decrease found on granite. Solum carbon pools ranged from 2 to 25 kg m-2 with the highest contents found in soils formed on andesitic lahar and in mid-elevation soils corresponding to the highest ecosystem net primary productivity. The degree of weathering and mineral assemblages exhibited a strong threshold change at the elevation of the permanent winter snowline (1200-1500 m). Measures of chemical weathering (e.g., clay and Fe oxide production) increased in a near-linear fashion to the winter snowline where they abruptly decreased by about 10-fold. The clay mineralogical assemblage in the rain-dominated weathering zone was dominated by kaolin minerals and was remarkably similar among all parent materials. Within the snow-dominated weathering zone, clay mineralogy was dominated by allophanic materials (allophone/imogolite) on the andesite and basalt compared to hydroxy-Al interlayered 2:1 layer silicates and gibbsite on the granite. Clay translocation resulting in the formation of argillic horizons was only found in the rain-dominated zone. With increasing elevation, soil development followed the order: Alfisols → Ultisols → Inceptisols (granite

  9. Smoothed particle hydrodynamic modeling of volcanic debris flows: Application to Huiloac Gorge lahars (Popocatépetl volcano, Mexico)

    NASA Astrophysics Data System (ADS)

    Haddad, Bouchra; Palacios, David; Pastor, Manuel; Zamorano, José Juan

    2016-09-01

    Lahars are among the most catastrophic volcanic processes, and the ability to model them is central to mitigating their effects. Several lahars recently generated by the Popocatépetl volcano (Mexico) moved downstream through the Huiloac Gorge towards the village of Santiago Xalitzintla. The most dangerous was the 2001 lahar, in which the destructive power of the debris flow was maintained throughout the extent of the flow. Identifying the zone of hazard can be based either on numerical or empirical models, but a calibration and validation process is required to ensure hazard map quality. The Geoflow-SPH depth integrated numerical model used in this study to reproduce the 2001 lahar was derived from the velocity-pressure version of the Biot-Zienkiewicz model, and was discretized using the smoothed particle hydrodynamics (SPH) method. The results of the calibrated SPH model were validated by comparing the simulated deposit depth with the field depth measured at 16 cross sections distributed strategically along the gorge channel. Moreover, the dependency of the results on topographic mesh resolution, initial lahar mass shape and dimensions is also investigated. The results indicate that to accurately reproduce the 2001 lahar flow dynamics the channel topography needed to be discretized using a mesh having a minimum 5 m resolution, and an initial lahar mass shape that adopted the source area morphology. Field validation of the calibrated model showed that there was a satisfactory relationship between the simulated and field depths, the error being less than 20% for 11 of the 16 cross sections. This study demonstrates that the Geoflow-SPH model was able to accurately reproduce the lahar path and the extent of the flow, but also reproduced other parameters including flow velocity and deposit depth.

  10. Secondary lahar hazard assessment for Villa la Angostura, Argentina, using Two-Phase-Titan modelling code during 2011 Cordón Caulle eruption

    NASA Astrophysics Data System (ADS)

    Córdoba, G.; Villarosa, G.; Sheridan, M. F.; Viramonte, J. G.; Beigt, D.; Salmuni, G.

    2015-04-01

    This paper presents the results of lahar modelling in the town of Villa La Angostura (Neuquén-Argentina) based on the Two-Phase-Titan modelling computer code. The purpose of this exercise is to provide decision makers with a useful tool to assess lahar hazard during the 2011 Puyehue-Cordón Caulle Volcanic Complex eruption. The possible occurrence of lahars mobilized from recent ash falls that could reach the city was analysed. The performance of the Two-Phase-Titan model using 15 m resolution digital elevation models (DEMs) developed from optical satellite images and from radar satellite images was evaluated. The output of these modellings showed inconsistencies that, based on field observations, were attributed to bad adjustment of the DEMs to real topography. Further testing of results using more accurate radar-based 10 m DEM, provided more realistic predictions. This procedure allowed us to simulate the path of flows from Florencia, Las Piedritas and Colorado creeks, which are the most hazardous streams for debris flows in Villa La Angostura. The output of the modelling is a valuable tool for city planning and risk management especially considering the glacial geomorphic features of the region, the strong urban development growth and the land occupation that has occurred in the last decade in Villa La Angostura and its surroundings.

  11. Simulating the Osceola Mudflow Lahar Event in the Pacific Northwest using a GPU Based 2-Dimensional Hydraulic Model

    NASA Astrophysics Data System (ADS)

    Katz, B. G.; Eppert, S.; Lohmann, D.; Li, S.; Goteti, G.; Kaheil, Y. H.

    2011-12-01

    At 4,400 meters, Mount Rainer has been the point of origin for several major lahar events. The largest event, termed the "Osceola Mudflow," occurred 5,500 years ago and covered an area of approximately 550km2 with a total volume of deposited material from 2 to 4km3. Particularly deadly, large lahars are estimated to have maximum flow velocities in of 100km/h with a density often described as "Flowing Concrete." While rare, these events typically cause total destruction within a lahar inundation zone. It is estimated that approximately 150,000 people live on top of previous deposits left by lahars which can be triggered by anything from earthquakes to glacial and chemical erosion of volcanic bedrock over time to liquefaction caused by extreme rainfall events. A novel methodology utilizing a 2 dimensional hydraulic model has been implemented allowing for high resolution (30m) lahar inundation maps to be generated. The utility of this model above or in addition to other methodologies such as that of Iverson (1998), lies in its portability to other lahar zones as well as its ability to model any total volume specified by the user. The process for generating lahar flood plains requires few inputs including: a Digital Terrain Map of any resolution (DTM), a mask defining the locations for lahar genesis, a raster of friction coefficients, and a time series depicting uniform material accumulation over the genesis mask which is allowed to flow down-slope. Finally, a significant improvement in speed has been made for solving the two dimensional model by utilizing the latest in graphics processing unit (GPU) technology which has resulted in a greater than 200 times speed up in model run time over previous CPU-based methods. The model runs for the Osceola Mudflow compare favorably with USGS derived inundation regions as derived using field measurements and GIS based approaches such as the LAHARZ program suit. Overall gradation of low to high risk match well, however the new

  12. Lahar flow simulations using LAHARZ program: Application for the Popocatépetl volcano, Mexico

    NASA Astrophysics Data System (ADS)

    Muñoz-Salinas, E.; Castillo-Rodríguez, M.; Manea, V.; Manea, M.; Palacios, D.

    2009-05-01

    Lahars are a type of catastrophic event that is generated in volcanic environments. Due to the recent technological advances computer simulated models have become widely used by the scientific community to predict the distribution of catastrophic events, such as floods, debris avalanches, earth landslides or lahars, and to mitigate possible damage to populations located in a risk area. Our study risk area is located in the Popocatépetl volcano, a large and active stratovolcano in the Trans Mexican Volcanic Belt. Here, lahars have been occurring during the last century, and some of them have seriously affected nearby towns. Since one of the main input parameters into a lahar simulation is the Digital Elevation Model (DEM), this paper focuses on studying lahar simulations on Popocatépetl using DEMs obtained from various sources with the aid of a popular package known as LAHARZ. The first simulations were performed using a 1:50,000 topographic map with 20-m contours interval. The second set of tests used an enhanced DEM with profiles measured during fieldwork in the lahar channel. In the last set of tests, we added lateral artificial barriers to the channel profiles used in the previous enhanced DEM. Our results show that the use of high resolution and detailed DEMs does not necessarily guarantee a realistic lahar simulation as was previously thought, and a detailed explanation as to why the LAHARZ code computation creates ragged edges is offered. Our results also offer useful information on the level of detail necessary for an input DEM into LAHARZ in order to obtain lahar simulations closer to real lahars. The minimum resolution for the DEM used in simulations must be less than the width of the narrowest cross-section of the channel to realistically reproduce the lahars. For DEMs with lower resolution, sometimes the modeled lahars are shorter in length, and therefore underestimate real hazards to populations. Also, the results show that lateral barriers are

  13. The use of FLO2D numerical code in lahar hazard evaluation at Popocatépetl volcano: a 2001-lahar scenario

    NASA Astrophysics Data System (ADS)

    Caballero, L.; Capra, L.

    2014-07-01

    Lahar modelling represents an excellent tool to design hazard maps. It allows the definition of potential inundation zones for different lahar magnitude scenarios and sediment concentrations. Here we present the results obtained for the 2001 syneruptive lahar at Popocatépetl volcano, based on simulations performed with FLO2D software. An accurate delineation of this event is needed since it is one of the possible scenarios considered during a volcanic crisis. One of the main issues for lahar simulation using FLO2D is the calibration of the input hydrograph and rheologic flow properties. Here we verified that geophone data can be properly calibrated by means of peak discharge calculations obtained by superelevation method. Simulation results clearly show the influence of concentration and rheologic properties on lahar depth and distribution. Modifying rheologic properties during lahar simulation strongly affect lahar distribution. More viscous lahars have a more restricted aerial distribution, thicker depths, and resulting velocities are noticeable smaller. FLO2D proved to be a very successful tool to delimitate lahar inundation zones as well as to generate different lahar scenarios not only related to lahar volume or magnitude but also to take into account different sediment concentrations and rheologies widely documented to influence lahar prone areas.

  14. Secondary lahar hazard assessment for Villa la Angostura, Argentina, using Two-Phase-Titan modelling code during 2011 Cordón Caulle eruption

    NASA Astrophysics Data System (ADS)

    Córdoba, G.; Villarosa, G.; Sheridan, M. F.; Viramonte, J. G.; Beigt, D.; Salmuni, G.

    2014-10-01

    This paper shows the results of secondary lahar modelling in Villa La Angostura town (Neuquén-Argentina) based on the Two-Phase-Titan modelling computer code, which aimed to provide decision makers a useful tool to assess lahar hazard during the 2011 Puyehue-Cordón Caulle Volcanic Complex eruption. Possible occurrence of secondary lahars that could reach the city was analized. The performance of the Two-Phase-Titan model using 15 m resolution DEMs developed from optical satellite images and from radar satellite images was evaluated. The output of these modellings showed inconsistencies that, based on field observations, were attributed to bad adjustment of DEMs to real topography. Further testing of results using more accurate radar based 10 m DEM, proved more realistic predictions. The procedure allowed to simulate the path of flows from Florencia, Las Piedritas and Colorado creeks, which are the most influencing streams in Villa La Angostura. The output of the modelling is a valuable tool for city planning and risk management especially considering the glacial geomorphology features of the region, the strong urban development growth and the land occupation tendencies observed in last decade in Villa La Angostura and its surroundings.

  15. Digital Elevation Models

    USGS Publications Warehouse

    U.S. Geological Survey

    1993-01-01

    The Earth Science Information Center (ESIC) distributes digital cartographic/geographic data files produced by the U.S. Geological Survey (USGS) as part of the National Mapping Program. Digital cartographic data files may be grouped into four basic types. The first of these, called a Digital Line Graph (DLG), is the line map information in digital form. These data files include information on base data categories, such as transportation, hypsography, hydrography, and boundaries. The second type, called a Digital Elevation Model (DEM), consists of a sampled array of elevations for a number of ground positions at regularly spaced intervals. The third type is Land Use and Land Cover digital data which provides information on nine major classes of land use such as urban, agricultural, or forest as well as associated map data such as political units and Federal land ownership. The fourth type, the Geographic Names Information System, provides primary information for all known places, features, and areas in the United States identified by a proper name.

  16. Lahars of Mount Pinatubo, Philippines

    USGS Publications Warehouse

    Newhall, Christopher G.; Stauffer, Peter H.; Hendley, James W., II

    1997-01-01

    On June 15, 1991, Mount Pinatubo in the Philippines exploded in the second largest volcanic eruption on Earth this century. This eruption deposited more than 1 cubic mile (5 cubic kilometers) of volcanic ash and rock fragments on the volcano's slopes. Within hours, heavy rains began to wash this material down into the surrounding lowlands in giant, fast-moving mudflows called lahars. In the next four rainy seasons, lahars carried about half of the deposits off the volcano, causing even more destruction in the lowlands than the eruption itself.

  17. Lahar hazard zones for eruption-generated lahars in the Lassen Volcanic Center, California

    USGS Publications Warehouse

    Robinson, Joel E.; Clynne, Michael A.

    2012-01-01

    Lahar deposits are found in drainages that head on or near Lassen Peak in northern California, demonstrating that these valleys are susceptible to future lahars. In general, lahars are uncommon in the Lassen region. Lassen Peak's lack of large perennial snowfields and glaciers limits its potential for lahar development, with the winter snowpack being the largest source of water for lahar generation. The most extensive lahar deposits are related to the May 1915 eruption of Lassen Peak, and evidence for pre-1915 lahars is sparse and spatially limited. The May 1915 eruption of Lassen Peak was a small-volume eruption that generated a snow and hot-rock avalanche, a pyroclastic flow, and two large and four smaller lahars. The two large lahars were generated on May 19 and 22 and inundated sections of Lost and Hat Creeks. We use 80 years of snow depth measurements from Lassen Peak to calculate average and maximum liquid water depths, 2.02 meters (m) and 3.90 m respectively, for the month of May as estimates of the 1915 lahars. These depths are multiplied by the areal extents of the eruptive deposits to calculate a water volume range, 7.05-13.6x106 cubic meters (m3). We assume the lahars were a 50/50 mix of water and sediment and double the water volumes to provide an estimate of the 1915 lahars, 13.2-19.8x106 m3. We use a representative volume of 15x106 m3 in the software program LAHARZ to calculate cross-sectional and planimetric areas for the 1915 lahars. The resultant lahar inundation zone reasonably portrays both of the May 1915 lahars. We use this same technique to calculate the potential for future lahars in basins that head on or near Lassen Peak. LAHARZ assumes that the total lahar volume does not change after leaving the potential energy, H/L, cone (the height of the edifice, H, down to the approximate break in slope at its base, L); therefore, all water available to initiate a lahar is contained inside this cone. Because snow is the primary source of water for

  18. A New Two-phase Flow Model Applied to the 2007 Crater Lake Break-out Lahar, Mt. Ruapehu, New Zealand

    NASA Astrophysics Data System (ADS)

    Sheridan, M. F.; Cordoba, G.; Pitman, E.; Cronin, S. J.; Procter, J.

    2010-12-01

    The 2007 Crater Lake break-out lahar, Mt. Ruapehu, New Zealand, is a complex but well-characterized natural debris flow that follows an intricate course over an array of topographic features (see Manville et al., this conference). Detailed digital terrain data (DEM) and accurate flow characterization allow us to test our computational model with an unusually high level of control for such a large natural flood wave. The new two-phase flow code is imbedded within the TITAN2D framework (Patra et al. 2005) that is widely used in hazard assessment for both dry (granular) and wet (debris flow) flows (Murcia et al., 2010). Because TITAN2D is actually valid for dry flows (avalanches) we developed a new two-phase model based on balance laws for mass and momentum for each phase. The granular material is assumed to obey a Coulomb constitutive relation and the fluid is assumed to be inviscid. The Darcy-Weisbach formulation is used to account for bed friction, and a phenomenological drag coefficient mediates the momentum exchange between phases. The resulting system of 6 partial differential equations are depth averaged and correspond to the Savage and Hutter model in the limit of no fluid, and to the typical shallow water solutions (Ortiz, et al., 2005) for pure water. This model is capable of simulating particle volumetric fractions as dilute as 0.001 and as concentrated as 0.55. To confirm the usefulness of the new code for complex flows we used data from four observation stations at Ruapehu located at runout distances of 2 km, 5 km, 7 km and 9 km. The specific flow data that we compare with the model outcomes include: 1) arrival time of the flood front, 2) maximum flood depth, and 3) flow velocity. The computed values for these flow characteristics are all within about ± 10% of the observed figures. References: Manville, V., et al., 2010, Anatomy of a basin break-out flood: The 2007 Crater Lake break-out lahar, Mt. Ruapehu, New Zealand, this conference. Murcia, H

  19. Application of a wide-ranging two-phase Debris Flow Model to the 2007 Crater Lake break-out lahar at Mt. Ruapehu, New Zealand

    NASA Astrophysics Data System (ADS)

    Sheridan, M. F.; Cordoba, G. A.; Pitman, B.; Cronin, S. J.; Procter, J.

    2011-12-01

    landscapes because this code is widely used in hazard assessment (see Murcia et al., 2010). This new model is capable of simulating particle volumetric fractions as dilute as 0.001 and as concentrated as 0.55. The model has been successfully tested on artificial topographic channels as well as on some volcano landscapes. The 2007 Crater Lake break-out lahar, Mt. Ruapehu, New Zealand, is a complex but well-characterized natural debris flow that follows an intricate course through a variety of topographic features. Detailed digital terrain data (DEM) and accurate flow observations allow us to test our computational model with an unusually high level of control for such a large natural flood wave. We compare model and actual data recorded at four Ruapehu observation stations located at runout distances of 2 km, 5 km, 7 km and 9 km. The specific flow data that we compare include: 1) arrival time of the flood front, 2) maximum flood depth, and 3) flow velocity. The computed values for these flow characteristics are all within about ± 10% of the observed figures.

  20. The Mount Rainier Lahar Detection System

    NASA Astrophysics Data System (ADS)

    Lockhart, A. B.; Murray, T. L.

    2003-12-01

    To mitigate the risk of unheralded lahars from Mount Rainier, the U.S. Geological Survey, in cooperation with Pierce County, Washington, installed a lahar-detection system on the Puyallup and Carbon rivers that originate on Mount Rainier's western slopes. The system, installed in 1998, is designed to automatically detect the passage of lahars large enough to potentially affect populated areas downstream (approximate volume threshold 40 million cubic meters), while ignoring small lahars, earthquakes, extreme weather and floods. Along each river valley upstream, arrays of independent lahar-monitoring stations equipped with geophones and short tripwires telemeter data to a pair of redundant computer base stations located in and near Tacoma at existing public safety facilities that are staffed around the clock. Monitored data consist of ground-vibration levels, tripwire status, and transmissions at regular intervals. The base stations automatically evaluate these data to determine if a dangerous lahar is passing through the station array. The detection algorithm requires significant ground vibration to occur at those stations in the array that are above the anticipated level of inundation, while lower level `deadman' stations, inundated by the flow, experience tripwire breakage or are destroyed. Once a base station detects a lahar, it alerts staff who execute a call-down of public-safety officials and schools, initiating evacuation of areas potentially at risk. Because the system's risk-mitigation task imposes high standards of reliability on all components, it has been under test for several years. To date, the system has operated reliably and without false alarms, including during the nearby M6.8 Nisqually Earthquake on February 28, 2001. The system is being turned over to Pierce County, and activated as part of their lahar warning system.

  1. Lahars simulation and field calibration in Popocatépetl Volcano (Mexico)

    NASA Astrophysics Data System (ADS)

    Haddad, B.; Zamorano, J. J.; Pastor, M.; Andrés, N.; Tanarro, L. M.; Palacios, D.

    2012-04-01

    The term "lahar" refers to the process generated in volcanoes by high sediment concentration flows that range from hyperconcentrated to Debris flows. This complex dynamic system represents a threat to people living near volcanoes. In order to delimitate hazardous area, mathematical models should be applied and tested. These models depend strongly on data collected in the field and, an additionally, good DEM is required to produce satisfactory results. Recent Popocatépetl lahars are well documented and as a result, they can be used to assess the accuracy of numerical models. In this work, SPH (Smoothed Particle Hydrodynamics) depth integrated model created by Pastor in 2005 is applied to reproduce Popocatépetl lahars. The mathematical model is derived from the velocity-pressure version of the Biot-Zienkiewicz model and the assumed rheology corresponds to the Bingham model. On the other hand, a systematic collection of field data it's carried out by GFAM group in Popocatépetl volcano and it's included updating channel topography; as well as the factors: i) run out area boundary; ii) estimation of the velocity of the flow and iii) depth distribution of lahar's deposit. All this field data it's used for back analyses and calibration of the rheological parameters. Besides the calibration of rheological parameters, it is also investigated the effect of the topographic mesh resolution. Moreover, flow depth obtained by SPH model is systematically compared with field evidences along the lahar's path. Research funded by CGL2009-7343 project, Government of Spain.

  2. Causes, Dynamics and Impacts of Lahar Mass Flows due to the April 2015 Eruption of Calbuco Volcano, Chile

    NASA Astrophysics Data System (ADS)

    Dussaillant, Alejandro; Russell, Andy; Meier, Claudio; Rivera, Andres; Mella, Mauricio; Garrido, Natalia; Hernandez, Jorge; Napoleoni, Felipe; Gonzalez, Cristian

    2016-04-01

    provides large volumes of sediment to distal portions of fluvial systems radiating from Calbuco, continuing impact on infrastructure and settlements, including secondary lahars due to rain and melt events. The database generated by this study hopes to contribute to further studies into lahars, including its use to test lahar numerical models.

  3. Lahar simulation with SPH and field calibration at the Colima Volcano (Mexico)

    NASA Astrophysics Data System (ADS)

    Calvo, Leticia; Haddad, Bouchra; Capra, Lucia; Palacios, David

    2015-04-01

    As a result of the frequent effusive activity of Volcán de Colima (10° 30'44''N, 103° 37'02'' W), the most active volcano in Mexico, plenty of rain triggered lahars are produced, especially during the rainy season. Along the recent period of activity, particularly from 2010, many of these lahars channelled through the main ravines of the volcano and reach large distances, representing high risk for more than 10,000 people at the surroundings. Modeling of lahars has become an important tool in the assessment of the related hazards, in order to undertake appropriate mitigation actions and reduce the associated risks. Recent lahars at the Colima Volcano are well documented, so they can be used to prove the accuracy of modelling. In this work, we used the SPH (Smoothed Particle Hydrodynamics) method, a depth integrated coupled model created by Pastor in 2005, to replicate the propagation stage of 3 recent Colima lahars occurred on Montegrande ravine in 1992, 2011 and 2012. The studied events include hyperconcentrated, debris and a mixture of the previous flow natures. The inputs used for the SPH simulations were the initial point, volume of each lahar and an adapted morphology of its mass. Field data used to verify the SPH results include the stopping point of the lahar, its path, velocity and height values, as the floodplain area. All this information was a result of fieldwork recognition (cross section profiles of the inner part of the ravine) and free satellite imagery analysis. The best results were obtained using Bingham rheology. The proposed parameters to simulate Colima lahars were 20 Pa of yield strength and 30 Pa.s of viscosity for the 1992 lahar (hyperconcentrated flow), 200 Pa and 50 Pa.s in case of the 2011 debris flow, and finally 20 Pa and 24 Pa.s for the 2012 event, whose nature evolved from debris to an hyperconcentrated flow. In all cases a 1900 kg/m3 density was used. Highly accurate results showed the relevant role played by rheological

  4. Sedimentological processes in lahars: Insights from optically stimulated luminescence analysis

    NASA Astrophysics Data System (ADS)

    Muñoz-Salinas, Esperanza; Bishop, Paul; Zamorano, Jose-Juan; Sanderson, David

    2012-01-01

    A lahar is a sediment-laden flow capable of major destructive impacts on infrastructure and human life. How lahars transport sediment is thus a key issue for understanding lahar sedimentology and behavior, especially in terms of the lahar's hydraulic and rheological properties, which can be substantially altered as the lahar gains or loses material during its travel. In this research, we analyze lahar entrainment processes by evaluating luminescence signals (total photon counts) from lahar sediments using blue luminescence signals (BLSL). A portable OSL reader that analyzes several grams of polymineral and polygrain-size samples was applied. We use data from three lahars in the Tenenepanco and Huiloac gorges on Popocatépetl volcano (Mexico) to elucidate the ways in which lahars may gain, lose and transport sediment during flow. Sediment samples for luminescence analysis were taken through the full thickness of the lahar deposits at eight different sites for two lahars (1997 and 2001) and at one site for a lahar that occurred > 500 yrs ago. Mean luminescence values obtained for the 1997 lahar help to evaluate the relationship between the lahar entrainment processes and the drainage at the different sites. For the 2001 lahar this relationship was also established taking into account the channel slope. The main conclusion is that luminescence signals (total photon counts) can be used to detect entrained material in the flow because the bulked materials modify the final OSL signals. The application of this type of luminescence analysis thus has the potential to provide a better understanding of sediment entrainment in these sediment-laden flows.

  5. 1-Meter Digital Elevation Model specification

    USGS Publications Warehouse

    Arundel, Samantha T.; Archuleta, Christy-Ann M.; Phillips, Lori A.; Roche, Brittany L.; Constance, Eric W.

    2015-01-01

    In January 2015, the U.S. Geological Survey National Geospatial Technical Operations Center began producing the 1-Meter Digital Elevation Model data product. This new product was developed to provide high resolution bare-earth digital elevation models from light detection and ranging (lidar) elevation data and other elevation data collected over the conterminous United States (lower 48 States), Hawaii, and potentially Alaska and the U.S. territories. The 1-Meter Digital Elevation Model consists of hydroflattened, topographic bare-earth raster digital elevation models, with a 1-meter x 1-meter cell size, and is available in 10,000-meter x 10,000-meter square blocks with a 6-meter overlap. This report details the specifications required for the production of the 1-Meter Digital Elevation Model.

  6. Multidisciplinary determination of Lahar Erosion Dynamics at the Colima Volcano (Mexico)

    NASA Astrophysics Data System (ADS)

    Calvo, Leticia; Renschler, Crish; Haddad, Bouchra; Palacios, David

    2015-04-01

    Volcán de Colima (10° 30'44''N, 103° 37'02'' W) is currently the most active volcano in Mexico and the North American plate. Associated to its frequent volcanic activity, renovated in 1998 and later in 2010, secondary processes like lahars, triggered by rain mixing with the loose pyroclastic debris produced, are common. Colima lahars channelled through the main water drainages (ravines), and reach large distances along their path (from 7 to 15 km long) burying farmland and all kind of human infrastructures at the surrounding area. The inner part of the ravines is greatly affected by lahars, especially by the bulking processes, so establish an appropriate method to determine its affection rate seems to be needed. In order to analyze 1-year lahar erosion dynamics inside one of the most active ravine (Montegrande 2011-2012 period), our team proposed a multidisciplinary perspective that combines numerical modeling (ArcGeoWEPP), fieldwork recognition and free satellite imagery, in the assessment of the related hazards. On the one hand, ArcGeoWEPP model allowed simulation of watersheds and hillslope profiles within ravines, taking into account climate parameters, land and vegetation covers. This tool was especially useful in areas where the terrain complexity prevented access. The results of this model were combined with 16 real cross-section topographies observed inside the Montegrande ravine and the floodplain delineation of lahars created from satellite imagery. The total 1-year volume of debris at Montegrande was finally reached, but also the erosive, sedimentary and balanced areas were identified, so as the lahar and its deposit dimensions. 750,000 tons per year were eroded inside the Montegrande ravine during 2011-2012 lahars, 805,000 tons if the hillslopes of the surrounding area were considered, and 580,000 tons were deposited along the path. The flood plain area was 1,100,000 m2. Numerical models combined with field data obtained from different sources seems

  7. Lahar hazard assessment using Titan2D for an alluvial fan with rapidly changing geomorphology: Whangaehu River, Mt. Ruapehu

    NASA Astrophysics Data System (ADS)

    Procter, J. N.; Cronin, S. J.; Fuller, I. C.; Sheridan, M.; Neall, V. E.; Keys, H.

    2010-03-01

    Rapid changes in small areas at the apex of alluvial fans may have devastating consequences by directing downstream flood or lahar impacts into catchments of widely varying population or infrastructure vulnerability. During a series of lahars in 1995 at Mt. Ruapehu, New Zealand, aggradation of the Whangaehu fan apex (draining the eastern edifice) caused the onset of avulsion of flows northward into the highly vulnerable Tongariro catchment. An earth training dike (or bund) was constructed to protect this catchment by retaining flows on the southern side and the normal lahar outlet path to the south. Surveys in 2001, late 2005, and following a major lahar in March 2007 now show net degradation of a channel in the Whangaehu fan apex, bordered by the bund. This indicates a net increase in the channel capacity at this site and shows that the bund remains at its effective design capacity. Past hydrological modelling used for the bund design provided a large range of discharge estimates but lacked precise constraints on the size and nature of lahars from eruption and lake-breakout events. New modelling has been carried out using Titan2D to examine the impacts of a 6 × 10 6 m 3 volume granular flow down this catchment. This simulates either an eruption or a lake breakout-induced lahar with a historically typical volumetric bulking factor of 4. These simulations predict minimum discharges between 1800 and 2100 m 3/s at the bund site. By comparison, the largest 1995 flow at this site was estimated at around 1200 m 3/s. Further, any single modelled flow from the normal outlet channel of Crater Lake could not be induced to overtop the bund because discharge appears to be limited by the narrow upper reaches of the Whangaehu Gorge. Theoretical discharge levels required to overtop the bund are estimated to be > 6800 m 3/s, assuming no aggradation of the channel by the decelerating flow. Maximum potential discharge at the bund site is additionally modified by potential

  8. Early succession on lahars spawned by Mount St. Helens.

    PubMed

    Del Moral, R

    1998-06-01

    The effects of isolation on primary succession are poorly documented. I monitored vegetation recovery on two Mount St. Helens lahars (mud flows) with different degrees of isolation using contiguous plots. Seventeen years after the eruption, species richness was stable, but cover continued to increase. That isolation affects community structure was confirmed in several ways. The dominance hierarchies of the lahars differed sharply. Detrended correspondence analysis on Lahar I showed a trend related to distance from an adjacent woodland, whereas vegetation on Lahar II was relatively homogeneous. Spectra of growth forms and dispersal types also differed. Lahar I was dominated by species with modest dispersal ability, while Lahar II was dominated by species with better dispersal. Variation between plots should decline through time, a prediction confirmed on Lahar II. Lahar I remained heterogeneous despite having developed significantly higher cover. Here, the increasing distance from the forest has prevented plots from becoming more homogeneous. At this stage of early primary succession, neither lahar is converging towards the species composition of adjacent vegetation. This study shows that isolation and differential dispersal ability combine to determine initial vegetation structure. Stochastic effects resulting from dispersal limitations may resist the more deterministic effects of competition that could lead to floristic convergence. PMID:21684966

  9. Nasadem Global Elevation Model: Methods and Progress

    NASA Astrophysics Data System (ADS)

    Crippen, R.; Buckley, S.; Agram, P.; Belz, E.; Gurrola, E.; Hensley, S.; Kobrick, M.; Lavalle, M.; Martin, J.; Neumann, M.; Nguyen, Q.; Rosen, P.; Shimada, J.; Simard, M.; Tung, W.

    2016-06-01

    NASADEM is a near-global elevation model that is being produced primarily by completely reprocessing the Shuttle Radar Topography Mission (SRTM) radar data and then merging it with refined ASTER GDEM elevations. The new and improved SRTM elevations in NASADEM result from better vertical control of each SRTM data swath via reference to ICESat elevations and from SRTM void reductions using advanced interferometric unwrapping algorithms. Remnant voids will be filled primarily by GDEM3, but with reduction of GDEM glitches (mostly related to clouds) and therefore with only minor need for secondary sources of fill.

  10. Magnitude and frequency of lahars and lahar-runout flows in the Toutle-Cowlitz River system

    SciTech Connect

    Scott, K.M.

    1989-01-01

    The recurrence interval of a lahar or lahar-runout flow at least large enough to inundate flood plains 50 kilometers from Mount St. Helens is less than 100 years. Lahars are volcanic debris flows and their deposits; lahar-runout flows are the hyperconcentrated streamflow evolved from distal lahars. The recurrence interval is conditional on eruptive state and is based on the most recent 4,500 years of the volcano's approximately 40,000- to 50,000-year history. The 100-year recurrence interval is within a normal time frame for long-term planning. Therefore engineering works in the Toutle River system should be designed for lahars, as well as floods, of a particular frequency. Unlike a water flood, a lahar that has a flow depth at least 1 meter on flood plains can cause a significant part of the maximum possible damage. Trees are killed, many structures are inundated and made unusable even if they are not crushed by timber floating in the lahar, and agriculture is not feasible for periods of as much as several years. The largest lahar in this history of the watershed was formed by the bulking of the sediment in a flood surge that originated from breaching of a natural dam of ancestral Spirit Lake. The flow had a peak discharge of 300,000 to 300,000 m{sup 3}/s at a distance of 30 to 50 km from the volcano, and was the first of four lake-breakout lahars that occurred during a span of several years near the end of Pine Creek time. This series of lahars is interpreted as an analog of the events that would have happened, without engineering intervention, after the 1980 eruption. In 1980, a debris avalanche catastrophically raised Spirit Lake more than 60 m and created new lakes in blocked tributaries.

  11. Rapid Loss of Andean Alpine Glaciers: A Reflection on Cotopaxi´s Long-Distance Historical Lahars and Future Lahar Scenarios

    NASA Astrophysics Data System (ADS)

    Mothes, P. A.; Hall, M. L.; Samaniego, P.; Francou, B.; Castro, M.; Hidalgo, X.

    2007-05-01

    Andean alpine glaciers are in rapid retreat, as witnessed by actual measurements, comparative imagery and popular memory. Overall glacier losses will diminish future water availability for human consumption as well as for lahar generation, the product of mixing incandescent eruptive materials with glacial ice and snow. The field study and modeling of long-distance historical lahars from Cotopaxi volcano, Ecuador has shown them to be some of the most voluminous and longest reported. Based on back calculations, peak discharges were commonly between 45,000-60,000 m3/sec, velocities reached 70 km/hr, and run outs attained 325 km. The last "super" debris flow was produced at Cotopaxi in 1877. Observations made after the 1877 eruption reported that the glacier had suffered about 10 meters of ice stripped off the top and the incision of deep gullies from melting and erosion by the scoria block-rich pyroclastic flows. Average reductions of 45% and 60%, respectively, of the area and volume of Cotopaxi´s 19 alpine glaciers during the last 30 years have left an ice cap of only 13 km2 and a volume of 0.60 km3. Descriptions by astute 18th and 19th century observers lead us to conclude that Cotopaxi glaciers were much more robust then, surpassing a total area of about 30 km2, a fact which contributed to generating large volume lahars and high discharges, during the waning "Little Ice Age". If an eruption similar to that of 1877 occurs at Cotopaxi in the future, reduced glacier sizes and the glaciers´ preferential distribution upon the cone will likely attenuate volcano-ice interactions and will lower the probability of "super" lahars being produced during eruptive periods. However, in the last 2000 years of eruptive activity, explosive eruptions display a large size span-- from weakly explosive events (VEI= 2) to highly explosive eruptive cycles (VEI= 4-5). Given the uncertainty of the size of the next explosive eruption of Cotopaxi, several scenarios for lahar generation must

  12. US GeoData Digital Elevation Models

    USGS Publications Warehouse

    U.S. Geological Survey

    2000-01-01

    Digital elevation model (DEM) data are arrays of regularly spaced elevation values referenced horizontally either to a Universal Transverse Mercator (UTM) projection or to a geographic coordinate system. The grid cells are spaced at regular intervals along south to north profiles that are ordered from west to east. The U.S. Geological Survey (USGS) produces five primary types of elevation data: 7.5-minute DEM, 30-minute DEM, 1-degree DEM, 7.5-minute Alaska DEM, and 15-minute Alaska DEM.

  13. US GeoData digital elevation models

    USGS Publications Warehouse

    U.S. Geological Survey

    1997-01-01

    Digital elevation model (DEM) data consist of a sampled array of regularly spaced elevation values referenced horizontally either to a Universal Transverse Mercator (UTM) projection or to a geographic coordinate system. The grid cells are spaced at regular intervals along south to north profiles that are ordered from west to east. the U.S> Geological Survey (USGS) produces five primary types of elevation data: 7.5-minute DEM, 30-minute DEM, 1-degree DEM, 7.5-minute Alaska DEM, and 15-minute Alaska DEM.

  14. A comparison of the Landsat image and LAHARZ-simulated lahar inundation hazard zone by the 2010 Merapi eruption

    NASA Astrophysics Data System (ADS)

    Lee, Seul-Ki; Lee, Chang-Wook; Lee, Saro

    2015-06-01

    Located above the Java subduction zone, Merapi Volcano is an active stratovolcano with a volcanic activity cycle of 1-5 years. Most Merapi eruptions are relatively small with volcanic explosivity index (VEI) of 1-3. However, the most recent eruption, which occurred in 2010, was quite violent with a VEI of 4 and 386 people were killed. In this study, lahars and pyroclastic flow zones were detected using optical Landsat images and the lahar and pyroclastic flow zone simulated using the LAHARZ program. To detect areal extents of lahar and pyroclastic flows using Landsat images, supervised classification was performed after atmospheric correction by using a cosine of the solar zenith correction (COST) model. As a result, the extracted dimensions of pyroclastic flows are nearly identical to the Calatrava Volcanic Province (CVP) monthly reports. Then, areas of potential lahar and pyroclastic flow inundation based on flow volume using the LAHARZ program were simulated and mapped. Finally, the detected lahars and pyroclastic flow zones were compared with the simulated potential zones using LAHARZ program and verified. Results showed satisfactory similarity (55.63 %) between the detected and simulated zone. The simulated zones using the LAHARZ program can be used as an essential volcanic hazard map for preventing life and property damages for Merapi Volcano and other hazardous volcanic areas. Also, the LAHARZ program can be used to map volcano hazards in other hazardous volcanic areas.

  15. Laharz_py: GIS tools for automated mapping of lahar inundation hazard zones

    USGS Publications Warehouse

    Schilling, Steve P.

    2014-01-01

    Laharz_py is written in the Python programming language as a suite of tools for use in ArcMap Geographic Information System (GIS). Primarily, Laharz_py is a computational model that uses statistical descriptions of areas inundated by past mass-flow events to forecast areas likely to be inundated by hypothetical future events. The forecasts use physically motivated and statistically calibrated power-law equations that each has a form A = cV2/3, relating mass-flow volume (V) to planimetric or cross-sectional areas (A) inundated by an average flow as it descends a given drainage. Calibration of the equations utilizes logarithmic transformation and linear regression to determine the best-fit values of c. The software uses values of V, an algorithm for idenitifying mass-flow source locations, and digital elevation models of topography to portray forecast hazard zones for lahars, debris flows, or rock avalanches on maps. Laharz_py offers two methods to construct areas of potential inundation for lahars: (1) Selection of a range of plausible V values results in a set of nested hazard zones showing areas likely to be inundated by a range of hypothetical flows; and (2) The user selects a single volume and a confidence interval for the prediction. In either case, Laharz_py calculates the mean expected A and B value from each user-selected value of V. However, for the second case, a single value of V yields two additional results representing the upper and lower values of the confidence interval of prediction. Calculation of these two bounding predictions require the statistically calibrated prediction equations, a user-specified level of confidence, and t-distribution statistics to calculate the standard error of regression, standard error of the mean, and standard error of prediction. The portrayal of results from these two methods on maps compares the range of inundation areas due to prediction uncertainties with uncertainties in selection of V values. The Open-File Report

  16. Lahar hazards at Agua volcano, Guatemala

    USGS Publications Warehouse

    Schilling, S.P.; Vallance, J.W.; Matías, O.; Howell, M.M.

    2001-01-01

    At 3760 m, Agua volcano towers more than 3500 m above the Pacific coastal plain to the south and 2000 m above the Guatemalan highlands to the north. The volcano is within 5 to 10 kilometers (km) of Antigua, Guatemala and several other large towns situated on its northern apron. These towns have a combined population of nearly 100,000. It is within about 20 km of Escuintla (population, ca. 100,000) to the south. Though the volcano has not been active in historical time, or about the last 500 years, it has the potential to produce debris flows (watery flows of mud, rock, and debris—also known as lahars when they occur on a volcano) that could inundate these nearby populated areas.

  17. Lahar hazards at Mombacho Volcano, Nicaragua

    USGS Publications Warehouse

    Vallance, J.W.; Schilling, S.P.; Devoli, G.

    2001-01-01

    Mombacho volcano, at 1,350 meters, is situated on the shores of Lake Nicaragua and about 12 kilometers south of Granada, a city of about 90,000 inhabitants. Many more people live a few kilometers southeast of Granada in 'las Isletas de Granada and the nearby 'Peninsula de Aseses. These areas are formed of deposits of a large debris avalanche (a fast moving avalanche of rock and debris) from Mombacho. Several smaller towns with population, in the range of 5,000 to 12,000 inhabitants are to the northwest and the southwest of Mombacho volcano. Though the volcano has apparently not been active in historical time, or about the last 500 years, it has the potential to produce landslides and debris flows (watery flows of mud, rock, and debris -- also known as lahars when they occur on a volcano) that could inundate these nearby populated areas. -- Vallance, et.al., 2001

  18. Instrumental lahar monitoring at Merapi Volcano, Central Java, Indonesia

    USGS Publications Warehouse

    Lavigne, F.; Thouret, J.-C.; Voight, B.; Young, K.; LaHusen, R.; Marso, J.; Suwa, H.; Sumaryono, A.; Sayudi, D.S.; Dejean, M.

    2000-01-01

    More than 50 volcanic debris flows or lahars were generated around Mt Merapi during the first rainy season following the nuees ardentes of 22 November 1994. The rainfalls that triggered the lahars were analyzed, using such instruments as weather radar and telemetered rain gauges. Lahar dynamics were also monitored, using new non-contact detection instrumentation installed on the slopes of the volcano. These devices include real-time seismic amplitude measurement (RSAM), seismic spectral amplitude measurement (SSAM) and acoustic flow monitoring (AFM) systems. Calibration of the various systems was accomplished by field measurements of flow velocities and discharge, contemporaneously with instrumental monitoring. The 1994–1995 lahars were relatively short events, their duration in the Boyong river commonly ranging between 30 min and 1 h 30 min. The great majority (90%) of the lahars was recognized at Kaliurang village between 13:00 and 17:30 h, due to the predominance of afternoon rainfalls. The observed mean velocity of lahar fronts ranged between 1.1 and 3.4 m/s, whereas the peak velocity of the flows varied from 11 to 15 m/s, under the Gardu Pandang viewpoint location at Kaliurang, to 8–10 m/s at a section 500 m downstream from this site. River slopes vary from 28 to 22 m/km at the two sites. Peak discharges recorded in various events ranged from 33 to 360 m3/s, with the maximum value of peak discharge 360 m3/s, on 20 May 1995. To improve the lahar warning system along Boyong river, some instrumental thresholds were proposed: large and potentially hazardous lahars may be detected by RSAM units exceeding 400, SSAM units exceeding 80 on the highest frequency band, or AFM values greater than 1500 mV on the low-gain, broad-band setting.

  19. Digital Elevation Model Mosaic of Mercury

    NASA Technical Reports Server (NTRS)

    Cook, A. C.; Watters, T. R.; Robinson, M. S.

    2001-01-01

    At CEPS (Center for Earth and Planetary Studies) work has been underway since 2000 to semi-automatically stereo match all Mariner 10 stereo pairs. The resulting matched image coordinates are converted into longitude, latitude, and height points and then combined to form a map projected Digital Elevation Model (DEM) mosaic of the planet's surface. Stereo images from Mariner 10 cover one quarter of the planet's surface, mostly in the southern hemisphere. Additional information is contained in the original extended abstract.

  20. Stochastic Downscaling of Digital Elevation Models

    NASA Astrophysics Data System (ADS)

    Rasera, Luiz Gustavo; Mariethoz, Gregoire; Lane, Stuart N.

    2016-04-01

    High-resolution digital elevation models (HR-DEMs) are extremely important for the understanding of small-scale geomorphic processes in Alpine environments. In the last decade, remote sensing techniques have experienced a major technological evolution, enabling fast and precise acquisition of HR-DEMs. However, sensors designed to measure elevation data still feature different spatial resolution and coverage capabilities. Terrestrial altimetry allows the acquisition of HR-DEMs with centimeter to millimeter-level precision, but only within small spatial extents and often with dead ground problems. Conversely, satellite radiometric sensors are able to gather elevation measurements over large areas but with limited spatial resolution. In the present study, we propose an algorithm to downscale low-resolution satellite-based DEMs using topographic patterns extracted from HR-DEMs derived for example from ground-based and airborne altimetry. The method consists of a multiple-point geostatistical simulation technique able to generate high-resolution elevation data from low-resolution digital elevation models (LR-DEMs). Initially, two collocated DEMs with different spatial resolutions serve as an input to construct a database of topographic patterns, which is also used to infer the statistical relationships between the two scales. High-resolution elevation patterns are then retrieved from the database to downscale a LR-DEM through a stochastic simulation process. The output of the simulations are multiple equally probable DEMs with higher spatial resolution that also depict the large-scale geomorphic structures present in the original LR-DEM. As these multiple models reflect the uncertainty related to the downscaling, they can be employed to quantify the uncertainty of phenomena that are dependent on fine topography, such as catchment hydrological processes. The proposed methodology is illustrated for a case study in the Swiss Alps. A swissALTI3D HR-DEM (with 5 m resolution

  1. Digital Elevation Models Aid the Analysis of Flows at Hrad Vallis, Mars

    NASA Astrophysics Data System (ADS)

    Mouginis-Mark, P. J.; Hamilton, C.; Garbeil, H.

    2015-12-01

    We have identified several landforms in the Hrad Vallis region of Mars (33.0o - 35.5oN, 216o - 218oW), which suggest that this area was covered by an ice sheet concurrent with volcanic eruptions. Using digital elevation models derived from High Resolution Imaging Science (HiRISE) and Context Camera (CTX) data, a reexamination of the area reveals a complex history including flow inflation and topographic control by transient topographic barriers. Among Amazonian-age outflow channels, Hrad Vallis is exceptional as it exhibits good evidence of magma/water interactions. It is inferred to have formed in association with a shallow igneous sill that melted part of the martian cryosphere and/or released water from an extensive aquifer to produce enormous lahar-like mud flows. Exposed ~30 m high dikes, 20 m high eroded mounds, and flow paths that are inconsistent with present-day topographic gradients, lead us to speculate that this area was covered by at least ~40 m of material (eolian deposits or ice) at the time of volcanic dike intrusion and flow emplacement. This material was subsequently removed leaving no clear morphologic signs (e.g., wind streaks, if eolian material; moraines, if ice). We favor the ice model because if this area was once ice-covered, it offers a plausible mode of formation (as pingoes) for some enigmatic 30 m high domes in the vicinity. At least one 120 km long flow from Hrad Vallis was emplaced as a pahoehoe-like flow that was confined by topographic obstacles and subsequently inflated to thickness of ~45 m. Although the direct relationship between this flow and Hrad Vallis remains to be determined, the inflated flow suggests a longer period of eruption/emplacement at a slower effusion rate than was previously believed.

  2. Lahar Infrasound Associated with Villarrica's March 3, 2015 Eruption

    NASA Astrophysics Data System (ADS)

    Johnson, J. B.; Palma, J.

    2015-12-01

    The paroxysmal eruption of Volcan Villarrica on March 3rd, 2015 produced a moderate-sized lahar, which descended more than 20 km within the Rio Correntoso/Turbio drainage. A three-element infrasound array 10 km from the summit, and 4 km from the lahar's closest approach, was used to track the flow's evolution as it progressed downstream. Array processing using cross correlation lag times as well as semblance techniques places important constraints on the lahar's dynamics, including the detection of an early flow pulse that traveled from 3 to 13 km at an average speed of 36 m/s. After the first six minutes of lahar advancement the signal evolves into a relatively stationary infrasonic tremor located ~11.5 km downstream and consistent with a notch in intervening topography. Diminishing tremor amplitude over the course of more than two hours constrain the flow duration and indicates progressively decreasing flow energy and/or confinement of the flow to more distant reaches. This study demonstrates the powerful capabilities of infrasound arrays for lahar study and suggests its potential implementation for hazard monitoring.

  3. Comparative lahar hazard mapping at Volcan Citlaltépetl, Mexico using SRTM, ASTER and DTED-1 digital topographic data

    NASA Astrophysics Data System (ADS)

    Hubbard, Bernard E.; Sheridan, Michael F.; Carrasco-Núñez, Gerardo; Díaz-Castellón, Rodolfo; Rodríguez, Sergio Raúl

    2007-02-01

    In this study, we evaluated and compared the utility of spaceborne SRTM and ASTER DEMs with baseline DTED-1 "bald-earth" topography for mapping lahar inundation hazards from volcan Citlaltépetl, Mexico, a volcano which has had a history of producing debris flows of various extents. In particular, we tested the utility of these topographic datasets for resolving ancient valley-filling deposits exposed around the flanks of the volcano, for determining their magnitude using paleohydrologic methods and for forecasting their inundation limits in the future. We also use the three datasets as inputs to a GIS stream inundation flow model, LAHARZ, and compare the results. In general all three datasets, with spatial resolution of 90 m or better, were capable of resolving debris flow and lahar deposits at least 3 × 10 6 m 3 in volume or larger. Canopy- and slope-related height errors in the ASTER and SRTM DEMs limit their utility for measuring valley-filling cross-sectional area and deriving flow magnitude for the smallest deposits using a cross-sectional area to volume scaling equation. Height errors in the ASTER and SRTM DEMs also causes problems in resolving stream valley hydrography which controls lahar flow paths and stream valley morphology which controls lahar filling capacity. However, both of the two spaceborne DEM datasets are better than DTED-1 at resolving fine details in stream hydrography and erosional morphologies of volcaniclastics preserved in the valleys around the more humid, eastern flanks of the volcanic range. The results of LAHARZ flow inundation modeling using all three DEMs as inputs are remarkably similar and co-validate one another. For example, at Citlaltépetl all lahar simulations show that the city of Orizaba is the most vulnerable to flows similar in magnitude to, or larger than, one that occurred in 1920. Many of the other cities and towns illustrated are built higher up on terrace deposits of older debris flows, and are safe from all but

  4. A global digital elevation model - GTOP030

    USGS Publications Warehouse

    1999-01-01

    GTOP030, the U.S. Geological Survey's (USGS) digital elevation model (DEM) of the Earth, provides the flrst global coverage of moderate resolution elevation data.  The original GTOP30 data set, which was developed over a 3-year period through a collaborative effort led by the USGS, was completed in 1996 at the USGS EROS Data Center in Sioux Falls, South Dakota.  The collaboration involved contributions of staffing, funding, or source data from cooperators including the National Aeronautics and Space Administration (NASA), the United Nations Environment Programme Global Resource Information Database (UNEP/GRID), the U.S. Agency for International Development (USAID), the Instituto Nacional de Estadistica Geografia e Informatica (INEGI) of Mexico, the Geographical Survey Institute (GSI) of Japan, Manaaki Whenua Landcare Research of New Zealand, and the Scientific Committee on Antarctic Research (SCAR). In 1999, work was begun on an update to the GTOP030 data set. Additional data sources are being incorporated into GTOP030 with an enhanced and improved data set planned for release in 2000.

  5. Simulating Lahars Using A Rotating Drum

    NASA Astrophysics Data System (ADS)

    Neather, Adam; Lube, Gert; Jones, Jim; Cronin, Shane

    2014-05-01

    A large (0.5 m in diameter, 0.15 m wide) rotating drum is used to investigate the erosion and deposition mechanics of lahars. To systematically simulate the conditions occurring in natural mass flows our experimental setup differs from the common rotating drum employed in industrial/engineering studies. Natural materials with their typical friction properties are used, as opposed to the frequently employed spherical glass beads; the drum is completely water-proof, so solid/air and solid/liquid mixtures can be investigated; the drum velocity and acceleration can be precisely controlled using a software interface to a micro-controller, allowing for the study of steady, unsteady and intermediate flow regimes. The drum has a toughened glass door, allowing high-resolution, high-speed video recording of the material inside. Vector maps of the velocities involved in the flows are obtained using particle image velocimetry (PIV). The changes in velocity direction and/or magnitude are used to locate the primary internal boundaries between layers of opposite flow direction, as well as secondary interfaces between shear layers. A range of variables can be measured: thickness and number of layers; the curvature of the free surface; frequency of avalanching; position of the centre of mass of the material; and the velocity profiles of the flowing material. Experiments to date have focussed on dry materials, and have had a fill factor of approximately 0.3. Combining these measured variables allows us to derive additional data of interest, such as mass and momentum flux. It is these fluxes that we propose will allow insight into the erosion/deposition mechanics of a lahar. A number of conclusions can be drawn to date. A primary interface separates flowing and passive region (this interface has been identified in previous studies). As well as the primary interface, the flowing layer separates into individual shear layers, with individual erosion/deposition and flow histories. This

  6. Factors controlling erosion/deposition phenomena related to lahars at Volcán de Colima, Mexico

    NASA Astrophysics Data System (ADS)

    Vázquez, Rosario; Capra, Lucia; Coviello, Velio

    2016-08-01

    comparing rainfalls associated with lahars that originated after the last main eruptive episode that occurred in 2004-2005, we observed that higher accumulated rainfall was needed to trigger lahars in the 2013 and 2014 seasons, which points to a progressive stabilization of the volcano slope during a post-eruptive period. These results can be used as a tool to foresee the channel response to future volcanic activity, to improve the input parameters for lahar modeling and to better constrain the hazard zonation at Volcán de Colima.

  7. Alteration, slope-classified alteration, and potential lahar inundation maps of volcanoes for the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Volcano Archive

    USGS Publications Warehouse

    Mars, John C.; Hubbard, Bernard E.; Pieri, David; Linick, Justin

    2015-01-01

    This study was undertaken during 2012–2013 in cooperation with the National Aeronautics and Space Administration (NASA). Since completion of this study, a new lahar modeling program (LAHAR_pz) has been released, which may produce slightly different modeling results from the LAHARZ model used in this study. The maps and data from this study should not be used in place of existing volcano hazard maps published by local authorities. For volcanoes without hazard maps and (or) published lahar-related hazard studies, this work will provide a starting point from which more accurate hazard maps can be produced. This is the first dataset to provide digital maps of altered volcanoes and adjacent watersheds that can be used for assessing volcanic hazards, hydrothermal alteration, and other volcanic processes in future studies.

  8. SHADED RELIEF, HILLSHADE, DIGITAL ELEVATION MODEL (DEM), NEVADA

    EPA Science Inventory

    Shaded relief of the state of Nevada developed from 1-degree US Geological Survey (USGS) Digital Elevation Models (DEMs). DEM is a terminology adopted by the USGS to describe terrain elevation data sets in a digital raster form.

  9. SHADED RELIEF, HILLSHADE, DIGITAL ELEVATION MODEL (DEM), ARIZONA

    EPA Science Inventory

    Shaded relief of the state of Arizona developed from 1-degree US Geological Survey (USGS) Digital Elevation Models (DEMs). DEM is a terminology adopted by the USGS to describe terrain elevation data sets in a digital raster form.

  10. Probabilistic analysis of rain-triggered lahar initiation at Tungurahua volcano

    NASA Astrophysics Data System (ADS)

    Jones, Robbie; Manville, Vern; Andrade, Daniel

    2015-08-01

    Semi-continuous production of pyroclastic material by intermittent strombolian, vulcanian and sub-plinian eruptions at Volcán Tungurahua, Ecuador has created a persistent rain-triggered lahar hazard during the 1999-present eruptive episode. Lahars threaten the city of Baños, which lies approximately 8 km from the crater, as well as other villages and vital infrastructure situated in close proximity to the dense radial drainage network of the volcano. This study analyses the initiation of rain-triggered lahars and the influence of antecedent rainfall on this process in two northern instrumented drainages, La Pampa and the Vazcun. Analysis of lahar-triggering rainfall intensity and duration between March 2012 and June 2013 yields a power-law relationship, whilst receiver operating characteristic (ROC) analysis indicates that peak rainfall intensity (10, 30 and 60 min) is the most effective single predictor of lahar occurrence. The probability of a lahar exceeding a pre-defined magnitude increases with peak rainfall intensity. Incorporation of antecedent rainfall (24 h and 3, 5 and 7 days) as a secondary variable significantly impacts lahar probabilities, particularly during moderate-high-intensity rainfall events. The resultant two- and three-dimensional lahar probability matrices are applied to rainfall data between 1st July and 31st December 2013 with the aim of predicting lahar occurrence. Composite lahar indicators comprised from the mean lahar probability estimates of individual matrices are shown to perform this task most effectively. ROC analysis indicates a probability > 80 % that these composite indicators will generate a higher estimated lahar probability for a randomly selected lahar event than a randomly selected non-lahar event. This method provides an average of 24 min of additional warning time compared with the current acoustic flow monitors (AFMs) used for lahar detection, effectively doubling warning times for key downstream infrastructure in the

  11. Estimation of lahar flow velocity on Popocatépetl volcano (Mexico)

    NASA Astrophysics Data System (ADS)

    Muñoz-Salinas, E.; Manea, V. C.; Palacios, D.; Castillo-Rodriguez, M.

    2007-11-01

    During 1997 and 2001 two lahars took place in the Tenenepanco and Huiloac gorges on the northeastern sector of the Popocatépetl volcano in Mexico. These lahars were the result of volcanic activity characterized by a pyroclastic fall and flow. The lahars were triggered by the mass failures produced hours after the eruption. The well-known superelevation technique commonly used in this kind of study was applied to infer the lahar's mean velocity. Estimations show a broad flow velocity range of 1.3-13.8 m/s depending on the following geomorphological parameters: channel depth and slope, peak discharge and distance from the initial source. These results were compared with lahar velocities from Mount St Helens (1980) and Nevado del Ruiz (1985), which were obtained by the same superelevation technique. According to the data, the lahars from Nevado del Ruiz and Popocatépetl were triggered by similar conditions, while at Mount St Helens they were the result of the conversion of a pyroclastic surge cloud. Actually, the Mount St Helens lahars were almost twice as fast as the ones at Nevado del Ruiz and Popocatépetl. This observation shows the importance of the lahars' genesis in that the pyroclastic surge clouds produce faster lahars than do pyroclastic falls. Such results should be taken into account in lahar risk assessment studies for communities located near active volcanoes.

  12. Digital elevation modeling via curvature interpolation for lidar data

    Technology Transfer Automated Retrieval System (TEKTRAN)

    Digital elevation model (DEM) is a three-dimensional (3D) representation of a terrain's surface - for a planet (including Earth), moon, or asteroid - created from point cloud data which measure terrain elevation. Its modeling requires surface reconstruction for the scattered data, which is an ill-p...

  13. Lahar Risk on the NE Flank of Popocatepetl Volcano

    NASA Technical Reports Server (NTRS)

    Delgado, H.; Huesca, E. A. Gonzalez; Abrams, M.

    1994-01-01

    Popocatepetl volcano has an altitude of 5452 meters and is capped by glaciers which represent a volume of less than 0.017 km(sup 3) of ice. These glaciers are distributed on the northwest-north face of the cone, starting at 4900 m.a.s.l. The ablation runoff of the glaciers is channelized to the north through the Barranca Central and Barranca Ventorillo which at the lower altitude join together to form a larger canyon bent to the northeast following a spur made of volcanic rock from the extinct Iztacchuatl volcano. At 3200 m.a.s.l. a sudden change in morphology marks the starting point of lahar deposits. Several of these mudflows have been mapped. The San Nicolas Lahar covers a poorly developed soil where pottery and other cultural remains have been found. There are now several towns in the area.

  14. Integrated research in constitutive modelling at elevated temperatures, part 1

    NASA Technical Reports Server (NTRS)

    Haisler, W. E.; Allen, D. H.

    1986-01-01

    Topics covered include: numerical integration techniques; thermodynamics and internal state variables; experimental lab development; comparison of models at room temperature; comparison of models at elevated temperature; and integrated software development.

  15. Model-based approach for elevator performance estimation

    NASA Astrophysics Data System (ADS)

    Esteban, E.; Salgado, O.; Iturrospe, A.; Isasa, I.

    2016-02-01

    In this paper, a dynamic model for an elevator installation is presented in the state space domain. The model comprises both the mechanical and the electrical subsystems, including the electrical machine and a closed-loop field oriented control. The proposed model is employed for monitoring the condition of the elevator installation. The adopted model-based approach for monitoring employs the Kalman filter as an observer. A Kalman observer estimates the elevator car acceleration, which determines the elevator ride quality, based solely on the machine control signature and the encoder signal. Finally, five elevator key performance indicators are calculated based on the estimated car acceleration. The proposed procedure is experimentally evaluated, by comparing the key performance indicators calculated based on the estimated car acceleration and the values obtained from actual acceleration measurements in a test bench. Finally, the proposed procedure is compared with the sliding mode observer.

  16. Digital Elevation Models of the Lunar Surface

    NASA Astrophysics Data System (ADS)

    Cook, A. C.; Robinson, M. S.

    1999-01-01

    Several digital elevation models (DEMs) have been produced at a scale of 1km/pixel and covering approximately one-fifth of the lunar surface. These were produced mostly by semiautomatically matching the stereo available between Clementine UV/VIS images, although some localized DEMs have been produced by applying this technique to Apollo Metric stereo pairs, or by digitizing an existing Apollo Metric contour map. The DEMS that result from Clementine UV/VIS images, although Of Poorer height accuracy (1300-600 in for a single matched point) than the Clementine laser altimeter point measurements (<+/-100 m), do provide considerably higher spatial resolution (e.g., every kilometer vs. every tens of kilometers) and allow topography in the polar regions to be determined. Nadir-pointing Clementine UV-VIS stereo pairs are automatically stereo matched using a patch-based matcher and fed through A stereo intersection camera model to yield a digital terrain model (DTM) of longitude, latitude, and height points. The DTM for each stereo pair is then replotted and interpolated to form map-projected DEM tiles. The DEM files can then be fitted to absolute height laser altimeter points, or iteratively to each other, to form a DEM mosaic. Uncertainties in UV-VIS camera pointing and the need to accumulate a sufficiently good topographic S/N ratio necessitates the use of 1 km pixels for the UV-VIS derived DEMs. For Apollo Metric stereo, an internal camera geometry correction and a full photogrammetric block adjustment must be performed using ground- control points to derive a DEM. The image scale of Apollo Metric, as well as the stereo angle, allow for a DEM with 100 m pixels and a height accuracy of +/- 25m. Apollo Metric imagery had previously been used to derive contour maps for much of the lunar equatorial regions; however, to recover this information in digital form these maps must be digitized. Most of the mare areas mapped contain noticeable topographic noise. This results from

  17. Implications of different digital elevation models and preprocessing techniques to delineate debris flow inundation hazard zones in El Salvador

    NASA Astrophysics Data System (ADS)

    Anderson, E. R.; Griffin, R.; Irwin, D.

    2013-12-01

    Heavy rains and steep, volcanic slopes in El Salvador cause numerous landslides every year, posing a persistent threat to the population, economy and environment. Although potential debris inundation hazard zones have been delineated using digital elevation models (DEMs), some disparities exist between the simulated zones and actual affected areas. Moreover, these hazard zones have only been identified for volcanic lahars and not the shallow landslides that occur nearly every year. This is despite the availability of tools to delineate a variety of landslide types (e.g., the USGS-developed LAHARZ software). Limitations in DEM spatial resolution, age of the data, and hydrological preprocessing techniques can contribute to inaccurate hazard zone definitions. This study investigates the impacts of using different elevation models and pit filling techniques in the final debris hazard zone delineations, in an effort to determine which combination of methods most closely agrees with observed landslide events. In particular, a national DEM digitized from topographic sheets from the 1970s and 1980s provide an elevation product at a 10 meter resolution. Both natural and anthropogenic modifications of the terrain limit the accuracy of current landslide hazard assessments derived from this source. Global products from the Shuttle Radar Topography Mission (SRTM) and the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global DEM (ASTER GDEM) offer more recent data but at the cost of spatial resolution. New data derived from the NASA Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) in 2013 provides the opportunity to update hazard zones at a higher spatial resolution (approximately 6 meters). Hydrological filling of sinks or pits for current hazard zone simulation has previously been achieved through ArcInfo spatial analyst. Such hydrological processing typically only fills pits and can lead to drastic modifications of original elevation values

  18. Digital Elevation Models of the Lunar Surface

    NASA Astrophysics Data System (ADS)

    Cook, A. C.; Robinson, M. S.

    1999-01-01

    Several digital elevation models (DEMs) have been produced at a scale of 1km/pixel and covering approximately one-fifth of the lunar surface. These were produced mostly by semiautomatically matching the stereo available between Clementine UV/VIS images, although some localized DEMs have been produced by applying this technique to Apollo Metric stereo pairs, or by digitizing an existing Apollo Metric contour map. The DEMS that result from Clementine UV/VIS images, although Of Poorer height accuracy (1300-600 in for a single matched point) than the Clementine laser altimeter point measurements (<+/-100 m), do provide considerably higher spatial resolution (e.g., every kilometer vs. every tens of kilometers) and allow topography in the polar regions to be determined. Nadir-pointing Clementine UV-VIS stereo pairs are automatically stereo matched using a patch-based matcher and fed through A stereo intersection camera model to yield a digital terrain model (DTM) of longitude, latitude, and height points. The DTM for each stereo pair is then replotted and interpolated to form map-projected DEM tiles. The DEM files can then be fitted to absolute height laser altimeter points, or iteratively to each other, to form a DEM mosaic. Uncertainties in UV-VIS camera pointing and the need to accumulate a sufficiently good topographic S/N ratio necessitates the use of 1 km pixels for the UV-VIS derived DEMs. For Apollo Metric stereo, an internal camera geometry correction and a full photogrammetric block adjustment must be performed using ground- control points to derive a DEM. The image scale of Apollo Metric, as well as the stereo angle, allow for a DEM with 100 m pixels and a height accuracy of +/- 25m. Apollo Metric imagery had previously been used to derive contour maps for much of the lunar equatorial regions; however, to recover this information in digital form these maps must be digitized. Most of the mare areas mapped contain noticeable topographic noise. This results from

  19. INTEGER ELEVATION MODEL GRIDS FOR US EPA REGION 9

    EPA Science Inventory

    Integer Digital Elevation Models in GRID format for the mainland US administrative boundary of the US EPA Region 9 developed from the United States Geological Survey (USGS) National Elevation Dataset (NED). The administrative boundary is represented by the state boundaries of Ca...

  20. Elevation control system model for the DSS 13 antenna

    NASA Technical Reports Server (NTRS)

    Gawronski, W.; Mellstrom, J. A.

    1991-01-01

    In order to meet the requirements for precision pointing of 34-m antennas, adequate control design and simulation software have to be developed along with a detailed description of the supporting analytical tools. This article describes a control system model for the elevation drive of the DSS 13 antenna. The model allows one to simulate elevation dynamics, cross-coupled dynamics in azimuth and elevation, and RF pointing error. A modal state-space model of the antenna structure was obtained from its finite-element model with a free rotating tipping structure. Model reduction techniques were applied separately for the antenna model and rate-loop model, thereby reducing the system order to one-third of the original one while preserving its dynamic properties. Extensive simulation results illustrate properties of the model.

  1. Gulf Stream model. [which considers surface elevation deviations

    NASA Technical Reports Server (NTRS)

    1974-01-01

    Surface elevation deviations in the Gulf Stream region off the eastern coast of the United States between Wallops Island, Virginia and Miami, Florida were investigated. The main causes of surface elevation deviations are geoid perturbations due to the continental shelf and the geostrophic adjustment of the density field due to the Gulf Stream. Quantitative surface elevation profiles were calculated based on geophysical measurements of gravity anomalies and hydrographic data. The results are presented graphically along with contemporaneous weather data. Comparisons are made between the profiles based on hydrographic data and a mean theoretical model. The theory of geostrophic flows including some classical Gulf Stream models is also presented briefly.

  2. Lahar infrasound associated with Volcán Villarrica's 3 March 2015 eruption

    NASA Astrophysics Data System (ADS)

    Johnson, Jeffrey B.; Palma, Jose L.

    2015-08-01

    The paroxysmal 2015 eruption of Volcán Villarrica (Chile) produced a 2.5 h long lahar, which descended more than 20 km within the Rio Correntoso/Turbio drainage and destroyed two small bridges. A three-element infrasound array 10 km from the summit, and 4 km from the lahar's closest approach, was used to study the flow's progression. Array processing using cross-correlation lag times and semblance places constraints on the lahar's dynamics, including detection of an initial flow pulse that traveled from 2 to 12 km at an average speed of 38 m/s. Subsequently, the lahar signal evolved to a relatively stationary infrasonic tremor located 10 to 12 km from the vent and adjacent to a topographic notch, through which sound may have preferentially diffracted toward the recording site. This study demonstrates the powerful capabilities of infrasound arrays for lahar study and suggests their potential application for future hazard monitoring.

  3. Causes, Dynamics and Impacts of Lahars Generated by the April, 2015 Calbuco Eruption, Chile.

    NASA Astrophysics Data System (ADS)

    Russell, A. J.; Dussaillant, A. R.; Meier, C. I.; Rivera, A.; Barra, M. M.; Urzua, N. G.; Hernandez, J. F.; Napoleoni, F.; Gonzalez, C.

    2015-12-01

    Calbuco is a 2015m high, glacier capped, stratovolcano in the heavily populated Los Lagos district of southern Chile with a history of large volcanic eruptions in 1893-95, 1906-7, 1911-12, 1917, 1932, 1945, 1961 and 1972. Calbuco experienced a powerful 90 minute eruption at 18:04h on 22 April, 2015 followed by additional major eruptions at 01:00h and 13:10h on 23 & 30 April, respectively, resulting in the evacuation of 6500 people and the imposition of a 20 km radius exclusion zone. Pyroclastic flows descended into several river catchments radiating from the volcano with lahars travelling distances of up to 14 km, reaching populated areas. We present preliminary findings regarding the causes, dynamics and impacts of lahars generated by the April 2015 eruption. Pyroclastic flows melted glacier ice and snow generating the largest lahars in the Rio Este and Rio Blanco Sur on the southern flanks of the volcano. Lahar deposits in the Rio Blanco Norte were buried by pyroclastic flow deposits with measured temperatures of up to 282°C three months after emplacement. Lahar erosional impacts included bedrock erosion, alluvial channel incision, erosion of surficial deposits and the felling of large areas of forest. Depositional landforms included boulder run-ups on the outsides of channel bends, boulder clusters and large woody debris jams. Lahars deposited up to 8m of sediment within distal reaches. Deposits on the southern flanks of Calbuco indicate the passage of multiple pulses of contrasting rheology. Lahar occurrence and magnitude was controlled by the pre-eruption distribution of snow and ice on the volcano. Pre-existing lahar channels controlled flows to lower piedmont zones where routing was determined by palaeo lahar geomorphology. Ongoing erosion of proximal pyroclastic flow and lahar deposits provides large volumes of sediment to distal portions of fluvial systems radiating from Calbuco.

  4. Carving and adaptive drainage enforcement of grid digital elevation models

    NASA Astrophysics Data System (ADS)

    Soille, Pierre; Vogt, Jürgen; Colombo, Roberto

    2003-12-01

    An effective and widely used method for removing spurious pits in digital elevation models consists of filling them until they overflow. However, this method sometimes creates large flat regions which in turn pose a problem for the determination of accurate flow directions. In this study, we propose to suppress each pit by creating a descending path from it to the nearest point having a lower elevation value. This is achieved by carving, i.e., lowering, the terrain elevations along the detected path. Carving paths are identified through a flooding simulation starting from the river outlets. The proposed approach allows for adaptive drainage enforcement whereby river networks coming from other data sources are imposed to the digital elevation model only in places where the automatic river network extraction deviates substantially from the known networks. An improvement to methods for routing flow over flat regions is also introduced. Detailed results are presented over test areas of the Danube basin.

  5. Ice Cloud and Land Elevation Satellite (ICESat) Altimetry and Digital Elevation Models in the Polar Regions

    NASA Astrophysics Data System (ADS)

    Carabajal, C. C.; Boy, J.

    2012-12-01

    A global set of Ground Control Points (GCPs) from altimetry measurements from the Ice, Cloud and land Elevation Satellite (ICESat) has been produced with the support of the NASA's Earth Surface and Interior Program. The highest quality altimetry measurements that can be used for ground control have been selected by applying rigorous editing criteria. This database represents a key means to establishing a much-needed global topography reference frame to aid solid Earth application studies, particularly useful at high latitudes, where other topographic control is scarce. ICESat GCPs were used to characterize and quantify spatially varying elevation biases in Digital Elevation Models (DEMs) in the polar regions, assessing the horizontal and vertical accuracy of valuable topographic datasets produced by sensors like ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) missions, and datasets like GMTED2010 (Global Multi-resolution Terrain Elevation Data), developed by the USGS (United States Geological Survey) and NGA (National Geospatial-Intelligence Agency), a large improvement over the global GTOPO30 dataset. We have analyzed error statistics globally and per continent, in conjunction with MODIS (Moderate Resolution Imaging Spectro-Radiometer) and MERIS (Medium Resolution Imaging Spectrometer) land cover products, relief, topography and other DEM and altimetry specific parameters, and will present the results of these evaluations.

  6. Evaluation Digital Elevation Model Generated by Synthetic Aperture Radar Data

    NASA Astrophysics Data System (ADS)

    Makineci, H. B.; Karabörk, H.

    2016-06-01

    Digital elevation model, showing the physical and topographical situation of the earth, is defined a tree-dimensional digital model obtained from the elevation of the surface by using of selected an appropriate interpolation method. DEMs are used in many areas such as management of natural resources, engineering and infrastructure projects, disaster and risk analysis, archaeology, security, aviation, forestry, energy, topographic mapping, landslide and flood analysis, Geographic Information Systems (GIS). Digital elevation models, which are the fundamental components of cartography, is calculated by many methods. Digital elevation models can be obtained terrestrial methods or data obtained by digitization of maps by processing the digital platform in general. Today, Digital elevation model data is generated by the processing of stereo optical satellite images, radar images (radargrammetry, interferometry) and lidar data using remote sensing and photogrammetric techniques with the help of improving technology. One of the fundamental components of remote sensing radar technology is very advanced nowadays. In response to this progress it began to be used more frequently in various fields. Determining the shape of topography and creating digital elevation model comes the beginning topics of these areas. It is aimed in this work , the differences of evaluation of quality between Sentinel-1A SAR image ,which is sent by European Space Agency ESA and Interferometry Wide Swath imaging mode and C band type , and DTED-2 (Digital Terrain Elevation Data) and application between them. The application includes RMS static method for detecting precision of data. Results show us to variance of points make a high decrease from mountain area to plane area.

  7. The effects of wavelet compression on Digital Elevation Models (DEMs)

    USGS Publications Warehouse

    Oimoen, M.J.

    2004-01-01

    This paper investigates the effects of lossy compression on floating-point digital elevation models using the discrete wavelet transform. The compression of elevation data poses a different set of problems and concerns than does the compression of images. Most notably, the usefulness of DEMs depends largely in the quality of their derivatives, such as slope and aspect. Three areas extracted from the U.S. Geological Survey's National Elevation Dataset were transformed to the wavelet domain using the third order filters of the Daubechies family (DAUB6), and were made sparse by setting 95 percent of the smallest wavelet coefficients to zero. The resulting raster is compressible to a corresponding degree. The effects of the nulled coefficients on the reconstructed DEM are noted as residuals in elevation, derived slope and aspect, and delineation of drainage basins and streamlines. A simple masking technique also is presented, that maintains the integrity and flatness of water bodies in the reconstructed DEM.

  8. A Unified Low-Elevation-Angle Scintillation Model

    NASA Astrophysics Data System (ADS)

    Lee, C. H.; Cheung, K.-M.; Ho, C.

    2011-05-01

    Enabling communications at very low elevation angles can lengthen the duration of a tracking pass between a satellite and a ground station, which in turn can increase the amount of data return and possibly reduce the number of required supporting ground station tracking passes. Link performance, especially at very low angles and high frequencies, depends heavily on terrain, atmosphere, and weather conditions. Among the different contributions to attenuation, scintillation fading plays a very significant role and can impair the performance of the link. It is therefore necessary to accurately model the overall impact to the link due to scintillation fading. The current International Telecommunication Union ITU-R P.618-10 Recommendation describes three scintillation loss models as a function of elevation angle and percentage of time for which the loss exceeds a certain threshold. Implementation of the recommendation resulted in the uncovering of several issues. Particularly, it was identified that (i) iterative solutions to an implicit nonlinear exponential model, in some cases, are not guaranteed to exist, (ii) there is a discontinuity in fading values between models at the cross-over elevation angle, (iii) at certain low elevation angles scintillation from the shallow fade model generates unrealistically small losses, and (iv) for elevation angles lying between 4 and 5 deg, there are two applicable scintillation models that yield conflicting values. In this article, we develop a new approach to unify the different fading models within the current ITU recommendation and fully remove the discrepancies. We further validated our models with ITU-adopted scintillation data measured at Goonhilly, Great Britain, and data from several recent NASA Space Shuttle launches. This improved model was provisionally approved at the ITU International Meeting in Italy, November 2010, and is being evaluated by the ITU members for adoption into the next-version ITU Recommendation.

  9. Improving merge methods for grid-based digital elevation models

    NASA Astrophysics Data System (ADS)

    Leitão, J. P.; Prodanović, D.; Maksimović, Č.

    2016-03-01

    Digital Elevation Models (DEMs) are used to represent the terrain in applications such as, for example, overland flow modelling or viewshed analysis. DEMs generated from digitising contour lines or obtained by LiDAR or satellite data are now widely available. However, in some cases, the area of study is covered by more than one of the available elevation data sets. In these cases the relevant DEMs may need to be merged. The merged DEM must retain the most accurate elevation information available while generating consistent slopes and aspects. In this paper we present a thorough analysis of three conventional grid-based DEM merging methods that are available in commercial GIS software. These methods are evaluated for their applicability in merging DEMs and, based on evaluation results, a method for improving the merging of grid-based DEMs is proposed. DEMs generated by the proposed method, called MBlend, showed significant improvements when compared to DEMs produced by the three conventional methods in terms of elevation, slope and aspect accuracy, ensuring also smooth elevation transitions between the original DEMs. The results produced by the improved method are highly relevant different applications in terrain analysis, e.g., visibility, or spotting irregularities in landforms and for modelling terrain phenomena, such as overland flow.

  10. Characteristics of the syneruptive-spouted type lahar generated by the September 2014 eruption of Mount Ontake, Japan

    NASA Astrophysics Data System (ADS)

    Sasaki, Hisashi; Chiba, Tatsuro; Kishimoto, Hiroshi; Naruke, Shino

    2016-08-01

    Mount Ontake erupted at 11:52 am on September 27, 2014, which generated pyroclastic density currents, ballistic projectiles, ash falls, and a small-scale lahar that spouted directly from craters formed by the eruption. Because this lahar may have been generated by water released from within these craters, we refer to this lahar as a "syneruptive-spouted type lahar" in this study. The lahar of the 2014 eruption was small relative to the other syneruptive type lahars reported in the past that were snowmelt type or crater lake breakout type lahars. Nevertheless, in the 2014 event, the syneruptive-spouted type lahar extended approximately 5 km downstream from the Jigokudani crater via the Akagawa River, with an estimated total volume of ~1.2 × 105 m3. We have reviewed other representative syneruptive-spouted type lahars that have been reported in Japan. The syneruptive-spouted type lahar attributed to the September 2014 eruption had the longest runout distance and largest volume of all cases studied. The mineral assemblage identified from samples of the lahar deposits is similar to that of ash-fall deposits from the same eruption. Previous workers deduced that the ash was derived mainly from shallow depths (within 2 km of the surface). The syneruptive-spouted type lahar deposits are therefore also considered to have originated from shallow depths. A syneruptive-spouted type lahar is a small-scale phenomenon that causes little direct damage to infrastructure, but has long-term influence on water quality. Increases in turbidity and decreases in pH are expected to occur in the Mount Ontake area downstream of Nigorisawa after heavy rainfall events in the future. Therefore, the potential indirect (but long term) damage of syneruptive-spouted type lahars should be considered for hazard mapping and planning volcanic disaster prevention measures.

  11. Variations in population exposure and sensitivity to lahar hazards from Mount Rainier, Washington

    NASA Astrophysics Data System (ADS)

    Wood, Nathan; Soulard, Christopher

    2009-12-01

    Although much has been done to understand, quantify, and delineate volcanic hazards, there are fewer efforts to assess societal vulnerability to these hazards, particularly demographic differences in exposed populations or spatial variations in exposure to regional hazards. To better understand population diversity in volcanic hazard zones, we assess the number and types of people in a single type of hazard zone (lahars) for 27 communities downstream of Mount Rainier, Washington (USA). Using various socioeconomic and hazard datasets, we estimate that there are more than 78 000 residents, 59 000 employees, several dependent-population facilities (e.g., child-day-care centers, nursing homes) and numerous public venues (e.g., churches, hotels, museums) in a Mount Rainier lahar-hazard zone. We find that communities vary in the primary category of individuals in lahar-prone areas—exposed populations are dominated by residents in some communities (e.g., Auburn), employees in others (e.g., Tacoma), and tourists likely outnumber both of these groups in yet other areas (e.g., unincorporated Lewis County). Population exposure to potential lahar inundation varies considerably—some communities (e.g., Auburn) have large numbers of people but low percentages of them in hazard zones, whereas others (e.g., Orting) have fewer people but they comprise the majority of a community. A composite lahar-exposure index is developed to help emergency managers understand spatial variations in community exposure to lahars and results suggest that Puyallup has the highest combination of high numbers and percentages of people and assets in lahar-prone areas. Risk education and preparedness needs will vary based on who is threatened by future lahars, such as residents, employees, tourists at a public venue, or special-needs populations at a dependent-care facility. Emergency managers must first understand the people whom they are trying to prepare before they can expect these people to take

  12. Modeling low elevation GPS signal propagation in maritime atmospheric ducts

    NASA Astrophysics Data System (ADS)

    Zhang, Jinpeng; Wu, Zhensen; Wang, Bo; Wang, Hongguang; Zhu, Qinglin

    2012-05-01

    Using the parabolic wave equation (PWE) method, we model low elevation GPS L1 signal propagation in maritime atmospheric ducts. To consider sea surface impedance, roughness, and the effects of earth's curvature, we propose a new initial field model for the GPS PWE split-step solution. On the basis of the comparison between the proposed model and the conventional initial field model for a smooth, perfectly conducting sea surface on a planar earth, we conclude that both the amplitude and phase of the initial field are influenced by surface impedance and roughness, and that the interference behavior between direct and reflected GPS rays is affected by earth's curvature. The performance of the proposed model is illustrated with examples of low elevation GPS L1 signal propagation in three types of ducts: an evaporation duct, a surface-based duct, and an elevated duct. The GPS PWE is numerically implemented using the split-step discrete mixed Fourier transform algorithm to enforce impedance-type boundary conditions at the rough sea surface. Because the GPS signal is right hand circularly polarized, we calculate its power strength by combining the propagation predictions of the horizontally and the vertically polarized components. The effects of the maritime atmospheric ducts on low elevation GPS signal propagation are demonstrated according to the presented examples, and the potential applications of the GPS signals affected by ducts are discussed.

  13. How processing digital elevation models can affect simulated water budgets

    USGS Publications Warehouse

    Kuniansky, E.L.; Lowery, M.A.; Campbell, B.G.

    2009-01-01

    For regional models, the shallow water table surface is often used as a source/sink boundary condition, as model grid scale precludes simulation of the water table aquifer. This approach is appropriate when the water table surface is relatively stationary. Since water table surface maps are not readily available, the elevation of the water table used in model cells is estimated via a two-step process. First, a regression equation is developed using existing land and water table elevations from wells in the area. This equation is then used to predict the water table surface for each model cell using land surface elevation available from digital elevation models (DEM). Two methods of processing DEM for estimating the land surface for each cell are commonly used (value nearest the cell centroid or mean value in the cell). This article demonstrates how these two methods of DEM processing can affect the simulated water budget. For the example presented, approximately 20% more total flow through the aquifer system is simulated if the centroid value rather than the mean value is used. This is due to the one-third greater average ground water gradients associated with the centroid value than the mean value. The results will vary depending on the particular model area topography and cell size. The use of the mean DEM value in each model cell will result in a more conservative water budget and is more appropriate because the model cell water table value should be representative of the entire cell area, not the centroid of the model cell.

  14. Extraction of terrain features from digital elevation models

    USGS Publications Warehouse

    Price, Curtis V.; Wolock, David M.; Ayers, Mark A.

    1989-01-01

    Digital elevation models (DEMs) are being used to determine variable inputs for hydrologic models in the Delaware River basin. Recently developed software for analysis of DEMs has been applied to watershed and streamline delineation. The results compare favorably with similar delineations taken from topographic maps. Additionally, output from this software has been used to extract other hydrologic information from the DEM, including flow direction, channel location, and an index describing the slope and shape of a watershed.

  15. Lahar Hazards at Concepción volcano, Nicaragua

    USGS Publications Warehouse

    Vallance, J.W.; Schilling, S.P.; Devoli, G.; Howell, M.M.

    2001-01-01

    Concepción is one of Nicaragua’s highest and most active volcanoes. The symmetrical cone occupies the northeastern half of a dumbbell shaped island called Isla Ometepa. The dormant volcano, Maderas, occupies the southwest half of the island. A narrow isthmus connects Concepción and Maderas volcanoes. Concepción volcano towers more than 1600 m above Lake Nicaragua and is within 5 to 10 km of several small towns situated on its aprons at or near the shoreline. These towns have a combined population of nearly 5,000. The volcano has frequently produced debris flows (watery flows of mud, rock, and debris—also known as lahars when they occur on a volcano) that could inundate these nearby populated areas. Concepción volcano has erupted more than 25 times in the last 120 years. Its first recorded activity was in AD 1883. Eruptions in the past century, most of which have originated from a small summit crater, comprise moderate explosions, ash that falls out of eruption plumes (called tephra), and occasional lava flows. Near the summit area, there are accumulations of rock that were emplaced hot (pyroclastic deposits), most of which were hot enough to stick together during deposition (a process called welding). These pyroclastic rocks are rather weak, and tend to break apart easily. The loose volcanic rock remobilizes during heavy rain to form lahars. Volcanic explosions have produced blankets of tephra that are distributed downwind, which on Isla Ometepe is mostly to the west. Older deposits at the west end of the island that are up to 1 m thick indicate larger explosive events have happened at Concepción volcano in prehistoric time. Like pyroclastic-flow deposits, loose tephra on the steep slopes of the volcano provides source material that heavy rainstorms and earthquakes can mobilize to trigger debris flow.

  16. Reducing risk from lahar hazards: concepts, case studies, and roles for scientists

    USGS Publications Warehouse

    Pierson, Thomas C.; Wood, Nathan J.; Driedger, Carolyn L.

    2014-01-01

    Lahars are rapid flows of mud-rock slurries that can occur without warning and catastrophically impact areas more than 100 km downstream of source volcanoes. Strategies to mitigate the potential for damage or loss from lahars fall into four basic categories: (1) avoidance of lahar hazards through land-use planning; (2) modification of lahar hazards through engineered protection structures; (3) lahar warning systems to enable evacuations; and (4) effective response to and recovery from lahars when they do occur. Successful application of any of these strategies requires an accurate understanding and assessment of the hazard, an understanding of the applicability and limitations of the strategy, and thorough planning. The human and institutional components leading to successful application can be even more important: engagement of all stakeholders in hazard education and risk-reduction planning; good communication of hazard and risk information among scientists, emergency managers, elected officials, and the at-risk public during crisis and non-crisis periods; sustained response training; and adequate funding for risk-reduction efforts. This paper reviews a number of methods for lahar-hazard risk reduction, examines the limitations and tradeoffs, and provides real-world examples of their application in the U.S. Pacific Northwest and in other volcanic regions of the world. An overriding theme is that lahar-hazard risk reduction cannot be effectively accomplished without the active, impartial involvement of volcano scientists, who are willing to assume educational, interpretive, and advisory roles to work in partnership with elected officials, emergency managers, and vulnerable communities.

  17. A study of the Taisho lahar generated by the 1926 eruption of Tokachidake Volcano, central Hokkaido, Japan, and implications for the generation of cohesive lahars

    NASA Astrophysics Data System (ADS)

    Uesawa, Shimpei

    2014-01-01

    Understanding the generation mechanisms of lahars is important for improving volcanic hazard assessments. The Taisho lahar (TL) was generated during the 1926 eruption of Tokachidake Volcano, Japan, and was considered a typical snowmelt lahar caused by the runout of hot debris onto a snow-covered slope. A similar mechanism produced a huge mud flow during the 1985 eruption of Nevado del Ruiz, Colombia. However, the origin of water in such lahars remains a controversial topic because the calculated water mass is based on the assumption that all of the snow on the runout area of the TL was melted, although this is much less than the estimated water volume in the TL estimated by previous studies. I have re-examined proximal deposits of the TL and their paleomagnetic characteristics in order to better understand the eruption sequence and formation of the TL. The TL produced two debris avalanche deposits and a surge-like deposit that had relatively high emplacement temperature (~ 350 °C). The deposits are composed of hydrothermally altered andesitic gravel, sand and mud. The high clay content (3-5 wt.% clay in the < 2 mm fraction) and sedimentary characteristics indicate that the flow was a cohesive lahar, most likely induced by collapse of a hydrothermally altered pyroclastic cone (hypocenter). The presence of the surge deposit indicates that the TL was not caused by simple collapse of a cinder cone but by a phreatic explosion that resulted in sector collapse. This suggests that the hydrothermal system was related to the 1926 eruption. The present-day volcano has a large hydrothermal system (1 × 106 m3 water) beneath the active crater. This study indicates that hydrothermal system explosions can trigger cohesive lahars that contain both snow melt and hydrothermal pore water, and this indicates the need to monitor hydrothermal systems.

  18. The Rational Polynomial Coefficients Modification Using Digital Elevation Models

    NASA Astrophysics Data System (ADS)

    Alidoost, F.; Azizi, A.; Arefi, H.

    2015-12-01

    The high-resolution satellite imageries (HRSI) are as primary dataset for different applications such as DEM generation, 3D city mapping, change detection, monitoring, and deformation detection. The geo-location information of HRSI are stored in metadata called Rational Polynomial Coefficients (RPCs). There are many methods to improve and modify the RPCs in order to have a precise mapping. In this paper, an automatic approach is presented for the RPC modification using global Digital Elevation Models. The main steps of this approach are: relative digital elevation model generation, shift parameters calculation, sparse point cloud generation and shift correction, and rational polynomial fitting. Using some ground control points, the accuracy of the proposed method is evaluated based on statistical descriptors in which the results show that the geo-location accuracy of HRSI can be improved without using Ground Control Points (GCPs).

  19. Failure prediction in automatically generated digital elevation models

    NASA Astrophysics Data System (ADS)

    Gooch, M. J.; Chandler, J. H.

    2001-10-01

    Developments in digital photogrammetry have provided the ability to generate digital elevation models (DEMs) automatically and are increasingly used by geoscientists. Using overlapping imagery, dense grids of digital elevations can be collected at high speeds (150 points per second) with a high level of accuracy. The trend towards using PC-based hardware, the widespread use of geographical information systems, and the forthcoming availability of high-resolution satellite imagery over the Internet at ever lower costs mean that the use of automated digital photogrammetry for elevation modelling is likely to become more widespread. Automation can reduce the need for an in-depth knowledge of the subject thus rendering the technology an option for more users. One criticism of the trend towards the automated "black box" approach is the lack of quality control procedures within the software, particularly with reference to identifying areas of the DEM with low accuracy. The traditional method of accuracy assessment is through the use of check point data (data collected by an independent method which has a higher level of accuracy against which the DEM can be compared). Check point data are, however, rarely available and it is typically recommended that the user manually check and edit the data using stereo viewing methods, a potentially lengthy process which can negate the obvious speed advantages brought about by automation. A data processing model has been developed that is capable of identifying areas where elevations are unreliable and to which the user should pay attention when editing and checking the data. The software model developed will be explained and described in detail in the paper. Results from tests on different scales of imagery, different types of imagery and other software packages will also be presented to demonstrate the efficacy and significantly the generality of the technique with other digital photogrammetric software systems.

  20. Elevator model based on a tiny PLC for teaching automation

    NASA Astrophysics Data System (ADS)

    Kim, Kee Hwan; Lee, Young Dae

    2005-12-01

    The development of control related applications requires knowledge of different subject matters like mechanical components, control equipment and physics. To understand the behavior of these heterogeneous applications is not easy especially the students who begin to study the electronic engineering. In order to introduce to them the most common components and skills necessary to put together a functioning automated system, we have designed a simple elevator model controlled by a PLC which was designed based on a microcontroller.

  1. Digital elevation model visibility including Earth's curvature and atmosphere refraction

    NASA Astrophysics Data System (ADS)

    Santossilva, Ewerton; Vieiradias, Luiz Alberto

    1990-03-01

    There are some instances in which the Earth's curvature and the atmospheric refraction, optical or electronic, are important factors when digital elevation models are used for visibility calculations. This work deals with this subject, suggesting a practical approach to solve this problem. Some examples, from real terrain data, are presented. The equipment used was an IBM-PC like computer with a SITIM graphic card.

  2. Tree-ring based reconstruction of past lahar activity at Popocat

    NASA Astrophysics Data System (ADS)

    Bollschweiler, M.; Stoffel, M.; Vázquez Selem, L.; Palacios, D.

    2009-04-01

    Lahars are rapid, saturated flows of water and rock fragments that occur on volcanoes and that can be triggered either by volcanic activity or by intense precipitation falling on unconsolidated volcanic deposits. As their occurrence is unpredictable, as the flow contains sometimes considerably large rock fragments and as the flow is able to travel long distances even on gentle gradients, lahars represent one of the most destructive natural disasters in terms of loss of human lives and property damage. In order to realistically assess hazards, knowledge on the occurrence and timing of past lahar activity is of crucial importance. However, archival data on past events is usually scarce or completely missing. Tree-ring records have repeatedly proved to be a reliable data source for the reconstruction of past geomorphic events. However, tree rings have hardly ever been applied for the identification of past lahars. Therefore, it was the aim of this study (i) to identify and describe disturbances in tree growth induced by well-documented lahar events and on this basis (ii) to recognize older, unknown lahar events with tree-ring analyses. Based on these goals, we collected 140 tree-ring series from 62 trees (Abies religiosa, Pinus hartwegii, Pinus ayacahuite) standing inside or adjacent to the lahar channel in the Huiloac gorge at Popocatépetl volcano, central Mexico. Most commonly, the known lahar events of 1997 and 2001 resulted in abrupt changes in tree-ring width as well as injuries. The same growth disturbances could be identified in the tree-ring series, indicating that five previously unknown lahar events would have occurred during the 20th century. Popocatépetl is one of the best surveyed volcanoes in the world and past eruptions are precisely noted in archives. As most of these unknown events occurred during periods with no volcanic activity, we believe that they were rainfall-induced rather than related to volcanic activity. This study revealed the potential

  3. A wind model for an elevated STOL-port configuration

    NASA Technical Reports Server (NTRS)

    Peterka, J. A.; Cermak, J. E.

    1974-01-01

    Measurements of mean velocity magnitude and direction as well as three-dimensional turbulence intensity were made in the flow over a model of an elevated STOL-port. A 1:300 scale model was placed in a wind tunnel flow simulating the mean velocity profile and turbulence characteristics of atmospheric winds over a typical city environment excluding detailed wake structures of possible nearby buildings. Hot-wire anemometer measurements of velocity and turbulence were made along approach and departure paths of aircraft operating on the runway centerline and at specified lateral distances from the centerline. Approach flow directions simulated were 0 and 30 degrees to the runway centerline.

  4. Online, On Demand Access to Coastal Digital Elevation Models

    NASA Astrophysics Data System (ADS)

    Long, J.; Bristol, S.; Long, D.; Thompson, S.

    2014-12-01

    Process-based numerical models for coastal waves, water levels, and sediment transport are initialized with digital elevation models (DEM) constructed by interpolating and merging bathymetric and topographic elevation data. These gridded surfaces must seamlessly span the land-water interface and may cover large regions where the individual raw data sources are collected at widely different spatial and temporal resolutions. In addition, the datasets are collected from different instrument platforms with varying accuracy and may or may not overlap in coverage. The lack of available tools and difficulties in constructing these DEMs lead scientists to 1) rely on previously merged, outdated, or over-smoothed DEMs; 2) discard more recent data that covers only a portion of the DEM domain; and 3) use inconsistent methodologies to generate DEMs. The objective of this work is to address the immediate need of integrating land and water-based elevation data sources and streamline the generation of a seamless data surface that spans the terrestrial-marine boundary. To achieve this, the U.S. Geological Survey (USGS) is developing a web processing service to format and initialize geoprocessing tasks designed to create coastal DEMs. The web processing service is maintained within the USGS ScienceBase data management system and has an associated user interface. Through the map-based interface, users define a geographic region that identifies the bounds of the desired DEM and a time period of interest. This initiates a query for elevation datasets within federal science agency data repositories. A geoprocessing service is then triggered to interpolate, merge, and smooth the data sources creating a DEM based on user-defined configuration parameters. Uncertainty and error estimates for the DEM are also returned by the geoprocessing service. Upon completion, the information management platform provides access to the final gridded data derivative and saves the configuration parameters

  5. A linear model of stationary elevator traveling and compensation cables

    NASA Astrophysics Data System (ADS)

    Zhu, W. D.; Ren, H.

    2013-06-01

    Based on a recent asymptotic analysis of a nonlinear model of a slack cable, a computationally efficient, linear model is developed for calculating the natural frequencies, mode shapes, and dynamic responses of stationary elevator traveling and compensation cables. The linear cable model consists of two vertical cable segments connected by a half-circular lower loop. The two vertical cable segments are modeled as a string with a variable tension due to the weight of the cable. The horizontal displacements of the cable segments consist of boundary-induced displacements and relative elastic displacements, where the boundary-induced displacements are interpolated from the displacements of the two lower ends of the cable segments, and the relative elastic displacements satisfy the corresponding homogeneous boundary conditions of the cable segments. The horizontal displacement of the lower loop is interpolated from those of the two lower ends of the two cable segments, and the bending stiffness of the lower loop is modeled by a spring with a constant stiffness, which can be calculated from the nonlinear model. Given a car position, the natural frequencies and mode shapes of an elevator traveling or compensation cable are calculated using the linear model and compared with those from the nonlinear model. The calculated natural frequencies are also compared with those from a full-scale experiment. In addition, the dynamic responses of a cable under a boundary excitation are calculated and compared with those from the nonlinear model. There is a good agreement between the predictions from the linear and nonlinear models and between the measured natural frequencies from the full-scale experiment and the corresponding calculated ones.

  6. Catastrophic precipitation-triggered lahar at Casita volcano, Nicaragua: Occurrence, bulking and transformation

    USGS Publications Warehouse

    Scott, K.M.; Vallance, J.W.; Kerle, N.; Macias, J.L.; Strauch, W.; Devoli, G.

    2005-01-01

    A catastrophic lahar began on 30 October 1998, as hurricane precipitation triggered a small flank collapse of Casita volcano, a complex and probably dormant stratovolcano. The initial rockslide-debris avalanche evolved on the flank to yield a watery debris flood with a sediment concentration less than 60 per cent by volume at the base of the volcano. Within 2-5 km, however, the watery flow entrained (bulked) enough sediment to transform entirely to a debris flow. The debris flow, 6 km downstream and 1??2 km wide and 3 to 6 m deep, killed 2500 people, nearly the entire populations of the communities of El Porvenir and Rolando Rodriguez. These 'new towns' were developed in a prehistoric lahar pathway: at least three flows of similar size since 8330 14C years BP are documented by stratigraphy in the same 30-degree sector. Travel time between perception of the flow and destruction of the towns was only 2??5-3??0 minutes. The evolution of the flow wave occurred with hydraulic continuity and without pause or any extraordinary addition of water. The precipitation trigger of the Casita lahar emphasizes the nee d, in volcano hazard assessments, for including the potential for non-eruption-related collapse lahars with the more predictable potential of their syneruption analogues. The flow behaviour emphasizes that volcano collapses can yield not only volcanic debris avalanches with restricted runouts, but also mobile lahars that enlarge by bulking as they flow. Volumes and hence inundation areas of collapse-runout lahars can increase greatly beyond their sources: the volume of the Casita lahar bulked to at least 2??6 times the contributing volume of the flank collapse and 4??2 times that of the debris flood. At least 78 per cent of the debris flow matrix (sediment < -1??0??; 2 mm) was entrained during flow. Copyright c 2004 John Wiley & Sons, Ltd.

  7. Elevated temperature alters carbon cycling in a model microbial community

    NASA Astrophysics Data System (ADS)

    Mosier, A.; Li, Z.; Thomas, B. C.; Hettich, R. L.; Pan, C.; Banfield, J. F.

    2013-12-01

    Earth's climate is regulated by biogeochemical carbon exchanges between the land, oceans and atmosphere that are chiefly driven by microorganisms. Microbial communities are therefore indispensible to the study of carbon cycling and its impacts on the global climate system. In spite of the critical role of microbial communities in carbon cycling processes, microbial activity is currently minimally represented or altogether absent from most Earth System Models. Method development and hypothesis-driven experimentation on tractable model ecosystems of reduced complexity, as presented here, are essential for building molecularly resolved, benchmarked carbon-climate models. Here, we use chemoautotropic acid mine drainage biofilms as a model community to determine how elevated temperature, a key parameter of global climate change, regulates the flow of carbon through microbial-based ecosystems. This study represents the first community proteomics analysis using tandem mass tags (TMT), which enable accurate, precise, and reproducible quantification of proteins. We compare protein expression levels of biofilms growing over a narrow temperature range expected to occur with predicted climate changes. We show that elevated temperature leads to up-regulation of proteins involved in amino acid metabolism and protein modification, and down-regulation of proteins involved in growth and reproduction. Closely related bacterial genotypes differ in their response to temperature: Elevated temperature represses carbon fixation by two Leptospirillum genotypes, whereas carbon fixation is significantly up-regulated at higher temperature by a third closely related genotypic group. Leptospirillum group III bacteria are more susceptible to viral stress at elevated temperature, which may lead to greater carbon turnover in the microbial food web through the release of viral lysate. Overall, this proteogenomics approach revealed the effects of climate change on carbon cycling pathways and other

  8. Sparse Representation and Multiscale Methods - Application to Digital Elevation Models

    NASA Astrophysics Data System (ADS)

    Stefanescu, R. E. R.; Patra, A. K.; Bursik, M. I.

    2014-12-01

    In general, a Digital Elevation Model (DEM) is produced either digitizing existing maps and elevation values are interpolated from the contours, or elevation information is collected from stereo imagery on digital photogrammetric workstations. Both methods produce a DEM to the required specification, but each method contains a variety of possible production scenarios, and each method results in DEM cells with totally different character. Common artifacts found in DEM are missing-values at different location which can influence the output of the application that uses this particular DEM. In this work we introduce a numerically-stable multiscale scheme to evaluate the missing-value DEM's quantity of interest (elevation, slope, etc.). This method is very efficient for the case when dealing with large high resolution DEMs that cover large area, resulting in O(106-1010) data points. Our scheme relies on graph-based algorithms and low-rank approximations of the entire adjacency matrix of the DEM's graph. When dealing with large data sets such as DEMs, the Laplacian or kernel matrix resulted from the interaction of the data points is stupendously big. One needs to identify a subspace that capture most of the action of the kernel matrix. By the application of a randomized projection on the graph affinity matrix, a well-conditioned basis is identified for it numerical range. This basis is later used in out-of-sample extension at missing-value location. In many cases, this method beats its classical competitors in terms of accuracy, speed, and robustness.

  9. Mount Baker lahars and debris flows, ancient, modern, and future

    USGS Publications Warehouse

    Tucker, David S; Scott, Kevin M.; Grossman, Eric E.; Linneman, Scott

    2014-01-01

    Holocene lahars and large debris flows (>106 m3) have left recognizable deposits in the Middle Fork Nooksack valley. A debris flow in 2013 resulting from a landslide in a Little Ice Age moraine had an estimated volume of 100,000 m3, yet affected turbidity for the entire length of the river, and produced a slug of sediment that is currently being reworked and remobilized in the river system. Deposits of smaller-volume debris flows, deposited as terraces in the upper valley, may be entirely eroded within a few years. Consequently, the geologic record of small debris flows such as those that occurred in 2013 is probably very fragmentary. Small debris flows may still have significant impacts on hydrology, biology, and human uses of rivers downstream. Impacts include the addition of waves of fine sediment to stream loads, scouring or burying salmon-spawning gravels, forcing unplanned and sudden closure of municipal water intakes, damaging or destroying trail crossings, extending river deltas into estuaries, and adding to silting of harbors near river mouths.

  10. Accuracy and reliability of the Hungarian digital elevation model (DEM)

    NASA Astrophysics Data System (ADS)

    Detrekoi, Akos; Melykuti, Gabor; Szabo, Gyorgy

    1994-08-01

    In the period 1991-92 a 50 X 50 meter grid digital elevation model (DEM) was created in Hungary. The design and the quality control of DEM are discussed in this paper. The paper has three parts: (1) the data acquisition methods for DEM by scanning and photogrammetry are discussed, (2) a general overview about the accuracy and reliability of DEMs is given, and (3) the algorithm for the checking of data and some general conclusions about the control activity of the Hungarian DEM are reviewed.

  11. Validation of Orthorectified Interferometric Radar Imagery and Digital Elevation Models

    NASA Technical Reports Server (NTRS)

    Smith Charles M.

    2004-01-01

    This work was performed under NASA's Verification and Validation (V&V) Program as an independent check of data supplied by EarthWatch, Incorporated, through the Earth Science Enterprise Scientific Data Purchase (SDP) Program. This document serves as the basis of reporting results associated with validation of orthorectified interferometric interferometric radar imagery and digital elevation models (DEM). This validation covers all datasets provided under the first campaign (Central America & Virginia Beach) plus three earlier missions (Indonesia, Red River: and Denver) for a total of 13 missions.

  12. Registering Thematic Mapper imagery to digital elevation models

    NASA Technical Reports Server (NTRS)

    Frew, J.

    1984-01-01

    The problems encountered when attempting to register Landsat Thematic Mapper (TM) data to U.S. geological survey digital elevation models (DEMs) are examined. It is shown that TM and DEM data are not available in the same map projection, necessitating geometric transformation of one of the data type, that the TM data are not accurately located in their nominal projection, and that TM data have higher resolution than most DEM data, but oversampling the DEM data to TM resolution introduces systematic noise. Further work needed in this area is discussed.

  13. Applications of hydrologic information automatically extracted from digital elevation models

    USGS Publications Warehouse

    Jenson, S.K.

    1991-01-01

    Digital elevation models (DEMs) can be used to derive a wealth of information about the morphology of a land surface. Traditional raster analysis methods can be used to derive slope, aspect, and shaded relief information; recently-developed computer programs can be used to delineate depressions, overland flow paths, and watershed boundaries. These methods were used to delineate watershed boundaries for a geochemical stream sediment survey, to compare the results of extracting slope and flow paths from DEMs of varying resolutions, and to examine the geomorphology of a Martian DEM. -Author

  14. Void-Filled SRTM Digital Elevation Model of Afghanistan

    USGS Publications Warehouse

    Chirico, Peter G.; Barrios, Boris

    2005-01-01

    EXPLANATION The purpose of this data set is to provide a single consistent elevation model to be used for national scale mapping, GIS, remote sensing applications, and natural resource assessments for Afghanistan's reconstruction. For 11 days in February of 2000, the National Aeronautics and Space Administration (NASA), the National Geospatial-Intelligence Agency ian Space Agency (ASI) flew X-band and C-band radar interferometry onboard the Space Shuttle Endeavor. The mission covered the Earth between 60?N and 57?S and will provide interferometric digital elevation models (DEMs) of approximately 80% of the Earth's land mass when processing is complete. The radar-pointing angle was approximately 55? at scene center. Ascending and descending orbital passes generated multiple interferometric data scenes for nearly all areas. Up to eight passes of data were merged to form the final processed Shuttle Radar Topography Mission (SRTM) DEMs. The effect of merging scenes averages elevation values recorded in coincident scenes and reduces, but does not completely eliminate, the amount of area with layover and terrain shadow effects. The most significant form of data processing for the Afghanistan DEM was gap-filling areas where the SRTM data contained a data void. These void areas are as a result of radar shadow, layover, standing water, and other effects of terrain as well as technical radar interferometry phase unwrapping issues. To fill these gaps, topographic contours were digitized from 1:200,000 - scale Soviet General Staff Topographic Maps which date from the middle to late 1980's. Digital contours were gridded to form elevation models for void areas and subsequently were merged with the SRTM data through GIS and image processing techniques. The data contained in this publication includes SRTM DEM quadrangles projected and clipped in geographic coordinates for the entire country. An index of all available SRTM DEM quadrangles is displayed here: Index_Geo_DD.pdf. Also

  15. Integrated firn elevation change model for glaciers and ice caps

    NASA Astrophysics Data System (ADS)

    Saß, Björn; Sauter, Tobias; Braun, Matthias

    2016-04-01

    We present the development of a firn compaction model in order to improve the volume to mass conversion of geodetic glacier mass balance measurements. The model is applied on the Arctic ice cap Vestfonna. Vestfonna is located on the island Nordaustlandet in the north east of Svalbard. Vestfonna covers about 2400 km² and has a dome like shape with well-defined outlet glaciers. Elevation and volume changes measured by e.g. satellite techniques are becoming more and more popular. They are carried out over observation periods of variable length and often covering different meteorological and snow hydrological regimes. The elevation change measurements compose of various components including dynamic adjustments, firn compaction and mass loss by downwasting. Currently, geodetic glacier mass balances are frequently converted from elevation change measurements using a constant conversion factor of 850 kg m‑³ or the density of ice (917 kg m‑³) for entire glacier basins. However, the natural conditions are rarely that static. Other studies used constant densities for the ablation (900 kg m‑³) and accumulation (600 kg m‑³) areas, whereby density variations with varying meteorological and climate conditions are not considered. Hence, each approach bears additional uncertainties from the volume to mass conversion that are strongly affected by the type and timing of the repeat measurements. We link and adapt existing models of surface energy balance, accumulation and snow and firn processes in order to improve the volume to mass conversion by considering the firn compaction component. Energy exchange at the surface is computed by a surface energy balance approach and driven by meteorological variables like incoming short-wave radiation, air temperature, relative humidity, air pressure, wind speed, all-phase precipitation, and cloud cover fraction. Snow and firn processes are addressed by a coupled subsurface model, implemented with a non-equidistant layer

  16. Preparation of the Digital Elevation Model for Orthophoto CR Production

    NASA Astrophysics Data System (ADS)

    Švec, Z.; Pavelka, K.

    2016-06-01

    The Orthophoto CR is produced in co-operation with the Land Survey Office and the Military Geographical and Hydrometeorological Office. The product serves to ensure a defence of the state, integrated crisis management, civilian tasks in support of the state administration and the local self-government of the Czech Republic as well. It covers the whole area of the Republic and for ensuring its up-to-datedness is reproduced in the biennial period. As the project is countrywide, it keeps the project within the same parameters in urban and rural areas as well. Due to economic reasons it cańt be produced as a true ortophoto because it requires large side and forward overlaps of the aerial photographs and a preparation of the digital surface model instead of the digital terrain model. Use of DTM without some objects of DSM for orthogonalization purposes cause undesirable image deformations in the Orthophoto. There are a few data sets available for forming a suitable elevation model. The principal source should represent DTMs made from data acquired by the airborne laser scanning of the entire area of the Czech Republic that was carried out in the years 2009-2013, the DMR4G in the grid form and the DMR5G in TIN form respectively. It can be replenished by some vector objects (bridges, dams, etc.) taken from the geographic base data of the Czech Republic or obtained by new stereo plotting. It has to be taken into account that the option of applying DSM made from image correlation is also available. The article focuses on the possibilities of DTM supplement for ortogonalization. It looks back to the recent transition from grid to hybrid elevation models, problems that occurred, its solution and getting some practical remarks. Afterwards it assesses the current state and deals with the options for updating the model. Some accuracy analysis are included.

  17. Post-glacial lahars of the Sandy River Basin, Mount Hood, Oregon.

    USGS Publications Warehouse

    Cameron, K.A.; Pringle, P.

    1986-01-01

    Within the last 10 000 years, three significant lahar-producing periods have occurred at Mount Hood, Oregon. The Timberline eruptive period occurred between 1400 and 1800 years BP. It was by far the most voluminous of the three periods, producing enough clastic debris to bury the glacial topography of the southwest face of the mountain beneath a smooth debris fan. Timberline-age lahars traveled the length of the Zigzag and Sandy Rivers, a distance in excess of 90 km. Between 400 and 600 years BP, the Zigzag eruptive period produced deposits along the middle reaches of the Zigzag River and the upper Sandy River. The Old Maid eruptive period occurred between 180 and 270 years BP. A single lahar flowed down the Sandy River at least as far as Brightwood, 30 km from Crater Rock. -from Authors

  18. Levee crest elevation profiles derived from airborne lidar-based high resolution digital elevation models in south Louisiana

    USGS Publications Warehouse

    Palaseanu-Lovejoy, Monica; Thatcher, Cindy A.; Barras, John A.

    2014-01-01

    This study explores the feasibility of using airborne lidar surveys to derive high-resolution digital elevation models (DEMs) and develop an automated procedure to extract levee longitudinal elevation profiles for both federal levees in Atchafalaya Basin and local levees in Lafourche Parish. Generally, the use of traditional manual surveying methods to map levees is a costly and time consuming process that typically produces cross-levee profiles every few hundred meters, at best. The purpose of our paper is to describe and test methods for extracting levee crest elevations in an efficient, comprehensive manner using high resolution lidar generated DEMs. In addition, the vertical uncertainty in the elevation data and its effect on the resultant estimate of levee crest heights is addressed in an assessment of whether the federal levees in our study meet the USACE minimum height design criteria.

  19. ACE2 Global Digital Elevation Model : User Analysis

    NASA Astrophysics Data System (ADS)

    Smith, R. G.; Berry, P. A. M.; Benveniste, J.

    2013-12-01

    Altimeter Corrected Elevations 2 (ACE2), first released in October 2009, is the Global Digital Elevation Model (GDEM) created by fusing the high accuracy of over 100 million altimeter retracked height estimates, derived primarily from the ERS-1 Geodetic Mission, with the high frequency content available within the near-global Shuttle Radar Topography Mission. This novel ACE2 GDEM is freely available at 3”, 9”, 30” and 5' and has been distributed via the web to over 680 subscribers. This paper presents the results of a detailed analysis of geographical distribution of subscribed users, along with fields of study and potential uses. Investigations have also been performed to determine the most popular spatial resolutions and the impact these have on the scope of data downloaded. The analysis has shown that, even though the majority of users have come from Europe and America, a significant number of website hits have been received from South America, Africa and Asia. Registered users also vary widely, from research institutions and major companies down to individual hobbyists looking at data for single projects.

  20. Record of late holocene debris avalanches and lahars at Iliamna Volcano, Alaska

    USGS Publications Warehouse

    Waythomas, C.F.; Miller, T.P.; Beget, J.E.

    2000-01-01

    Iliamna Volcano is a 3053-meter high, glaciated stratovolcano in the southern Cook Inlet region of Alaska and is one of seven volcanoes in this region that have erupted multiple times during the past 10,000 yr. Prior to our studies of Iliamna Volcano, little was known about the frequency, magnitude, and character of Holocene volcanic activity. Here we present geologic evidence of the most recent eruptive activity of the volcano and provide the first outline of Late Holocene debris-avalanche and lahar formation. Iliamna has had no documented historical eruptions but our recent field investigations indicate that the volcano has erupted at least twice in the last 300 yr. Clay-rich lahar deposits dated by radiocarbon to ???1300 and ???90 yr BP are present in two major valleys that head on the volcano. These deposits indicate that at least two large, possibly deep-seated, flank failures of the volcanic edifice have occurred in the last 1300 yr. Noncohesive lahar deposits likely associated with explosive pyroclastic eruptions date to 2400-1300,>1500,???300, and <305 yr BP. Debris-avalanche deposits from recent and historical small-volume slope failures of the hydrothermally altered volcanic edifice cover most of the major glaciers on the volcano. Although these deposits consist almost entirely of hydrothermally altered rock debris and snow and ice, none of the recently generated debris avalanches evolved to lahars. A clay-rich lahar deposit that formed <90??60 radiocarbon yr BP and entered the Johnson River Valley southeast of the volcano cannot be confidently related to an eruption of Iliamna Volcano, which has had no known historical eruptions. This deposit may record an unheralded debris avalanche and lahar. ?? 2000 Elsevier Science B.V. All rights reserved.

  1. The October 16, 2013 rainfall-induced landslides and associated lahars at Izu Oshima Volcano, Japan

    NASA Astrophysics Data System (ADS)

    Miyabuchi, Yasuo; Maeno, Fukashi; Nakada, Setsuya

    2015-09-01

    Intense rainfall related to the typhoon T1326 on October 15-16, 2013 (total 824 mm; maximum hourly rainfall 118.5 mm) triggered numerous landslides and associated lahars at Izu Oshima Volcano, the northernmost part of Izu Mariana volcanic arc, Japan. The landslides were concentrated mainly in a 2-km2 area located on the western slope of Izu Oshima Volcano. Most of the landslides were shallow soil slips (< 2 m thick) in unconsolidated fallout tephra layers overlying lava flows and pyroclastic rocks. The rupture surfaces of them were located near the base of Y1 tephra (AD 1777-1778) and/or the base of Y4 tephra (AD 1421). The Y1 and Y4 tephras differ from the underlying paleosols in permeability, grain size and degree of compaction. The saturated hydraulic conductivities of the paleosols were one to two orders of magnitude smaller than those of the overlying Y1 and Y4 tephras. Most landslides mobilized completely into lahars, traveling along stream channels or flat slopes and flooding at the foot of the volcano. The associated lahars severely damaged inhabited areas and caused thirty five fatalities. Although the lahars eroded slopes and transported boulders up to 1 m in diameter and a large amount of woody debris, they contained more than 90% of sand-to-silt-size particles, similar in composition to the original sliding materials. Sediment discharge volumes from three basins were estimated at 1.8-4.1 × 104 m3/km2, based on debris volumes trapped by sediment retention dams. The characteristics of rainfall-induced landslides and associated lahars at Izu Oshima Volcano in 2013 provide an important lesson about future non-eruption-related landslide and lahar hazards at tephra-rich volcanoes.

  2. Elevated Temperature Testing and Modeling of Advanced Toughened Ceramic Materials

    NASA Technical Reports Server (NTRS)

    Keith, Theo G.

    2005-01-01

    The purpose of this report is to provide a final report for the period of 12/1/03 through 11/30/04 for NASA Cooperative Agreement NCC3-776, entitled "Elevated Temperature Testing and Modeling of Advanced Toughened Ceramic Materials." During this final period, major efforts were focused on both the determination of mechanical properties of advanced ceramic materials and the development of mechanical test methodologies under several different programs of the NASA-Glenn. The important research activities made during this period are: 1. Mechanical properties evaluation of two gas-turbine grade silicon nitrides. 2) Mechanical testing for fuel-cell seal materials. 3) Mechanical properties evaluation of thermal barrier coatings and CFCCs and 4) Foreign object damage (FOD) testing.

  3. Registratiom of TM data to digital elevation models

    NASA Technical Reports Server (NTRS)

    1984-01-01

    Several problems arise when attempting to register LANDSAT thematic mapper data to U.S. B Geological Survey digital elevation models (DEMs). The TM data are currently available only in a rotated variant of the Space Oblique Mercator (SOM) map projection. Geometric transforms are thus; required to access TM data in the geodetic coordinates used by the DEMs. Due to positional errors in the TM data, these transforms require some sort of external control. The spatial resolution of TM data exceeds that of the most commonly DEM data. Oversampling DEM data to TM resolution introduces systematic noise. Common terrain processing algorithms (e.g., close computation) compound this problem by acting as high-pass filters.

  4. Fast Ray Tracing of Lunar Digital Elevation Models

    NASA Technical Reports Server (NTRS)

    McClanahan, Timothy P.; Evans, L. G.; Starr, R. D.; Mitrofanov, I.

    2009-01-01

    Ray-tracing (RT) of Lunar Digital Elevation Models (DEM)'s is performed to virtually derive the degree of radiation incident to terrain as a function of time, orbital and ephemeris constraints [I- 4]. This process is an integral modeling process in lunar polar research and exploration due to the present paucity of terrain information at the poles and mission planning activities for the anticipated spring 2009 launch of the Lunar Reconnaissance Orbiter (LRO). As part of the Lunar Exploration Neutron Detector (LEND) and Lunar Crater Observation and Sensing Satellite (LCROSS) preparations RI methods are used to estimate the critical conditions presented by the combined effects of high latitude, terrain and the moons low obliquity [5-7]. These factors yield low incident solar illumination and subsequently extreme thermal, and radiation conditions. The presented research uses RT methods both for radiation transport modeling in space and regolith related research as well as to derive permanently shadowed regions (PSR)'s in high latitude topographic minima, e.g craters. These regions are of scientific and human exploration interest due to the near constant low temperatures in PSRs, inferred to be < 100 K. Hydrogen is thought to have accumulated in PSR's through the combined effects of periodic cometary bombardment and/or solar wind processes, and the extreme cold which minimizes hydrogen sublimation [8-9]. RT methods are also of use in surface position optimization for future illumination dependent on surface resources e.g. power and communications equipment.

  5. Advances in the data compression of digital elevation models

    NASA Astrophysics Data System (ADS)

    Kidner, David B.; Smith, Derek H.

    2003-10-01

    The maintenance and dissemination of spatial databases requires efficient strategies for handling the large volumes of data that are now publicly available. In particular, satellite and aerial imagery, radar, LiDAR, and digital elevation models (DEMs) are being utilised by a sizeable user-base, for predominantly environmental applications. The efficient dissemination of such datasets has become a key issue in the development of web-based and distributed computing environments. However, the physical size of these datasets is a major bottleneck in their storage and transmission. The problem is often exaggerated when the data is supplied in less efficient, proprietary or national data formats. This paper presents a methodology for the lossless compression of DEMs, based on the statistical correlation of terrain data in local neighbourhoods. Most data and image compression algorithms fail to capitalise fully on the inherent redundancy in spatial data. At the same time, users often prefer a uniform solution to all their data compression requirements, but these solutions may be far from optimal. The approach presented here can be thought of as a simple pre-processing of the elevation data before the use of traditional data compression software frequently applied to spatial data sets, such as GZIP. Identification and removal of the spatial redundancy in terrain data, with the use of optimal predictors for DEMs and optimal statistical encoders such as Arithmetic Coding, gives even higher compression ratios. Both GZIP and our earlier approach of combining a simple linear prediction algorithm with Huffman Coding are shown to be far from optimal in identifying and removing the spatial redundancy in DEMs. The new approaches presented here typically halve the file sizes of our earlier approach, and give a 40-62% improvement on GZIP-compressed DEMs.

  6. Modeling Ka-band low elevation angle propagation statistics

    NASA Technical Reports Server (NTRS)

    Russell, Thomas A.; Weinfield, John; Pearson, Chris; Ippolito, Louis J.

    1995-01-01

    The statistical variability of the secondary atmospheric propagation effects on satellite communications cannot be ignored at frequencies of 20 GHz or higher, particularly if the propagation margin allocation is such that link availability falls below 99 percent. The secondary effects considered in this paper are gaseous absorption, cloud absorption, and tropospheric scintillation; rain attenuation is the primary effect. Techniques and example results are presented for estimation of the overall combined impact of the atmosphere on satellite communications reliability. Statistical methods are employed throughout and the most widely accepted models for the individual effects are used wherever possible. The degree of correlation between the effects is addressed and some bounds on the expected variability in the combined effects statistics are derived from the expected variability in correlation. Example estimates are presented of combined effects statistics in the Washington D.C. area of 20 GHz and 5 deg elevation angle. The statistics of water vapor are shown to be sufficient for estimation of the statistics of gaseous absorption at 20 GHz. A computer model based on monthly surface weather is described and tested. Significant improvement in prediction of absorption extremes is demonstrated with the use of path weather data instead of surface data.

  7. Preliminary observations of voluminous ice-rich and water-rich lahars generated during the 2009 eruption of Redoubt, Alaska

    USGS Publications Warehouse

    Waythomas, Christopher F.; Pierson, Thomas C.; Major, Jon J.; Scott, William E.

    2012-01-01

    Redoubt Volcano in south-central Alaska began erupting on March 15, 2009, and by April 4, 2009, had produced at least 20 explosive events that generated plumes of ash and lahars. The 3,108-m high, snow- and -ice-clad stratovolcano has an ice-filled summit crater that is breached to the north. The volcano supports about 4 km3 of ice and snow and about 1 km3 of this makes up the Drift glacier on the northern side of the volcano. Explosive eruptions between March 22 and April 4, which included the destruction of at least two lava domes, triggered significant lahars in the Drift River valley on March 23 and April 4 and several smaller lahars between March 24 and March 31. High-flow marks, character of deposits, areas of inundation, and estimates of flow velocity revealed that the lahars on March 23 and April 4 were the largest of the eruption. In the 2-km-wide upper Drift River valley, average flow depths were about 3–5 m. Average peak-flow velocities were likely between 10 and 15 ms-1, and peak discharges were on the order of 104–105 m3s-1. The area inundated by lahars on March 23 was at least 100 km2 and on April 4 about 125 km2. The lahars emplaced on March 23 and April 4 had volumes on the order of 107–108 m3 and were similar in size to the largest lahar of the 1989–90 eruption. The March 23 lahars were primarily flowing slurries of snow and ice entrained from the Drift glacier and seasonal snow and tabular blocks of river ice from the Drift River valley. Only a single, undifferentiated deposit up to 5 m thick was found and contained about 80–95 percent of poorly sorted, massive to imbricate assemblages of snow and ice. The deposit was frozen soon after it was emplaced and later eroded and buried by the April 4 lahar. The lahar of April 4, in contrast, was primarily a hyperconcentrated flow, as interpreted from 1- to 6-m thick deposits of massive to horizontally stratified sand-to-fine-gravel. Rock material in the April 4 lahar deposit is predominantly

  8. Incremental terrain processing for large digital elevation models

    NASA Astrophysics Data System (ADS)

    Ye, Z.

    2012-12-01

    Incremental terrain processing for large digital elevation models Zichuan Ye, Dean Djokic, Lori Armstrong Esri, 380 New York Street, Redlands, CA 92373, USA (E-mail: zye@esri.com, ddjokic@esri.com , larmstrong@esri.com) Efficient analyses of large digital elevation models (DEM) require generation of additional DEM artifacts such as flow direction, flow accumulation and other DEM derivatives. When the DEMs to analyze have a large number of grid cells (usually > 1,000,000,000) the generation of these DEM derivatives is either impractical (it takes too long) or impossible (software is incapable of processing such a large number of cells). Different strategies and algorithms can be put in place to alleviate this situation. This paper describes an approach where the overall DEM is partitioned in smaller processing units that can be efficiently processed. The processed DEM derivatives for each partition can then be either mosaicked back into a single large entity or managed on partition level. For dendritic terrain morphologies, the way in which partitions are to be derived and the order in which they are to be processed depend on the river and catchment patterns. These patterns are not available until flow pattern of the whole region is created, which in turn cannot be established upfront due to the size issues. This paper describes a procedure that solves this problem: (1) Resample the original large DEM grid so that the total number of cells is reduced to a level for which the drainage pattern can be established. (2) Run standard terrain preprocessing operations on the resampled DEM to generate the river and catchment system. (3) Define the processing units and their processing order based on the river and catchment system created in step (2). (4) Based on the processing order, apply the analysis, i.e., flow accumulation operation to each of the processing units, at the full resolution DEM. (5) As each processing unit is processed based on the processing order defined

  9. Developing building-damage scales for lahars: application to Merapi volcano, Indonesia

    NASA Astrophysics Data System (ADS)

    Jenkins, Susanna F.; Phillips, Jeremy C.; Price, Rebecca; Feloy, Kate; Baxter, Peter J.; Hadmoko, Danang Sri; de Bélizal, Edouard

    2015-09-01

    Lahar damage to buildings can include burial by sediment and/or failure of walls, infiltration into the building and subsequent damage to contents. The extent to which a building is damaged will be dictated by the dynamic characteristics of the lahar, i.e. the velocity, depth, sediment concentration and grain size, as well as the structural characteristics and setting of the building in question. The focus of this paper is on quantifying how buildings may respond to impact by lahar. We consider the potential for lahar damage to buildings on Merapi volcano, Indonesia, as a result of the voluminous deposits produced during the large (VEI 4) eruption in 2010. A building-damage scale has been developed that categorises likely lahar damage levels and, through theoretical calculations of expected building resistance to impact, approximate ranges of impact pressures. We found that most weak masonry buildings on Merapi would be destroyed by dilute lahars with relatively low velocities (ca. 3 m/s) and pressures (ca. 5 kPa); however, the majority of stronger rubble stone buildings may be expected to withstand higher velocities (to 6 m/s) and pressures (to 20 kPa). We applied this preliminary damage scale to a large lahar in the Putih River on 9 January 2011, which inundated and caused extensive building damage in the village of Gempol, 16 km southwest of Merapi. The scale was applied remotely through the use of public satellite images and through field studies to categorise damage and estimate impact pressures and velocities within the village. Results were compared with those calculated independently from Manning's calculations for flow velocity and depth within Gempol village using an estimate of flow velocity at one upstream site as input. The results of this calculation showed reasonable agreement with an average channel velocity derived from travel time observations. The calculated distribution of flow velocities across the area of damaged buildings was consistent with

  10. High-resolution digital elevation model of lower Cowlitz and Toutle Rivers, adjacent to Mount St. Helens, Washington, based on an airborne lidar survey of October 2007

    USGS Publications Warehouse

    Mosbrucker, Adam

    2015-01-01

    The lateral blast, debris avalanche, and lahars of the May 18th, 1980, eruption of Mount St. Helens, Washington, dramatically altered the surrounding landscape. Lava domes were extruded during the subsequent eruptive periods of 1980–1986 and 2004–2008. More than three decades after the emplacement of the 1980 debris avalanche, high sediment production persists in the Toutle River basin, which drains the northern and western flanks of the volcano. Because this sediment increases the risk of flooding to downstream communities on the Toutle and lower Cowlitz Rivers, the U.S. Army Corps of Engineers (USACE), under the direction of Congress to maintain an authorized level of flood protection, continues to monitor and mitigate excess sediment in North and South Fork Toutle River basins to help reduce this risk and to prevent sediment from clogging the shipping channel of the Columbia River. From October 22–27, 2007, Watershed Sciences, Inc., under contract to USACE, collected high-precision airborne lidar (light detection and ranging) data that cover 273 square kilometers (105 square miles) of lower Cowlitz and Toutle River tributaries from the Columbia River at Kelso, Washington, to upper North Fork Toutle River (below the volcano's edifice), including lower South Fork Toutle River. These data provide a digital dataset of the ground surface, including beneath forest cover. Such remotely sensed data can be used to develop sediment budgets and models of sediment erosion, transport, and deposition. The U.S. Geological Survey (USGS) used these lidar data to develop digital elevation models (DEMs) of the study area. DEMs are fundamental to monitoring natural hazards and studying volcanic landforms, fluvial and glacial geomorphology, and surface geology. Watershed Sciences, Inc., provided files in the LASer (LAS) format containing laser returns that had been filtered, classified, and georeferenced. The USGS produced a hydro-flattened DEM from ground-classified points at

  11. High-resolution digital elevation model of Mount St. Helens crater and upper North Fork Toutle River basin, Washington, based on an airborne lidar survey of September 2009

    USGS Publications Warehouse

    Mosbrucker, Adam

    2014-01-01

    The lateral blast, debris avalanche, and lahars of the May 18th, 1980, eruption of Mount St. Helens, Washington, dramatically altered the surrounding landscape. Lava domes were extruded during the subsequent eruptive periods of 1980–1986 and 2004–2008. More than three decades after the emplacement of the 1980 debris avalanche, high sediment production persists in the North Fork Toutle River basin, which drains the northern flank of the volcano. Because this sediment increases the risk of flooding to downstream communities on the Toutle and Cowlitz Rivers, the U.S. Army Corps of Engineers (USACE), under the direction of Congress to maintain an authorized level of flood protection, built a sediment retention structure on the North Fork Toutle River in 1989 to help reduce this risk and to prevent sediment from clogging the shipping channel of the Columbia River. From September 16–20, 2009, Watershed Sciences, Inc., under contract to USACE, collected high-precision airborne lidar (light detection and ranging) data that cover 214 square kilometers (83 square miles) of Mount St. Helens and the upper North Fork Toutle River basin from the sediment retention structure to the volcano's crater. These data provide a digital dataset of the ground surface, including beneath forest cover. Such remotely sensed data can be used to develop sediment budgets and models of sediment erosion, transport, and deposition. The U.S. Geological Survey (USGS) used these lidar data to develop digital elevation models (DEMs) of the study area. DEMs are fundamental to monitoring natural hazards and studying volcanic landforms, fluvial and glacial geomorphology, and surface geology. Watershed Sciences, Inc., provided files in the LASer (LAS) format containing laser returns that had been filtered, classified, and georeferenced. The USGS produced a hydro-flattened DEM from ground-classified points at Castle, Coldwater, and Spirit Lakes. Final results averaged about five laser last

  12. Budapest, Hungary, Perspective View, SRTM Elevation Model with Landsat Overlay

    NASA Technical Reports Server (NTRS)

    2004-01-01

    After draining the northern flank of the Alps Mountains in Germany and Austria, the Danube River flows east as it enters this west-looking scene (upper right) and forms the border between Slovakia and Hungary. The river then leaves the border as it enters Hungary and transects the Transdanubian Mountains, which trend southwest to northeast. Upon exiting the mountains, the river turns southward, flowing past Budapest (purplish blue area) and along the western margin of the Great Hungarian Plain.

    South and west of the Danube, the Transdanubian Mountains have at most only about 400 meters (about 1300 feet) of relief but they exhibit varied landforms, which include volcanic, tectonic, fluvial (river), and eolian (wind) features. A thick deposit of loess (dust deposits likely blown from ancient glacial outwash) covers much of this area, and winds from the northwest, funneled between the Alps and the Carpathian Mountains, are apparently responsible for a radial pattern of erosional streaks across the entire region.

    This image was generated from a Landsat satellite image draped over an elevation model produced by the Shuttle Radar Topography Mission (SRTM). The view uses a 3-times vertical exaggeration to enhance topographic expression. The false colors of the scene result from displaying Landsat bands 1, 4, and 7 in blue, green, and red, respectively. Band 1 is visible blue light, but bands 4 and 7 are reflected infrared light. This band combination maximizes color contrasts between the major land cover types, namely vegetation (green), bare ground (red), and water (blue). Shading of the elevation model was used to further highlight the topographic features.

    Elevation data used in this image was acquired by the Shuttle Radar Topography Mission (SRTM) aboard the Space Shuttle Endeavour, launched on February 11, 2000. SRTM used the same radar instrument that comprised the Spaceborne Imaging Radar-C/X-Band Synthetic Aperture Radar (SIR-C/X-SAR) that flew twice on

  13. Identifying pyroclastic and lahar deposits and assessing erosion and lahar hazards at active volcanoes using multi-temporal HSR image analysis and techniques for change detection

    NASA Astrophysics Data System (ADS)

    Kassouk, Zeineb; Thouret, Jean-Claude; Oehler, Jean-François; Solikhin, Akhmad

    2014-05-01

    The increasing availability of high-spatial resolution (HSR) remote sensing images leads to new opportunities for hazard assessment in the case of active volcanoes. Object-oriented analysis (OOA) of HSR images helps to simultaneously exploit spatial, spectral and contextual information. Here, we identify and delineate pyroclastic density current (PDC) and post-eruption lahar deposits on the south flank of Merapi volcano, Indonesia, after the large 2010 eruption. GeoEye-1 (2010 and 2011) and Pleiades (2012) images were analyzed with an adjusted object-oriented method. The PDC deposits include valley-confined block-and-ash flows (BAFs), unconfined, overbank pyroclastic flows (OPFs), and high-energy surges or ash-cloud surges. We follow up the evolution of the pyroclastic and lahar deposits through changes in the spectral indices calculated in segmented features, which represent the principal units of deposits and devastated areas. The object-oriented analysis has been applied to the pseudo image comprising of three spectral indices (NDWI water index; NDVI vegetation index; and NDRSI Red Soil Index). This pseudo image has enabled us to delineate fifteen units of PDC and lahar deposits, and damaged forests and settlements in the Gendol-Opak catchment (c.80 sqkm). The units represent 75% of classes obtained by photointerpretation of the same image and supported by field observations. A combination of NDWI and NDVI helps to separate areas affected by surges (NDWI <0.2 and 0.1lahar deposits (NDRSI> 0.3 and NDWI<0.1). NDRSI values close to 0 are assigned to scoria-rich PFs darker than other PF deposits. Bivariate analyses of three spectral indices, NDWI, NDVI and NDRSI, show the temporal evolution of the delineated deposits and areas between 2010 and 2012. The NDVI/NDWI 2010 plot shows two clusters: NDVI and NDWI close to 0

  14. Object representations at multiple scales from digital elevation models.

    PubMed

    Drăguţ, Lucian; Eisank, Clemens

    2011-06-15

    In the last decade landform classification and mapping has developed as one of the most active areas of geomorphometry. However, translation from continuous models of elevation and its derivatives (slope, aspect, and curvatures) to landform divisions (landforms and landform elements) is filtered by two important concepts: scale and object ontology. Although acknowledged as being important, these two issues have received surprisingly little attention.This contribution provides an overview and prospects of object representation from DEMs as a function of scale. Relationships between object delineation and classification or regionalization are explored, in the context of differences between general and specific geomorphometry. A review of scales issues in geomorphometry-ranging from scale effects to scale optimization techniques-is followed by an analysis of pros and cons of using cells and objects in DEM analysis. Prospects for coupling multi-scale analysis and object delineation are then discussed. Within this context, we propose discrete geomorphometry as a possible approach between general and specific geomorphometry. Discrete geomorphometry would apply to and describe land-surface divisions defined solely by the criteria of homogeneity in respect to a given land-surface parameter or a combination of several parameters. Homogeneity, in its turn, should always be relative to scale. PMID:21760655

  15. Object representations at multiple scales from digital elevation models

    PubMed Central

    Drăguţ, Lucian; Eisank, Clemens

    2011-01-01

    In the last decade landform classification and mapping has developed as one of the most active areas of geomorphometry. However, translation from continuous models of elevation and its derivatives (slope, aspect, and curvatures) to landform divisions (landforms and landform elements) is filtered by two important concepts: scale and object ontology. Although acknowledged as being important, these two issues have received surprisingly little attention. This contribution provides an overview and prospects of object representation from DEMs as a function of scale. Relationships between object delineation and classification or regionalization are explored, in the context of differences between general and specific geomorphometry. A review of scales issues in geomorphometry—ranging from scale effects to scale optimization techniques—is followed by an analysis of pros and cons of using cells and objects in DEM analysis. Prospects for coupling multi-scale analysis and object delineation are then discussed. Within this context, we propose discrete geomorphometry as a possible approach between general and specific geomorphometry. Discrete geomorphometry would apply to and describe land-surface divisions defined solely by the criteria of homogeneity in respect to a given land-surface parameter or a combination of several parameters. Homogeneity, in its turn, should always be relative to scale. PMID:21760655

  16. Dune Morphometry in the Age of Digital Elevation Models

    NASA Astrophysics Data System (ADS)

    Lancaster, N.

    2014-12-01

    Dune patterns can be characterized in many different ways. Relationships between dune height, width and spacing, and the spatial variation in these parameters have been widely employed to provide quantitative information that can be used to describe dune patterns and make comparisons between dunes in widely separated areas, as well as to identify different generations of dunes. Digital elevation models (e.g. ASTER GDEM) provide a rich resource of data for analyses of dune patterns at landscape scales in several ways, including: (1) more extensive analyses using traditional measures, such as relationships between dune height and spacing, and the spatial variation in these parameters; and (2) estimation of sediment thickness on a regional scale. Analyses of data for Arabian and Namibian sand seas and dune fields show that dune height and spacing relationships are much more variable than previously reported and call into question existing models. Regional patterns of sediment thickness reveal areas of erosion, bypass, and accumulation that can be related to transport pathways and wind regimes. The widespread occurrence of complex dune patterns as well as the magnitude of the newly available data sets however requires more sophisticated analyses than simple extraction of dune morphometric parameters using GIS approaches. Geostatistical analyses using spatial autocorrelation, Fourier, and Wavelet methods have been employed in analyses of sub-aqueous bedforms and show promise for dune systems. Automated or semi-automated identification of dune length, width, spacing, and trends using advanced image analysis techniques such as linear segment detection is a potentially transformative approach. The strengths and weaknesses of these methods to provide pertinent geomorphic information are currently being evaluated, but they have the potential to provide new insights into the nature of dune patterns.

  17. Online Mars Digital Elevation Model Derived from Profiles

    NASA Astrophysics Data System (ADS)

    Delacourt, C.; Gros, N.; Allemand, P.; Baratoux, D.

    The topography of Mars is a key parameter for understanding the geological evolution of the planet. Since 1997, the Mars Orbital Laser Altimeter (MOLA), launched in the frame of Mars Global Surveyor, has acquired more than 600 million topographic measurements. However, despite the high vertical accuracy of those profiles, the main limitation of this technique appears when topographic maps are required. To create a Digital Elevation Model (DEM) or a topographic map, an interpolation on individual MOLA measurements on regular grids is required. Calculating the global full-resolution Martian DEM requires very intensive computation and large disk capacities. Only a few teams throughout the world have computed a full resolution DEM from MOLA data. Even if a scientist is interested in a small area of Mars, numerous profiles have to be processed and extracted from 44 CD-ROMs. To facilitate the exploitation of the high potential of MOLA data, we propose an Internet application that allows any user to extract the individual MOLA measurements from all profiles over a given area and to process local DEMs of the surface of Mars with adjustable parameters of computation. This facility, developed in Interactive Data Language by Research Systems, Inc., allows users to select the zone of interest and the resolution of the output DEM. After online processing, various products in standard formats can be downloaded on the user's computer, including DEMs, individual MOLA points, list and map of the MOLA ground tracks used for the DEM generation, and a quality map. This map is computed by evaluating the distance between each point of the DEM and the closest measurements of the altimeter. Furthermore, IDL tools are supplied to facilitate data integration and use.

  18. The ASTER Global Digital Elevation Model (GDEM) -for societal benefit -

    NASA Astrophysics Data System (ADS)

    Hato, M.; Tsu, H.; Tachikawa, T.; Abrams, M.; Bailey, B.

    2009-12-01

    The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model (GDEM) was developed jointly by the Ministry of Economy, Trade and Industry (METI) of Japan and the United States National Aeronautics and Space Administration (NASA) under the agreement of contribution to GEOSS and a public release was started on June 29th. ASTER GDEM can be downloaded to users from the Earth Remote Sensing Data Analysis Center (ERSDAC) of Japan and NASA’s Land Processes Distributed Active Archive Center (LP DAAC) free of charge. The ASTER instrument was built by METI and launched onboard NASA’s Terra spacecraft in December 1999. It has an along-track stereoscopic capability using its near infrared spectral band (NIR) and its nadir-viewing and backward-viewing telescopes to acquire stereo image data with a base-to-height ratio of 0.6. The ASTER GDEM was produced by applying newly-developed automated algorithm to more than 1.2 million NIR data Produced DEMs of all scene data was stacked after cloud masking and finally partitioned into 1° x 1°unit (called ‘tile’) data for convenience of distribution and handling by users. Before start of public distribution, ERSDAC and USGS/NASA together with many volunteers did validation and characterization by using a preliminary product of the ASTER GDEM. As a result of validation, METI and NASA evaluated that Version 1 of the ASTER GDEM has enough quality to be used as “experimental” or “research grade” data and consequently decided to release it. The ASTER GDEM covering almost all land area of from 83N to 83S on the earth represents as an important contribution to the global earth observation community. We will show our effort of development of ASTER GDEM and its accuracy and character.

  19. An evaluation of onshore digital elevation models for modelling tsunami inundation zones

    NASA Astrophysics Data System (ADS)

    Griffin, Jonathan; Latief, Hamzah; Kongko, Widjo; Harig, Sven; Horspool, Nick; Hanung, Raditya; Rojali, Aditia; Maher, Nicola; Fuchs, Annika; Hossen, Jakir; Upi, Supriyati; Edi, Dewanto; Rakowsky, Natalja; Cummins, Phil

    2015-06-01

    A sensitivity study is undertaken to assess the utility of different onshore digital elevation models (DEM) for simulating the extent of tsunami inundation using case studies from two locations in Indonesia. We compare airborne IFSAR, ASTER and SRTM against high resolution LiDAR and stereo-camera data in locations with different coastal morphologies. Tsunami inundation extents modelled with airborne IFSAR DEMs are comparable with those modelled with the higher resolution datasets and are also consistent with historical run-up data, where available. Large vertical errors and poor resolution of the coastline in the ASTER and SRTM elevation datasets cause the modelled inundation extent to be much less compared with the other datasets and observations. Therefore ASTER and SRTM should not be used to underpin tsunami inundation models. a model mesh resolution of 25 m was sufficient for estimating the inundated area when using elevation data with high vertical accuracy in the case studies presented here. Differences in modelled inundation between digital terrain models (DTM) and digital surface models (DSM) for LiDAR and IFSAR are greater than differences between the two data types. Models using DTM may overestimate inundation while those using DSM may underestimate inundation when a constant Manning’s roughness value is used. We recommend using DTM for modelling tsunami inundation extent with further work needed to resolve the scale at which surface roughness should be parameterised.

  20. Quality assessment of Digital Elevation Model (DEM) in view of the Altiplano hydrological modeling

    NASA Astrophysics Data System (ADS)

    Satgé, F.; Arsen, A.; Bonnet, M.; Timouk, F.; Calmant, S.; Pilco, R.; Molina, J.; Lavado, W.; Crétaux, J.; HASM

    2013-05-01

    Topography is crucial data input for hydrological modeling but in many regions of the world, the only way to characterize topography is the use of satellite-based Digital Elevation Models (DEM). In some regions, the quality of these DEMs remains poor and induces modeling errors that may or not be compensated by model parameters tuning. In such regions, the evaluation of these data uncertainties is an important step in the modeling procedure. In this study, which focuses on the Altiplano region, we present the evaluation of the two freely available DEM. The shuttle radar topographic mission (SRTM), a product of the National Aeronautics and Space Administration (NASA) and the Advanced Space Born Thermal Emission and Reflection Global Digital Elevation Map (ASTER GDEM), data provided by the Ministry of Economy, Trade and Industry of Japan (MESI) in collaboration with the NASA, are widely used. While the first represents a resolution of 3 arc seconds (90m) the latter is 1 arc second (30m). In order to select the most reliable DEM, we compared the DEM elevation with high qualities control points elevation. Because of its large spatial coverture (track spaced of 30 km with a measure of each 172 m) and its high vertical accuracy which is less than 15 cm in good weather conditions, the Geoscience Laser Altimeter System (GLAS) on board on the Ice, Cloud and Land elevation Satellite of NASA (ICESat) represent the better solution to establish a high quality elevation database. After a quality check, more than 150 000 ICESat/GLAS measurements are suitable in terms of accuracy for the Altiplano watershed. This data base has been used to evaluate the vertical accuracy for each DEM. Regarding to the full spatial coverture; the comparison has been done for both, all kind of land coverture, range altitude and mean slope.

  1. High Resolution Digital Elevation Models of Pristine Explosion Craters

    NASA Technical Reports Server (NTRS)

    Farr, T. G.; Krabill, W.; Garvin, J. B.

    2004-01-01

    In order to effectively capture a realistic terrain applicable to studies of cratering processes and landing hazards on Mars, we have obtained high resolution digital elevation models of several pristine explosion craters at the Nevada Test Site. We used the Airborne Terrain Mapper (ATM), operated by NASA's Wallops Flight Facility to obtain DEMs with 1 m spacing and 10 cm vertical errors of 4 main craters and many other craters and collapse pits. The main craters that were mapped are Sedan, Scooter, Schooner, and Danny Boy. The 370 m diameter Sedan crater, located on Yucca Flat, is the largest and freshest explosion crater on Earth that was formed under conditions similar to hypervelocity impact cratering. As such, it is effectively pristine, having been formed in 1962 as a result of a controlled detonation of a 100 kiloton thermonuclear device, buried at the appropriate equivalent depth of burst required to make a simple crater. Sedan was formed in alluvium of mixed lithology and subsequently studied using a variety of field-based methods. Nearby secondary craters were also formed at the time and were also mapped by ATM. Adjacent to Sedan and also in alluvium is Scooter, about 90 m in diameter and formed by a high-explosive event. Schooner (240 m) and Danny Boy (80 m) craters were also important targets for ATM as they were excavated in hard basalt and therefore have much rougher ejecta. This will allow study of ejecta patterns in hard rock as well as engineering tests of crater and rock avoidance and rover trafficability. In addition to the high resolution DEMs, crater geometric characteristics, RMS roughness maps, and other higher-order derived data products will be generated using these data. These will provide constraints for models of landing hazards on Mars and for rover trafficability. Other planned studies will include ejecta size-frequency distribution at the resolution of the DEM and at finer resolution through air photography and field measurements

  2. Orthographic terrain views using data derived from digital elevation models

    NASA Technical Reports Server (NTRS)

    Dubayah, R. O.; Dozier, J.

    1986-01-01

    A fast algorithm for producing three-dimensional orthographic terrain views uses digital elevation data and co-registered imagery. These views are created using projective geometry and are designed for display on high-resolution raster graphics devices. The algorithm's effectiveness is achieved by (1) the implementation of two efficient gray-level interpolation routines that offer the user a choice between speed and smoothness, and (2) a unique visible surface determination procedure based on horizon angles derived from the elevation data set.

  3. Seismic characterisation of lahars at Volcán de Colima, Mexico

    NASA Astrophysics Data System (ADS)

    Vázquez, R.; Suriñach, E.; Capra, L.; Arámbula-Mendoza, R.; Reyes-Dávila, G.

    2016-02-01

    Volcán de Colima is one of the most active volcanoes in Mexico, not only for its eruptive history, but also for its annual occurrence of lahars. This makes the volcano a natural laboratory for monitoring and studying lahar processes. Since 2011, monitoring instruments have been deployed along the highly active Montegrande ravine, with at least three lahar events per year. Here, we report the datasets of three events collected during the 2012 and 2013 seasons, then interpret the acquired data. An event classification scheme based on lahar magnitude, duration and seismic characteristics is presented to distinguish "single-pulse" events (SPEs) from "multi-pulse" events (MPE). SPEs lasted approximately 60 min, had average velocities of ~2 m/s and mean peak discharges of ~24 m3/s. MPEs endured for more than 3 h, reached mean velocities of ~4.5 m/s and peak discharges of ~60 m3/s (for block-rich surges). The seismic signal-analysis also allowed us to discriminate physical flow fluctuations within single lahars, i.e. between the arrival of block-rich fronts and subsequent variations in flow discharge. The exponential regression analysis showed a best fit, with correlation coefficients around 0.92 and exponential coefficients of ~0.01 s, for the block-rich fronts, with seismic amplitudes increasing from 4.8 × 10-4 to 2.3 × 10-3 m/s and frequency ranges from 10 to 20 Hz. The variations in flow discharge were distinguished by lower amplitudes of ~5.7 × 10-4 m/s than those of the block-rich fronts and with frequency ranges of 10-40 Hz. The results presented in this paper demonstrate that the seismic data of events allowed us to describe and discriminate among different flow types; these records are thus a useful investigation tool for lahar events that have a seismic record but are not observed directly. We propose that a seismic early warning system can be developed to help civil protection authorities in designing risk mitigation strategies.

  4. Mitigation of hazards from future lahars from Mount Merapi in the Krasak River channel near Yogyakarta, central Java

    USGS Publications Warehouse

    Ege, John R.; Sutikno

    1983-01-01

    Procedures for reducing hazards from future lahars and debris flows in the Krasak River channel near Yogyakarta, Central Java, Indonesia, include (1) determining the history of the location, size, and effects of previous lahars and debris flows, and (2) decreasing flow velocities. The first may be accomplished by geologic field mapping along with acquiring information by interviewing local residents, and the second by increasing the cross sectional area of the river channel and constructing barriers in the flow path.

  5. Comparative lahar hazard mapping at Volcan Citlaltépetl, Mexico using SRTM, ASTER and DTED-1 digital topographic data

    USGS Publications Warehouse

    Hubbard, Bernard E.; Sheridan, Michael F.; Carrasco-Nunez, Gerardo; Diaz-Castellon, Rodolfo; Rodriguez, Sergio R.

    2007-01-01

    Finally, ASTERs 60 km swath width and 8% duty cycle presents a challenge for mapping lahar inundation hazards at E–W oriented stream valleys in low-latitude areas with persistent cloud cover. However, its continued operations enhances its utility as a means for updating the continuous but one-time coverage of SRTM, and for filling voids in the SRTM dataset such as those that occur along steep-sided valleys prone to hazards from future lahars.

  6. Monitoring the Dynamic of a Fluvial Channel after Lahar Disturbance: Huiloac Gorge (Popocatepetl Volcano, Mexico)

    NASA Astrophysics Data System (ADS)

    Andres, N.; Palacios, D.; Zamorano, J. J.; Tanarro, L. M.; Renschler, C.; Sanjosé, J. J.; Atkinson, A.

    2009-04-01

    Volcanic eruptions generate disturbances that affect hydrological systems (Major, 2003) by depositing large volumes of sediments in watersheds that exceed amounts common to non-volcanic river systems (Montgomery, 2005). If the eruption releases abundant melt water, the river system may respond immediately by forming hazardous flows called lahars. River system recovery following eruptive and laharic impact is an important process, but it has received little attention (Gran and Montgomery, 2005) despite the fact that Major et al. (2000) and Hayes et al. (2002) have shown that these disruptions cause long term instability and their effects persist for decades. Lahar deposits resulting from interaction between volcanic activity and the glacier located above the Huiloac Gorge on the northern slope of Popocatepetl volcano (19°02´ N, 98°62´ W, 5,424 m), have infilled the gorge (Palacios, 1995; Palacios et al., 1998 and 2001; Capra et al., 2004; Muñoz, 2007). All of the major lahars that occurred on the volcano in 1995 (4 km), 1997 (21 km), and 2001 (14 km) have channelled through Huiloac Gorge, and have dramatically altered its morphology and dynamics through erosion and deposition. The present study traces these changes in the aftermath of the laharic events that occurred from 1997-2001. A sector of the channel, located at 3200m-3240m altitude, of 500 m long and 15 to 20 m wide, in the mid-section of the gorge, was chosen as the control site. Precipitation is heaviest there and is most apt to trigger secondary post-eruptive lahars. ArcGis software was used to draw 6 geomorphic maps of the site showing spatial variations in the landforms for the period February 2002 - February 2008. In addition, 29 cross-profiles were made of the gorge for the same time interval, excluding February 2004. The volume of sediment eroded and deposited was calculated for each date by comparing variations in the height of the floor and banks of the gorge depicted in the cross-profile, and

  7. Uncertainty Analysis of LROC NAC Derived Elevation Models

    NASA Astrophysics Data System (ADS)

    Burns, K.; Yates, D. G.; Speyerer, E.; Robinson, M. S.

    2012-12-01

    One of the primary objectives of the Lunar Reconnaissance Orbiter Camera (LROC) [1] is to gather stereo observations with the Narrow Angle Camera (NAC) to generate digital elevation models (DEMs). From an altitude of 50 km, the NAC acquires images with a pixel scale of 0.5 meters, and a dual NAC observation covers approximately 5 km cross-track by 25 km down-track. This low altitude was common from September 2009 to December 2011. Images acquired during the commissioning phase and those acquired from the fixed orbit (after 11 December 2011) have pixel scales that range from 0.35 meters at the south pole to 2 meters at the north pole. Alimetric observations obtained by the Lunar Orbiter Laser Altimeter (LOLA) provide measurements of ±0.1 m between the spacecraft and the surface [2]. However, uncertainties in the spacecraft positioning can result in offsets (±20m) between altimeter tracks over many orbits. The LROC team is currently developing a tool to automatically register alimetric observations to NAC DEMs [3]. Using a generalized pattern search (GPS) algorithm, the new automatic registration adjusts the spacecraft position and pointing information during times when NAC images, as well as LOLA measurements, of the same region are acquired to provide an absolute reference frame for the DEM. This information is then imported into SOCET SET to aide in creating controlled NAC DEMs. For every DEM, a figure of merit (FOM) map is generated using SOCET SET software. This is a valuable tool for determining the relative accuracy of a specific pixel in a DEM. Each pixel in a FOM map is given a value to determine its "quality" by determining if the specific pixel was shadowed, saturated, suspicious, interpolated/extrapolated, or successfully correlated. The overall quality of a NAC DEM is a function of both the absolute and relative accuracies. LOLA altimetry provides the most accurate absolute geodetic reference frame with which the NAC DEMs can be compared. Offsets

  8. Eruption-triggered avalanche, flood, and lahar at mount st. Helens--effects of winter snowpack.

    PubMed

    Waitt, R B; Pierson, T C; Macleod, N S; Janda, R J; Voight, B; Holcomb, R T

    1983-09-30

    An explosive eruption of Mount St. Helens on 19 March 1982 had substantial impact beyond the vent because hot eruption products interacted with a thick snowpack. A blast of hot pumice, dome rocks, and gas dislodged crater-wall snow that avalanched through the crater and down the north flank. Snow in the crater swiftly melted to form a transient lake, from which a destructive flood and lahar swept down the north flank and the North Fork Toutle River. PMID:17759014

  9. Eruption-triggered avalanche, flood, and lahar at Mount St. Helens - Effects of winter snowpack

    USGS Publications Warehouse

    Waitt, R.B., Jr.; Pierson, T.C.; MacLeod, N.S.; Janda, R.J.; Voight, B.; Holcomb, R.T.

    1983-01-01

    An explosive eruption of Mount St. Helens on 19 March 1982 had substantial impact beyond the vent because hot eruption products interacted with a thick snowpack. A blast of hot pumice, dome rocks, and gas dislodged crater-wall snow that avalanched through the crater and down the north flank. Snow in the crater swiftly melted to form a transient lake, from which a destructive flood and lahar swept down the north flank and the North Fork Toutle River.

  10. Unusual Ice-Rich Lahars From the 23 September 1995 Eruption of Mt. Ruapehu, New Zealand

    NASA Astrophysics Data System (ADS)

    Fagents, S. A.; Houghton, B. F.

    2007-05-01

    Interactions between explosive eruptions and summit ice or snowpack are among the most hazardous lahar- forming events at snow-capped volcanoes, as demonstrated by the 1985 disaster at Nevado del Ruiz, Colombia. A case of contrasting impact is the 23 September 1995 eruption at Mt. Ruapehu, New Zealand, in which an ice- rich lahar traveled <3 km before coming to rest on relatively steep slopes (10-16°). Addressing the diverse thermal and fluid dynamic consequences of such events is therefore critical to understanding the threat to surrounding areas. In the 1995 Ruapehu eruption, a series of explosions through the Crater Lake ejected jets of hot tephra and warm lake water upon the snow cover, generating two bifurcating lahars that took paths directly down a ski-field. The late hour and limited runout precluded any loss of life or property. Fieldwork immediately after the event acquired samples and multiple thickness profiles across the deposit. Componentry and geochemical analyses showed that the deposit consisted of 10-30 wt% tephra and lake sediments, 1-4 wt% frozen lake water, and 65-85 wt% snow and ice. We infer that the evolving thermal balance within the lahar induced it to come to rest as the liquid component (lake water and melted snow/ice) progressively froze with distance from the source. We will present progress in development of a treatment of the heat transfer and fluid dynamics tephra-water-ice mixtures to constrain conditions leading to the widely varying runout distances and associated hazards of such events.

  11. Integrated research in constitutive modelling at elevated temperatures, part 2

    NASA Technical Reports Server (NTRS)

    Haisler, W. E.; Allen, D. H.

    1986-01-01

    Four current viscoplastic models are compared experimentally with Inconel 718 at 1100 F. A series of tests were performed to create a sufficient data base from which to evaluate material constants. The models used include Bodner's anisotropic model; Krieg, Swearengen, and Rhode's model; Schmidt and Miller's model; and Walker's exponential model.

  12. An approach to source characterization of tremor signals associated with eruptions and lahars

    NASA Astrophysics Data System (ADS)

    Kumagai, Hiroyuki; Mothes, Patricia; Ruiz, Mario; Maeda, Yuta

    2015-11-01

    Tremor signals are observed in association with eruption activity and lahar descents. Reduced displacement ( D R) derived from tremor signals has been used to quantify tremor sources. However, tremor duration is not considered in D R, which makes it difficult to compare D R values estimated for different tremor episodes. We propose application of the amplitude source location (ASL) method to characterize the sources of tremor signals. We used this method to estimate the tremor source location and source amplitude from high-frequency (5-10 Hz) seismic amplitudes under the assumption of isotropic S-wave radiation. We considered the source amplitude to be the maximum value during tremor. We estimated the cumulative source amplitude ( I s) as the offset value of the time-integrated envelope of the vertical seismogram of tremor corrected for geometrical spreading and medium attenuation in the 5-10-Hz band. For eruption tremor signals, we also estimated the cumulative source pressure ( I p) from an infrasonic envelope waveform corrected for geometrical spreading. We studied these parameters of tremor signals associated with eruptions and lahars and explosion events at Tungurahua volcano, Ecuador. We identified two types of eruption tremor at Tungurahua: noise-like inharmonic waveforms and harmonic oscillatory signals. We found that I s increased linearly with increasing source amplitude for lahar tremor signals and explosion events, but I s increased exponentially with increasing source amplitude for inharmonic eruption tremor signals. The source characteristics of harmonic eruption tremor signals differed from those of inharmonic tremor signals. We found a linear relation between I s and I p for both explosion events and eruption tremor. Because I p may be proportional to the total mass involved during an eruption episode, this linear relation suggests that I s may be useful to quantify eruption size. The I s values we estimated for inharmonic eruption tremor were

  13. Monitoring the Dynamic of a Fluvial Channel after Lahar Disturbance: Huiloac Gorge (Popocatepetl Volcano, Mexico)

    NASA Astrophysics Data System (ADS)

    Andres, N.; Palacios, D.; Zamorano, J. J.; Tanarro, L. M.; Renschler, C.; Sanjosé, J. J.; Atkinson, A.

    2009-04-01

    Volcanic eruptions generate disturbances that affect hydrological systems (Major, 2003) by depositing large volumes of sediments in watersheds that exceed amounts common to non-volcanic river systems (Montgomery, 2005). If the eruption releases abundant melt water, the river system may respond immediately by forming hazardous flows called lahars. River system recovery following eruptive and laharic impact is an important process, but it has received little attention (Gran and Montgomery, 2005) despite the fact that Major et al. (2000) and Hayes et al. (2002) have shown that these disruptions cause long term instability and their effects persist for decades. Lahar deposits resulting from interaction between volcanic activity and the glacier located above the Huiloac Gorge on the northern slope of Popocatepetl volcano (19°02´ N, 98°62´ W, 5,424 m), have infilled the gorge (Palacios, 1995; Palacios et al., 1998 and 2001; Capra et al., 2004; Muñoz, 2007). All of the major lahars that occurred on the volcano in 1995 (4 km), 1997 (21 km), and 2001 (14 km) have channelled through Huiloac Gorge, and have dramatically altered its morphology and dynamics through erosion and deposition. The present study traces these changes in the aftermath of the laharic events that occurred from 1997-2001. A sector of the channel, located at 3200m-3240m altitude, of 500 m long and 15 to 20 m wide, in the mid-section of the gorge, was chosen as the control site. Precipitation is heaviest there and is most apt to trigger secondary post-eruptive lahars. ArcGis software was used to draw 6 geomorphic maps of the site showing spatial variations in the landforms for the period February 2002 - February 2008. In addition, 29 cross-profiles were made of the gorge for the same time interval, excluding February 2004. The volume of sediment eroded and deposited was calculated for each date by comparing variations in the height of the floor and banks of the gorge depicted in the cross-profile, and

  14. Comparison of 7.5-minute and 1-degree digital elevation models

    NASA Technical Reports Server (NTRS)

    Isaacson, Dennis L.; Ripple, William J.

    1995-01-01

    We compared two digital elevation models (DEM's) for the Echo Mountain SE quadrangle in the Cascade Mountains of Oregon. Comparisons were made between 7.5-minute (1:24,000-scale) and 1-degree (1:250,000-scale) images using the variables of elevation, slope aspect, and slope gradient. Both visual and statistical differences are presented.

  15. Comparison of 7.5-minute and 1-degree digital elevation models

    NASA Technical Reports Server (NTRS)

    Isaacson, Dennis L.; Ripple, William J.

    1990-01-01

    Two digital elevation models are compared for the Echo Mountain SE quadrangle in the Cascade Mountains of Oregon. Comparisons were made between 7.5-minute (1:24,000-scale) and 1-degree (1:250,000-scale) images using the variables of elevation, slope aspect, and slope gradient. Both visual and statistical differences are presented.

  16. Voluminous ice-rich and water-rich lahars generated during the 2009 eruption of Redoubt Volcano, Alaska

    NASA Astrophysics Data System (ADS)

    Waythomas, Christopher F.; Pierson, Thomas C.; Major, Jon J.; Scott, William E.

    2013-06-01

    Redoubt Volcano in south-central Alaska began erupting on March 15, 2009, and by April 4, 2009, had produced at least 20 explosive events that generated multiple plumes of ash and numerous lahars. The 3108-m-high, snow- and ice-clad stratovolcano has an ice-filled summit crater that is breached to the north. The volcano supports about 4 km3 of ice and snow and about 1 km3 of this makes up the Drift glacier on the north side of the volcano. Explosive eruptions between March 23 and April 4, which included the destruction of at least two lava domes, triggered significant lahars in the Drift River valley on March 23 and April 4, and several smaller lahars between March 24 and March 31. Mud-line high-water marks, character of deposits, areas of inundation, and estimates of flow velocity revealed that the lahars on March 23 and April 4 were the largest of the eruption. In the 2-km-wide upper Drift River valley, average flow depths were at least 2-5 m. Average peak-flow velocities were likely between 10 and 15 ms- 1, and peak discharges were on the order of 104-105 m3 s- 1. The area inundated by lahars on March 23 was at least 100 km2 and on April 4 about 125 km2. Two substantial lahars emplaced on March 23 and one on April 4 had volumes on the order of 107-108 m3 and were similar in size to the largest lahar of the 1989-90 eruption. The two principal March 23 lahars were primarily flowing slurries of snow and ice derived from Drift glacier and the Drift River valley where seasonal snow and tabular blocks of river ice were entrained and incorporated into the lahars. Despite morphologic evidence of two lahars, only a single deposit up to 5 m thick was found in most places and it contained about 80-95% of poorly sorted, massive to imbricate assemblages of snow and ice clasts. The deposit was frozen soon after it was emplaced and later eroded and buried by the April 4 lahar. The lahar of April 4, in contrast, was primarily a hyperconcentrated flow, as interpreted from 1- to

  17. Influence of elevation modelling on hydrodynamic simulations of a tidally-dominated estuary

    NASA Astrophysics Data System (ADS)

    Falcão, Ana Paula; Mazzolari, Andrea; Gonçalves, Alexandre B.; Araújo, Maria Amélia V. C.; Trigo-Teixeira, António

    2013-08-01

    Hydrodynamic simulation of estuaries requires a single digital elevation model (DEM) resulting from merging of both topographic and bathymetric data. These two datasets are usually produced using different technologies, co-ordinate systems and datums. Intertidal data in particular are often lacking due to the difficulty of data acquisition using conventional survey techniques. This paper presents a fast, accurate and low-cost methodology to fill this gap and highlights the effect of the digital elevation model characteristics, such as the interpolation method and spatial resolution, on modelled water levels and flooded areas. The Lima river estuary, located in North-western Portugal, is used as a case study. Validation tests for commonly available spatial interpolators showed ordinary kriging to be the most adequate interpolator. Digital elevation models with regular grids of 5 m and 50 m resolution were used, together with the original (not interpolated) elevation dataset, as input to a finite element hydrodynamic model for astronomic tide simulation. Results indicate that the larger differences between using different elevation models occur at low tide during spring tide, marginally impacting the flood modelling. The effect of a vertical offset of the chart datum with respect to a part of the digital elevation model was finally investigated, showing a limited influence in the determination of the water levels.

  18. An algorithm for treating flat areas and depressions in digital elevation models using linear interpolation

    EPA Science Inventory

    Digital elevation model (DEM) data are essential to hydrological applications and have been widely used to calculate a variety of useful topographic characteristics, e.g., slope, flow direction, flow accumulation area, stream channel network, topographic index, and others. Excep...

  19. Quantifying Slopes with Digital Elevation Models of the Verdugo Hills, California: Effects of Resolution

    NASA Technical Reports Server (NTRS)

    Fielding, E. J.; Burbank, D. W.; Duncan, C. C.

    1996-01-01

    Quantification of surface slope angles is valuable in a wide variety of earth sciences. Slopes measured from digital elevation models (DEMs) or other topographic data sets depend strongly on the length scale or window size used in the slope calculations.

  20. Lahar Hazards at Casita and San Cristóbal Volcanoes, Nicaragua

    USGS Publications Warehouse

    Vallance, J.W.; Schilling, S.P.; Devoli, G.; Reid, M.E.; Howell, M.M.; Brien, D.L.

    2004-01-01

    to form ash-fall deposits (tephra), debris avalanches, lava flows, and hot flowing mixtures of ash and rock (called pyroclastic flows). The chronology of activity at Casita is rather poorly known. Its last documented eruption occurred 8300 years ago, and included a pyroclastic flow. Tephra deposits exposed in the east crater suggest the possibility of subsequent eruptions. Work prior to Hurricane Mitch suggested that a part of the volcano’s apron that included the area inundated during the 1998 event south of Casita was a lahar pathway. Erosion during Hurricane Mitch revealed that at least three large lahars descended this pathway to distances of up to 10 km. This report describes the hazards of landslides and lahars in general, and discusses potential hazards from future landslides and lahars at San Cristóbal and Casita volcanoes in particular. The report also shows, in the accompanying lahar hazard-zonation maps, which areas are likely to be at risk from future landslides and lahars at Casita and San Cristóbal.

  1. An uncoupled viscoplastic constitutive model for metals at elevated temperature

    NASA Technical Reports Server (NTRS)

    Haisler, W. E.; Cronenworth, J.

    1983-01-01

    An uncoupled constitutive model for predicting the transient response of thermal and rate dependent, inelastic material behavior is presented. The uncoupled model assumes that there is a temperature below which the total strain consists essentially of elastic and rate insensitive inelastic strains only. Above this temperature, the rate dependent inelastic strain (creep) dominates. The rate insensitive inelastic strain component is modeled in an incremental form with a yield function, flow rule and hardening law. Revisions to the hardening rule permit the model to predict temperature-dependent kinematic-isotropic hardening behavior, cyclic saturation, asymmetric stress-strain response upon stress reversal, and variable Bauschinger effect. The rate dependent inelastic strain component is modeled using a rate equation in terms of back stress, drag stress and exponent n as functions of temperature and strain. A sequence of hysteresis loops and relaxation tests are utilized to define the rate dependent inelastic strain rate. Evaluation of the model is performed by comparison with experiments involving various thermal and mechanical load histories on 5086 aluminum alloy, 304 stainless steel and Hastelloy-X.

  2. Eruption of the nevado del ruiz volcano, Colombia, on 13 november 1985: tephra fall and lahars.

    PubMed

    Naranjo, J L; Sigurdsson, H; Carey, S N; Fritz, W

    1986-08-29

    A small Plinian eruption of the Nevado del Ruiz volcano in Colombia ejected 3.5 x 10(10) kilograms of mixed dacite and andesite tephra on 13 November 1985, with a maximum column height of 31 kilometers above sea level. Small pyroclastic flows and surges, generated during the initial stage of the eruption, caused surface melting of approximately 10% of the volcano's ice cap, leading to meltwater floods. The erosive floods incorporated soils and loose sediments from the volcano's flanks and developed into lahars, which claimed at least 25,000 lives. PMID:17732038

  3. Downstream dilution of a lahar: transition from debris flow to hyperconcentrated streamflow.

    USGS Publications Warehouse

    Pierson, T.C.; Scott, K.M.

    1985-01-01

    Nearly instantaneous melting of snow and ice by the March 19, 1982, eruption of Mount St. Helens, released a 4 X 106 m3 flood of water from the crater that was converted to a lahar (volcanic debris flow) through erosion and incorporation of sediment by the time it reached the base of the volcano. Over the next 81 km that it traveled down the Toutle River, the flood wave was progressively diluted through several mechanisms. A transformation from debris flow to hyperconcentrated streamflow began to occur about 27 km downstream from the crater, when the total sediment concentration had decreased to about 78% by weight (57% by volume).-from Authors

  4. Validating Flood Mapping Products Using a Digital Elevation Model Comparison Technique

    NASA Astrophysics Data System (ADS)

    Fayne, J.

    2014-12-01

    This preliminary study assessed the validity of a pixel analysis elevation comparison technique and determined necessary steps for improvement. The pixel analysis sought to assess the probability of a flood occurring in a particular area by comparing the spatial extent of flood mapping products to the local elevation. The method was developed to determine if the physical relationship between elevation and floods as shown in satellite images is accurately represented in a flood mapping product. The data incorporated in this study are raster digital elevation model (DEM) tiles, a scene from Landsat 5 during a flood period, and a scene from the NASA DEVELOP Flood Disasters Team Flood Product. Pixels representing flooded areas were compared to the elevation height pixels using horizontal transect lines to create pixel value profiles across a 727 km transect of Vietnam and Cambodia. The elevation model comparison validates the Flood Product by depicting water presence in alignment with areas of low elevation. Initial findings indicate that the technique can be used to improve the assessment of flood mapping products in transects less than 10 km. Future research will focus on streamlining the pixel analysis process to yield comprehensive results for larger areas.

  5. Frozen Martian lahars? Evaluation of morphology, degradation and geologic development in the Utopia-Elysium transition zone

    NASA Astrophysics Data System (ADS)

    Pedersen, G. B. M.

    2013-09-01

    Regional coverage of high-resolution data from the CTX camera has permitted new, detailed morphologic analysis of the enigmatic Utopia-Elysium flows which dominate the transition zone between Elysium volcanic province and Utopia Planitia. Based on topographic and morphologic analysis of the Galaxias region, this study supports the lahar hypothesis put forth by previous works and suggests that the center and the margins of the outflow deposits have very diverse morphologies that can be explained by varying degrees of water drainage and freezing. Regular channel and flood plain deposits are found in the central part of the outflow deposits, whereas the marginal deposits are interpreted to contain significant amount of ice because of their distinct morphological properties (smooth, lobate flow-fronts with upward convex snouts, unusual crater morphologies, raised rim fractures and localized flow fronts indicating rheomorphism). Thus, this study suggest that, unlike terrestrial lahars, lahar emplacement under Martian conditions only drain in the central parts, whereas the water in the margins of the outflow deposit (∼75% of the total outflow deposit in the Galaxias region) freezes up resulting in a double-layered deposit consisting of ice-rich core with an ice-poor surface layer. It is here furthermore suggested that continued intrusive volcanic activity was highly affected by the presence of the ice-rich lahar deposits, generating ground-ice-volcano interactions resulting in a secondary suite of morphologies. These morphologies include seventeen ridges that are interpreted to be möberg ridges (due to their NW-SE orientation, distinct ridge-crests and association with fractures and linear ridges) and depressions with nested faults interpreted to be similar to terrestrial ice-cauldrons, which form by enhanced subglacial geothermal activity including subglacial volcanic eruptions. These sub-lahar intrusions caused significant volatile loss in the ice-rich core of the

  6. Process-scale modeling of elevated wintertime ozone in Wyoming.

    SciTech Connect

    Kotamarthi, V. R.; Holdridge, D. J.; Environmental Science Division

    2007-12-31

    Measurements of meteorological variables and trace gas concentrations, provided by the Wyoming Department of Environmental Quality for Daniel, Jonah, and Boulder Counties in the state of Wyoming, were analyzed for this project. The data indicate that highest ozone concentrations were observed at temperatures of -10 C to 0 C, at low wind speeds of about 5 mph. The median values for nitrogen oxides (NOx) during these episodes ranged between 10 ppbv and 20 ppbv (parts per billion by volume). Measurements of volatile organic compounds (VOCs) during these periods were insufficient for quantitative analysis. The few available VOCs measurements indicated unusually high levels of alkanes and aromatics and low levels of alkenes. In addition, the column ozone concentration during one of the high-ozone episodes was low, on the order of 250 DU (Dobson unit) as compared to a normal column ozone concentration of approximately 300-325 DU during spring for this region. Analysis of this observation was outside the scope of this project. The data analysis reported here was used to establish criteria for making a large number of sensitivity calculations through use of a box photochemical model. Two different VOCs lumping schemes, RACM and SAPRC-98, were used for the calculations. Calculations based on this data analysis indicated that the ozone mixing ratios are sensitive to (a) surface albedo, (b) column ozone, (c) NOx mixing ratios, and (d) available terminal olefins. The RACM model showed a large response to an increase in lumped species containing propane that was not reproduced by the SAPRC scheme, which models propane as a nearly independent species. The rest of the VOCs produced similar changes in ozone in both schemes. In general, if one assumes that measured VOCs are fairly representative of the conditions at these locations, sufficient precursors might be available to produce ozone in the range of 60-80 ppbv under the conditions modeled.

  7. Calculation and Error Analysis of a Digital Elevation Model of Hofsjokull, Iceland from SAR Interferometry

    NASA Technical Reports Server (NTRS)

    Barton, Jonathan S.; Hall, Dorothy K.; Sigurosson, Oddur; Williams, Richard S., Jr.; Smith, Laurence C.; Garvin, James B.

    1999-01-01

    Two ascending European Space Agency (ESA) Earth Resources Satellites (ERS)-1/-2 tandem-mode, synthetic aperture radar (SAR) pairs are used to calculate the surface elevation of Hofsjokull, an ice cap in central Iceland. The motion component of the interferometric phase is calculated using the 30 arc-second resolution USGS GTOPO30 global digital elevation product and one of the ERS tandem pairs. The topography is then derived by subtracting the motion component from the other tandem pair. In order to assess the accuracy of the resultant digital elevation model (DEM), a geodetic airborne laser-altimetry swath is compared with the elevations derived from the interferometry. The DEM is also compared with elevations derived from a digitized topographic map of the ice cap from the University of Iceland Science Institute. Results show that low temporal correlation is a significant problem for the application of interferometry to small, low-elevation ice caps, even over a one-day repeat interval, and especially at the higher elevations. Results also show that an uncompensated error in the phase, ramping from northwest to southeast, present after tying the DEM to ground-control points, has resulted in a systematic error across the DEM.

  8. Utility of ethological analysis to overcome locomotor confounds in elevated maze models of anxiety.

    PubMed

    Weiss, S M; Wadsworth, G; Fletcher, A; Dourish, C T

    1998-01-01

    The elevated plus-maze is a commonly used model to identify putative anxiolytic and anxiogenic drugs. However, the validity of elevated plus-maze and other recently developed variants such as the elevated zero-maze has recently been questioned on the grounds that both the reference anxiolytic drug chlordiazepoxide and the psychostimulant d-amphetamine increase open arm exploration and stimulate locomotor activity. These findings suggest that measures of "anxiety" in the elevated maze cannot be adequately dissociated from simple changes in locomotor activity, which may confound the interpretation of results obtained using these models. A variety of approaches to assess drug effects on locomotor activity in the elevated maze have been suggested, including the use of total and closed arm entries, as well as supplementary tests such as exploration of the holeboard apparatus. However, all these approaches utilise the measurement of exploration in a novel environment, and as such, could potentially be influenced by either changes in anxiety or locomotor activity. Recently, it has been shown that ethological measures of "risk assessment", such as stretched-attend postures and head-dipping, are sensitive indicators of drug-effects in the elevated maze. The present study assessed the utility of ethological analysis in dissociating locomotor activity from "anxiety" by comparing the effects of d-amphetamine to those of chlordiazepoxide in the rat elevated zero-maze. The results showed that both chlordiazepoxide and d-amphetamine increase the amount of time spent in the open arms and reduce "risk assessment" without increasing line crossing or rearing. These results confirm that under certain test conditions, psychostimulants are capable of producing "false-positives" in elevated maze models, and that both traditional methods and the ethological measures used in this study fail to unequivocally dissociate drug effects on anxiety from effects on locomotor activity. Further

  9. Assessment of Required Accuracy of Digital Elevation Data for Hydrologic Modeling

    NASA Technical Reports Server (NTRS)

    Kenward, T.; Lettenmaier, D. P.

    1997-01-01

    The effect of vertical accuracy of Digital Elevation Models (DEMs) on hydrologic models is evaluated by comparing three DEMs and resulting hydrologic model predictions applied to a 7.2 sq km USDA - ARS watershed at Mahantango Creek, PA. The high resolution (5 m) DEM was resempled to a 30 m resolution using method that constrained the spatial structure of the elevations to be comparable with the USGS and SIR-C DEMs. This resulting 30 m DEM was used as the reference product for subsequent comparisons. Spatial fields of directly derived quantities, such as elevation differences, slope, and contributing area, were compared to the reference product, as were hydrologic model output fields derived using each of the three DEMs at the common 30 m spatial resolution.

  10. Rapid calculation of terrain parameters for radiation modeling from digital elevation data

    NASA Technical Reports Server (NTRS)

    Dozier, Jeff; Frew, James

    1990-01-01

    Digital elevation models are now widely used to calculate terrain parameters to determine incoming solar and longwave radiation for use in surface climate models, interpretation of remote-sensing data, and parameters in hydrologic models. Because of the large number of points in an elevation grid, fast algorithms are useful to save computation time. A description is given of rapid methods for calculating slope and azimuth, solar illumination angle, horizons, and view factors for radiation from sky and terrain. Calculation time is reduced by fast algorithms and lookup tables.

  11. Lahars in and around the Taipei basin: Implications for the activity of the Shanchiao fault

    NASA Astrophysics Data System (ADS)

    Song, Sheng-Rong; Chen, Tsu-Mo; Tsao, Shuhjong; Chen, Huei-Fen; Liu, Huan-Chi

    2007-11-01

    In the last decade, more than 21 deep geological cores have been drilled in the Taipei basin to obtain a firmer grasp of its basic geology and engineering properties prior to the construction of new infrastructure. Thirteen of those cores contain lahar deposits, with the number of layers varying from one to three and the thickness of each layer varying from several to over 100 m. Based on their occurrence, petrology and geochemistry, it has been determined that the deposits originated from the southern slope of the Tatun Volcano Group (TVG). K-Ar age dating has shown that the lower layer of lahars was deposited less than 0.4 Ma, and this is clearly correlated to outcrops in the Kauntu, Chengtzeliao and Shihtzutao areas. These findings may well suggest that the Taipei basin has been formed in last 0.4 Ma and that the Shanchiao normal fault commenced its activity within this period. The surface trace and the activity of the Shanchiao normal fault have also been inferred and subsequently defined from stratigraphic data derived from these cores.

  12. Relationship between geomorphology and lithotypes of lahar deposit from Chokai volcano, Japan

    NASA Astrophysics Data System (ADS)

    Minami, Y.; Ohba, T.; Hayashi, S.; Kataoka, K.

    2013-12-01

    Chokai volcano, located in the northern Honshu arc in Japan, is an andesitic stratovolcano that collapsed partly at ca. 2500 years ago. A post collapse lahar deposit (Shirayukigawa lahar deposit) is distributed in the northern foot of the volcanic edifice. The deposit consists of 16 units of debris flow, hyperconcentrated flow and streamflow deposits. The Shirayukigawa lahar deposit has a total thickness of 30 m and overlies the 2.5-ka Kisakata debris avalanche deposit. Shirayukigawa lahar deposit forms volcanic fan and volcanic apron. The volcanic fan is subdivided into four areas on the basis of slope angles and of geomorphological features: 1) steeply sloped area, 2) moderately sloped area, 3) gently sloped area and 4) horizontal area. From sedimentary facies and structures, each unit of the Shirayukigawa lahar deposit is classified into one of four lithotypes: clast-supported debris flow deposit (Cc), matrix-supported debris flow deposit (Cm1), hyperconcentrated flow deposit (Cm2) and streamflow deposit (Sl). Each type has the following lithological characteristics. The lithotypes are well correlated with the geomorphology of the volcanic fan. The steeply-sloped and the moderately-sloped areas are dominated by Cc, Cm1, and Cm2, and The horizontal area are dominated by Sl. Debris flow deposit (Cc) is massive, very poorly sorted, partly graded, and clast-supported with polymictic clasts dominated by subrounded to rounded volcanic clasts. Matrix is sandy to muddy. Preferred clast orientation are present. Debris flow deposit (Cm1) is massive, very poorly sorted, and matrix-supported with polymictic clasts dominated by subrounded to rounded volcanic clasts. Matrix is sandy to muddy. Some layers exhibit coarse-tail normal/inverse grading. Most clasts are oriented. Hyperconcentrated flow deposit (Cm2) is massive to diffusely laminated, very poorly sorted and matrix-supported with polymictic clasts dominated by subrounded to rounded volcanic rocks. Matrix is sandy. The

  13. Elevator control system modeling and PLC program generating based on SIPN

    NASA Astrophysics Data System (ADS)

    Zhong, Yanru; Gao, Xiaofeng; Huang, Meifa

    2006-11-01

    This paper researches into the modeling method of elevator control system and the algorithm of generating PLC program based on Signal Interpreted Petri Net (SIPN). We also analyze the properties of the SIPN model in the system. SIPN is obtained by adding input and output signals into the ordinary Petri Net. Input signals are related to every transition in SIPN, which express the firing conditions of the transition, while output signals related to the corresponding place, which express the control information of the place. In the SIPN model of the system, the relationships are established from one to one correspondence between the Input and Output (I/O) signals of SIPN and the I/O points of PLC. The places of the SIPN stand for the statuses of the elevator, and the transitions represent the changes of the statuses. The input signals of system are the actions of pressing buttons and some sensor signals, which represent the conditions of changing the states of elevator. The output signals are the actions of elevator control functions, such as stop, open, and close the door, which describe the output control information of the elevator in every place. A case study demonstrates the validity both of the SIPN model and the algorithm.

  14. Tectonic development of the Northwest Bonaparte Basin, Australia by using Digital Elevation Model (DEM)

    NASA Astrophysics Data System (ADS)

    Wahid, Ali; Salim, Ahmed Mohamed Ahmed; Ragab Gaafar, Gamal; Yusoff, AP Wan Ismail Wan

    2016-02-01

    The Bonaparte Basin consist of majorly offshore part is situated at Australia's NW continental margin, covers an area of approx. 270,000km2. Bonaparte Basin having a number of sub-basins and platform areas of Paleozoic and Mesozoic is structurally complex. This research established the geologic and geomorphologic studies using Digital Elevation Model (DEM) as a substitute approach in morphostructural analysis to unravel the geological complexities. Although DEMs have been in practice since 1990s, they still have not become common tool for mapping studies. The research work comprised of regional structural analysis with the help of integrated elevation data, satellite imageries, available open topograhic images and internal geological maps with interpreted seismic. The structural maps of the study area have been geo-referenced which further overlaid onto SRTM data and satellite images for combined interpretation which facilitate to attain Digital Elevation Model of the study area. The methodology adopts is to evaluate and redefine development of geodynamic processes involved in formation of Bonaparte Basin. The main objectives is to establish the geological histories by using digital elevation model. The research work will be useful to incorporate different tectonic events occurred at different Geological times in a digital elevation model. The integrated tectonic analysis of different digital data sets benefitted substantially from combining them into a common digital database. Whereas, the visualization software facilitates the overlay and combined interpretation of different data sets which is helpful to reveal hidden information not obvious or accessible otherwise for regional analysis.

  15. A Model for Predicting Persistent Elevation of Factor VIII among Patients with Acute Ischemic Stroke

    PubMed Central

    Samai, Alyana A.; Boehme, Amelia K.; Shaban, Amir; George, Alexander J.; Dowell, Lauren; Monlezun, Dominique J.; Leissinger, Cindy; Schluter, Laurie; El Khoury, Ramy; Martin-Schild, Sheryl

    2016-01-01

    Background and Purpose Elevated levels of coagulation factor VIII (FVIII) may persist independent of the acute-phase response; however, this relationship has not been investigated relative to acute ischemic stroke (AIS). We examined the frequency and predictors of persistently elevated FVIII in AIS patients. Methods AIS patients admitted between July 2008 and May 2014 with elevated baseline FVIII levels and repeat FVIII levels drawn for more than 7 days postdischarge were included. The patients were dichotomized by repeat FVIII level for univariate analysis at 150% and 200% activity thresholds. An adjusted model was developed to predict the likelihood of persistently elevated FVIII levels. Results Among 1616 AIS cases, 98 patients with elevated baseline FVIII had repeat FVIII levels. Persistent FVIII elevation was found in more than 75% of patients. At the 150% threshold, the prediction score ranged from 0 to 7 and included black race, female sex, prior stroke, hyperlipidemia, smoking, baseline FVIII > 200%, and baseline von Willebrand factor (vWF) level greater than 200%. At the 200% threshold, the prediction score ranged from 0–5 and included female sex, prior stroke, diabetes mellitus, baseline FVIII level greater 200%, and baseline vWF level greater than 200%. For each 1-point increase in score, the odds of persistent FVIII at both the 150% threshold (odds ratio [OR] = 10.4, 95% confidence interval [CI] 1.63–66.9, P = .0134) and 200% threshold (OR = 10.2, 95% CI 1.82–57.5, P = .0083) increased 10 times. Conclusion Because an elevated FVIII level confers increased stroke risk, our model for anticipating a persistently elevated FVIII level may identify patients at high risk for recurrent stroke. FVIII may be a target for secondary stroke prevention. PMID:26777556

  16. The Importance of Precise Digital Elevation Models (DEM) in Modelling Floods

    NASA Astrophysics Data System (ADS)

    Demir, Gokben; Akyurek, Zuhal

    2016-04-01

    Digital elevation Models (DEM) are important inputs for topography for the accurate modelling of floodplain hydrodynamics. Floodplains have a key role as natural retarding pools which attenuate flood waves and suppress flood peaks. GPS, LIDAR and bathymetric surveys are well known surveying methods to acquire topographic data. It is not only time consuming and expensive to obtain topographic data through surveying but also sometimes impossible for remote areas. In this study it is aimed to present the importance of accurate modelling of topography for flood modelling. The flood modelling for Samsun-Terme in Blacksea region of Turkey is done. One of the DEM is obtained from the point observations retrieved from 1/5000 scaled orthophotos and 1/1000 scaled point elevation data from field surveys at x-sections. The river banks are corrected by using the orthophotos and elevation values. This DEM is named as scaled DEM. The other DEM is obtained from bathymetric surveys. 296 538 number of points and the left/right bank slopes were used to construct the DEM having 1 m spatial resolution and this DEM is named as base DEM. Two DEMs were compared by using 27 x-sections. The maximum difference at thalweg of the river bed is 2m and the minimum difference is 20 cm between two DEMs. The channel conveyance capacity in base DEM is larger than the one in scaled DEM and floodplain is modelled in detail in base DEM. MIKE21 with flexible grid is used in 2- dimensional shallow water flow modelling. The model by using two DEMs were calibrated for a flood event (July 9, 2012). The roughness is considered as the calibration parameter. From comparison of input hydrograph at the upstream of the river and output hydrograph at the downstream of the river, the attenuation is obtained as 91% and 84% for the base DEM and scaled DEM, respectively. The time lag in hydrographs does not show any difference for two DEMs and it is obtained as 3 hours. Maximum flood extents differ for the two DEMs

  17. TERRAIN: A computer program to process digital elevation models for modeling surface flow

    SciTech Connect

    Schwartz, P.M.; Levine, D.A.; Hunsaker, C.T.; Timmins, S.P.

    1995-08-01

    This document provides a step by step procedure, TERRAIN, for processing digital elevation models to calculate overland flow paths, watershed boundaries, slope, and aspect. The algorithms incorporated into TERRAIN have been used at two different geographic scales: first for small research watersheds where surface wetness measurements are made, and second for regional water modeling for entire counties. For small areas methods based on flow distribution may be more desirable, especially if time-dependent flow models are to be used. The main improvement in TERRAIN compared with earlier programs on which it is based is that it combines the conditioning routines, which remove depressions to avoid water storage, into a single process. Efficiency has also been improved, reducing run times as much as 10:1 and enabling the processing of very large grids in strips for regional modeling. Additionally, the ability to calculate the nutrient load delivered any cell in a watershed has been added. These improvements make TERRAIN a powerful tool for modeling surface flow.

  18. Vapor pressure and boiling point elevation of slash pine black liquors: Predictive models with statistical approach

    SciTech Connect

    Zaman, A.A.; McNally, T.W.; Fricke, A.L.

    1998-01-01

    Vapor-liquid equilibria and boiling point elevation of slash pine kraft black liquors over a wide range of solid concentrations (up to 85% solids) has been studied. The liquors are from a statistically designed pulping experiment for pulping slash pine in a pilot scale digester with four cooking variables of effective alkali, sulfidity, cooking time, and cooking temperature. It was found that boiling point elevation of black liquors is pressure dependent, and this dependency is more significant at higher solids concentrations. The boiling point elevation data at different solids contents (at a fixed pressure) were correlated to the dissolved solids (S/(1 {minus} S)) in black liquor. Due to the solubility limit of some of the salts in black liquor, a change in the slope of the boiling point elevation as a function of the dissolved solids was observed at a concentration of around 65% solids. An empirical method was developed to describe the boiling point elevation of each liquor as a function of pressure and solids mass fraction. The boiling point elevation of slash pine black liquors was correlated quantitatively to the pulping variables, using different statistical procedures. These predictive models can be applied to determine the boiling point rise (and boiling point) of slash pine black liquors at processing conditions from the knowledge of pulping variables. The results are presented, and their utility is discussed.

  19. A theoretical model of cytosolic calcium elevation following wounding in urothelial cell monolayers

    NASA Astrophysics Data System (ADS)

    Appleby, Peter A.; Shabir, Saqib; Southgate, Jennifer; Walker, Dawn

    2013-02-01

    Scratch wounding of a urothelial cell monolayer triggers a number of events including the release of soluble, diffusible signalling factors and mechanical stimulation of cells at the wound edge. These events cause a sustained elevation in cytosolic calcium concentration in the cells surrounding the wound and a transient rise in those further away. The precise form of this calcium transient is believed to play a central role in determining the subsequent response of individual cells and ultimately leads to a co-ordinated, population-level response that rapidly closes the wound. Here we present a framework for modelling the initial phases of this process. We combine a PDE model of diffusion in the extracellular medium and an ODE model of calcium signalling that has been tailored to represent urothelial cells. The ODE model is capable of generating a wide range of calcium transients, including spikes, bursts, oscillations and sustained elevations in the cytosolic calcium concentration. In multi-cell simulations of scratch wounding in a perfusion flow we find that the spatial position of the cells relative to the wound site leads to distinct classes of calcium response, with cells proximal to the wound exhibiting a sustained elevation and cells distal to the wound exhibiting a more transient elevation. We compare these results to existing experimental data and generate a number of novel predictions that could be used to test the model experimentally.

  20. Test of the depression distress amplification model in young adults with elevated risk of current suicidality.

    PubMed

    Capron, Daniel W; Lamis, Dorian A; Schmidt, Norman B

    2014-11-30

    Suicide is a leading cause of death among young adults and the rate of suicide has been increasing for decades. A depression distress amplification model posits that young adults with comorbid depression and anxiety have elevated suicide rates due to the intensification of their depressive symptoms by anxiety sensitivity cognitive concerns. The current study tested the effects of anxiety sensitivity subfactors as well as the depression distress amplification model in a very large sample of college students with elevated suicide risk. Participants were 721 college students who were at elevated risk of suicidality (scored>0 on the Beck Scale for Suicide Ideation). Consistent with prior work, anxiety sensitivity cognitive concerns, but not physical or social concerns, were associated with suicidal ideation. Consistent with the depression distress amplification model, in individuals high in depression, anxiety sensitivity cognitive concerns predicted elevated suicidal ideation but not among those with low depression. The results of this study corroborate the role of anxiety sensitivity cognitive concerns and the depression distress amplification model in suicidal ideation among a large potentially high-risk group of college students. The depression distress amplification model suggests a specific mechanism, anxiety sensitivity cognitive concerns, that may be responsible for increased suicide rates among those with comorbid anxiety and depression. PMID:25063018

  1. Segregating two simultaneous sounds in elevation using temporal envelope: Human psychophysics and a physiological model.

    PubMed

    Johnson, Jeffrey S; O'Connor, Kevin N; Sutter, Mitchell L

    2015-07-01

    The ability to segregate simultaneous sound sources based on their spatial locations is an important aspect of auditory scene analysis. While the role of sound azimuth in segregation is well studied, the contribution of sound elevation remains unknown. Although previous studies in humans suggest that elevation cues alone are not sufficient to segregate simultaneous broadband sources, the current study demonstrates they can suffice. Listeners segregating a temporally modulated noise target from a simultaneous unmodulated noise distracter differing in elevation fall into two statistically distinct groups: one that identifies target direction accurately across a wide range of modulation frequencies (MF) and one that cannot identify target direction accurately and, on average, reports the opposite direction of the target for low MF. A non-spiking model of inferior colliculus neurons that process single-source elevation cues suggests that the performance of both listener groups at the population level can be accounted for by the balance of excitatory and inhibitory inputs in the model. These results establish the potential for broadband elevation cues to contribute to the computations underlying sound source segregation and suggest a potential mechanism underlying this contribution. PMID:26233004

  2. Segregating two simultaneous sounds in elevation using temporal envelope: Human psychophysics and a physiological model

    PubMed Central

    Johnson, Jeffrey S.; O'Connor, Kevin N.; Sutter, Mitchell L.

    2015-01-01

    The ability to segregate simultaneous sound sources based on their spatial locations is an important aspect of auditory scene analysis. While the role of sound azimuth in segregation is well studied, the contribution of sound elevation remains unknown. Although previous studies in humans suggest that elevation cues alone are not sufficient to segregate simultaneous broadband sources, the current study demonstrates they can suffice. Listeners segregating a temporally modulated noise target from a simultaneous unmodulated noise distracter differing in elevation fall into two statistically distinct groups: one that identifies target direction accurately across a wide range of modulation frequencies (MF) and one that cannot identify target direction accurately and, on average, reports the opposite direction of the target for low MF. A non-spiking model of inferior colliculus neurons that process single-source elevation cues suggests that the performance of both listener groups at the population level can be accounted for by the balance of excitatory and inhibitory inputs in the model. These results establish the potential for broadband elevation cues to contribute to the computations underlying sound source segregation and suggest a potential mechanism underlying this contribution. PMID:26233004

  3. Digital elevation model (DEM) of Cascadia, latitude 39N-53N, longitude 116W-133W

    USGS Publications Warehouse

    Haugerud, Ralph A.

    1999-01-01

    This report contains a 250-meter digital elevation model (DEM) for Cascadia (latitude 39N - 53N, longitude 116W - 133W), a region that encompasses the Cascade volcanic arc, the Cascadia subduction zone, and the Juan de Fuca Ridge system. The DEM is distributed as file cascdem.tar.gz (39 MB; 78MB uncompressed).

  4. 90-METER DIGITAL ELEVATION MODEL (DEM) FOR THE MID-ATLANTIC INTEGRATED ASSESSMENT (MAIA) STUDY AREA

    EPA Science Inventory

    This data set is a geographic information system (GIS) coverage of the 90-meter digital elevation model (DEM) for the United States Environmental Protection Agency (USEPA) Mid-Atlantic Integrated Assessment (MAIA) Project region. The coverage was produced using US Geological Su...

  5. VERTICAL DISPERSION FROM SURFACE AND ELEVATED RELEASES: AN INVESTIGATION OF A NON-GAUSSIAN PLUME MODEL

    EPA Science Inventory

    The vertical diffusion of a passive tracer released from surface and elevated sources in a neutrally-stratified boundary layer has been studied by comparing field and laboratory experiments with a non-Gaussian K-theory model that assumes power-law profiles for the mean velocity a...

  6. Lunar Topography and Basins Mapped Using a Clementine Stereo Digital Elevation Model

    NASA Technical Reports Server (NTRS)

    Cook, A. C.; Spudis, P. D.; Robinson, M. S.; Watters, T. R.

    2002-01-01

    Planet-wide (1 km/pixel and 5 km/pixel) Digital Elevation Models (DEM) of the Moon have been produced using Clementine UVVIS (Ultraviolet-Visible) stereo. Six new basins have been discovered, two suspected basins have been confirmed, and the dimensions of existing basins better defined. Additional information is contained in the original extended abstract.

  7. The effects of elevation data representation on mesoscale atmospheric model simulations

    SciTech Connect

    Walker, H.; Leone, J.M. Jr.; Kim, Jinwon

    1996-01-01

    Mesoscale atmospheric model simulations rely on descriptions of the land surface characteristics, which must be developed from geographic databases. Certain features of the geographic data, such as its resolution and accuracy, as well as the method of processing for use in the model, can be very important in producing accurate model simulations. The work described here is part of research effort into the relationship between these aspects of geographic data and the performance of mesoscale atmospheric models and is particularly focused on elevation data and how it is prepared for use in such models. A source for digital elevation data will typically not be at the resolution required for a given model simulation and so a resampling step is required. In addition, predictive non-linear model often cannot accept forcing at high spatial frequencies due to the terrain, thus smoothing is also required. The effect of different means of resampling and smoothing elevation data on two types of model simulations is investigated. At smaller spatial scales, nocturnal drainage winds in mountain valleys in Colorado are examined for effects on the general characteristics as well as the details of the flows. At the larger end of the mesoscale, extended simulations of California weather are examined for effects on orographic lifting, low-level convergence and divergence and ultimately rain and snow distribution.

  8. Modeling the population dynamics of Culex quinquefasciatus (Diptera: Culicidae), along an elevational gradient in Hawaii.

    PubMed

    Ahumada, Jorge A; Lapointe, Dennis; Samuel, Michael D

    2004-11-01

    We present a population model to understand the effects of temperature and rainfall on the population dynamics of the southern house mosquito, Culex quinquefasciatus Say, along an elevational gradient in Hawaii. We use a novel approach to model the effects of temperature on population growth by dynamically incorporating developmental rate into the transition matrix, by using physiological ages of immatures instead of chronological age or stages. We also model the effects of rainfall on survival of immatures as the cumulative number of days below a certain rain threshold. Finally, we incorporate density dependence into the model as competition between immatures within breeding sites. Our model predicts the upper altitudinal distributions of Cx. quinquefasciatus on the Big Island of Hawaii for self-sustaining mosquito and migrating summer sink populations at 1,475 and 1,715 m above sea level, respectively. Our model predicts that mosquitoes at lower elevations can grow under a broader range of rainfall parameters than middle and high elevation populations. Density dependence in conjunction with the seasonal forcing imposed by temperature and rain creates cycles in the dynamics of the population that peak in the summer and early fall. The model provides a reasonable fit to the available data on mosquito abundance for the east side of Mauna Loa, Hawaii. The predictions of our model indicate the importance of abiotic conditions on mosquito dynamics and have important implications for the management of diseases transmitted by Cx. quinquefasciatus in Hawaii and elsewhere. PMID:15605655

  9. Controls on interior West Antarctic Ice Sheet Elevations: inferences from geologic constraints and ice sheet modeling

    NASA Astrophysics Data System (ADS)

    Ackert, Robert P.; Putnam, Aaron E.; Mukhopadhyay, Sujoy; Pollard, David; DeConto, Robert M.; Kurz, Mark D.; Borns, Harold W.

    2013-04-01

    Knowledge of the West Antarctic Ice Sheet (WAIS) response to past sea level and climate forcing is necessary to predict its response to warmer temperatures in the future. The timing and extent of past interior WAIS elevation changes provides insight to WAIS behavior and constraints for ice sheet models. Constraints prior to the Last Glacial Maximum (LGM) however, are rare. Surface exposure ages of glacial erratics near the WAIS divide at Mt. Waesche in Marie Byrd Land, and at the Ohio Range in the Transantarctic Mountains, range from ∼10 ka to >500 ka without a dependence on elevation. The probability distribution functions (PDF) of the exposure ages at both locations, are remarkably similar. During the last glaciation, maximum interior ice elevations as recorded by moraines and erratics were reached between 10 ka and 12 ka. However, most exposure ages are older than the LGM and cluster around ∼40 ka and ∼80 ka. The peak in the exposure age distributions at ∼40 ka includes ages of alpine moraine boulders at Mercer Ridge in the Ohio Range. Comparison of the PDF of exposures ages from the Ohio Range and Mt. Waesche with the temperature record from the Fuji Dome ice core indicates that the youngest peak in the exposure age distributions corresponds to the abrupt warming during the Last Glacial termination. A prominent peak in the Ohio Range PDF corresponds to the penultimate termination (stage 5e). During the intervening glacial period, there is not a consistent relationship between the peaks in the PDF at each location and temperature. A combined ice sheet/ice shelf model with forcing scaled to marine δ18O predicts that interior WAIS elevations near the ice divide have varied ∼300 m over the Last Glacial cycle. Peaks in the PDF correspond to model highstands over the last 200 ka. In the simulated elevation history, maximum ice elevations at Ohio Range (+100 m) and Mt. Waesche (+60 m) occur at ∼10 ka, in agreement with observations from these sites

  10. A seamless, high-resolution digital elevation model (DEM) of the north-central California coast

    USGS Publications Warehouse

    Foxgrover, Amy C.; Barnard, Patrick L.

    2012-01-01

    A seamless, 2-meter resolution digital elevation model (DEM) of the north-central California coast has been created from the most recent high-resolution bathymetric and topographic datasets available. The DEM extends approximately 150 kilometers along the California coastline, from Half Moon Bay north to Bodega Head. Coverage extends inland to an elevation of +20 meters and offshore to at least the 3 nautical mile limit of state waters. This report describes the procedures of DEM construction, details the input data sources, and provides the DEM for download in both ESRI Arc ASCII and GeoTIFF file formats with accompanying metadata.

  11. ICESat Lidar and Global Digital Elevation Models: Application to DESDynI

    NASA Technical Reports Server (NTRS)

    Carabajal, Claudia C.; Harding, David J.; Suchdeo, Vijay P.

    2010-01-01

    Geodetic control is extremely important in the production and quality control of topographic data sets, enabling elevation results to be referenced to an absolute vertical datum. Global topographic data with improved geodetic accuracy achieved using global Ground Control Point (GCP) databases enable more accurate characterization of land topography and its change related to solid Earth processes, natural hazards and climate change. The multiple-beam lidar instrument that will be part of the NASA Deformation, Ecosystem Structure and Dynamics of Ice (DESDynI) mission will provide a comprehensive, global data set that can be used for geodetic control purposes. Here we illustrate that potential using data acquired by NASA's Ice, Cloud and land Elevation Satellite (ICEsat) that has acquired single-beam, globally distributed laser altimeter profiles (+/-86deg) since February of 2003 [1, 2]. The profiles provide a consistently referenced elevation data set with unprecedented accuracy and quantified measurement errors that can be used to generate GCPs with sub-decimeter vertical accuracy and better than 10 m horizontal accuracy. Like the planned capability for DESDynI, ICESat records a waveform that is the elevation distribution of energy reflected within the laser footprint from vegetation, where present, and the ground where illuminated through gaps in any vegetation cover [3]. The waveform enables assessment of Digital Elevation Models (DEMs) with respect to the highest, centroid, and lowest elevations observed by ICESat and in some cases with respect to the ground identified beneath vegetation cover. Using the ICESat altimetry data we are developing a comprehensive database of consistent, global, geodetic ground control that will enhance the quality of a variety of regional to global DEMs. Here we illustrate the accuracy assessment of the Shuttle Radar Topography Mission (SRTM) DEM produced for Australia, documenting spatially varying elevation biases of several meters

  12. Gravity modeling reveals that the "Miocene Pyrenean peneplain" developed at high elevation

    NASA Astrophysics Data System (ADS)

    Bosch, Gemma V.; Van Den Driessche, Jean; Robert, Alexandra; Babault, Julien; Le Carlier, Christian

    2016-04-01

    Geodynamics that shaped the present morphology of the western Mediterranean are mostly linked to the African-Eurasia collision and the extension related to the Mediterranean opening. The Pyrenean chain formed by the collision between the Iberian microplate and the Eurasian plate from the Eocene to the late Oligocene. This resulted in lithosphere thickening especially below the Central Pyrenees that becomes thinner eastwards. Whether the later thinning of the lithosphere in the easternmost Pyrenees involves the removal of the lithospheric mantle or not is debated. This issue joins the problematics about the origin of the high-elevation of the "Miocene Pyrenean peneplain" remnants. Indeed the most striking feature of the Pyrenean morphology is the occurrence of high-elevation, low relief erosional surfaces that are interpreted as the remnants of a Miocene single planation surface, dissected and reworked by Quaternary fluvial and glacial erosion. Two end-member interpretations have proposed to explain the high elevation of this original surface. The first considers that the Miocene Pyrenean peneplain develops near sea-level and was later uplifted, the second claims that the planation surface developed at high elevation in response to the inhibition of erosion consecutively to the progressive rise of the base-level of the Pyrenean drainage network. The first interpretation implies the return to normal crustal thickness by erosion and later uplift by removal of the lithospheric mantle. The second interpretation considers that the mean elevation of the original planation surface matches the thickness of the lithosphere below the chain, taking into account some hundred meters of isostatic rebound due to Quaternary erosion. To test these interpretations, we first restore the Miocene original planation surface by mapping and interpolating the high-elevation, low relief surfaces across the Pyrenees. We then performed 1D and 2D gravity models that we compare with recent

  13. Coupling of digital elevation model and rainfall-runoff model in storm drainage network design

    NASA Astrophysics Data System (ADS)

    Gumbo, Bekithemba; Munyamba, Nelson; Sithole, George; Savenije, Hubert H. G.

    Often planners and engineers are faced with various options and questions in storm drainage network design e.g. flow pattern, direction, runoff quantity and therefore size of drain, or scenario after a road, airfield or building has been constructed. In most instances planning without drainage in mind has caused failure or extensive damage to property including the storm water drains which channel the water away. With the advent of various modelling and geographic information systems (GIS) tools this problem can be averted. The University of Zimbabwe’s (UZ) main campus had its storm drainage network reconstructed at a cost of about US$100 000, because of persistent flooding. This paper describes a method of assessing the effectiveness of storm drainage networks by combining a digital elevation model (DEM) with a rainfall-runoff model based on the Soil Conservation Service South African manual (SCS-SA). The UZ campus was used as the test site. The DEM was generated from aerial photographs and the data imported into ArcView. The 3.0 km 2 basin was then delineated into sub-catchments using ArcView Hydro extension tools. The land-use, watershed and soil map of the UZ were merged in ArcView and initial curve numbers (CN) assigned. Using three years of daily rainfall data, runoff and peak flows were calculated for each sub-catchment. By overlaying the natural flow lines derived from the DEM with the reconstructed physical drains a comparison of the flow direction and the orientation of the drains was achieved. Peak flows where calculated for each delineated watershed and the results used to check the adequacy of the trapezoidal concrete lined drains. A combination of a DEM and rainfall-runoff model within a GIS platform proves to be useful in estimating runoff on partly urbanised watersheds and in determining the size and orientation of storm drains. It is particularly useful for new areas where development is being contemplated.

  14. A computational-grid based system for continental drainage network extraction using SRTM digital elevation models

    NASA Technical Reports Server (NTRS)

    Curkendall, David W.; Fielding, Eric J.; Pohl, Josef M.; Cheng, Tsan-Huei

    2003-01-01

    We describe a new effort for the computation of elevation derivatives using the Shuttle Radar Topography Mission (SRTM) results. Jet Propulsion Laboratory's (JPL) SRTM has produced a near global database of highly accurate elevation data. The scope of this database enables computing precise stream drainage maps and other derivatives on Continental scales. We describe a computing architecture for this computationally very complex task based on NASA's Information Power Grid (IPG), a distributed high performance computing network based on the GLOBUS infrastructure. The SRTM data characteristics and unique problems they present are discussed. A new algorithm for organizing the conventional extraction algorithms [1] into a cooperating parallel grid is presented as an essential component to adapt to the IPG computing structure. Preliminary results are presented for a Southern California test area, established for comparing SRTM and its results against those produced using the USGS National Elevation Data (NED) model.

  15. Comparative Analysis of Global Digital Elevation Models and Ultra-Prominent Mountain Peaks

    NASA Astrophysics Data System (ADS)

    Grohmann, Carlos H.

    2016-06-01

    Global Digital Elevation Models (GDEMs) are datasets of vital importance for regional-scale analysis in areas such as geomorphology, [paleo]climatology, oceanography and biodiversity. In this work I present a comparative assessment of the datasets ETOPO1 (1' resolution), GTOPO30, GLOBE, SRTM30 PLUS, GMTED2010 and ACE2 (30") against the altitude of the world's ultra prominent peaks. GDEMs' elevations show an expected tendency of underestimating the peak's altitude, but differences reach 3,500 m. None of the GDEMs captures the full range of elevation on Earth and they do not represent well the altitude of the most prominent peaks. Some of these problems could be addressed with the release of NASADEM, but the smoothing effect caused by moving-window resampling can only be tackled by using new techniques, such as scale-adaptative kernels and curvature-based terrain generalisation.

  16. Lossless data compression of grid-based digital elevation models: A png image format evaluation

    NASA Astrophysics Data System (ADS)

    Scarmana, G.

    2014-05-01

    At present, computers, lasers, radars, planes and satellite technologies make possible very fast and accurate topographic data acquisition for the production of maps. However, the problem of managing and manipulating this data efficiently remains. One particular type of map is the elevation map. When stored on a computer, it is often referred to as a Digital Elevation Model (DEM). A DEM is usually a square matrix of elevations. It is like an image, except that it contains a single channel of information (that is, elevation) and can be compressed in a lossy or lossless manner by way of existing image compression protocols. Compression has the effect of reducing memory requirements and speed of transmission over digital links, while maintaining the integrity of data as required. In this context, this paper investigates the effects of the PNG (Portable Network Graphics) lossless image compression protocol on floating-point elevation values for 16-bit DEMs of dissimilar terrain characteristics. The PNG is a robust, universally supported, extensible, lossless, general-purpose and patent-free image format. Tests demonstrate that the compression ratios and run decompression times achieved with the PNG lossless compression protocol can be comparable to, or better than, proprietary lossless JPEG variants, other image formats and available lossless compression algorithms.

  17. Appending High-Resolution Elevation Data to GPS Speed Traces for Vehicle Energy Modeling and Simulation

    SciTech Connect

    Wood, E.; Burton, E.; Duran, A.; Gonder, J.

    2014-06-01

    Accurate and reliable global positioning system (GPS)-based vehicle use data are highly valuable for many transportation, analysis, and automotive considerations. Model-based design, real-world fuel economy analysis, and the growing field of autonomous and connected technologies (including predictive powertrain control and self-driving cars) all have a vested interest in high-fidelity estimation of powertrain loads and vehicle usage profiles. Unfortunately, road grade can be a difficult property to extract from GPS data with consistency. In this report, we present a methodology for appending high-resolution elevation data to GPS speed traces via a static digital elevation model. Anomalous data points in the digital elevation model are addressed during a filtration/smoothing routine, resulting in an elevation profile that can be used to calculate road grade. This process is evaluated against a large, commercially available height/slope dataset from the Navteq/Nokia/HERE Advanced Driver Assistance Systems product. Results will show good agreement with the Advanced Driver Assistance Systems data in the ability to estimate road grade between any two consecutive points in the contiguous United States.

  18. Subreflector model depending on elevation for the Tianma 65m Radio Telescope

    NASA Astrophysics Data System (ADS)

    Sun, Zheng-Xiong; Wang, Jin-Qing; Chen, Lan

    2016-08-01

    A subreflector adjustment system for the Tianma 65 m radio telescope, administered by Shanghai Astronomical Observatory, has been installed to compensate for gravitational deformation of the main reflector and the structure supporting the subreflector. The position and attitude of the subreflector are variable in order to improve the efficiency at different elevations. The subreflector model has the goal of improving the antenna's performance. A new fitting formulation which is different from the traditional formulation is proposed to reduce the fitting error in the Y direction. The only difference in the subreflector models of the 65m radio telescope is the bias of a constant term in the Z direction. We have investigated the effect of movements of the subreflector on the pointing of the antenna. The results of these performance measurements made by moving the antenna in elevation show that the subreflector model can effectively improve the efficiency of the 65 m radio telescope at each elevation. An antenna efficiency of about 60% at the Ku band is reached in the whole angular range of elevation.

  19. Initial Everglades Depth Estimation Network (EDEN) Digital Elevation Model Research and Development

    USGS Publications Warehouse

    Jones, John W.; Price, Susan D.

    2007-01-01

    Introduction The Everglades Depth Estimation Network (EDEN) offers a consistent and documented dataset that can be used to guide large-scale field operations, to integrate hydrologic and ecological responses, and to support biological and ecological assessments that measure ecosystem responses to the Comprehensive Everglades Restoration Plan (Telis, 2006). To produce historic and near-real time maps of water depths, the EDEN requires a system-wide digital elevation model (DEM) of the ground surface. Accurate Everglades wetland ground surface elevation data were non-existent before the U.S. Geological Survey (USGS) undertook the collection of highly accurate surface elevations at the regional scale. These form the foundation for EDEN DEM development. This development process is iterative as additional high accuracy elevation data (HAED) are collected, water surfacing algorithms improve, and additional ground-based ancillary data become available. Models are tested using withheld HAED and independently measured water depth data, and by using DEM data in EDEN adaptive management applications. Here the collection of HAED is briefly described before the approach to DEM development and the current EDEN DEM are detailed. Finally future research directions for continued model development, testing, and refinement are provided.

  20. Anatomy of a sudden onset flood: The 18 March 2007 Crater Lake break-out lahar, Mt. Ruapehu, New Zealand

    NASA Astrophysics Data System (ADS)

    Manville, Vern

    2010-05-01

    Sudden onset floods are significant hazards in many environments and regions around the world, resulting in loss of life, damage to infrastructure, and dramatic geomorphic changes due to the very high rate of energy expenditure associated with high flow velocities and depths. In these transient, highly dynamic phenomena, the evolving floodwave typically undergoes a complex series of interactions with the flowpath, involving multiple feedback loops, to produce spatially and temporarily varying sediment loads that affect their density, viscosity, mobility, peak discharge and absolute volume, and hence their hazardousness. However, due to their unpredictability, ephemeral nature, and high energy they are challenging to characterise adequately, making mitigation difficult. At Mt. Ruapehu, New Zealand, an outburst flood caused by collapse of a tephra barrier impounding a summit Crater Lake generated an opportunity to characterise the downstream evolution of a single discrete floodwave through: (i) multi-parameter measurement of time-series hydraulic parameters at key locations along the channel; (ii) visual observations at manual sampling and gauging sites; (iii) capture of geomorphic changes using pre- and post-event high resolution topographic surveys and vertical aerial and oblique imagery; (iv) stratigraphic logging and granulometric analyses of the flood deposits; and (v) testing and calibration of numerical models using the newly acquired data. On 18 March 2007, the natural dam failed, releasing c. 1 Mm3 of warm, highly mineralised water in less than two hours. The flood rapidly bulked by entraining debris along the steep gorge of the upper Whangaehu valley, more than tripling in size, before debouching onto the Whangaehu Fan where it braided into multiple distributary channels. The flood recollected into a single channel to continue downstream to the coast 215 km from source. Stage, discharge, frontal and flow velocity, sediment-load, and geochemical data show

  1. Perturbation and melting of snow and ice by the 13 November 1985 eruption of Nevado del Ruiz, Colombia, and consequent mobilization, flow and deposition of lahars

    USGS Publications Warehouse

    Pierson, T.C.; Janda, R.J.; Thouret, J.-C.; Borrero, C.A.

    1990-01-01

    A complex sequence of pyroclastic flows and surges erupted by Nevado del Ruiz volcano on 13 November 1985 interacted with snow and ice on the summit ice cap to trigger catastrophic lahars (volcanic debris flows), which killed more than 23,000 people living at or beyond the base of the volcano. The rapid transfer of heat from the hot eruptive products to about 10 km2 of the snowpack, combined with seismic shaking, produced large volumes of meltwater that flowed downslope, liquefied some of the new volcanic deposits, and generated avalanches of saturated snow, ice and rock debris within minutes of the 21:08 (local time) eruption. About 2 ?? 107 m3 of water was discharged into the upper reaches of the Molinos, Nereidas, Guali, Azufrado and Lagunillas valleys, where rapid entrainment of valley-fill sediment transformed the dilute flows and avalanches to debris flows. Computed mean velocities of the lahars at peak flow ranged up to 17 m s-1. Flows were rapid in the steep, narrow upper canyons and slowed with distance away from the volcano as flow depth and channel slope diminished. Computed peak discharges ranged up to 48,000 m3 s-1 and were greatest in reaches 10 to 20 km downstream from the summit. A total of about 9 ?? 107 m3 of lahar slurry was transported to depositional areas up to 104 km from the source area. Initial volumes of individual lahars increased up to 4 times with distance away from the summit. The sedimentology and stratigraphy of the lahar deposits provide compelling evidence that: (1) multiple initial meltwater pulses tended to coalesce into single flood waves; (2) lahars remained fully developed debris flows until they reached confluences with major rivers; and (3) debris-flow slurry composition and rheology varied to produce gradationally density-stratified flows. Key lessons and reminders from the 1985 Nevado del Ruiz volcanic eruption are: (1) catastrophic lahars can be generated on ice- and snow-capped volcanoes by relatively small eruptions; (2

  2. In-flow evolution of lahar deposits from video-imagery with implications for post-event deposit interpretation, Mount Semeru, Indonesia

    NASA Astrophysics Data System (ADS)

    Starheim, Colette C. A.; Gomez, Christopher; Davies, Tim; Lavigne, Franck; Wassmer, Patrick

    2013-04-01

    The hazardous and unpredictable nature of lahars makes them challenging to study, yet the in-flow processes characterizing these events are important to understand. As a result, much of the previous research on lahar sedimentation and flow processes has been derived from experimental flows or stratigraphic surveys of post-event deposits. By comparison, little is known on the time-dependent sediment and flow dynamics of lahars in natural environments. Using video-footage of seven lahars on the flanks of Semeru Volcano (East Java, Indonesia), the present study offers new insights on the in-flow evolution of sediment in natural lahars. Video analysis revealed several distinctive patterns of sediment entrainment and deposition that varied with time-related fluctuations in flow. These patterns were used to generate a conceptual framework describing possible processes of formation for subsurface architectural features identified in an earlier lateral survey of lahar deposits on Semeru Volcano (Gomez and Lavigne, 2010a). The formation of lateral discontinuities was related to the partial erosion of transitional bank deposits followed by fresh deposition along the erosional contact. This pattern was observed over the course of several lahar events and within individual flows. Observations similarly offer potential explanations for the formation of lenticular features. Depending on flow characteristics, these features appeared to form by preferential erosion or deposition around large stationary blocks, and by deposition along channel banks during episodes of channel migration or channel constriction. Finally, conditions conducive to the deposition of fine laminated beds were observed during periods of attenuating and surging flow. These results emphasize the difficulties associated with identifying process-structure relationships solely from post-event deposit interpretation and illustrate that an improved understanding of the time-dependent sediment dynamics in lahars may

  3. Elevating your elevator talk

    Technology Transfer Automated Retrieval System (TEKTRAN)

    An important and often overlooked item that every early career researcher needs to do is compose an elevator talk. The elevator talk, named because the talk should not last longer than an average elevator ride (30 to 60 seconds), is an effective method to present your research and yourself in a clea...

  4. Comparison of Surface Flow Features from Lidar-Derived Digital Elevation Models with Historical Elevation and Hydrography Data for Minnehaha County, South Dakota

    USGS Publications Warehouse

    Poppenga, Sandra K.; Worstell, Bruce B.; Stoker, Jason M.; Greenlee, Susan K.

    2009-01-01

    The U.S. Geological Survey (USGS) has taken the lead in the creation of a valuable remote sensing product by incorporating digital elevation models (DEMs) derived from Light Detection and Ranging (lidar) into the National Elevation Dataset (NED), the elevation layer of 'The National Map'. High-resolution lidar-derived DEMs provide the accuracy needed to systematically quantify and fully integrate surface flow including flow direction, flow accumulation, sinks, slope, and a dense drainage network. In 2008, 1-meter resolution lidar data were acquired in Minnehaha County, South Dakota. The acquisition was a collaborative effort between Minnehaha County, the city of Sioux Falls, and the USGS Earth Resources Observation and Science (EROS) Center. With the newly acquired lidar data, USGS scientists generated high-resolution DEMs and surface flow features. This report compares lidar-derived surface flow features in Minnehaha County to 30- and 10-meter elevation data previously incorporated in the NED and ancillary hydrography datasets. Surface flow features generated from lidar-derived DEMs are consistently integrated with elevation and are important in understanding surface-water movement to better detect surface-water runoff, flood inundation, and erosion. Many topographic and hydrologic applications will benefit from the increased availability of accurate, high-quality, and high-resolution surface-water data. The remotely sensed data provide topographic information and data integration capabilities needed for meeting current and future human and environmental needs.

  5. A Seamless, High-Resolution, Coastal Digital Elevation Model (DEM) for Southern California

    USGS Publications Warehouse

    Barnard, Patrick L.; Hoover, Daniel

    2010-01-01

    A seamless, 3-meter digital elevation model (DEM) was constructed for the entire Southern California coastal zone, extending 473 km from Point Conception to the Mexican border. The goal was to integrate the most recent, high-resolution datasets available (for example, Light Detection and Ranging (Lidar) topography, multibeam and single beam sonar bathymetry, and Interferometric Synthetic Aperture Radar (IfSAR) topography) into a continuous surface from at least the 20-m isobath to the 20-m elevation contour. This dataset was produced to provide critical boundary conditions (bathymetry and topography) for a modeling effort designed to predict the impacts of severe winter storms on the Southern California coast (Barnard and others, 2009). The hazards model, run in real-time or with prescribed scenarios, incorporates atmospheric information (wind and pressure fields) with a suite of state-of-the-art physical process models (tide, surge, and wave) to enable detailed prediction of water levels, run-up, wave heights, and currents. Research-grade predictions of coastal flooding, inundation, erosion, and cliff failure are also included. The DEM was constructed to define the general shape of nearshore, beach and cliff surfaces as accurately as possible, with less emphasis on the detailed variations in elevation inland of the coast and on bathymetry inside harbors. As a result this DEM should not be used for navigation purposes.

  6. Evaluating LMA and CLAMP: Using information criteria to choose a model for estimating elevation

    NASA Astrophysics Data System (ADS)

    Miller, I.; Green, W.; Zaitchik, B.; Brandon, M.; Hickey, L.

    2005-12-01

    The morphology of leaves and composition of the flora respond strongly to the moisture and temperature of their environment. Elevation and latitude correlate, at first order, to these atmospheric parameters. An obvious modern example of this relationship between leaf morphology and environment is the tree line, where boreal forests give way to artic (high latitude) or alpine (high elevation) tundra. Several quantitative methods, all of which rely on uniformitarianism, have been developed to estimate paleoelevation using fossil leaf morphology. These include 1) the univariate leaf-margin analysis (LMA), which estimates mean annual temperature (MAT) by the positive linear correlation between MAT and P, the proportion of entire or smooth to non-entire or toothed margined woody dicot angiosperm leaves within a flora and 2) the Climate Leaf Analysis Multivariate Program (CLAMP) which uses Canonical Correspondence Analysis (CCA) to estimate MAT, moist enthalpy, and other atmospheric parameters using 31 explanatory leaf characters from woody dicot angiosperms. Given a difference in leaf-estimated MAT or moist enthalpy between contemporaneous, synlatitudinal fossil floras-one at sea-level, the other at an unknown paleoelevation-paleoelevation may be estimated. These methods have been widely applied to orogenic settings and concentrate particularly in the Western US. We introduce the use of information criteria to compare different models for estimating elevation and show how the additional complexity of the CLAMP analytical methodology does not necessarily improve on the elevation estimates produced by simpler regression models. In addition, we discuss the signal-to-noise ratio in the data, give confidence intervals for detecting elevations, and address the problem of spatial autocorrelation and irregular sampling in the data.

  7. A comparative study of spherical and flat-Earth geopotential modeling at satellite elevations

    NASA Technical Reports Server (NTRS)

    Parrott, M. H.; Hinze, W. J.; Braile, L. W.

    1985-01-01

    Flat-Earth and spherical-Earth geopotential modeling of crustal anomaly sources at satellite elevations are compared by computing gravity and scalar magnetic anomalies perpendicular to the strike of variably dimensioned rectangular prisms at altitudes of 150, 300, and 450 km. Results indicate that the error caused by the flat-Earth approximation is less than 10% in most geometric conditions. Generally, error increase with larger and wider anomaly sources at higher altitudes. For most crustal source modeling applications at conventional satellite altitudes, flat-Earth modeling can be justified and is numerically efficient.

  8. A geomorphology based approach for digital elevation model fusion - case study in Danang City, Vietnam

    NASA Astrophysics Data System (ADS)

    Tran, T. A.; Raghavan, V.; Masumoto, S.; Vinayaraj, P.; Yonezawa, G.

    2014-04-01

    Global Digital Elevation Model (DEM) is considered as vital spatial information and finds wide use in several applications. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global DEM (GDEM) and Shuttle Radar Topographic Mission (SRTM) DEM offer almost global coverage and provide elevation data for geospatial analysis. However, GDEM and SRTM still contain some height errors that affect the quality of elevation data significantly. This study aims to examine methods to improve the resolution as well as accuracy of available free DEMs by data fusion technique and evaluating the results with high quality reference DEM. The DEM fusion method is based on the accuracy assessment of each global DEM and geomorphological characteristics of the study area. Land cover units were also considered to correct the elevation of GDEM and SRTM with respect to the bare earth surface. Weighted averaging method was used to fuse the input DEMs based on landform classification map. According to the landform types, the different weights were used for GDEM and SRTM. Finally, a denoising algorithm (Sun et al., 2007) was applied to filter the output fused DEM. This fused DEM shows excellent correlation to the reference DEM having correlation coefficient R2 = 0.9986 and the accuracy was also improved from Root Mean Square Error (RMSE) 14.9 m in GDEM and 14.8 m in SRTM into 11.6 m in fused DEM.

  9. Digital Elevation Models of TYCHO Crater and the Lunar Polar Regions

    NASA Astrophysics Data System (ADS)

    Margot, J. L.; Campbell, D. B.; Jurgens, R. F.; Slade, M. A.

    1998-09-01

    Earth-based radar interferometry [1] has been used to map the lunar polar regions and Tycho Crater at high spatial ( ~ 100 m) and height ( ~ 50 m) resolutions. Compared to existing topographic data sets, the radar observations offer digital elevation models with dense horizontal spacing and improved height resolution. Earth-based radars can also provide measurements of the largely unknown topography in the polar regions. Elevation data and radar imagery obtained with the Goldstone X-band system (lambda = 3.5 cm) are presented for the Tycho Crater area, with a spatial resolution of 200 m and a height resolution of 30 m. A careful comparison of the radar-derived topography with Clementine altimetry points [2] reveals a very good agreement between the two techniques. Rms deviations between the radar-derived heights and 87 Clementine points available over the 200 x 200 km scene are ~ 100 m. The digital elevation model allows detailed morphometry of the 85 km diameter crater: the floor of Tycho lies 3970 m below a 1738 km radius sphere, and the crater's central peak rises 2400 m above the floor. The average rim crest elevation is 730 m above the 1738 km datum, giving a mean rim to floor depth of 4700 m. The floor has two distinct units with the western section being higher in elevation by ~ 200 m. This dichotomy is consistent with an asymmetry in the crater shape which reveals that maximum wall slumping occured in the western and southwestern regions of the crater. Digital elevation models of the polar regions are being used to estimate the location of permanently shadowed areas which may harbor ice deposits [3]. The range of illumination conditions over the lunar polar regions could be sampled by an imaging instrument in a polar orbit during a full terrestrial year. Alternatively, topographic maps obtained with Earth-based radar can be used to model the illumination conditions over the entire solar illumination cycle. [1] I. I. Shapiro et al. (1972). Science, 178, 939

  10. Modelling Periglacial Processes on Low-Relief High-Elevation Surfaces

    NASA Astrophysics Data System (ADS)

    Andersen, J. L.; Knudsen, M. F.; Egholm, D. L.

    2014-12-01

    Are low-relief high-elevation surfaces generally a result of uplift of flat surfaces formed close to sea-level or can they be formed "in situ" by climate dependent surface processes such as those associated with glaciation? This question is important to resolve in order to understand the geological history in many regions of the world. The glacial buzzsaw concept suggests that intense glacial erosion focused at the equilibrium-line altitude (ELA) leads to a concentration in surface area close to the ELA. However, even in predominantly glacial landscapes, such as the Scandinavian Mountains, the high surfaces often have nonglacial characteristics and show pre-glacial inheritance of cosmogenic nuclide concentrations. This suggests that subglacial erosion is not the dominant erosion mechanism. Anderson (2002) showed that periglacial frost-driven processes working over timescales of millions of years form low-curvature parabolic surfaces at accordant elevation across intervening valleys. The rate of these erosion processes is limited by the slow diffusive transport of regolith in the periglacial domain. We elaborate on this idea by quantifying frost cracking and frost creep in a numerical model as a function of mean annual air temperature and sediment thickness. This allows us to incorporate periglacial processes into a long-term landscape evolution model where surface elevation, sediment thickness, and climate evolve over time. With this model we are able to explore the slow feedbacks between periglacial erosion, sediment transport, and the evolving topography. We show that smooth peaks, convex hillslopes, and a few meters thick regolith cover at high elevation are emergent properties of the landscape evolution model. By varying climate and other model parameters, we discuss how the landscape evolution model can be used for obtaining more insight into the conditions needed for formation of low-relief surfaces at high elevation. Anderson, R. S. Modeling the tor

  11. Modeling the population dynamics of Culex quinquefasciatus (Diptera: Culcidae), along an elevational gradient in Hawaii

    USGS Publications Warehouse

    Ahumada, Jorge A.; LaPointe, Dennis; Samuel, Michael D.

    2004-01-01

    We present a population model to understand the effects of temperature and rainfall on the population dynamics of the southern house mosquito, Culex quinquefasciatus Say, along an elevational gradient in Hawaii. We use a novel approach to model the effects of temperature on population growth by dynamically incorporating developmental rate into the transition matrix, by using physiological ages of immatures instead of chronological age or stages. We also model the effects of rainfall on survival of immatures as the cumulative number of days below a certain rain threshold. Finally, we incorporate density dependence into the model as competition between immatures within breeding sites. Our model predicts the upper altitudinal distributions of Cx. quinquefasciatus on the Big Island of Hawaii for self-sustaining mosquito and migrating summer sink populations at 1,475 and 1,715 m above sea level, respectively. Our model predicts that mosquitoes at lower elevations can grow under a broader range of rainfall parameters than middle and high elevation populations. Density dependence in conjunction with the seasonal forcing imposed by temperature and rain creates cycles in the dynamics of the population that peak in the summer and early fall. The model provides a reasonable fit to the available data on mosquito abundance for the east side of Mauna Loa, Hawaii. The predictions of our model indicate the importance of abiotic conditions on mosquito dynamics and have important implications for the management of diseases transmitted by Cx. quinquefasciatus in Hawaii and elsewhere.

  12. Modeling glaucoma in rats by sclerosing aqueous outflow pathways to elevate intraocular pressure.

    PubMed

    Morrison, John C; Cepurna, William O; Johnson, Elaine C

    2015-12-01

    Injection of hypertonic saline via episcleral veins toward the limbus in laboratory rats can produce elevated intraocular pressure (IOP) by sclerosis of aqueous humor outflow pathways. This article describes important anatomic characteristics of the rat optic nerve head (ONH) that make it an attractive animal model for human glaucoma, along with the anatomy of rat aqueous humor outflow on which this technique is based. The injection technique itself is also described, with the aid of a supplemental movie, including necessary equipment and specific tips to acquire this skill. Outcomes of a successful injection are presented, including IOP elevation and patterns of optic nerve injury. These concepts are then specifically considered in light of the use of this model to assess potential neuroprotective therapies. Advantages of the hypertonic saline model include a delayed and relatively gradual IOP elevation, likely reproduction of scleral and ONH stresses and strains that may be important in producing axonal injury, and its ability to be applied to any rat (and potentially mouse) strain, leaving the unmanipulated fellow eye as an internal control. Challenges include the demanding surgical skill required by the technique itself, a wide range of IOP response, and mild corneal clouding in some animals. However, meticulous application of the principles detailed in this article and practice will allow most researchers to attain this useful skill for studying cellular events of glaucomatous optic nerve damage. PMID:26003399

  13. Impact of land-surface elevation and riparian evapotranspiration seasonality on groundwater budget in MODFLOW models

    NASA Astrophysics Data System (ADS)

    Ajami, Hoori; Meixner, Thomas; Maddock, Thomas; Hogan, James F.; Guertin, D. Phillip

    2011-09-01

    Riparian groundwater evapotranspiration (ETg) constitutes a major component of the water balance especially in many arid and semi-arid environments. Although spatial and temporal variability of riparian ETg are controlled by climate, vegetation and subsurface characteristics, depth to water table (DTWT) is often considered the major controlling factor. Relationships between ETg rates and DTWT, referred to as ETg curves, are implemented in MODFLOW ETg packages (EVT, ETS1 and RIP-ET) with different functional forms. Here, the sensitivity of the groundwater budget in MODFLOW groundwater models to ETg parameters (including ETg curves, land-surface elevation and ETg seasonality) are investigated. A MODFLOW model of the hypothetical Dry Alkaline Valley in the Southwestern USA is used to show how spatial representation of riparian vegetation and digital elevation model (DEM) processing methods impact the water budget when RIPGIS-NET (a GIS-based ETg program) is used with MODFLOW's RIP-ET package, and results are compared with the EVT and ETS1 packages. Results show considerable impact on ETg and other groundwater budget components caused by spatial representation of riparian vegetation, vegetation type, fractional coverage areas and land-surface elevation. RIPGIS-NET enhances ETg estimation in MODFLOW by incorporating vegetation and land-surface parameters, providing a tool for ecohydrology studies, riparian ecosystem management and stream restoration.

  14. Application of a two-pool model to soil carbon dynamics under elevated CO2.

    PubMed

    van Groenigen, Kees Jan; Xia, Jianyang; Osenberg, Craig W; Luo, Yiqi; Hungate, Bruce A

    2015-12-01

    Elevated atmospheric CO2 concentrations increase plant productivity and affect soil microbial communities, with possible consequences for the turnover rate of soil carbon (C) pools and feedbacks to the atmosphere. In a previous analysis (Van Groenigen et al., 2014), we used experimental data to inform a one-pool model and showed that elevated CO2 increases the decomposition rate of soil organic C, negating the storage potential of soil. However, a two-pool soil model can potentially explain patterns of soil C dynamics without invoking effects of CO2 on decomposition rates. To address this issue, we refit our data to a two-pool soil C model. We found that CO2 enrichment increases decomposition rates of both fast and slow C pools. In addition, elevated CO2 decreased the carbon use efficiency of soil microbes (CUE), thereby further reducing soil C storage. These findings are consistent with numerous empirical studies and corroborate the results from our previous analysis. To facilitate understanding of C dynamics, we suggest that empirical and theoretical studies incorporate multiple soil C pools with potentially variable decomposition rates. PMID:26313640

  15. Perspectives on open access high resolution digital elevation models to produce global flood hazard layers

    NASA Astrophysics Data System (ADS)

    Sampson, Christopher; Smith, Andrew; Bates, Paul; Neal, Jeffrey; Trigg, Mark

    2015-12-01

    Global flood hazard models have recently become a reality thanks to the release of open access global digital elevation models, the development of simplified and highly efficient flow algorithms, and the steady increase in computational power. In this commentary we argue that although the availability of open access global terrain data has been critical in enabling the development of such models, the relatively poor resolution and precision of these data now limit significantly our ability to estimate flood inundation and risk for the majority of the planet's surface. The difficulty of deriving an accurate 'bare-earth' terrain model due to the interaction of vegetation and urban structures with the satellite-based remote sensors means that global terrain data are often poorest in the areas where people, property (and thus vulnerability) are most concentrated. Furthermore, the current generation of open access global terrain models are over a decade old and many large floodplains, particularly those in developing countries, have undergone significant change in this time. There is therefore a pressing need for a new generation of high resolution and high vertical precision open access global digital elevation models to allow significantly improved global flood hazard models to be developed.

  16. Elevated MicroRNA-33 in Sarcoidosis and a Carbon Nanotube Model of Chronic Granulomatous Disease.

    PubMed

    Barna, Barbara P; McPeek, Matthew; Malur, Anagha; Fessler, Michael B; Wingard, Christopher J; Dobbs, Larry; Verbanac, Kathryn M; Bowling, Mark; Judson, Marc A; Thomassen, Mary Jane

    2016-06-01

    We established a murine model of multiwall carbon nanotube (MWCNT)-induced chronic granulomatous disease, which resembles human sarcoidosis pathology. At 60 days after oropharyngeal MWCNT instillation, bronchoalveolar lavage (BAL) cells from wild-type mice exhibit an M1 phenotype with elevated proinflammatory cytokines and reduced peroxisome proliferator-activated receptor γ (PPARγ)-characteristics also present in human sarcoidosis. Based upon MWCNT-associated PPARγ deficiency, we hypothesized that the PPARγ target gene, ATP-binding cassette (ABC) G1, a lipid transporter with antiinflammatory properties, might also be repressed. Results after MWCNT instillation indicated significantly repressed ABCG1, but, surprisingly, lipid transporter ABCA1 was also repressed, suggesting a possible second pathway. Exploration of potential regulators revealed that microRNA (miR)-33, a lipid transporter regulator, was strikingly elevated (13.9 fold) in BAL cells from MWCNT-instilled mice but not sham control mice. Elevated miR-33 was also detected in murine granulomatous lung tissue. In vitro studies confirmed that lentivirus-miR-33 overexpression repressed both ABCA1 and ABCG1 (but not PPARγ) in cultured murine alveolar macrophages. BAL cells of patients with sarcoidosis also displayed elevated miR-33 together with reduced ABCA1 and ABCG1 messenger RNA and protein compared with healthy control subjects. Moreover, miR-33 was elevated within sarcoidosis granulomatous tissue. The findings suggest that alveolar macrophage miR-33 is up-regulated by proinflammatory cytokines and may perpetuate chronic inflammatory granulomatous disease by repressing antiinflammatory functions of ABCA1 and ABCG1 lipid transporters. The results also suggest two possible pathways for transporter dysregulation in granulomatous disease-one associated with intrinsic PPARγ status and the other with miR-33 up-regulation triggered by environmental challenges, such as MWCNT. PMID:26641802

  17. a New High-Resolution Elevation Model of Greenland Derived from Tandem-X

    NASA Astrophysics Data System (ADS)

    Wessel, B.; Bertram, A.; Gruber, A.; Bemm, S.; Dech, S.

    2016-06-01

    In this paper we present for the first time the new digital elevation model (DEM) for Greenland produced by the TanDEM-X (TerraSAR add-on for digital elevation measurement) mission. The new, full coverage DEM of Greenland has a resolution of 0.4 arc seconds corresponding to 12 m. It is composed of more than 7.000 interferometric synthetic aperture radar (InSAR) DEM scenes. X-Band SAR penetrates the snow and ice pack by several meters depending on the structures within the snow, the acquisition parameters, and the dielectricity constant of the medium. Hence, the resulting SAR measurements do not represent the surface but the elevation of the mean phase center of the backscattered signal. Special adaptations on the nominal TanDEM-X DEM generation are conducted to maintain these characteristics and not to raise or even deform the DEM to surface reference data. For the block adjustment, only on the outer coastal regions ICESat (Ice, Cloud, and land Elevation Satellite) elevations as ground control points (GCPs) are used where mostly rock and surface scattering predominates. Comparisons with ICESat data and snow facies are performed. In the inner ice and snow pack, the final X-Band InSAR DEM of Greenland lies up to 10 m below the ICESat measurements. At the outer coastal regions it corresponds well with the GCPs. The resulting DEM is outstanding due to its resolution, accuracy and full coverage. It provides a high resolution dataset as basis for research on climate change in the arctic.

  18. Elastic model for the gravity and elevation changes before the 2001 eruption of Etna volcano

    NASA Astrophysics Data System (ADS)

    Carbone, Daniele; Currenti, Gilda; Del Negro, Ciro

    2007-03-01

    For 5 months before the 2001 Mt. Etna eruption, a progressive gravity decrease was measured along a profile of stations on the southern slope of the volcano. Between January and July 2001, the amplitude of the change reached 80 μGal, while the wavelength of the anomaly was of the order of 15 km. Elevation changes observed through GPS measurements during a period encompassing the 5-month gravity decrease, remained within 4 6 cm over the entire volcano and within 2 4 cm in the zone covered by the microgravity profile. We review both gravity and elevation changes by a model assuming the formation of new cracks, uniformly distributed in a rectangular prism. The inversion problem was formulated following a global optimization approach based on the use of Genetic Algorithms. Although it is possible to explain the observed gravity changes by means of the proposed analytical formulation, the results show that calculated elevation changes are significantly higher than those observed. Two alternative hypotheses are proposed to account for this apparent discrepancy: (1) that the assumptions behind the analytical formulation, used to invert the data, are fallacious at Etna, and thus, numerical models should be utilized; (2) that a second process, enabling a considerable mass decrease to occur without deformation, acted together with the formation of new cracks in the source volume.

  19. Spatial disaggregation of satellite-derived irradiance using a high-resolution digital elevation model

    SciTech Connect

    Ruiz-Arias, Jose A.; Tovar-Pescador, Joaquin; Cebecauer, Tomas; Suri, Marcel

    2010-09-15

    Downscaling of the Meteosat-derived solar radiation ({proportional_to}5 km grid resolution) is based on decomposing the global irradiance and correcting the systematic bias of its components using the elevation and horizon shadowing that are derived from the SRTM-3 digital elevation model (3 arc sec resolution). The procedure first applies the elevation correction based on the difference between coarse and high spatial resolution. Global irradiance is split into direct, diffuse circumsolar and diffuse isotropic components using statistical models, and then corrections due to terrain shading and sky-view fraction are applied. The effect of reflected irradiance is analysed only in the theoretical section. The method was applied in the eastern Andalusia, Spain, and the validation was carried out for 22 days on April, July and December 2006 comparing 15-min estimates of the satellite-derived solar irradiance and observations from nine ground stations. Overall, the corrections of the satellite estimates in the studied region strongly reduced the mean bias of the estimates for clear and cloudy days from roughly 2.3% to 0.4%. (author)

  20. 2009 ERUPTION OF REDOUBT VOLCANO: Lahars, Oil, and the Role of Science in Hazards Mitigation (Invited)

    NASA Astrophysics Data System (ADS)

    Swenson, R.; Nye, C. J.

    2009-12-01

    In March, 2009, Redoubt Volcano erupted for the third time in 45 years. More than 19 explosions produced ash plumes to 60,000 ft asl, lahar flows of mud and ice down the Drift river ~30 miles to the coast, and tephra fall up to 1.5 mm onto surrounding communities. The eruption had severe impact on many operations. Airlines were forced to cancel or divert hundreds of international and domestic passenger and cargo flights, and Anchorage International airport closed for over 12 hours. Mudflows and floods down the Drift River to the coast impacted operations at the Drift River Oil Terminal (DROT) which was forced to shut down and ultimately be evacuated. Prior mitigation efforts to protect the DROT oil tank farm from potential impacts associated with a major eruptive event were successful, and none of the 148,000 barrels of oil stored at the facility was spilled or released. Nevertheless, the threat of continued eruptive activity at Redoubt, with the possibility of continued lahar flows down the Drift River alluvial fan, required an incident command post be established so that the US Coast Guard, Alaska Dept. of Environmental Conservation, and the Cook Inlet Pipeline Company could coordinate a response to the potential hazards. Ultimately, the incident command team relied heavily on continuous real-time data updates from the Alaska Volcano Observatory, as well as continuous geologic interpretations and risk analysis by the USGS Volcanic Hazards group, the State Division of Geological and Geophysical Surveys and the University of Alaska Geophysical Institute, all members of the collaborative effort of the Alaska Volcano Observatory. The great success story that unfolded attests to the efforts of the incident command team, and their reliance on real-time scientific analysis from scientific experts. The positive results also highlight how pre-disaster mitigation and monitoring efforts, in concert with hazards response planning, can be used in a cooperative industry

  1. The enormous Chillos Valley Lahar: An ash-flow-generated debris flow from Cotopaxi Volcano, Ecuador

    USGS Publications Warehouse

    Mothes, P.A.; Hall, M.L.; Janda, R.J.

    1998-01-01

    The Chillos Valley Lahar (CVL), the largest Holocene debris flow in area and volume as yet recognized in the northern Andes, formed on Cotopaxi volcano's north and northeast slopes and descended river systems that took it 326 km north-northwest to the Pacific Ocean and 130+ km east into the Amazon basin. In the Chillos Valley, 40 km downstream from the volcano, depths of 80-160 m and valley cross sections up to 337000m2 are observed, implying peak flow discharges of 2.6-6.0 million m3/s. The overall volume of the CVL is estimated to be ???3.8 km3. The CVL was generated approximately 4500 years BP by a rhyolitic ash flow that followed a small sector collapse on the north and northeast sides of Cotopaxi, which melted part of the volcano's icecap and transformed rapidly into the debris flow. The ash flow and resulting CVL have identical components, except for foreign fragments picked up along the flow path. Juvenile materials, including vitric ash, crystals, and pumice, comprise 80-90% of the lahar's deposit, whereas rhyolitic, dacitic, and andesitic lithics make up the remainder. The sand-size fraction and the 2- to 10-mm fraction together dominate the deposit, constituting ???63 and ???15 wt.% of the matrix, respectively, whereas the silt-size fraction averages less than ???10 wt.% and the clay-size fraction less than 0.5 wt.%. Along the 326-km runout, these particle-size fractions vary little, as does the sorting coefficient (average = 2.6). There is no tendency toward grading or improved sorting. Limited bulking is recognized. The CVL was an enormous non-cohesive debris flow, notable for its ash-flow origin and immense volume and peak discharge which gave it characteristics and a behavior akin to large cohesive mudflows. Significantly, then, ash-flow-generated debris flows can also achieve large volumes and cover great areas; thus, they can conceivably affect large populated regions far from their source. Especially dangerous, therefore, are snowclad volcanoes

  2. A flooding algorithm for extracting drainage networks from unprocessed digital elevation models

    NASA Astrophysics Data System (ADS)

    Rueda, Antonio; Noguera, José M.; Martínez-Cruz, Carmen

    2013-09-01

    A new method for extracting the drainage network from a digital elevation model (DEM) is presented. It is based on the well-known D8 approach that simulates the overland flow but uses a more elaborate water transfer model that is inspired by the natural behaviour of water. The proposed solution has several advantages: it works on unprocessed DEMs avoiding the problems caused by pits and flats, can generate watercourses with a width greater than one cell and detects fluvial landforms like lakes, marshes or river islands that are not directly handled by most previous solutions.

  3. Analysis of accuracy of digital elevation models created from captured data by digital photogrammetry method

    NASA Astrophysics Data System (ADS)

    Hudec, P.

    2011-12-01

    A digital elevation model (DEM) is an important part of many geoinformatic applications. For the creation of DEM, spatial data collected by geodetic measurements in the field, photogrammetric processing of aerial survey photographs, laser scanning and secondary sources (analogue maps) are used. It is very important from a user's point of view to know the vertical accuracy of a DEM. The article describes the verification of the vertical accuracy of a DEM for the region of Medzibodrožie, which was created using digital photogrammetry for the purposes of water resources management and modeling and resolving flood cases based on geodetic measurements in the field.

  4. Modeling Laser Altimeter Return Waveforms Over Complex Vegetation Using High-Resolution Elevation Data

    NASA Technical Reports Server (NTRS)

    Blair, J. Bryan; Hofton, Michelle A.

    1999-01-01

    The upcoming generation of laser altimeters record the interaction of emitted laser radiation with terrestrial surfaces in the form of a digitized waveform. We model these laser altimeter return waveforms as the sum of the reflections from individual surfaces within laser footprints, accounting for instrument-specific properties. We compare over 1000 modeled and recorded waveform pairs using the Pearson correlation. We show that we reliably synthesize the vertical structure information for vegetation canopies contained in a medium-large diameter laser footprint from a high-resolution elevation data set.

  5. Constitutive Modeling and Testing of Polymer Matrix Composites Incorporating Physical Aging at Elevated Temperatures

    NASA Technical Reports Server (NTRS)

    Veazie, David R.

    1998-01-01

    Advanced polymer matrix composites (PMC's) are desirable for structural materials in diverse applications such as aircraft, civil infrastructure and biomedical implants because of their improved strength-to-weight and stiffness-to-weight ratios. For example, the next generation military and commercial aircraft requires applications for high strength, low weight structural components subjected to elevated temperatures. A possible disadvantage of polymer-based composites is that the physical and mechanical properties of the matrix often change significantly over time due to the exposure of elevated temperatures and environmental factors. For design, long term exposure (i.e. aging) of PMC's must be accounted for through constitutive models in order to accurately assess the effects of aging on performance, crack initiation and remaining life. One particular aspect of this aging process, physical aging, is considered in this research.

  6. Determination of Martian Northern Polar Insolation Levels Using a Geodetic Elevation Model

    NASA Technical Reports Server (NTRS)

    Arrell, J. R.; Zuber, M. T.

    2000-01-01

    Solar insolation levels at the Martian polar caps bear significantly on the seasonal and climatic cycling of volatiles on that planet. In the northern hemisphere, the Martian surface slopes downhill from the equator to the pole such that the north polar cap is situated in a 5-km-deep hemispheric-scale depression. This large-scale topographic setting plays an important role in the insolation of the northern polar cap. Elevations measured by the Mars Orbiter Laser Altimeter (MOLA) provide comprehensive, high-accuracy topographical information required to precisely determine polar insolation. In this study, we employ a geodetic elevation model to quantify the north polar insolation and consider implications for seasonal and climatic changes. Additional information is contained in original extended abstract.

  7. Geomorphological feature extraction from a digital elevation model through fuzzy knowledge-based classification

    NASA Astrophysics Data System (ADS)

    Argialas, Demetre P.; Tzotsos, Angelos

    2003-03-01

    The objective of this research was the investigation of advanced image analysis methods for geomorphological mapping. Methods employed included multiresolution segmentation of the Digital Elevation Model (DEM) GTOPO30 and fuzzy knowledge based classification of the segmented DEM into three geomorphological classes: mountain ranges, piedmonts and basins. The study area was a segment of the Basin and Range Physiographic Province in Nevada, USA. The implementation was made in eCognition. In particular, the segmentation of GTOPO30 resulted into primitive objects. The knowledge-based classification of the primitive objects based on their elevation and shape parameters, resulted in the extraction of the geomorphological features. The resulted boundaries in comparison to those by previous studies were found satisfactory. It is concluded that geomorphological feature extraction can be carried out through fuzzy knowledge based classification as implemented in eCognition.

  8. Revealing topographic lineaments through IHS enhancement of DEM data. [Digital Elevation Model

    NASA Technical Reports Server (NTRS)

    Murdock, Gary

    1990-01-01

    Intensity-hue-saturation (IHS) processing of slope (dip), aspect (dip direction), and elevation to reveal subtle topographic lineaments which may not be obvious in the unprocessed data are used to enhance digital elevation model (DEM) data from northwestern Nevada. This IHS method of lineament identification was applied to a mosiac of 12 square degrees using a Cray Y-MP8/864. Square arrays from 3 x 3 to 31 x 31 points were tested as well as several different slope enhancements. When relatively few points are used to fit the plane, lineaments of various lengths are observed and a mechanism for lineament classification is described. An area encompassing the gold deposits of the Carlin trend and including the Rain in the southeast to Midas in the northwest is investigated in greater detail. The orientation and density of lineaments may be determined on the gently sloping pediment surface as well as in the more steeply sloping ranges.

  9. Satellite Elevation Magnetic and Gravity Models of Major South American Plate Tectonic Features

    NASA Technical Reports Server (NTRS)

    Vonfrese, R. R. B.; Hinze, W. J.; Braile, L. W.; Lidiak, E. G.; Keller, G. R. (Principal Investigator); Longacre, M. B.

    1984-01-01

    Some MAGSAT scalar and vector magnetic anomaly data together with regional gravity anomaly data are being used to investigate the regional tectonic features of the South American Plate. An initial step in this analysis is three dimensional modeling of magnetic and gravity anomalies of major structures such as the Andean subduction zone and the Amazon River Aulacogen at satellite elevations over an appropriate range of physical properties using Gaus-Legendre quadrature integration method. In addition, one degree average free-air gravity anomalies of South America and adjacent marine areas are projected to satellite elevations assuming a spherical Earth and available MAGSAT data are processed to obtain compatible data sets for correlation. Correlation of these data sets is enhanced by reduction of the MAGSAT data to radial polarization because of the profound effect of the variation of the magnetic inclination over South America.

  10. Development of a seamless multisource topographic/bathymetric elevation model of Tampa Bay

    USGS Publications Warehouse

    Gesch, D.; Wilson, R.

    2001-01-01

    Many applications of geospatial data in coastal environments require knowledge of the nearshore topography and bathymetry. However, because existing topographic and bathymetric data have been collected independently for different purposes, it has been difficult to use them together at the land/water interface owing to differences in format, projection, resolution, accuracy, and datums. As a first step toward solving the problems of integrating diverse coastal datasets, the U.S. Geological Survey (USGS) and the National Oceanic and Atmospheric Administration (NOAA) are collaborating on a joint demonstration project to merge their data for the Tampa Bay region of Florida. The best available topographic and bathymetric data were extracted from the USGS National Elevation Dataset and the NOAA hydrographic survey database, respectively. Before being merged, the topographic and bathymetric datasets were processed with standard geographic information system tools to place them in a common horizontal reference frame. Also, a key part of the preprocessing was transformation to a common vertical reference through the use of VDatum, a new tool created by NOAA's National Geodetic Survey for vertical datum conversions. The final merged product is a seamless topographic/bathymetric model covering the Tampa Bay region at a grid spacing of 1 arc-second. Topographic LIDAR data were processed and merged with the bathymetry to demonstrate the incorporation of recent third party data sources for several test areas. A primary application of a merged topographic/bathymetric elevation model is for user-defined shoreline delineation, in which the user decides on the tidal condition (for example, low or high water) to be superimposed on the elevation data to determine the spatial position of the water line. Such a use of merged topographic/bathymetric data could lead to the development of a shoreline zone, which could reduce redundant mapping efforts by federal, state, and local agencies

  11. Google Earth's derived digital elevation model: A comparative assessment with Aster and SRTM data

    NASA Astrophysics Data System (ADS)

    Rusli, N.; Majid, M. R.; Din, A. H. M.

    2014-02-01

    This paper presents a statistical analysis showing additional evidence that Digital Elevation Model (DEM) derived from Google Earth is commendable and has a good correlation with ASTER (Advanced Space-borne Thermal Emission and Reflection Radiometer) and SRTM (Shuttle Radar Topography Mission) elevation data. The accuracy of DEM elevation points from Google Earth was compared against that of DEMs from ASTER and SRTM for flat, hilly and mountainous sections of a pre-selected rural watershed. For each section, a total of 5,000 DEM elevation points were extracted as samples from each type of DEM data. The DEM data from Google Earth and SRTM for flat and hilly sections are strongly correlated with the R2 of 0.791 and 0.891 respectively. Even stronger correlation is shown for the mountainous section where the R2 values between Google Earth's DEM and ASTER's and between Google Earth's DEM and SRTM's DEMs are respectively 0.917 and 0.865. Further accuracy testing was carried out by utilising the DEM dataset to delineate Muar River's watershed boundary using ArcSWAT2009, a hydrological modelling software. The result shows that the percentage differences of the watershed size delineated from Google Earth's DEM compared to those derived from Department of Irrigation and Drainage's data (using 20m-contour topographic map), ASTER and SRTM data are 9.6%, 10.6%, and 7.6% respectively. It is therefore justified to conclude that the DEM derived from Google Earth is relatively as acceptable as DEMs from other sources.

  12. Use of upscaled elevation and surface roughness data in two-dimensional surface water models

    USGS Publications Warehouse

    Hughes, J.D.; Decker, J.D.; Langevin, C.D.

    2011-01-01

    In this paper, we present an approach that uses a combination of cell-block- and cell-face-averaging of high-resolution cell elevation and roughness data to upscale hydraulic parameters and accurately simulate surface water flow in relatively low-resolution numerical models. The method developed allows channelized features that preferentially connect large-scale grid cells at cell interfaces to be represented in models where these features are significantly smaller than the selected grid size. The developed upscaling approach has been implemented in a two-dimensional finite difference model that solves a diffusive wave approximation of the depth-integrated shallow surface water equations using preconditioned Newton–Krylov methods. Computational results are presented to show the effectiveness of the mixed cell-block and cell-face averaging upscaling approach in maintaining model accuracy, reducing model run-times, and how decreased grid resolution affects errors. Application examples demonstrate that sub-grid roughness coefficient variations have a larger effect on simulated error than sub-grid elevation variations.

  13. Geomorphology and sedimentary features, and temporal component-change of lahar deposits at the northern foot of Chokai volcano, NE Japan

    NASA Astrophysics Data System (ADS)

    Minami, Y.; Ohba, T.; Kataoka, K.; Hayashi, S.

    2014-12-01

    Chokai volcano is an andesitic stratovolcano that collapsed to the north ca.2500 years ago. The post-collapse fan deposits are distributed in the northern foot of the volcano, and to reveal their depositional process in terms of modern sedimentology, we carried out the geological study includung digging survey, as well as geomorphological analysis, mineralogy, and 14C chronology. Consequently, the geological study revealed that the fan deposits consist of more than 16 units, which are debris flow, hyperconcentrated flow and streamflow deposits. We give hare general name lahar deposits for these deposits. The lahar deposits have a total thickness of 30 m, and overlie the 2.5-ka Kisakata debris avalanche deposit. The lahar deposits form a part of volcanic fan and volcaniclastic apron of Chokai volcano. In proximal areas (steep or moderate sloped areas), the lahar flowed down as debris flows, and in the distal area (horizontal area) the lahars transformed into hyperconcentrated flow or stream flows but partly arrived the area as debris flow. The hyperconcentrated flows or stream flows reached the horizontal area at least four times, supposed by AMS dating (the ages of some lahar deposits are 2200, 1500-1600, 1000-1200, and 100-200 yBP). The lahar deposits contain clasts of altered andesite, fresh andesite, mudstone and sandstone. Proportions of altered andesite clasts to total clasts decrease upwards in stratigraphic sequence. Matrices of the lower eight units are composed of grayish-blue clay, and are different from those of the upper eight units, composed of brownish yellow volcanic sand. The stratigraphic variation in matrix component is consistent with the change in matrix mineral assemblage, possibly reflecting changes in the source materials from Chokai volcano.

  14. Quantification of soil losses from tourist trails - use of Digital Elevation Models

    NASA Astrophysics Data System (ADS)

    Tomczyk, Aleksandra

    2010-05-01

    Tourism impacts in protected mountain areas are one of the main concerns for land managers. Impact to environment is most visible at locations of highly concentrated activities like tourist trails, campsites etc. The main indicators of the tourist trail degradation are: vegetation loss (trampling of vegetation cover), change of vegetation type and composition, widening of the trails, muddiness and soil erosion. The last one is especially significant, since it can cause serious transformation of the land surface. Such undesirable changes cannot be repaired without high-cost management activities, and, in some cases they can made the trails difficult and unsafe to use. Scientific understanding of soil erosion related to human impact can be useful for more effective management of the natural protected areas. The aim of this study was to use of digital elevation models (DEMs) to precisely quantify of soil losses from tourist trails. In the study precise elevation data were gathered in several test fields of 4 by 5 m spatial dimension. Measurements were taken in 13 test fields, located in two protected natural areas in south Poland: Gorce National Park and Popradzki Landscape Park. The measuring places were located on trails characterized by different slope, type of vegetation and type of use. Each test field was established by four special marks, firmly dug into the ground. Elevation data were measured with the electronic total station. Irregular elevation points were surveying with essential elements of surrounding terrain surface being included. Moreover, surveys in fixed profile lines were done. For each test field a set of 30 measurements in control points has been collected and these data provide the base for verification of digital elevation models. Average density of the surveying was 70 points per square meter (1000 - 1500 elevation points per each test fields). Surveys in each test field were carried out in August and September of 2008, June 2009 and August

  15. The ASTER Global Digital Elevation Model version 2.0 - Early Validation Results

    NASA Astrophysics Data System (ADS)

    Meyer, D. J.; Tachikawa, T.; Abrams, M.; Tsu, H.; Hato, M.; Gesch, D. B.; Crippen, R. E.

    2011-12-01

    The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) on NASA's Terra spacecraft is capable of collecting in-track stereo using nadir- and aft looking near infrared cameras. Since 2001, these stereo pairs have been used to produce single-scene (60 x 60 km) digital elevation models having vertical (root-mean-squared-error) accuracies generally between 10 m and 25 m. On June 29, 2009, NASA and the Ministry of Economy, Trade and Industry (METI) of Japan released a Global Digital Elevation Model (GDEM) to users worldwide at no charge as a contribution to the Global Earth Observing System of Systems (GEOSS). This "version 1.0" ASTER GDEM (GDEM v1.0) was compiled from over 1.2 million scene-based DEMs covering land surface between 83°N and 83°S latitudes. This GDEM is a 1 arc-second elevation grid divided and distributed as 1°-by-1° tiles. ** A joint US-Japan validation team assessed the accuracy of the GDEM v1.0, augmented by a team of 20 cooperators selected through an Announcement of Opportunity (AO). In summary, the GDEM v1.0 was found to have an overall accuracy of around 20 meters at the 95% confidence interval. The team also noted several artifacts associated with poor coverage, cloud contamination, water masking issues and the stacking process used to produce the GDEM from individual scene-based DEMs. An independent horizontal resolution study estimated the effective spatial resolution of the GDEM v1.0 to be on the order of 120 meters. ** NASA & METI will release a second version of the ASTER GDEM (v2.0) in mid-October, 2011. The GDEM v2.0 has the same gridding and tile structure as v1.0, but benefits from the inclusion of 300,000 additional scenes to improve coverage, a smaller correlation kernel (5x5 versus 9x9 for v1.0) yielding higher spatial resolution, and improved water masking. This abstract presents early validation results available at the time of submission. Early results indicate: (1) the overall accuracy (both horizontal and

  16. Predicting Bed Grain Size in Threshold Channels Using Lidar Digital Elevation Models

    NASA Astrophysics Data System (ADS)

    Snyder, N. P.; Nesheim, A. O.; Wilkins, B. C.; Edmonds, D. A.

    2011-12-01

    Over the past 20 years, researchers have developed GIS-based algorithms to extract channel networks and measure longitudinal profiles from digital elevation models (DEMs), and have used these to study stream morphology in relation to tectonics, climate and ecology. The accuracy of stream elevations from traditional DEMs (10-50 m pixels) is typically limited by the contour interval (3-20 m) of the rasterized topographic map source. This is a particularly severe limitation in low-relief watersheds, where 3 m of channel elevation change may occur over several km. Lidar DEMs (~1 m pixels) allow researchers to resolve channel elevation changes of ~0.5 m, enabling reach-scale calculations of gradient, which is the most important parameter for understanding channel processes at that scale. Lidar DEMs have the additional advantage of allowing users to make estimates of channel width. We present a process-based model that predicts median bed grain size in threshold gravel-bed channels from lidar slope and width measurements using the Shields and Manning equations. We compare these predictions to field grain size measurements in segments of three Maine rivers. Like many paraglacial rivers, these have longitudinal profiles characterized by relatively steep (gradient >0.002) and flat (gradient <0.0005) segments, with length scales of several km. This heterogeneity corresponds to strong variations in channel form, sediment supply, bed grain size, and aquatic habitat characteristics. The model correctly predicts bed sediment size within a factor of two in ~70% of the study sites. The model works best in single-thread channels with relatively low sediment supply, and poorly in depositional, multi-thread and/or fine (median grain size <20 mm) reaches. We evaluate the river morphology (using field and lidar measurements) in the context of the Parker et al. (2007) hydraulic geometry relations for single-thread gravel-bed rivers, and find correspondence in the locations where both

  17. Back to the Future: Have Remotely Sensed Digital Elevation Models Improved Hydrological Parameter Extraction?

    NASA Astrophysics Data System (ADS)

    Jarihani, B.

    2015-12-01

    Digital Elevation Models (DEMs) that accurately replicate both landscape form and processes are critical to support modeling of environmental processes. Pre-processing analysis of DEMs and extracting characteristics of the watershed (e.g., stream networks, catchment delineation, surface and subsurface flow paths) is essential for hydrological and geomorphic analysis and sediment transport. This study investigates the status of the current remotely-sensed DEMs in providing advanced morphometric information of drainage basins particularly in data sparse regions. Here we assess the accuracy of three available DEMs: (i) hydrologically corrected "H-DEM" of Geoscience Australia derived from the Shuttle Radar Topography Mission (SRTM) data; (ii) the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) version2 1-arc-second (~30 m) data; and (iii) the 9-arc-second national GEODATA DEM-9S ver3 from Geoscience Australia and the Australian National University. We used ESRI's geospatial data model, Arc Hydro and HEC-GeoHMS, designed for building hydrologic information systems to synthesize geospatial and temporal water resources data that support hydrologic modeling and analysis. A coastal catchment in northeast Australia was selected as the study site where very high resolution LiDAR data are available for parts of the area as reference data to assess the accuracy of other lower resolution datasets. This study provides morphometric information for drainage basins as part of the broad research on sediment flux from coastal basins to Great Barrier Reef, Australia. After applying geo-referencing and elevation corrections, stream and sub basins were delineated for each DEM. Then physical characteristics for streams (i.e., length, upstream and downstream elevation, and slope) and sub-basins (i.e., longest flow lengths, area, relief and slopes) were extracted and compared with reference datasets from LiDAR. Results showed that

  18. Enhancing a RADARSAT/ICESat Digital Elevation Model of West Antarctica Using MODIS Imagery

    NASA Astrophysics Data System (ADS)

    Haran, T. M.; Scambos, T. A.

    2007-12-01

    An image enhancement approach is used to develop a new digital elevation map of West Antarctica, combining multiple MODIS images and both radar altimetry and ICESat laser altimetry Digital Elevation Model (DEM) data. The method combines the wide image coverage of MODIS, and its high radiometric sensitivity (which equates to high sunward slope sensitivity), with the high precision and accuracy of ICESat and combined ICESat and radar altimetry DEMs. We calibrate brightness-to-slope relationships for several MODIS images of the central West Antarctic using smoothed DEMs derived from both sources. Using the calibrations, we then created, first, a slope map of the ice sheet surface from the image data (regressing slope information from many images), and then integrated this absolute slope map to yield complete DEMs for the region. ICESat (as of September 2007) has acquired a series of eleven near-repeat tracks over the Antarctic during the period September 2003 to April 2007, covering the continent to 86 deg S. ICESat data are acquired as a series of spot elevations, averaging a ~60m diameter surface region every ~172m. However, ICESat track paths have spacings wide enough (2 km at 85 deg; 20 - 50 km at 75 deg) that some surface ice dynamical features (e.g. flowlines, undulations, ice rises) are missed by the track data used to construct the ICESat DEM. Radar altimetry can provide some of the missing data north of 81.5 deg, but only to a maximum resolution of about 5 km. A set of cloud-cleared MODIS band 1 data from both the Aqua and Terra platforms acquired during the 2003-2004 austral summer, used in generating the Mosaic of Antarctica, MOA, surface morphology image map, were used for the image enhancement. Past analyses of the slope-brightness relationship for MODIS have shown ice surface slope precisions of +/- 0.00015. ICESat spot elevations have nominal precisions of ~5 cm under ideal conditions, although thin-cloud effects and mislocation errors can magnify these

  19. Filling the voids in the SRTM elevation model — A TIN-based delta surface approach

    NASA Astrophysics Data System (ADS)

    Luedeling, Eike; Siebert, Stefan; Buerkert, Andreas

    The Digital Elevation Model (DEM) derived from NASA's Shuttle Radar Topography Mission is the most accurate near-global elevation model that is publicly available. However, it contains many data voids, mostly in mountainous terrain. This problem is particularly severe in the rugged Oman Mountains. This study presents a method to fill these voids using a fill surface derived from Russian military maps. For this we developed a new method, which is based on Triangular Irregular Networks (TINs). For each void, we extracted points around the edge of the void from the SRTM DEM and the fill surface. TINs were calculated from these points and converted to a base surface for each dataset. The fill base surface was subtracted from the fill surface, and the result added to the SRTM base surface. The fill surface could then seamlessly be merged with the SRTM DEM. For validation, we compared the resulting DEM to the original SRTM surface, to the fill DEM and to a surface calculated by the International Center for Tropical Agriculture (CIAT) from the SRTM data. We calculated the differences between measured GPS positions and the respective surfaces for 187,500 points throughout the mountain range (ΔGPS). Comparison of the means and standard deviations of these values showed that for the void areas, the fill surface was most accurate, with a standard deviation of the ΔGPS from the mean ΔGPS of 69 m, and only little accuracy was lost by merging it to the SRTM surface (standard deviation of 76 m). The CIAT model was much less accurate in these areas (standard deviation of 128 m). The results show that our method is capable of transferring the relative vertical accuracy of a fill surface to the void areas in the SRTM model, without introducing uncertainties about the absolute elevation of the fill surface. It is well suited for datasets with varying altitude biases, which is a common problem of older topographic information.

  20. An anisotropic thermomechanical damage model for concrete at transient elevated temperatures.

    PubMed

    Baker, Graham; de Borst, René

    2005-11-15

    The behaviour of concrete at elevated temperatures is important for an assessment of integrity (strength and durability) of structures exposed to a high-temperature environment, in applications such as fire exposure, smelting plants and nuclear installations. In modelling terms, a coupled thermomechanical analysis represents a generalization of the computational mechanics of fracture and damage. Here, we develop a fully coupled anisotropic thermomechanical damage model for concrete under high stress and transient temperature, with emphasis on the adherence of the model to the laws of thermodynamics. Specific analytical results are given, deduced from thermodynamics, of a novel interpretation on specific heat, evolution of entropy and the identification of the complete anisotropic, thermomechanical damage surface. The model is also shown to be stable in a computational sense, and to satisfy the laws of thermodynamics. PMID:16243703

  1. Analysis of 30 Years Sea Surface Elevation (sse) Data Obtained From A Global Ocean Model

    NASA Astrophysics Data System (ADS)

    Wenzel, M.

    In this presentation the sea surface elevation (SSE) as obtained from a global OGCM will be analysed. The model used is the Hamburg LSG model with a 2 degree hori- zontal resolution, 23 layers in the vertical and a 10 day timestep. The model has a free surface and includes the thermo- and halosteric effects. It is integrated for 50 years forced by monthly NCEP reanalyzes data (1950-1999). The last 30 years of model output are analyzed on regional to global scale to judge the role of the different contributions to the interannual sea level variations as there are: horizontal redistribution of volume (mass), surface freshwater flux (precipitaion- evaporation) and steric effects. On the global scale the surface freshwater flux can clearly be identified as the main contributor, while on regional to local scale it is the steric effect.

  2. Ventricular dilation and elevated aqueductal pulsations in a new experimental model of communicating hydrocephalus

    SciTech Connect

    Wagshul, M.; Smith, S.; Wagshul, M.; McAllister, J.P.; Rashid, S.; Li, J.; Egnor, M.R.; Walker, M.L.; Yu, M.; Smith, S.D.; Zhang, G.; Chen, J.J.; Beneveniste, H.

    2009-03-01

    In communicating hydrocephalus (CH), explanations for the symptoms and clear-cut effective treatments remain elusive. Pulsatile flow through the cerebral aqueduct is often significantly elevated, but a clear link between abnormal pulsations and ventriculomegaly has yet to be identified. We sought to demonstrate measurement of pulsatile aqueductal flow of CSF in the rat, and to characterize the temporal changes in CSF pulsations in a new model of CH. Hydrocephalus was induced by injection of kaolin into the basal cisterns of adult rats (n = 18). Ventricular volume and aqueductal pulsations were measured on a 9.4 T MRI over a one month period. Half of the animals developed ventricular dilation, with increased ventricular volume and pulsations as early as one day post-induction, and marked chronic elevations compared to intact controls (volume: 130.15 {+-} 83.21 {mu}l vs. 15.52 {+-} 2.00 {mu}l; pulsations: 114.51 nl {+-} 106.29 vs. 0.72 {+-} 0.13 nl). Similar to the clinical presentation, the relationship between ventricular size and pulsations was quite variable. However, the pulsation time-course revealed two distinct sub-types of hydrocephalic animals: those with markedly elevated pulsations which persisted over time, and those with mildly elevated pulsations which returned to near normal levels after one week. These groups were associated with severe and mild ventriculomegaly respectively. Thus, aqueductal flow can be measured in the rat using high-field MRI and basal cistern-induced CH is associated with an immediate change in CSF pulsatility. At the same time, our results highlight the complex nature of aqueductal pulsation and its relationship to ventricular dilation.

  3. A geomorphology-based approach for digital elevation model fusion - case study in Danang city, Vietnam

    NASA Astrophysics Data System (ADS)

    Tran, T. A.; Raghavan, V.; Masumoto, S.; Vinayaraj, P.; Yonezawa, G.

    2014-07-01

    Global digital elevation models (DEM) are considered a source of vital spatial information and find wide use in several applications. The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global DEM (GDEM) and Shuttle Radar Topographic Mission (SRTM) DEM offer almost global coverage and provide elevation data for geospatial analysis. However, GDEM and SRTM still contain some height errors that affect the quality of elevation data significantly. This study aims to examine methods to improve the resolution as well as accuracy of available free DEMs by data fusion techniques and evaluating the results with a high-quality reference DEM. The DEM fusion method is based on the accuracy assessment of each global DEM and geomorphological characteristics of the study area. Land cover units were also considered to correct the elevation of GDEM and SRTM with respect to the bare-earth surface. The weighted averaging method was used to fuse the input DEMs based on a landform classification map. According to the landform types, the different weights were used for GDEM and SRTM. Finally, a denoising algorithm (Sun et al., 2007) was applied to filter the output-fused DEM. This fused DEM shows excellent correlation to the reference DEM, having a correlation coefficient R2 = 0.9986, and the accuracy was also improved from a root mean square error (RMSE) of 14.9 m in GDEM and 14.8 m in SRTM to 11.6 m in the fused DEM. The results of terrain-related parameters extracted from this fused DEM such as slope, curvature, terrain roughness index and normal vector of topographic surface are also very comparable to reference data.

  4. Investigation of potential sea level rise impact on the Nile Delta, Egypt using digital elevation models.

    PubMed

    Hasan, Emad; Khan, Sadiq Ibrahim; Hong, Yang

    2015-10-01

    In this study, the future impact of Sea Level Rise (SLR) on the Nile Delta region in Egypt is assessed by evaluating the elevations of two freely available Digital Elevation Models (DEMs): the SRTM and the ASTER-GDEM-V2. The SLR is a significant worldwide dilemma that has been triggered by recent climatic changes. In Egypt, the Nile Delta is projected to face SLR of 1 m by the end of the 21th century. In order to provide a more accurate assessment of the future SLR impact on Nile Delta's land and population, this study corrected the DEM's elevations by using linear regression model with ground elevations from GPS survey. The information for the land cover types and future population numbers were derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) land cover and the Gridded Population of the Worlds (GPWv3) datasets respectively. The DEM's vertical accuracies were assessed using GPS measurements and the uncertainty analysis revealed that the SRTM-DEM has positive bias of 2.5 m, while the ASTER-GDEM-V2 showed a positive bias of 0.8 m. The future inundated land cover areas and the affected population were illustrated based on two SLR scenarios of 0.5 m and 1 m. The SRTM DEM data indicated that 1 m SLR will affect about 3900 km(2) of cropland, 1280 km(2) of vegetation, 205 km(2) of wetland, 146 km(2) of urban areas and cause more than 6 million people to lose their houses. The overall vulnerability assessment using ASTER-GDEM-V2 indicated that the influence of SLR will be intense and confined along the coastal areas. For instance, the data indicated that 1 m SLR will inundate about 580 Km(2) (6%) of the total land cover areas and approximately 887 thousand people will be relocated. Accordingly, the uncertainty analysis of the DEM's elevations revealed that the ASTER-GDEM-V2 dataset product was considered the best to determine the future impact of SLR on the Nile Delta region. PMID:26410824

  5. Watershed boundaries and digital elevation model of Oklahoma derived from 1:100,000-scale digital topographic maps

    USGS Publications Warehouse

    Cederstrand, J.R.; Rea, A.H.

    1995-01-01

    This document provides a general description of the procedures used to develop the data sets included on this compact disc. This compact disc contains watershed boundaries for Oklahoma, a digital elevation model, and other data sets derived from the digital elevation model. The digital elevation model was produced using the ANUDEM software package, written by Michael Hutchinson and licensed from the Centre for Resource and Environmental Studies at The Australian National University. Elevation data (hypsography) and streams (hydrography) from digital versions of the U.S. Geological Survey 1:100,000-scale topographic maps were used by the ANUDEM package to produce a hydrologically conditioned digital elevation model with a 60-meter cell size. This digital elevation model is well suited for drainage-basin delineation using automated techniques. Additional data sets include flow-direction, flow-accumulation, and shaded-relief grids, all derived from the digital elevation model, and the hydrography data set used in producing the digital elevation model. The watershed boundaries derived from the digital elevation model have been edited to be consistent with contours and streams from the U.S. Geological Survey 1:100,000-scale topographic maps. The watershed data set includes boundaries for 11-digit Hydrologic Unit Codes (watersheds) within Oklahoma, and 8-digit Hydrologic Unit Codes (cataloging units) outside Oklahoma. Cataloging-unit boundaries based on 1:250,000-scale maps outside Oklahoma for the Arkansas, Red, and White River basins are included. The other data sets cover Oklahoma, and where available, portions of 1:100,000-scale quadrangles adjoining Oklahoma.

  6. Modeling Water-Surface Elevations and Virtual Shorelines for the Colorado River in Grand Canyon, Arizona

    USGS Publications Warehouse

    Magirl, Christopher S.; Breedlove, Michael J.; Webb, Robert H.; Griffiths, Peter G.

    2008-01-01

    Using widely-available software intended for modeling rivers, a new one-dimensional hydraulic model was developed for the Colorado River through Grand Canyon from Lees Ferry to Diamond Creek. Solving one-dimensional equations of energy and continuity, the model predicts stage for a known steady-state discharge at specific locations, or cross sections, along the river corridor. This model uses 2,680 cross sections built with high-resolution digital topography of ground locations away from the river flowing at a discharge of 227 m3/s; synthetic bathymetry was created for topography submerged below the 227 m3/s water surface. The synthetic bathymetry was created by adjusting the water depth at each cross section up or down until the model?s predicted water-surface elevation closely matched a known water surface. This approach is unorthodox and offers a technique to construct one-dimensional hydraulic models of bedrock-controlled rivers where bathymetric data have not been collected. An analysis of this modeling approach shows that while effective in enabling a useful model, the synthetic bathymetry can differ from the actual bathymetry. The known water-surface profile was measured using elevation data collected in 2000 and 2002, and the model can simulate discharges up to 5,900 m3/s. In addition to the hydraulic model, GIS-based techniques were used to estimate virtual shorelines and construct inundation maps. The error of the hydraulic model in predicting stage is within 0.4 m for discharges less than 1,300 m3/s. Between 1,300-2,500 m3/s, the model accuracy is about 1.0 m, and for discharges between 2,500-5,900 m3/s, the model accuracy is on the order of 1.5 m. In the absence of large floods on the flow-regulated Colorado River in Grand Canyon, the new hydraulic model and the accompanying inundation maps are a useful resource for researchers interested in water depths, shorelines, and stage-discharge curves for flows within the river corridor with 2002 topographic

  7. High Resolution Photogrammetric Digital Elevation Models Across Calving Fronts and Meltwater Channels in Greenland

    NASA Astrophysics Data System (ADS)

    Le Bel, D. A.; Brown, S.; Zappa, C. J.; Bell, R. E.; Frearson, N.; Tinto, K. J.

    2014-12-01

    Photogrammetric digital elevation models (DEMs) are a powerful approach for understanding elevation change and dynamics along the margins of the large ice sheets. The IcePod system, mounted on a New York Air National Guard LC-130, can measure high-resolution surface elevations with a Riegl VQ580 scanning laser altimeter and Imperx Bobcat IGV-B6620 color visible-wavelength camera (6600x4400 resolution); the surface temperature with a Sofradir IRE-640L infrared camera (spectral response 7.7-9.5 μm, 640x512 resolution); and the structure of snow and ice with two radar systems. We show the use of IcePod imagery to develop DEMs across calving fronts and meltwater channels in Greenland. Multiple over-flights of the Kangerlussaq Airport ramp have provided a test of the technique at a location with accurate, independently-determined elevation. Here the photogrammetric DEM of the airport, constrained by ground control measurements, is compared with the Lidar results. In July 2014 the IcePod ice-ocean imaging system surveyed the calving fronts of five outlet glaciers north of Jakobshavn Isbrae. We used Agisoft PhotoScan to develop a DEM of each calving front using imagery captured by the IcePod systems. Adjacent to the ice sheet, meltwater plumes foster mixing in the fjord, moving warm ocean water into contact with the front of the ice sheet where it can undercut the ice front and trigger calving. The five glaciers provide an opportunity to examine the calving front structure in relation to ocean temperature, fjord circulation, and spatial scale of the meltwater plumes. The combination of the accurate DEM of the calving front and the thermal imagery used to constrain the temperature and dynamics of the adjacent plume provides new insights into the ice-ocean interactions. Ice sheet margins provide insights into the connections between the surface meltwater and the fate of the water at the ice sheet base. Surface meltwater channels are visualized here for the first time using

  8. Statistical Evaluation of Fitting Accuracy of Global and Local Digital Elevation Models in Iran

    NASA Astrophysics Data System (ADS)

    Alidoost, F.; Samadzadegan, F.

    2013-09-01

    Digital Elevation Models (DEMs) are one of the most important data for various applications such as hydrological studies, topography mapping and ortho image generation. There are well-known DEMs of the whole world that represent the terrain's surface at variable resolution and they are also freely available for 99% of the globe. However, it is necessary to assess the quality of the global DEMs for the regional scale applications.These models are evaluated by differencing with other reference DEMs or ground control points (GCPs) in order to estimate the quality and accuracy parameters over different land cover types. In this paper, a comparison of ASTER GDEM ver2, SRTM DEM with more than 800 reference GCPs and also with a local elevation model over the area of Iran is presented. This study investigates DEM's characteristics such as systematic error (bias), vertical accuracy and outliers for DEMs using both the usual (Mean error, Root Mean Square Error, Standard Deviation) and the robust (Median, Normalized Median Absolute Deviation, Sample Quantiles) descriptors. Also, the visual assessment tools are used to illustrate the quality of DEMs, such as normalized histograms and Q-Q plots. The results of the study confirmed that there is a negative elevation bias of approximately 5 meters of GDEM ver2. The measured RMSE and NMAD for elevation differences of GDEM-GCPs are 7.1 m and 3.2 m, respectively, while these values for SRTM and GCPs are 9.0 m and 4.4 m. On the other hand, in comparison with the local DEM, GDEM ver2 exhibits the RMSE of about 6.7 m, a little higher than the RMSE of SRTM (5.1 m).The results of height difference classification and other statistical analysis of GDEM ver2-local DEM and SRTM-local DEM reveal that SRTM is slightly more accurate than GDEM ver2. Accordingly, SRTM has no noticeable bias and shift from Local DEM and they have more consistency to each other, while GDEM ver2 has always a negative bias.

  9. Application of Low-Cost Digital Elevation Models to Detect Change in Forest Carbon Sequestration Projects

    SciTech Connect

    Kenneth Glenn MacDicken

    2007-07-31

    This two-year study evaluated advanced multispectral digital imagery applications for assessment of forest carbon stock change. A series of bench and field studies in North Carolina and Ohio tested aerial assessments of forest change between two time periods using two software packages (ERDAS and TERREST) for Digital Elevation Model (DEM) creation, automated classification software (eCognition) for canopy segmentation and a multiple ranging laser designed to improve quality of elevation data. Results of the DEM software comparison showed that while TERREST has the potential to produce much higher resolution DEM than ERDAS, it is unable to resolve crucial canopy features adequately. Lab tests demonstrated that additional laser data improves image registration and Z-axis DEM quality. Data collected in the field revealed difficult challenges in correctly modeling the location of laser strike and subsequently determining elevations in both software packages. Automated software segmentation of tree canopies provided stem diameter and biomass carbon estimates that were within 3% of comparable ground based estimates in the Ohio site and produced similar biomass estimates for a limited number of plots in the Duke forest. Tree height change between time periods and canopy segmentation from multispectral imagery allowed calculation of forest carbon stock change at costs that are comparable to those for ground-based methods. This work demonstrates the potential of lower cost imagery systems enhanced with laser data to collect high quality imagery and paired laser data for forestry and environmental applications. Additional research on automated canopy segmentation and multi-temporal image registration is needed to refine these methods for commercial use.

  10. Analytical basis for determining slope lines in grid digital elevation models

    NASA Astrophysics Data System (ADS)

    Orlandini, Stefano; Moretti, Giovanni; Gavioli, Andrea

    2014-01-01

    An analytical basis for the determination of slope lines in grid digital elevation models is provided by using the D8-LTD method (eight slope directions, least transverse deviation). The D8-LTD method's capability to predict consistently exact slope lines as the grid cell size goes to zero is shown analytically by applying mathematical analysis methods. The use of cumulative, least transverse deviations is found to be the key factor allowing for globally unbiased approximations of slope lines. The D8-LTD method's properties are also demonstrated numerically by using digital elevation models of a synthetic sloping surface obtained from the Himmelblau function. It is shown that slope lines obtained from the D8-LTD method can approximate the exact slope lines as close as desired by selecting a grid cell size that is small enough. In contrast, the standard D8 method is found to produce significantly biased results even when small grid cells are used. The D8-LTD method outperforms the D8 method over a wide range of grid cell sizes (up to 20 m in this application), beyond which grid cell size becomes too large to validly represent the underlying sloping surface. It is therefore concluded that the D8-LTD method should be used in preference to the standard D8 method in order to obtain slope lines that are only limited in reliability by the detail of topographic data, and not by the accuracy of the slope direction method applied.

  11. Delineating small karst watersheds based on digital elevation model and eco-hydrogeological principles

    NASA Astrophysics Data System (ADS)

    Jie Luo, Guang; Jie Wang, Shi; Bai, Xiao Yong; Liu, Xiu Ming; Cheng, An Yun

    2016-03-01

    Dominated by specific eco-hydrogeological backgrounds, a small watershed delineated by using the traditional method is always inauthentic in karst regions because it cannot accurately reflect the eco-hydrological process of the dual structure of the surface and subsurface. This study proposes a new method for the delineation of small watersheds based on digital elevation models (DEMs) and eco-hydrogeological principles in karst regions. This method is applied to one section of the tributary area (Sancha River) of the Yangtze River in China. By comparing the quantity, shape, superimposition, and characteristics of the internal hydrological process of a small watershed extracted by using the digital elevation model with that extracted by using the proposed method of this study, we conclude that the small karst watersheds extracted by the new method accurately reflect the hydrological process of the river basin. Furthermore, we propose that the minimum unit of the river basin in karst regions should be the watershed, whose exit is the corrosion and corrasion baselevel and a further division of watershed may cause a significant inconsistency with the true eco-hydrological process.

  12. Hydrologic analysis of a flood based on a new Digital Elevation Model

    NASA Astrophysics Data System (ADS)

    Nishio, M.; Mori, M.

    2015-06-01

    These The present study aims to simulate the hydrologic processes of a flood, based on a new, highly accurate Digital Elevation Model (DEM). The DEM is provided by the Ministry of Land, Infrastructure, Transport and Tourism (MLIT) of Japan, and has a spatial resolution of five meters. It was generated by the new National Project in 2012. The Hydrologic Engineering Center - Hydrologic Modeling System (HEC-HMS) is used to simulate the hydrologic process of a flood of the Onga River in Iizuka City, Japan. A large flood event in the typhoon season in 2003 caused serious damage around the Iizuka City area. Precise records of rainfall data from the Automated Meteorological Data Acquisition System (AMeDAS) were input into the HEC-HMS. The estimated flood area of the simulation results by HEC-HMS was identical to the observed flood area. A watershed aggregation map is also generated by HEC-HMS around the Onga River.

  13. A model for the biaxial post-yield behavior of extruded powder aluminum at elevated temperature

    SciTech Connect

    Woods, T.O.; Berghaus, D.G. ); Peacock, H.B. )

    1990-01-01

    A model has been developed which describes the post-yield behavior of extruded powder aluminum tested biaxially in tension and torsion at elevated temperature. Plots of shear stress versus shear strain for the powder aluminum loaded in simple torsion show that the shear stress increases linearly to the yield point, then remains relatively constant in a pure plastic type of behavior. For the tension-torsion tests, there is an initial linear region up to the yield point followed by a fairly linear decrease in shear stress. A similar linear decrease in axial stress with increasing axial strain is observed in uniaxial tension tests. The model for post-yield behavior of extruded powder aluminum gives a quantified description of the macroscopic material behavior in terms of changes in the laminar powder aluminum structure.

  14. A model for the biaxial post-yield behavior of extruded powder aluminum at elevated temperature

    SciTech Connect

    Woods, T.O.; Berghaus, D.G.; Peacock, H.B.

    1990-12-31

    A model has been developed which describes the post-yield behavior of extruded powder aluminum tested biaxially in tension and torsion at elevated temperature. Plots of shear stress versus shear strain for the powder aluminum loaded in simple torsion show that the shear stress increases linearly to the yield point, then remains relatively constant in a pure plastic type of behavior. For the tension-torsion tests, there is an initial linear region up to the yield point followed by a fairly linear decrease in shear stress. A similar linear decrease in axial stress with increasing axial strain is observed in uniaxial tension tests. The model for post-yield behavior of extruded powder aluminum gives a quantified description of the macroscopic material behavior in terms of changes in the laminar powder aluminum structure.

  15. Modeling the effects of elevated temperatures on action potential propagation in unmyelinated axons

    NASA Astrophysics Data System (ADS)

    Ganguly, Mohit; Jenkins, Michael W.; Chiel, Hillel J.; Jansen, E. Duco

    2016-03-01

    Infrared lasers (λ=1.87 μm) are capable of inducing a thermally mediated nerve block in Aplysia and rat nerves. While this block is spatially precise and reversible in sensory and motor neurons, the mechanism of block is not clearly understood. Model predictions show that, at elevated temperatures, the rates of opening and closing of the voltage gated ion channels are disrupted and normal functioning of the gates is hindered. A model combining NEURON with Python is presented here that can simulate the behavior of unmyelinated nerve axons in the presence of spatially and temporally varying temperature distributions. Axon behavior and underlying mechanism leading to conduction block is investigated. The ability to understand the photothermal interaction of laser light and temperature dependence of membrane ion channels in-silico will help speed explorations of parameter space and guide future experiments testing the feasibility of selectively blocking pain conduction fibers (Photonic Analgesia of Nerves (PAIN)) in humans.

  16. Extracting low-resolution river networks from high-resolution digital elevation models

    USGS Publications Warehouse

    Olivera, F.; Lear, M.S.; Famiglietti, J.S.; Asante, Kwasi

    2002-01-01

    Including a global river network in the land component of global climate models (GCMs) is necessary in order to provide a more complete representation of the hydrologic cycle. The process of creating these networks is called river network upscaling and consists of lowering the resolution of already available fine networks to make them compatible with GCMs. Fine-resolution river networks have a level of detail appropriate for analysis at the watershed scale but are too intensive for global hydrologic studies. A river network upscaling algorithm, which processes fine-resolution digital elevation models to determine the flow directions that best describe the flow patterns in a coarser user-defined scale, is presented. The objectives of this study were to develop an algorithm that advances the previous work in the field by being applicable at a global scale, allowing for the upscaling to be performed in a projected environment, and generating evenly distributed flow directions.

  17. A stochastic, evolutionary model for range shifts and richness on tropical elevational gradients under Quaternary glacial cycles

    PubMed Central

    Colwell, Robert K.; Rangel, Thiago F.

    2010-01-01

    Quaternary glacial–interglacial cycles repeatedly forced thermal zones up and down the slopes of mountains, at all latitudes. Although no one doubts that these temperature cycles have left their signature on contemporary patterns of geography and phylogeny, the relative roles of ecology and evolution are not well understood, especially for the tropics. To explore key mechanisms and their interactions in the context of chance events, we constructed a geographical range-based, stochastic simulation model that incorporates speciation, anagenetic evolution, niche conservatism, range shifts and extinctions under late Quaternary temperature cycles along tropical elevational gradients. In the model, elevational patterns of species richness arise from the differential survival of founder lineages, consolidated by speciation and the inheritance of thermal niche characteristics. The model yields a surprisingly rich variety of realistic patterns of phylogeny and biogeography, including close matches to a variety of contemporary elevational richness profiles from an elevational transect in Costa Rica. Mountaintop extinctions during interglacials and lowland extinctions at glacial maxima favour mid-elevation lineages, especially under the constraints of niche conservatism. Asymmetry in temperature (greater duration of glacial than of interglacial episodes) and in lateral area (greater land area at low than at high elevations) have opposing effects on lowland extinctions and the elevational pattern of species richness in the model—and perhaps in nature, as well. PMID:20980317

  18. Assessing the quality for hydrological applications of digital elevation models derived from contours

    NASA Astrophysics Data System (ADS)

    Wise, Stephen

    2000-07-01

    Digital elevation models (DEMs) are becoming increasingly important tools in hydrological research and in water resources management. The quality of DEMs, however, normally is reported simply as the root mean square error of elevation, a statistic that fails to capture the numerous sources of error in DEMs or to predict their effect on the result of using the DEM. This paper presents a review of other approaches to assessing DEM quality, and argues that a full assessment of DEM quality must focus on the accuracy and reliability of the final product of the DEM analysis. A number of DEMs for the Slapton Ley catchments in Devon derived from digitized contour data are compared in an initial assessment of their sustainability for use in hydrological work. Two are available for purchase from data suppliers, and five more were created using a variety of interpolation techniques in widely available geographical information system software. The different interpretation methods produce DEMs with different artefacts, although analyses of the distribution of elevation values, and visual techniques, suggested that none of these were of a particularly pronounced nature. The results of using the DEMs to derive drainage networks and catchment areas showed that at the broad scale there was a high level of agreement between the DEMs. There were, however, important differences of detail. For example, some DEMs predicted drainage lines that occasionally crossed the original contours. The results of calculating the TOPMODEL topographic index showed far more variation, because the index is calculated for each pixel in the area, rather than being an aggregate result derived from numerous pixels. The main conclusion was that care should always be taken to assess the quality of a DEM before attempting to use it, and that results should always be checked to ensure that they appear to be reasonable.

  19. Validation of the ASTER Global Digital Elevation Model Version 2 over the conterminous United States

    USGS Publications Warehouse

    Gesch, Dean B.; Oimoen, Michael J.; Zhang, Zheng; Meyer, David J.; Danielson, Jeffrey J.

    2012-01-01

    The ASTER Global Digital Elevation Model Version 2 (GDEM v2) was evaluated over the conterminous United States in a manner similar to the validation conducted for the original GDEM Version 1 (v1) in 2009. The absolute vertical accuracy of GDEM v2 was calculated by comparison with more than 18,000 independent reference geodetic ground control points from the National Geodetic Survey. The root mean square error (RMSE) measured for GDEM v2 is 8.68 meters. This compares with the RMSE of 9.34 meters for GDEM v1. Another important descriptor of vertical accuracy is the mean error, or bias, which indicates if a DEM has an overall vertical offset from true ground level. The GDEM v2 mean error of -0.20 meters is a significant improvement over the GDEM v1 mean error of -3.69 meters. The absolute vertical accuracy assessment results, both mean error and RMSE, were segmented by land cover to examine the effects of cover types on measured errors. The GDEM v2 mean errors by land cover class verify that the presence of aboveground features (tree canopies and built structures) cause a positive elevation bias, as would be expected for an imaging system like ASTER. In open ground classes (little or no vegetation with significant aboveground height), GDEM v2 exhibits a negative bias on the order of 1 meter. GDEM v2 was also evaluated by differencing with the Shuttle Radar Topography Mission (SRTM) dataset. In many forested areas, GDEM v2 has elevations that are higher in the canopy than SRTM.

  20. Morphological convexity measures for terrestrial basins derived from digital elevation models

    NASA Astrophysics Data System (ADS)

    Lim, Sin Liang; Daya Sagar, B. S.; Chet Koo, Voon; Tien Tay, Lea

    2011-09-01

    Geophysical basins of terrestrial surfaces have been quantitatively characterized through a host of indices such as topological quantities (e.g. channel bifurcation and length ratios), allometric scaling exponents (e.g. fractal dimensions), and other geomorphometric parameters (channel density, Hack's and Hurst exponents). Channel density, estimated by taking the ratio between the length of channel network ( L) and the area of basin ( A) in planar form, provides a quantitative index that has hitherto been related to various geomorphologically significant processes. This index, computed by taking the planar forms of channel network and its corresponding basin, is a kind of convexity measure in the two-dimensional case. Such a measure - estimated in general as a function of basin area and channel network length, where the important elevation values of the topological region within a basin and channel network are ignored - fails to capture the spatial variability between homotopic basins possessing different altitude-ranges. Two types of convexity measures that have potential to capture the terrain elevation variability are defined as the ratio of (i) length of channel network function and area of basin function and (ii) areas of basin and its convex hull functions. These two convexity measures are estimated in three data sets that include (a) synthetic basin functions, (b) fractal basin functions, and (c) realistic digital elevation models (DEMs) of two regions of peninsular Malaysia. It is proven that the proposed convexity measures are altitude-dependent and that they could capture the spatial variability across the homotopic basins of different altitudes. It is also demonstrated on terrestrial DEMs that these convexity measures possess relationships with other quantitative indexes such as fractal dimensions and complexity measures (roughness indexes).

  1. A time series of TanDEM-X digital elevation models to monitor a glacier surge

    NASA Astrophysics Data System (ADS)

    Wendt, Anja; Mayer, Christoph; Lambrecht, Astrid; Floricioiu, Dana

    2016-04-01

    Bivachny Glacier, a tributary of the more than 70 km long Fedchenko Glacier in the Pamir Mountains, Central Asia, is a surge-type glacier with three known surges during the 20th century. In 2011, the most recent surge started which, in contrast to the previous ones, evolved down the whole glacier and reached the confluence with Fedchenko Glacier. Spatial and temporal glacier volume changes can be derived from high-resolution digital elevation models (DEMs) based on bistatic InSAR data from the TanDEM-X mission. There are nine DEMs available between 2011 and 2015 covering the entire surge period in time steps from few months up to one year. During the surge, the glacier surface elevation increased by up to 130 m in the lower part of the glacier; and change rates of up to 0.6 m per day were observed. The surface height dataset was complemented with glacier surface velocity information from TerraSAR-X/ TanDEM-X data as well as optical Landsat imagery. While the glacier was practically stagnant in 2000 after the end of the previous surge in the 1990s, the velocity increase started in 2011 in the upper reaches of the ablation area and successively moved downwards and intensified, reaching up to 4.0 m per day. The combination of surface elevation changes and glacier velocities, both of high temporal and spatial resolution, provides the unique opportunity to describe and analyse the evolution of the surge in unprecedented detail. Especially the relation between the mobilization front and the local mass transport provides insight into the surge dynamics.

  2. Bathymetric survey and digital elevation model of Little Holland Tract, Sacramento-San Joaquin Delta, California

    USGS Publications Warehouse

    Snyder, Alexander G.; Lacy, Jessica R.; Stevens, Andrew W.; Carlson, Emily M.

    2016-01-01

    The U.S. Geological Survey conducted a bathymetric survey in Little Holland Tract, a flooded agricultural tract, in the northern Sacramento-San Joaquin Delta (the “Delta”) during the summer of 2015. The new bathymetric data were combined with existing data to generate a digital elevation model (DEM) at 1-meter resolution. Little Holland Tract (LHT) was historically diked off for agricultural uses and has been tidally inundated since an accidental levee breach in 1983. Shallow tidal regions such as LHT have the potential to improve habitat quality in the Delta. The DEM of LHT was developed to support ongoing studies of habitat quality in the area and to provide a baseline for evaluating future geomorphic change. The new data comprise 138,407 linear meters of real-time-kinematic (RTK) Global Positioning System (GPS) elevation data, including both bathymetric data collected from personal watercraft and topographic elevations collected on foot at low tide. A benchmark (LHT15_b1) was established for geodetic control of the survey. Data quality was evaluated both by comparing results among surveying platforms, which showed systematic offsets of 1.6 centimeters (cm) or less, and by error propagation, which yielded a mean vertical uncertainty of 6.7 cm. Based on the DEM and time-series measurements of water depth, the mean tidal prism of LHT was determined to be 2,826,000 cubic meters. The bathymetric data and DEM are available at http://dx.doi.org/10.5066/F7RX9954. 

  3. The conservation value of elevation data accuracy and model sophistication in reserve design under sea-level rise.

    PubMed

    Zhu, Mingjian; Hoctor, Tom; Volk, Mike; Frank, Kathryn; Linhoss, Anna

    2015-10-01

    Many studies have explored the value of using more sophisticated coastal impact models and higher resolution elevation data in sea-level rise (SLR) adaptation planning. However, we know little about to what extent the improved models and data could actually lead to better conservation outcomes under SLR. This is important to know because high-resolution data are likely to not be available in some data-poor coastal areas in the world and running more complicated coastal impact models is relatively time-consuming, expensive, and requires assistance by qualified experts and technicians. We address this research question in the context of identifying conservation priorities in response to SLR. Specifically, we investigated the conservation value of using more accurate light detection and ranging (Lidar)-based digital elevation data and process-based coastal land-cover change models (Sea Level Affecting Marshes Model, SLAMM) to identify conservation priorities versus simple "bathtub" models based on the relatively coarse National Elevation Dataset (NED) in a coastal region of northeast Florida. We compared conservation outcomes identified by reserve design software (Zonation) using three different model dataset combinations (Bathtub-NED, Bathtub-Lidar, and SLAMM-Lidar). The comparisons show that the conservation priorities are significantly different with different combinations of coastal impact models and elevation dataset inputs. The research suggests that it is valuable to invest in more accurate coastal impact models and elevation datasets in SLR adaptive conservation planning because this model-dataset combination could improve conservation outcomes under SLR. Less accurate coastal impact models, including ones created using coarser Digital Elevation Model (DEM) data can still be useful when better data and models are not available or feasible, but results need to be appropriately assessed and communicated. A future research priority is to investigate how

  4. Preparing for Volcanic Hazards: An Examination of Lahar Knowledge, Risk Perception, and Preparedness around Mount Baker and Glacier Peak, WA

    NASA Astrophysics Data System (ADS)

    Corwin, K.; Brand, B. D.

    2015-12-01

    As the number of people living at risk from volcanic hazards in the U.S. Pacific Northwest continues to rise, so does the need for improved hazard science, mitigation, and response planning. The effectiveness of these efforts relies not only on scientists and policymakers, but on individuals and their risk perception and preparedness levels. This study examines the individual knowledge, perception, and preparedness of over 500 survey respondents living or working within the lahar zones of Mount Baker and Glacier Peak volcanoes. We (1) explore the common disconnect between accurate risk perception and adequate preparedness; (2) determine how participation in hazard response planning influences knowledge, risk perception, and preparedness; and (3) assess the effectiveness of current lahar hazard maps for public risk communication. Results indicate that a disconnect exists between perception and preparedness for the majority of respondents. While 82% of respondents accurately anticipate that future volcanic hazards will impact the Skagit Valley, this knowledge fails to motivate increased preparedness. A majority of respondents also feel "very responsible" for their own protection and provision of resources during a hazardous event (83%) and believe they have the knowledge and skills necessary to respond effectively to such an event (56%); however, many of these individuals still do not adequately prepare. When asked what barriers prevent them from preparing, respondents primarily cite a lack of knowledge about relevant local hazards. Results show that participation in response-related activities—a commonly recommended solution to this disconnect—minimally influences preparedness. Additionally, although local hazard maps successfully communicate the primary hazard—97% of respondents recognize the lahar hazard—many individuals incorrectly interpret other important facets of the maps. Those who participate in response-related activities fail to understand these

  5. Swift snowmelt and floods (lahars) caused by great pyroclastic surge at Mount St Helens volcano, Washington, 18 May 1980

    USGS Publications Warehouse

    Waitt, R.B.

    1989-01-01

    The initial explosions at Mount St. Helens, Washington, on the moring of 18 May 1980 developed into a huge pyroclastic surge that generated catastrophic floods off the east and west flanks of the volcano. Near-source surge deposits on the east and west were lithic, sorted, lacking in accretionary lapilli and vesiculated ash, not plastered against upright obstacles, and hot enough to char wood - all attributes of dry pyroclastic surge. Material deposited at the surge base on steep slopes near the volcano transformed into high-concentration lithic pyroclastic flows whose deposits contain charred wood and other features indicating that these flows were hot and dry. Stratigraphy shows that even the tail of the surge had passed the east and west volcano flanks before the geomorphically distinct floods (lahars) arrived. This field evidence undermines hypotheses that the turbulent surge was itself wet and that its heavy components segregated out to transform directly into lahars. Nor is there evidence that meters-thick snow-slab avalanches intimately mixed with the surge to form the floods. The floods must have instead originated by swift snowmelt at the base of a hot and relatively dry turbulent surge. Impacting hot pyroclasts probably transferred downslope momentum to the snow surface and churned snow grains into the surge base. Melting snow and accumulating hot surge debris may have moved initially as thousands of small thin slushflows. As these flows removed the surface snow and pyroclasts, newly uncovered snow was partly melted by the turbulent surge base; this and accumulating hot surge debris in turn began flowing, a self-sustaining process feeding the initial flows. The flows thus grew swiftly over tens of seconds and united downslope into great slushy ejecta-laden sheetfloods. Gravity accelerated the floods to more than 100 km/h as they swept down and off the volcano flanks while the snow component melted to form great debris-rich floods (lahars) channeled into

  6. The 26 May 2006 magnitude 6.4 Yogyakarta earthquake south of Mt. Merapi volcano: Did lahar deposits amplify ground shaking and thus lead to the disaster?

    NASA Astrophysics Data System (ADS)

    Walter, T. R.; Wang, R.; Luehr, B.-G.; Wassermann, J.; Behr, Y.; Parolai, S.; Anggraini, A.; Günther, E.; Sobiesiak, M.; Grosser, H.; Wetzel, H.-U.; Milkereit, C.; Sri Brotopuspito, P. J. K.; Harjadi, P.; Zschau, J.

    2008-05-01

    Indonesia is repeatedly unsettled by severe volcano- and earthquake-related disasters, which are geologically coupled to the 5-7 cm/a tectonic convergence of the Australian plate beneath the Sunda Plate. On Saturday, 26 May 2006, the southern coast of central Java was struck by an earthquake at 2254 UTC in the Sultanate Yogyakarta. Although the magnitude reached only M w = 6.4, it left more than 6,000 fatalities and up to 1,000,000 homeless. The main disaster area was south of Mt. Merapi Volcano, located within a narrow topographic and structural depression along the Opak River. The earthquake disaster area within the depression is underlain by thick volcaniclastic deposits commonly derived in the form of lahars from Mt. Merapi Volcano, which had a major influence leading to the disaster. In order to more precisely understand this earthquake and its consequences, a 3-month aftershock measurement campaign was performed from May to August 2006. We here present the first location results, which suggest that the Yogyakarta earthquake occurred at 10-20 km distance east of the disaster area, outside of the topographic depression. Using simple model calculations taking material heterogeneity into account we illustrate how soft volcaniclastic deposits may locally amplify ground shaking at distance. As the high degree of observed damage may have been augmented by the seismic response of the volcaniclastic Mt. Merapi deposits, this work implies that the volcano had an indirect effect on the level of earthquake destruction.

  7. The 1997 and 2001 lahars of Popocatépetl volcano (Central Mexico): textural and sedimentological constraints on their origin and hazards

    NASA Astrophysics Data System (ADS)

    Capra, L.; Poblete, M. A.; Alvarado, R.

    2004-03-01

    Popocatépetl volcano is the most active volcano in central Mexico, and represents a high risk for more than 40 million people, including Mexico City. In 1994, volcanic activity at Popocatépetl renewed with the formation of ash-rich plumes up to 7 km high. In April 1996, lava emissions filled the crater and were accompanied by a series of explosions that produced eruptive columns up to 8 km high. Associated with explosive events in 1997 and 2001, two major lahars events occurred, leaving debris flow deposits along Huiloac Gorge for as far as 15 km, to the town of Santiago Xalitzinta. The 1997 debris flow deposit originated after a prolonged emission of ash which caused glacier melt and a rapid release of water (1×10 7 m 3). The amount of melting water was sufficient to gradually erode the river bed causing a flood that gradually transformed from a debris flow to a hyperconcentrated flow. In contrast, the 2001 debris flow that originated from a post-depositional remobilization of a pumice flow deposit, did not experience any flow transformation and carried 25% water at maximum. The different behavior of these two lahars has important hazard implications. Both lahars reached Xalitzintla town, but at that point, the 1997 lahar had already transformed into a sediment-loaded stream flow. The 2001 lahar, by contrast, maintained the characteristics of a debris flow, being more competent, and with greater destructive power. What happened with these lahars demonstrates how important it is to take into consideration secondary volcanic phenomena. Even though they were not large flows, they were capable of threatening populated areas, even during periods of volcanic quiescence or reduced magmatic activity.

  8. Assessing the quality of digital elevation models obtained from mini unmanned aerial vehicles for overland flow modelling in urban areas

    NASA Astrophysics Data System (ADS)

    Leitão, João P.; Moy de Vitry, Matthew; Scheidegger, Andreas; Rieckermann, Jörg

    2016-04-01

    Precise and detailed digital elevation models (DEMs) are essential to accurately predict overland flow in urban areas. Unfortunately, traditional sources of DEM, such as airplane light detection and ranging (lidar) DEMs and point and contour maps, remain a bottleneck for detailed and reliable overland flow models, because the resulting DEMs are too coarse to provide DEMs of sufficient detail to inform urban overland flows. Interestingly, technological developments of unmanned aerial vehicles (UAVs) suggest that they have matured enough to be a competitive alternative to satellites or airplanes. However, this has not been tested so far. In this study we therefore evaluated whether DEMs generated from UAV imagery are suitable for urban drainage overland flow modelling. Specifically, 14 UAV flights were conducted to assess the influence of four different flight parameters on the quality of generated DEMs: (i) flight altitude, (ii) image overlapping, (iii) camera pitch, and (iv) weather conditions. In addition, we compared the best-quality UAV DEM to a conventional lidar-based DEM. To evaluate both the quality of the UAV DEMs and the comparison to lidar-based DEMs, we performed regression analysis on several qualitative and quantitative metrics, such as elevation accuracy, quality of object representation (e.g. buildings, walls and trees) in the DEM, which were specifically tailored to assess overland flow modelling performance, using the flight parameters as explanatory variables. Our results suggested that, first, as expected, flight altitude influenced the DEM quality most, where lower flights produce better DEMs; in a similar fashion, overcast weather conditions are preferable, but weather conditions and other factors influence DEM quality much less. Second, we found that for urban overland flow modelling, the UAV DEMs performed competitively in comparison to a traditional lidar-based DEM. An important advantage of using UAVs to generate DEMs in urban areas is

  9. ASTER-Derived 30-Meter-Resolution Digital Elevation Models of Afghanistan

    USGS Publications Warehouse

    Chirico, Peter G.; Warner, Michael B.

    2007-01-01

    INTRODUCTION The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is an imaging instrument aboard the Terra satellite, launched on December 19, 1999, as part of the National Aeronautics and Space Administration's (NASA) Earth Observing System (EOS). The ASTER sensor consists of three subsystems: the visible and near infrared (VNIR), the shortwave infrared (SWIR), and the thermal infrared (TIR), each with a different spatial resolution (VNIR, 15 meters; SWIR, 30 meters, TIR 90 meters). The VNIR system has the capability to generate along-track stereo images that can be used to create digital elevation models (DEMs) at 30-meter resolution. Currently, the only available DEM dataset for Afghanistan is the 90-meter-resolution Shuttle Radar Topography Mission (SRTM) data. This dataset is appropriate for macroscale DEM analysis and mapping. However, ASTER provides a low cost opportunity to generate higher resolution data. For this publication, study areas were identified around populated areas and areas where higher resolution elevation data were desired to assist in natural resource assessments. The higher resolution fidelity of these DEMs can also be used for other terrain analysis including landform classification and geologic structure analysis. For this publication, ASTER scenes were processed and mosaicked to generate 36 DEMs which were created and extracted using PCI Geomatics' OrthoEngine 3D Stereo software. The ASTER images were geographically registered to Landsat data with at least 15 accurate and well distributed ground control points with a root mean square error (RMSE) of less that one pixel (15 meters). An elevation value was then assigned to each ground control point by extracting the elevation from the 90-meter SRTM data. The 36 derived DEMs demonstrate that the software correlated on nearly flat surfaces and smooth slopes accurately. Larger errors occur in cloudy and snow-covered areas, lakes, areas with steep slopes, and

  10. Combining criteria for delineating lahar- and flash-flood-prone hazard and risk zones for the city of Arequipa, Peru

    NASA Astrophysics Data System (ADS)

    Thouret, J.-C.; Enjolras, G.; Martelli, K.; Santoni, O.; Luque, J. A.; Nagata, M.; Arguedas, A.; Macedo, L.

    2013-02-01

    Arequipa, the second largest city in Peru, is exposed to many natural hazards, most notably earthquakes, volcanic eruptions, landslides, lahars (volcanic debris flows), and flash floods. Of these, lahars and flash floods, triggered by occasional torrential rainfall, pose the most frequently occurring hazards that can affect the city and its environs, in particular the areas containing low-income neighbourhoods. This paper presents and discusses criteria for delineating areas prone to flash flood and lahar hazards, which are localized along the usually dry (except for the rainy season) ravines and channels of the Río Chili and its tributaries that dissect the city. Our risk-evaluation study is based mostly on field surveys and mapping, but we also took into account quality and structural integrity of buildings, available socio-economic data, and information gained from interviews with risk-managers officials. In our evaluation of the vulnerability of various parts of the city, in addition to geological and physical parameters, we also took into account selected socio-economic parameters, such as the educational and poverty level of the population, unemployment figures, and population density. In addition, we utilized a criterion of the "isolation factor", based on distances to access emergency resources (hospitals, shelters or safety areas, and water) in each city block. By combining the hazard, vulnerability and exposure criteria, we produced detailed risk-zone maps at the city-block scale, covering the whole city of Arequipa and adjacent suburbs. Not surprisingly, these maps show that the areas at high risk coincide with blocks or districts with populations at low socio-economic levels. Inhabitants at greatest risk are the poor recent immigrants from rural areas who live in unauthorized settlements in the outskirts of the city in the upper parts of the valleys. Such settlements are highly exposed to natural hazards and have little access to vital resources. Our