Geographic Resources Analysis Support System Research Articles

Artificial Neural Networks for Mapping Coastal Lagoon of Chilika Lake, India, Using Earth Observation Data

This study presents the environmental mapping of the Chilika Lake coastal lagoon, India, using satellite images Landsat 8-9 OLI/TIRS processed using machine learning (ML) methods. The largest brackish water coastal lagoon in Asia, Chilika Lake, is a wetland of international importance included in the Ramsar site due to its rich biodiversity, productivity, and precious habitat for migrating birds and rare species. The vulnerable ecosystems of the Chilika Lagoon are subject to climate effects (monsoon effects) and anthropogenic activities (overexploitation through fishing and pollution by microplastics). Such environmental pressure results in the eutrophication of the lake, coastal erosion, fluctuations in size, and changes in land cover types in the surrounding landscapes. The habitat monitoring of the coastal lagoons is complex and difficult to implement with conventional Geographic Information System (GIS) methods. In particular, landscape variability, patch fragmentation, and landscape dynamics play a crucial role in environmental dynamics along the eastern coasts of the Bay of Bengal, which is strongly affected by the Indian monsoon system, which controls the precipitation pattern and ecosystem structure. To improve methods of environmental monitoring of coastal areas, this study employs the methods of ML and Artificial Neural Networks (ANNs), which present a powerful tool for computer vision, image classification, and analysis of Earth Observation (EO) data. Multispectral satellite data were processed by several ML image classification methods, including Random Forest (RF), Support Vector Machine (SVM), and the ANN-based MultiLayer Perceptron (MLP) Classifier. The results are compared and discussed. The ANN-based approach outperformed the other methods in terms of accuracy and precision of mapping. Ten land cover classes around the Chilika coastal lagoon were identified via spatio-temporal variations in land cover types from 2019 until 2024. This study provides ML-based maps implemented using Geographic Resources Analysis Support System (GRASS) GIS image analysis software and aims to support ML-based mapping approach of environmental processes over the Chilika Lake coastal lagoon, India.

Random Forest Classifier Algorithm of Geographic Resources Analysis Support System Geographic Information System for Satellite Image Processing: Case Study of Bight of Sofala, Mozambique

Mapping coastal regions is important for environmental assessment and for monitoring spatio-temporal changes. Although traditional cartographic methods using a geographic information system (GIS) are applicable in image classification, machine learning (ML) methods present more advantageous solutions for pattern-finding tasks such as the automated detection of landscape patches in heterogeneous landscapes. This study aimed to discriminate landscape patterns along the eastern coasts of Mozambique using the ML modules of a Geographic Resources Analysis Support System (GRASS) GIS. The random forest (RF) algorithm of the module ‘r.learn.train’ was used to map the coastal landscapes of the eastern shoreline of the Bight of Sofala, using remote sensing (RS) data at multiple temporal scales. The dataset included Landsat 8-9 OLI/TIRS imagery collected in the dry period during 2015, 2018, and 2023, which enabled the evaluation of temporal dynamics. The supervised classification of RS rasters was supported by the Scikit-Learn ML package of Python embedded in the GRASS GIS. The Bight of Sofala is characterized by diverse marine ecosystems dominated by swamp wetlands and mangrove forests located in the mixed saline–fresh waters along the eastern coast of Mozambique. This paper demonstrates the advantages of using ML for RS data classification in the environmental monitoring of coastal areas. The integration of Earth Observation data, processed using a decision tree classifier by ML methods and land cover characteristics enabled the detection of recent changes in the coastal ecosystem of Mozambique, East Africa.

Monitoring Seasonal Fluctuations in Saline Lakes of Tunisia Using Earth Observation Data Processed by GRASS GIS

