The use of high-resolution terrain data in gravity field prediction

Different types of gravity prediction methods for local and regional gravity evaluation are developed, tested, and compared. Four different test areas were particularly selected in view of different prediction requirements. Also different parts of the spectrum of the gravity field were considered

Evaluating the Differences of Gridding Techniques for Digital Elevation Models Generation and Their Influence on the Modeling of Stony Debris Flows Routing: A Case Study From Rovina di Cancia Basin (North-Eastern Italian Alps)

Debris \ufb02ows are among the most hazardous phenomena in mountain areas. To cope with debris \ufb02ow hazard, it is common to delineate the risk-prone areas through routing models. The most important input to debris \ufb02ow routing models are the topographic data, usually in the form of Digital Elevation Models (DEMs). The quality of DEMs depends on the accuracy, density, and spatial distribution of the sampled points; on the characteristics of the surface; and on the applied gridding methodology. Therefore, the choice of the interpolation method affects the realistic representation of the channel and fan morphology, and thus potentially the debris \ufb02ow routing modeling outcomes. In this paper, we initially investigate the performance of common interpolation methods (i.e., linear triangulation, natural neighbor, nearest neighbor, Inverse Distance to a Power, ANUDEM, Radial Basis Functions, and ordinary kriging) in building DEMs with the complex topography of a debris \ufb02ow channel located in the Venetian Dolomites (North-eastern Italian Alps), by using small footprint full- waveform Light Detection And Ranging (LiDAR) data. The investigation is carried out through a combination of statistical analysis of vertical accuracy, algorithm robustness, and spatial clustering of vertical errors, and multi-criteria shape reliability assessment. After that, we examine the in\ufb02uence of the tested interpolation algorithms on the performance of a Geographic Information System (GIS)-based cell model for simulating stony debris \ufb02ows routing. In detail, we investigate both the correlation between the DEMs heights uncertainty resulting from the gridding procedure and that on the corresponding simulated erosion/deposition depths, both the effect of interpolation algorithms on simulated areas, erosion and deposition volumes, solid-liquid discharges, and channel morphology after the event. The comparison among the tested interpolation methods highlights that the ANUDEM and ordinary kriging algorithms are not suitable for building DEMs with complex topography. Conversely, the linear triangulation, the natural neighbor algorithm, and the thin-plate spline plus tension and completely regularized spline functions ensure the best trade-off among accuracy and shape reliability. Anyway, the evaluation of the effects of gridding techniques on debris \ufb02ow routing modeling reveals that the choice of the interpolation algorithm does not signi\ufb01cantly affect the model outcomes

Geothermal Play Fairway Analysis, Part 2: GIS Methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development

River-bed armoring as a granular segregation phenomenon

Gravel-river beds typically have an "armored" layer of coarse grains on the surface, which acts to protect finer particles underneath from erosion. River bed-load transport is a kind of dense granular flow, and such flows are known to vertically segregate grains. The contribution of granular physics to river-bed armoring, however, has not been investigated. Here we examine these connections in a laboratory river with bimodal sediment size, by tracking the motion of particles from the surface to deep inside the bed, and find that armor develops by two distinct mechanisms. Bed-load transport in the near-surface layer drives rapid segregation, with a vertical advection rate proportional to the granular shear rate. Creeping grains beneath the bed-load layer give rise to slow but persistent segregation, which is diffusion dominated and insensitive to shear rate. We verify these findings with a continuum phenomenological model and discrete element method simulations. Our results suggest that river beds armor by granular segregation from below --- rather than fluid-driven sorting from above --- while also providing new insights on the mechanics of segregation that are relevant to a wide range of granular flows

Kinematic reconstruction of a deep-seated gravitational slope deformation by geomorphic analyses

On 4 November 2010, a deep-seated gravitational slope deformation (North Italy) reactivated with sudden ground movement. A 450,000 m2 mountainous area moved some metres downslope, but the undeniable signs were only connected to the triggering of a debris flow from the bulging area’s detrital cover and the presence of a continuous perimeter fracture near the crown area. Based on two detailed LiDAR surveys (2 m × 2 m) performed just a few days before and after the event, a quantitative topographic analysis was performed in a GIS environment, integrating morphometric terrain parameters (slope, aspect, surface roughness, hill shade, and curvature). The DEMs analysis highlighted some morphological changes related to deeper as well as shallow movements. Both global and sectorial displacements were widely verified and discussed, finally inferring that the geometry, persistence, and layout of all movements properly justify each current morphostructure, which has the shape of a typical Sackung-type structure with impulsive kinematics. Moreover, a targeted field survey allowed specific clues to be found that confirmed the global deduced dynamics of the slope deformation. Finally, thanks to a ground-based interferometric radar system (GB-InSAR) that was installed a few days after the reactivation, the residual deep-seated gravitational slope deformation (DSGSD) movements were also monitored. In the landslide lower bulging area, a localized material progression of small entities was observed for some months after the parossistic event, indicating a slow dissipation of forces in sectors more distant from the crown area

