A model for multiphase flow and transport in porous media including a phenomenological approach to account for deformation—a model concept and its validation within a code intercomparison study

Multiphase flow processes in unsaturated cohesive soils are often affected by deformation due to swelling and shrinking as a result of varying water contents. This paper presents a model concept which is denoted ‘phenomenological’ in terms of the processes responsible for soil deformation, since the effects of deformation on flow and transport are only considered by constitutive relations that allow an adaptation of the hydraulic properties. This new model is validated in a detailed intercomparison study with two state-of-the-art models that are capable of explicitly describing the processes relevant for the deformation. A ‘numerical experiment’ with a state-of-the-art reference model is used to produce ‘measurement data’ for an inverse-modelling-based estimation of the model input parameters for the phenomenological concept.     

Intrusion versus inversion—a 3D density model of the southern rim of the Northwest German Basin

An unsolved problem of regional importance for both the evolution and structure of the Northwest German Basin is the existence or non-existence of the so-called Bramsche Massif. Explaining the nature of this massif and the cause of a related strong, positive Bouguer anomaly (Bramsche Anomaly) is critical. In the study described here, we tested an existing “intrusion model” against a newer “inversion model” in the southern Northwest German Basin. In the intrusion model, the strongly-positive Bouguer anomaly represents the gravity effect of an intrusion at depths between 6 and 10 km. More recent interpretations invoke tectonic inversion rather than intrusion to explain increased burial and the low level of hydrocarbon maturity found in boreholes. We tested these different interpretations by constructing 3D forward density models to 15 km depth. The intrusion model was updated and adjusted to incorporate recent data and we also modelled pre-Zechstein structures using different scenarios. The final model has a very good fit between measured and modelled gravity fields. Based on currently available seismic and structural models, as well as borehole density measurements, we show that the positive Bouguer anomaly cannot be modeled without a high-density, intrusive-like body at depth. However, further in-sight into the crustal structures of the Bramsche region requires more detailed investigations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Forthcoming events

Tropical cyclones (hurricanes and typhoons) produce high winds that can generate waves capable of damaging coral reefs. As cyclones frequently pass through northeast Australia's Great Barrier Reef (GBR), it is important to understand how the spatial distribution of reef damage changes over time. However, direct measurements of wave damage, or even wave heights or wind speeds, are rare within the GBR. An important factor in estimating whether cyclone damage was possible is the magnitude and duration of high‐energy wind and waves. Thus, before the spatio‐temporal dynamics of past cyclone damage can be modelled, it is necessary to reconstruct the spread, intensity, and duration of high‐energy conditions during individual cyclones. This was done every hour along the track taken by each of 85 cyclones that passed near the GBR from 1969 to 2003, by implementing a cyclone wind hindcasting model directly within a raster GIS using cyclone data available from the Australian Bureau of Meteorology. Three measures of cyclone energy (maximum wind speed—MAX, duration of gales—GALES, and continuous duration of gales—CGALES) were derived from these data. For three cyclones, where field data documenting actual reef damage from cyclone‐generated waves were available, the predictive ability of each measure was assessed statistically. All three performed better in predicting reef damage at sites surveyed along the high‐energy reef front than those surveyed along the more protected reef back. MAX performed best for cyclone Joy (r ² = 0.5), while CGALES performed best for cyclones Ivor (r ² = 0.23) and Justin (r ² = 0.48). Using thresholds for MAX and GALES obtained via comparison with field data of damage, it was possible to produce a preliminary prediction of the risk of wave damage across the GBR from each of the 85 cyclones. The results suggest that while up to two‐thirds of the GBR was at risk from some damage for 30–50% of the time series (～18 out of 35 years), only scattered areas of the region were at risk more frequently than that.     

Comparison of existing methods for building triangular irregular network,models of terrain from grid digital elevation models

Abstract

Triangulated irregular networks (TINs) are increasingly popular for their efficiency in data storage and their ability to accommodate irregularly spaced elevation points for many applications of geographical information systems. This paper reviews and evaluates various methods for extracting TINs from dense digital elevation models (DEMs) on a sample DEM. Both structural and statistical comparisons show that the methods perform with different rates of success in different settings. Users of DEM to TIN conversion methods should be aware of the strengths and weaknesses of the methods in addition to their own purposes before conducting the conversion.     

Book Review

In this paper, we examine the applicability of spatial optimization as a generative modelling technique for sustainable land‐use allocation. Specifically, we test whether spatial optimization can be used to generate a number of compromise spatial alternatives that are both feasible and different from each other. We present a new spatial multiobjective optimization model, which encourages efficient utilization of urban space through infill development, compatibility of adjacent land uses, and defensible redevelopment. The model uses a density‐based design constraint developed by the authors. The constraint imposes a predefined level of consistent neighbourhood development to promote contiguity and compactness of urban areas. First, the model is tested on a hypothetical example. Further, we demonstrate a real‐world application of the model to land‐use planning in Chelan, a small environmental amenity town in the north‐central region of the State of Washington, USA. The results indicate that spatial optimization is a promising method for generating land‐use alternatives for further consideration in spatial decision‐making.     

Possibilities and limits of prospective GIS land cover modelling—a compared case study: Garrotxes (France) and Alta Alpujarra Granadina (Spain)