This study documents the changes in the Land Use/Land Cover (LULC) in the region of saline lakes in north Tunisia, Sahara Desert. Remote sensing data are a valuable data source in monitoring LULC in lacustrine landscapes, because variations in the extent of lakes are visible from space and can be detected on the images. In this study, changes in LULC of the salt pans of Tunisia were evaluated using a series of 12 Landsat 8-9 Operational Land Imager (OLI) and Thermal Infrared (TIRS) images. The images were processed with the Geographic Resources Analysis Support System (GRASS) Geographic Information System (GIS) software. The study area included four salt lakes of north Tunisia in the two regions of the Gulf of Hammamet and Gulf of Gabès: (1) Sebkhet de Sidi el Hani (Sousse Governorate), (2) Sebkha de Moknine (Mahdia Governorate), (3) Sebkhet El Rharra and (4) Sebkhet en Noual (Sfax). A quantitative estimate of the areal extent analysed in this study is 182 km × 185 km for each Landsat scene in two study areas: Gulf of Hammamet and Gulf of Gabès. The images were analysed for the period 2017–2023 on months February, April and July for each year. Spatio-temporal changes in LULC and their climate–environmental driving forces were analysed. The results were interpreted and the highest changes were detected by accuracy assessment, computing the class separability matrices, evaluating the means and standard deviation for each band and plotting the reject probability maps. Multi-temporal changes in LULC classes are reported for each image. The results demonstrated that changes in salt lakes were determined for winter/spring/summer months as detected changes in water/land/salt/sand/vegetation areas. The accuracy of the classified images was evaluated using pixel rejection probability values, which were filtered out using the ‘r.mapcalc’ module of GRASS GIS. The confidence levels were computed and visualised with a series of maps along with the error matrix and measured convergence level of classified pixels. This paper contributes to the environmental monitoring of Tunisian landscapes and analysis of climate effects on LULC in landscapes of north Africa.

Mapping Wetlands of Kenya Using Geographic Resources Analysis Support System (GRASS GIS) with Remote Sensing Data

Abstract Monitoring wetlands of Kenya is critical for analysis of environmental changes since they present unique ecosystems with special hydrological balance and biodiversity. In this study, the Landsat 8-9 OLI/TIRS satellite images for 2015-2022 were classified using GRASS GIS scripts to evaluate changes in the Lorian Swamp wetland, north-eastern Kenya. The results of the image analysis presented maps of land cover changes including wetlands. The study demonstrated technical effectiveness of the GRASS GIS for image analysis, and contributed to the environmental monitoring of African wetlands.

GDAL and PROJ Libraries Integrated with GRASS GIS for Terrain Modelling of the Georeferenced Raster Image

Libraries with pre-written codes optimize the workflow in cartography and reduce labour intensive data processing by iteratively applying scripts to implementing mapping tasks. Most existing Geographic Information System (GIS) approaches are based on traditional software with a graphical user’s interface which significantly limits their performance. Although plugins are proposed to improve the functionality of many GIS programs, they are usually ad hoc in finding specific mapping solutions, e.g., cartographic projections and data conversion. We address this limitation by applying the principled approach of Geospatial Data Abstraction Library (GDAL), library for conversions between cartographic projections (PROJ) and Geographic Resources Analysis Support System (GRASS) GIS for geospatial data processing and morphometric analysis. This research presents topographic analysis of the dataset using scripting methods which include several tools: (1) GDAL, a translator library for raster and vector geospatial data formats used for converting Earth Global Relief Model (ETOPO1) GeoTIFF in XY Cartesian coordinates into World Geodetic System 1984 (WGS84) by the ‘gdalwarp’ utility; (2) PROJ projection transformation library used for converting ETOPO1 WGS84 grid to cartographic projections (Cassini–Soldner equirectangular, Equal Area Cylindrical, Two-Point Equidistant Azimuthal, and Oblique Mercator); and (3) GRASS GIS by sequential use of the following modules: r.info, d.mon, d.rast, r.colors, d.rast.leg, d.legend, d.northarrow, d.grid, d.text, g.region, and r.contour. The depth frequency was analysed by the module ‘d.histogram’. The proposed approach provided a systematic way for morphometric measuring of topographic data and combine the advantages of the GDAL, PROJ, and GRASS GIS tools that include the informativeness, effectiveness, and representativeness in spatial data processing. The morphometric analysis included the computed slope, aspect, profile, and tangential curvature of the study area. The data analysis revealed the distribution pattern in topographic data: 24% of data with elevations below 400 m, 13% of data with depths −5000 to −6000 m, 4% of depths have values −3000 to −4000 m, the least frequent data (−6000 to 7000 m) <1%, 2% of depths have values −2000 to 3000 m in the basin, while other values are distributed proportionally. Further, by incorporating the generic coordinate transformation software library PROJ, the raster grid was transformed into various cartographic projections to demonstrate distortions in shape and area. Scripting techniques of GRASS GIS are demonstrated for applications in topographic modelling and raster data processing. The GRASS GIS shows the effectiveness for mapping and visualization, compatibility with libraries (GDAL, PROJ), technical flexibility in combining Graphical User Interface (GUI), and command-line data processing. The research contributes to the technical cartographic development.