Towards safer mining: the role of modelling software to find missing persons after a mine collapse

Purpose. The purpose of the study is to apply science and technology to determine the most likely location of a container in which three miners were trapped after the Lily mine disaster. Following the collapse of the Crown Pillar at Lily Mine in South Africa on the 5th of February 2016, there was a national outcry to find the three miners who were trapped in a surface container lamp room that disappeared in the sinkhole that formed during the surface col-lapse. Methods. At a visit to Lily Mine on the 9th of March, the Witwatersrand Mining Institute suggested a two-way strategy going forward to find the container in which the miners are trapped and buried. The first approach, which is the subject of this paper, is to test temporal 3D modeling software technology to locate the container, and second, to use scientific measurement and testing technologies. The overall methodology used was to first, request academia and research entities within the University to supply the WMI with ideas, which ideas list was compiled as responses came in. These were scrutinized and literature gathered for a conceptual study on which these ideas are likely to work. The software screening and preliminary testing of such software are discussed in this article. Findings. The findings are that software modeling is likely to locate the present position of the container, but accurate data and a combination of different advanced software packages will be required, but at tremendous cost. Originality. This paper presents original work on how software technology can be used to locate missing miners. Practical implications. The two approaches were not likely to recover the miners alive because of the considerable time interval, but will alert the rescue team and mine workers when they come in close proximity to them.Мета. Визначення можливого місця локалізації лампового приміщення контейнера, в якому опинилися три шахтаря після аварії на шахті Лілі (Барбертон, Мпумаланга) методом комп’ютерного моделювання. Після обвалення стельового цілика на шахті Лілі 5 лютого 2016 року почалася національна кампанія з порятунку трьох шахтарів, які залишилися у ламповому приміщенні поверхневого транспортного контейнера, що провалився в утворену після вибуху воронку. Методика. Співробітниками Гірничого Інституту (Уітуотерс) запропонована двостадійна стратегія пошуку контейнера, в якому існує ймовірність знаходження шахтарів. В рамках першого підходу (який розглядається у даній статті) для виявлення контейнера здійснювалось випробування комп’ютерної технології 3D-моделювання в часі. Другий підхід передбачав технологію проведення наукового вимірювання та експерименту. В цілому, методологія включала, насамперед, підключення викладацького та наукового складу університету до вирішення проблеми шляхом комплексної генерації ідей, які були об’єднані в загальний список, вивчені із залученням відповідних літературних джерел, і найбільш реалістичні ідеї були виділені із загального переліку. Дана стаття розглядає результати комп’ютерної експертизи цих ідей та перевірки надійності відповідного програмного забезпечення. Результати. Для зручності моделювання процес обвалення був розділений на три окремі фази: руйнування воронки, руйнування західного схилу та небезпека ковзання на південних схилах. Ідентифіковано програмні технології, які можуть імітувати рух контейнера у перших двох фазах обвалення. В результаті моделювання у програмному забезпеченні ParaView виявлено місце розташування даного контейнера. Виконано аналіз південного схилу за допомогою ArcGIS і складені карти небезпеки схилу для району, а також підземні карти порятунку з маршрутами евакуації. Встановлено, що комп’ютерне моделювання може визначити місцезнаходження контейнера, але для цього потрібні точні вихідні дані й комплекс дорогих високоефективних програмних пакетів. Наукова новизна. Вперше застосовано комплекс комп’ютерних технологій та програмного забезпечення для пошуку зниклих шахтарів після аварійних ситуацій у підземному просторі шахт. Практична значимість. При застосуванні двостадійної стратегії пошуку шахтарів, що опинилися під завалом порід, команда рятувальників отримає сигнал про наближення до їх місцезнаходження.Цель. Определение возможного места локализации лампового помещения контейнера, в котором оказались три шахтера после аварии на шахте Лили (Барбертон, Мпумаланга) методом компьютерного моделирования. После обрушения потолочного целика на шахте Лили 5 февраля 2016 года началась национальная кампания по спасению трех шахтеров, оставшихся в ламповом помещении поверхностного транспортного контейнера, который провалился в воронку, образовавшуюся после взрыва. Методика. Сотрудниками Горного Института (Уитуотерс) предложена двухстадийная стратегия поиска контейнера, в котором существует вероятность нахождения шахтеров. В рамках первого подхода (который рассматривается в данной статье) для обнаружения контейнера производилось испытание компьютерной технологии 3D-моделирования во времени. Второй подход предполагал технологию проведения научного измерения и эксперимента. В целом, методология включала, прежде всего, подключение преподавательского и научного состава университета к решению проблемы путем комплексной генерации идей, которые были объединены в общий список, изучены с привлечением соответствующих литературных источников, и наиболее реалистичные идеи были выделены из общего списка. Настоящая статья рассматривает результаты компьютерной экспертизы данных идей и проверки надежности соответствующего программного обеспечения. Результаты. Для удобства моделирования процесс обрушения был разделен на три отдельные фазы: разрушение воронки, разрушение западного склона и опасность скольжения на южных склонах. Идентифицированы программные технологии, которые могут имитировать движение контейнера в первых двух фазах обрушения. В результате моделирования в программном обеспечении ParaView выявлено местоположение данного контейнера. Выполнен анализа южного склона с помощью ArcGIS и составлены карты опасности склона для района, а также подземные карты спасения с маршрутами эвакуации. Установлено, что компьютерное моделирование может определить местонахождение контейнера, но для этого нужны точные исходные данные и комплекс дорогостоящих высокоэффективных программных пакетов. Научная новизна. Впервые применен комплекс компьютерных технологий и программного обеспечения для поиска пропавших шахтеров после аварийных ситуаций в подземном пространстве шахт. Практическая значимость. При применении двухстадийной стратегии поиска шахтеров, оказавшихся под завалом пород, команда горноспасателей получит сигнал о приближении к их местонахождению.The results of the article were obtained without the support of any of the projects or funding