This study focuses on the possibilities and the limits of a prospective GIS land cover modelling applied to two case studies (France and Spain). The methodology, based on available GIS tools, consists of using earlier land cover maps and relevant environmental factors (calibration data) to model actual, known land cover to validate the model. The model aggregates Markov chain analysis for time prediction and multi‐critera evaluation, multi‐objective evaluation and cellular automata to perform spatial contiguity of modelled land cover scores. The first results give an accurate, pixel by pixel prediction rate of approximately 75%. An important issue of this study consists of analysing prediction residues to improve the model.     

A land-use-based modelling chain to assess the impacts of Natural Water Retention Measures on Europe’s Green Infrastructure

This article details the process of integrating models to answer a specific policy-driven question: ‘What could be the impact of proposed Natural Water Retention Measures (NWRMs) on Europe’s Green Infrastructure (GI)?’ It describes the new Land Use Modelling Platform (LUMP), now enabling a high spatial scale (100-m) and large coverage (pan-European), whereby several sector-specific models contribute to assessing the impact of regional-level policy on a given spatial topic of concern. The configuration (land claims and land allocations modules) and calibration (accessibility and biophysical suitability) of the LUMP are explained. Four NWRM scenarios (riparian areas, afforestation, grassland and baseline scenario) are configured to run the simulations. For the reference: year 2006, the spatial representation of GI is based on land-use features of a refined version of the CORINE Land Cover (CLC), and resumed as connected components made of nodes and links.

Mathematical morphological image processing and network graph theory model, available from the free software package GUIDOS (the Joint Research Center of the European Commission), enabled the measurement of the GI connectivity and identified most critical links. Results show that the competition for land claimed by different economic sectors, combined with policy-driven rule-sets for the implementation of different NWRMs, yields very different results for the 2030 land-use projections, and subsequently for the morphology of GI. Three indicators associated with the morphology of GI are computed in order to assess the model outputs for 2030. The indicators are computed to answer the following questions: (1) How is the quantity of GI affected by each of the NWRM, and what proportion of that GI is most valuable? (2) What is the location of the most critical nodes and connectors of GI, and what land-use conversions occur under these? (3) Are the average components getting larger or smaller?

Whereas the grassland measure results in the largest net increase of GI, the afforestation measure results in the overall largest number of hectares of key nodes and links within the network. Land conversions occur under the critical GI nodes and links, with a large increase in agricultural areas, especially for the riparian measure under critical nodes and the grassland measure under critical links. Also predominant is the swapping of land from pasture to forest under critical links with the afforestation measure. The riparian measure most increases the average size of GI components, and all three measures contribute to bridging two large components which were divided in the 2006 land-use map, thus increasing the size of the largest component by more than 50%.     

Incorporating vegetation into visual exposure modelling in urban environments

Visual exposure modelling establishes the extent to which a nominated feature may be seen from a specified location. The advent of high-resolution light detection and ranging (LiDAR)-sourced elevation models has enabled visual exposure modelling to be applied in urban regions, for example, to calculate the field of view occupied by a landmark building when observed from a nearby street. Currently, visual exposure models access a single surface elevation model to establish the lines of sight (LoSs) between the observer and the landmark feature. This is a cause for concern in vegetated areas where trees are represented as solid protrusions in the surface model totally blocking the LoSs. Additionally, the observer's elevation, as read from the surface model, would be incorrectly set to the tree top height in those regions. The research presented here overcomes these issues by introducing a new visual exposure model, which accesses a bare earth terrain model, to establish the observer's true elevation even when passing through vegetated regions, a surface model for the city profile and an additional vegetation map. Where there is a difference between terrain and surface elevations, the vegetation map is consulted. In vegetated areas the LoS is permitted to continue its journey, either passing under the canopy with clear views or partially through it depending on foliage density, otherwise the LoS is terminated. This approach enables landmark visual exposure to be modelled more realistically, with consideration given to urban trees. The model's improvements are demonstrated through a number of real-world trials and compared to current visual exposure methods.     

A heuristic cellular automata approach for modelling urban land-use change based on simulated annealing

This article presents a novel cellular automata (CA) approach to simulate the spatio-temporal process of urban land-use change based on the simulated annealing (SA) algorithm. The SA algorithm enables dynamic optimisation of the CA's transition rules that would otherwise be difficult to configure using conventional mathematical methods. In this heuristic approach, an objective function is constructed based on a theoretical accumulative disagreement between the simulated land-use pattern and the actual land-use pattern derived from remotely sensed imagery. The function value that measures the mismatch between the actual and the simulated land-use patterns would be minimised randomly through the SA process. Hence, a set of attribution parameters that can be used in the CA model is achieved. An SA optimisation tool was developed using Matlab and incorporated into the cellular simulation in GIS to form an integrated SACA model. An application of the SACA model to simulate the spatio-temporal process of land-use change in Jinshan District of Shanghai Municipality, PR China, from 1992 to 2008 shows that this modelling approach is efficient and robust and can be used to reconstruct historical urban land-use patterns to assist with urban planning policy-making and actions. Comparison of the SACA model with a typical CA model based on a logistic regression method without the SA optimisation (also known as LogCA) shows that the SACA model generates better simulation results than the LogCA model, and the improvement of the SACA over the LogCA model is largely attributed to higher locational accuracy, a feature desirable in most spatially explicit simulations of geographical processes.