Accuracy assessment and performance analysis of raster to vector conversions on LULC data – India

Purpose The purpose of this paper is to convert real-world raster data into vector format and evaluate loss of accuracy in the conversion process. Open-source Geographic Information System (GIS) is used in this process and system resource utilizations were measured for conversion and accuracy analysis methods. Shape complexity attributes were analyzed in co-relation to the observed conversion errors. Design/methodology/approach The paper empirically evaluated the challenges and overheads involved in the format conversion algorithms available in open-source GIS with real-world land use and land cover (LULC) map data of India. Across the different LULC categories, geometric errors of varying density were observed in Quantum GIS (QGIS) algorithm. Area extents of original raster data were compared to the vector forms and the shape attributes such as average number of vertices and shape irregularity were evaluated to explore the possible correlation. Findings The results indicate that Geographic Resources Analysis Support System provides near error-free conversion algorithm. At the same time, the overall time taken for the conversion and the system resource utilizations were optimum as compared to the QGIS algorithm. Higher vector file sizes were generalized and accuracy loss was tested. Research limitations/implications Complete shape complexity analysis could not be achieved, as the weight factor for the irregularity of the shapes is to be varied based on the demography as well as on the LULC category. Practical implications Because of the higher system resource requirements of topological checker tool, positional accuracy checks for the converted objects could not be completed. Originality/value This paper addresses the need of accuracy analysis of real-world spatial data conversions from raster to vector format along with experimental setups challenges and impact of shape complexity.

Assessment of automated evapotranspiration estimates obtained using the GP-SEBAL algorithm for dry forest vegetation (Caatinga) and agricultural areas in the Brazilian semiarid region

Real evapotranspiration (ETr) plays a key role in water balance, especially for dry forest vegetation (Caatinga) and irrigated agricultural areas, as is the case in the Petrolina region located in the Brazilian semiarid region (BSAR). Population growth increases the need for food and the use of water resources, which are scarcer in semiarid regions. Thus, knowledge of the energy balance (EB) and radiative balance are indispensable for the management of water resources, especially in irrigated agricultural areas. Currently, one of the biggest challenges in the application of algorithms for estimating ETr based on satellite images is the calibration of the EB. The goal of this study was to evaluate the applicability of Geographic Resources Analysis Support System (GRASS) Python surface EB algorithm for land (GP-SEBAL) to calculate ETr for Caatinga vegetation and agricultural areas in the BSAR. GP-SEBAL was applied to two images from Landsat 8 obtained in 2013 and compared with manual SEBAL and observed data. In this study, kappa coefficient accuracy was used to compare the ground truth data and predicted land use and land cover. The performance of algorithms in determining the EB components, land-surface temperature (LST), net radiation (Rn), soil heat flux (G), latent heat flux (LE), sensible heat flux (H), and ETr were compared with statistical indices and uncertainty percentages (PUs) were calculated. In addition, a sensitivity analysis using a deviation of 30 W/m2 with an interval of 5 W/m2 was used to analyze the sensitivity of ETr. The results showed a strong correlation between GP-SEBAL and manual SEBAL for ETr, LST, Rn, LH, and H. The values suggest the regularity of estimating the EB in an automated manner. The results indicate a small PU difference in the estimated ETr values calculated using GP-SEBAL and SEBAL for two images of 10.87% and 20.11% for Caatinga vegetation and 15.01% and 35.16% for agriculture, respectively. We conclude that GP-SEBAL is an efficient automated algorithm for estimating ETr that uses a less complex process with a lower incidence of errors and considerably faster execution time than classical applications of SEBAL.

Seasonal hydrological loading from GPS observed data across contiguous USA using integrated R and Hadoop-GIS framework