Far-Field Ground Strain Failure Mode Assessment for Mineral Extraction Near Dams

Steeply sloped, high relief landforms with fractured sedimentary geology in proximity to underground mine voids have the potential to produce far-field ground strains which initiate strain failure modes for cross valley dams. The Ryerson Station Dam breach in 2005 brought attention to the possibility of increased hazard to dams due to far-field ground strain phenomena from underground mining.;The mapping study found that three United States Army Corps of Engineers (U.S.ACE) owned and operated dams were within the 1-mile buffer radius, and four within the 2.7-mile buffer radius of mine permits in West Virginia. Some U.S.ACE dams can be very large in their length, changing the offset distances used herein significantly. There were zero U.S.ACE dams found to have any underground mine permit directly underneath their location. The total frequency of occurrence for U.S.ACE dams considered nearby to underground mine permits in WV was 17.4%.;There are many more dams considered in the National Inventory of Dams (NID) than U.S.ACE dams alone. The same buffer radii were chosen to assess the offset distances of NID dams to underground mine permits as were used for the U.S.ACE dam assessment. There were found to be 115 NID dams within the 2.7-mile buffer radius with 79 of those within 1 mile, and 45 directly undermined. The total frequency of occurrence for NID dams existing nearby to permitted underground mines was found to be 18.9%

Two dimensional sediment transport model using parallel computers

Management and development of water bodies is vital for meeting domestic, agricultural, energy and industrial needs. To that end, dams, artificial channels, lakes and other water structures have been constructed. Management and development of these structures encounter problems of land erosion, reservoir silting, and degradation and aggradation of channel beds, which need to be addressed. Fundamental to these problems are sediment transport, erosion and deposition. Numerical modeling of sediment transport is the best tool to simulate sediment transport in a water body. This study develops a vertically integrated two-dimensional numerical sediment transport model. Sediment transport is simulated in two parts in this model: suspended load and bed load. A fractional step approach is used to solve the two-dimensional advection diffusion equation, which splits the advection-diffusion equation in to two separate parts: advection and diffusion. High resolution conservative algorithm is used to solve the advection part and a semi implicit finite difference scheme is used to solve the diffusion part. Different parallel numerical solvers are developed to solve linear system of equations resulting from diffusion part. Non-uniformity in sediment mixture which is quite common in real world problems is considered. The model is tested for different analytical and laboratory test cases. The model is coded for parallel computers so that enormous power of parallel computers can be exploited