Hydrological water mass redistribution on the Earth’s surface causes elastic loading deformation. In this study, vertical crustal deformation is assessed by an integrated spatio-temporal approach based on R statistical software and Hadoop-GIS framework at both regional and local scales. Specifically, the study examines the relationship between seasonal vertical loading deformation and seasonal hydrological loading from rain and snow. The vertical loading deformation is characterized by time-series estimated from continuous Global Positioning System (GPS) network across the contiguous USA for a timeframe of 48 months (January 1, 2013 to December 31, 2016) which coincides with California’s most extreme recent drought. The data processing used custom-built R programming scripts and spatial libraries that were integrated with Geographic Resources Analysis Support System (GRASS), Geographic Information System (GIS), and Hive framework which is a data warehouse extension of Apache Hadoop used as a database query interface. The relationships of vertical displacement are explored by visualization techniques such as spatial maps and wavelet coherence plots. Additional analysis of the relationships can be performed through interactive web-based interface.

Fractal surfaces of synthetical DEM generated by GRASS GIS module r.surf.fractal from ETOPO1 raster grid

The research problem is about to generate artificial fractal landscape surfaces from the Digital Elevation Model (DEM) using a stochastic algorithm by Geographic Resources Analysis Support System Geographic Information System (GRASS GIS) software. Fractal surfaces resemble appearance of natural topographic terrain and its structure using random surface modelling. Study area covers Kuril-Kamchatka region, Sea of Okhotsk, North Pacific Ocean. Techniques were included into GRASS GIS modules (r.relief, d.rast, r.slope.aspect, r.mapcalc) for raster calculation, processing and visualization. Module 'r.surf.fractal' was applied for generating synthetic fractal surface from ETOPO1 DEM GeoTIFF using algorithm of fractal analysis. Three tested dimensions of the fractal surfaces were automatically mapped and visualized. Algorithm of the automated fractal DEM modelling visualized variations in steepness and aspect of the artificially generated slopes in the mountains. Controllable topographic variation of the fractal surfaces was applied for three dimensions: dim=2.0001, 2.0050, 2.0100. Auxiliary modules were used for the visualization of DEMs (d.rast, r.colors, d.vect, r.contour, d.redraw, d.mon). Modules 'r.surf.gauss' and 'r.surf.random' were applied for artificial modelling as Gauss and random based mathematical surfaces, respectively. Univariate statistics for fractal surfaces were computed for comparative analysis of maps representing continuous fields by module 'r.univar': number of cells, min/max, range, mean, variance, standard deviation, variation coefficient and sum. The paper includes 9 maps and GRASS GIS codes used for visualization.

Land use Land Cover Mapping using Modified Ant Colony Optimization Technique

Land use Land cover classification is an important aspect for managing natural resources and monitoring environmental changes. Urban expansion becomes one of the major challenges for the administrator. The LANDSAT 8 images are processed using the open source GRASS (Geographic Resource Analysis Support System). Unsupervised classification technique based on Ant Colony Optimization (ACO) algorithm has been modified and proposed as Modified Ant Colony Optimization (MACO) for LULC classification. In order to improve the classification accuracy of the proposed algorithm, we have combined spatial, spectral and texture features to extract more information of homogeneous land surface. The classification accuracy of the proposed algorithm has been compared with other unsupervised classification methods such as k-means, ISODATA and ACO algorithms. The overall classification accuracy of proposed unsupervised MACO algorithm has been increased by 11.24 %, 8.24% for open space and water bodies class, respectively as compared to ACO algorithm.

ANALYSIS OF THE IMPACT OF VEGETATION DISTRIBUTION, URBANIZATION, AND SOLAR RADIATION ON THE SEASONAL VARIATION OF THE URBAN HEAT ISLAND EFFECT IN CEBU CITY USING LANDSAT AND GLOBAL HORIZONTAL IRRADIANCE DATA

Abstract. The Urban Heat Island (UHI) phenomenon is a manifestation of the abnormal amount of heat generated in urban areas and anthropogenic land surface modifications. While urbanization can improve material comfort and be a boon to the economy, the accompanying problems associated with urbanization like the UHI effect has implications on health, demand for water and energy, and impacts the microclimate. Land surface temperature (LST), the Normalized Difference Vegetation Index (NDVI), and the Normalized Difference Built-up Index (NDBI) were calculated from historical remotely-sensed Landsat data from 2013 to present. The global horizontal irradiance (GHI) was computed from the lidar-derived elevation model of Cebu City using the Geographical Resources Analysis Support System (GRASS). It is shown that annual variation in average temperatures in Cebu is generally less than 5 °C. Mean UHI temperatures in Cebu City do not show a clear trend over time, but categorizing data by season, namely the rainy season (June–November), the cool dry season (December–February), and the hot dry season (March–May), permits the emergence of a pattern. Surface temperatures for the cool dry season and hot dry season show a linearly increasing trend with R2 values of 0.916 and 0.514, respectively. This study further investigates the temporal change in the degree and extent of the UHI in Cebu City by analyzing LST maps. Regression analysis is done to determine how LST is affected by the distribution of vegetation (NDVI) and built-up (NDBI), and the seasonal variation in solar radiation through the GHI.

Classification of Vegetation to Estimate Forest Fire Danger Using Landsat 8 Images: Case Study

The vegetation cover of the Earth plays an important role in the life of mankind, whether it is natural forest or agricultural crop. The study of the variability of the vegetation cover, as well as observation of its condition, allows timely actions to make a forecast and monitor and estimate the forest fire condition. The objectives of the research were (i) to process the satellite image of the Gilbirinskiy forestry located in the basin of Lake Baikal; (ii) to select homogeneous areas of forest vegetation on the basis of their spectral characteristics; (iii) to estimate the level of forest fire danger of the area by vegetation types. The paper presents an approach for estimation of forest fire danger depending on vegetation type and radiant heat flux influence using geographic information systems (GIS) and remote sensing data. The Environment for Visualizing Images (ENVI) and the Geographic Resources Analysis Support System (GRASS) software were used to process satellite images. The area’s forest fire danger estimation and visual presentation of the results were carried out in ArcGIS Desktop software. Information on the vegetation was obtained using the analysis of the Landsat 8 Operational Land Imager (OLI) images for a typical forestry of the Lake Baikal natural area. The maps (schemes) of the Gilbirinskiy forestry were also used in the present article. The unsupervised k-means classification was used. Principal component analysis (PCA) was applied to increase the accuracy of decoding. The classification of forest areas according to the level of fire danger caused by the typical ignition source was carried out using the developed method. The final information product was the map displaying vector polygonal feature class, containing the type of vegetation and the level of fire danger for each forest quarter in the attribute table. The fire danger estimation method developed by the authors was applied to each separate quarter and showed realistic results. The method used may be applicable for other areas with prevailing forest vegetation.

Edaphic-Climatic Zoning of Eucalyptus Species in the Rio de Janeiro State, Brazil

ABSTRACT The aim of this study was to perform the edaphic-climatic zoning of six eucalyptus species (E. camaldulensis, C. citriodora, E. cloeziana , E. grandis, E. saligna, and E. urophylla) in the state of Rio de Janeiro (RJ). The edaphic-climatic zoning was determined by overlayer of thermal and water (water deficit [DEF] and rainfall), type of soils (pedology), and relief (altitude) suitability maps. Suitability maps were obtained using Geographic Resources Analysis Support System (GRASS-GIS). The species studied showed potential for cultivation in the entire RJ due to different tolerances and requirements to the environmental conditions. The exception was the E. saligna, which proved to be unsuitable for the entire state, while C. citriodora showed the greatest suitability. Among edaphic-climatic factors analyzed, air temperature and DEF were the most important limiting factor for the suitability of species.

ELECTRICITY INFRASTRUCTURE PLANNING: EVALUATING SOLAR POTENTIAL ASSESSMENTS FOR INFORMAL SETTLEMENTS USING GRASS AND FREELY AVAILABLE DATA

Abstract. Informal settlements, also known as slums or shanty towns, are characterised by rapid and unstructured expansion, poorly constructed buildings, and in some cases, they are on disputed land. Such settlements often lack basic services, such as electricity. As a result, informal settlement dwellers turn to hazardous alternative sources of energy, such as illegal electricity connections and paraffin. Solar power is a clean and safe alternative. However, informal settlements are often located on undesirable land on the urban fringe where the topography may hinder the use of solar energy. The high density of dwellings could also be a hindrance. Therefore, the solar potential needs to be assessed before any implementations are planned. Solar potential assessment functionality is generally available in geographic information system (GIS) products. The nature, cost and accessibility of datasets required for the assessment vary significantly. In this paper, we evaluate the results of solar potential assessments using GRASS (Geographic Resources Analysis Support System) for a number of different datasets. The assessments were done for two informal settlements in the City of Tshwane (South Africa): Alaska, which is nestled on a hill; and Phomolong, a densely populated settlement with a rather flat topography. The results show that solar potential assessments with open source GIS software and freely available data are feasible. This eliminates the need for lengthy and bureaucratic procurement processes and reduces the financial costs of assessing solar potential for informal settlements.

Development of a QGIS Plugin to Obtain Parameters and Elements of Plantation Trees and Vineyards with Aerial Photographs

Unmanned Aerial Vehicle (UAV) imagery allows for a new way of obtaining geographic information. In this work, a Geographical Information System (GIS) open source application was developed in QGIS software that estimates several parameters and metrics on tree crown through image analysis techniques (image segmentation and image classification) and fractal analysis. The metrics that have been estimated were: area, perimeter, number of trees, distance between trees, and a missing tree check. This methodology was tested on three different plantations: olive, eucalyptus, and vineyard. The application developed is free, open source and takes advantage of QGIS integration with external software. Several tools available from Orfeo Toolbox and Geographic Resources Analysis Support System (GRASS) GIS were employed to generate a classified raster image which allows calculating the metrics referred before. The application was developed in the Python 2.7 language. Also, some functions, modules, and classes from the QGIS Application Programming Interface (API) and PyQt4 API were used. This new plugin is a valuable tool, which allowed for automatizing several parameters and metrics on tree crown using GIS analysis tools, while considering data acquired by UAV.

The GRASS GIS temporal framework

ABSTRACTThe availability of continental and global-scale spatio-temporal geographical data sets and the requirement to efficiently process, analyse and manage them led to the development of the temporally enabled Geographic Resources Analysis Support System (GRASS GIS). We present the temporal framework that extends GRASS GIS with spatio-temporal capabilities. The framework provides comprehensive functionality to implement a full-featured temporal geographic information system (GIS) based on a combined field and object-based approach. A significantly improved snapshot approach is used to manage spatial fields of raster, three-dimensional raster and vector type in time. The resulting timestamped spatial fields are organised in spatio-temporal fields referred to as space-time data sets. Both types of fields are handled as objects in our framework. The spatio-temporal extent of the objects and related metadata is stored in relational databases, thus providing additional functionalities to perform SQL-based analysis. We present our combined field and object-based approach in detail and show the management, analysis and processing of spatio-temporal data sets with complex spatio-temporal topologies. A key feature is the hierarchical processing of spatio-temporal data ranging from topological analysis of spatio-temporal fields over boolean operations on spatio-temporal extents, to single pixel, voxel and vector feature access. The linear scalability of our approach is demonstrated by handling up to 1,000,000 raster layers in a single space-time data set. We provide several code examples to show the capabilities of the GRASS GIS Temporal Framework and present the spatio-temporal intersection of trajectory data which demonstrates the object-based ability of our framework.

Implementation of Algorithm for Satellite-Derived Bathymetry using Open Source GIS and Evaluation for Tsunami Simulation

Accurate and high resolution bathymetric data is a necessity for a wide range of coastal oceanographic research topics. Active sensing methods, such as ship-based soundings and Light Detection and Ranging (LiDAR), are expensive and time consuming solutions. Therefore, the significance of Satellite-Derived Bathymetry (SDB) has increased in the last ten years due to the availability of multi-constellation, multi-temporal, and multi-resolution remote sensing data as Open Data. Effective SDB algorithms have been proposed by many authors, but there is no ready-to-use software module available in the Geographical Information System (GIS) environment as yet. Hence, this study implements a Geographically Weighted Regression (GWR) based SDB workflow as a Geographic Resources Analysis Support System (GRASS) GIS module (i.image.bathymetry). Several case studies were carried out to examine the performance of the module in multi-constellation and multi-resolution satellite imageries for different study areas. The results indicate a strong correlation between SDB and reference depth. For instance, case study 1 (Puerto Rico, Northeastern Caribbean Sea) has shown an coefficient of determination (R2) of 0.98 and an Root Mean Square Error (RMSE) of 0.61 m, case study 2 (Iwate, Japan) has shown an R2 of 0.94 and an RMSE of 1.50 m, and case study 3 (Miyagi, Japan) has shown an R2 of 0.93 and an RMSE of 1.65 m. The reference depths were acquired by using LiDAR for case study 1 and an echo-sounder for case studies 2 and 3. Further, the estimated SDB has been used as one of the inputs for the Australian National University and Geoscience Australia (ANUGA) tsunami simulation model. The tsunami simulation results also show close agreement with post-tsunami survey data. The i.mage.bathymetry module developed as a part of this study is made available as an extension for the Open Source GRASS GIS to facilitate wide use and future improvements.

Identifying Urban Features from LiDAR for a High‐Resolution Urban Hydrologic Model

AbstractLight Detection and Ranging (LiDAR), is relatively inexpensive, provides high spatial resolution sampling at great accuracy, and can be used to generate surface terrain and land cover datasets for urban areas. These datasets are used to develop high‐resolution hydrologic models necessary to resolve complex drainage networks in urban areas. This work develops a five‐step algorithm to generate indicator fields for tree canopies, buildings, and artificial structures using Geographic Resources Analysis Support System (GRASS‐GIS), and a common computing language, Matrix Laboratory. The 54 km2 study area in Parker, Colorado consists of twenty‐four 1,500 × 1,500 m LiDAR subsets at 1 m resolution with varying degrees of urbanization. The algorithm correctly identifies 96% of the artificial structures within the study area; however, application success is dependent upon urban extent. Urban land use fractions below 0.2 experienced an increase in falsely identified building locations. ParFlow, a three‐dimensional, grid‐based hydrological model, uses these building and artificial structure indicator fields and digital elevation model for a hydrologic simulation. The simulation successfully develops the complex drainage network and simulates overland flow on the impervious surfaces (i.e., along the gutters and off rooftops) made possible through this spatial analysis process.

Building a hybrid DEVS and GRASS model using a composable cellular automaton

Modeling and simulation is pervasive throughout many different disciplines. As computing technology has provided more capability, the systems being modeled and simulated have grown larger and more complex. Often times, these large systems are managed as interacting subsystems. When it is necessary for the simulation to allow disparate subsystems to maintain their independence, then a hybrid model of the subsystems should be used. Furthermore, to ease the burden of verification and validation of simulation results, a proven system theoretical modeling specification should be used. However, many communities have already adopted nonsystem theoretical software solutions and established a group of domain experts familiar with these tools. This paper provides two things: a formal approach to building a hybrid model, and a discussion of how to incorporate a nonsystem theoretical software implementation into a proven framework. The first is done through the implementation of a Knowledge Interchange Broker (KIB) as an Interaction Model (IM). The second is accomplished by exemplifying the use of the IM in an agent-environment hybrid model. In the hybrid model, the agent is implemented in the Discrete-event System (DEVS) specification and the environment is implemented in the Geographical Resources Analysis Support System (GRASS) using a Composable Cellular Automaton (CCA) specification. This concept has been successfully applied to both example models and an interdisciplinary research project where the interactions between human activities and landscape processes are studied.

R.randomwalk v1, a multi-functional conceptual tool for mass movement routing

Abstract. We introduce r.randomwalk, a flexible and multi-functional open-source tool for backward and forward analyses of mass movement propagation. r.randomwalk builds on GRASS GIS (Geographic Resources Analysis Support System – Geographic Information System), the R software for statistical computing and the programming languages Python and C. Using constrained random walks, mass points are routed from defined release pixels of one to many mass movements through a digital elevation model until a defined break criterion is reached. Compared to existing tools, the major innovative features of r.randomwalk are (i) multiple break criteria can be combined to compute an impact indicator score; (ii) the uncertainties of break criteria can be included by performing multiple parallel computations with randomized parameter sets, resulting in an impact indicator index in the range 0–1; (iii) built-in functions for validation and visualization of the results are provided; (iv) observed landslides can be back analysed to derive the density distribution of the observed angles of reach. This distribution can be employed to compute impact probabilities for each pixel. Further, impact indicator scores and probabilities can be combined with release indicator scores or probabilities, and with exposure indicator scores. We demonstrate the key functionalities of r.randomwalk for (i) a single event, the Acheron rock avalanche in New Zealand; (ii) landslides in a 61.5 km2 study area in the Kao Ping Watershed, Taiwan; and (iii) lake outburst floods in a 2106 km2 area in the Gunt Valley, Tajikistan.