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.     

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.     

Identifying linear and non-linear behaviour in reduced complexity modelling output using surrogate data methods

Validation of model output is an important issue in the environmental sciences. This is particularly the case for reduced complexity modelling approaches where verification of the underlying equations is often problematic. Hence, methods that go beyond validation based on average values or the histogram of model output, are clearly advantageous. In this paper we show that a validation method based on the serial properties of time–space data, and developed in geomorphology, is a useful alternative but is nevertheless still insensitive to certain asymmetric characteristics of the signal. Hence, use of this method for validation may be enhanced by the adoption of an accompanying test for asymmetry in the model output.     

3D gravity modelling for Anyongbok Seamount in the East Sea

A seamount chain with an approximately WNW trend is observed in the northeastern Ulleung Basin. It has been argued that these seamounts, including two islands called Ulleung and Dok islands, were formed by a hotspot process or by ridge related volcanism. Many geological and geophysical studies have been done for all the seamounts and islands in the chain except Anyongbok Seamount, which is close to the proposed spreading ridge. We first report morphological characteristics, sediment distribution patterns, and the crustal thickness of Anyongbok Seamount using multibeam bathymetry data, seismic reflection profiles, and 3D gravity modeling. The morphology of Anyongbok Seamount shows a cone shaped feature and is characterized by the development of many flank cones and flank rift zones. The estimated surface volume is about 60 km³, and implies that the seamount is smaller than the other seamounts in the chain. No sediments have been observed on the seamount except the lower slope, which is covered by more than 1,000 m of strata. The crustal structure obtained from a 3D gravity modeling (GFR = 3.11, SD 3.82 = mGal) suggests that the seamount was formed around the boundary of the Ulleung Plateau and the Ulleung Basin, and the estimated crustal thickness is about 20 km, which is a little thicker than other nearby seamounts distributed along the northeastern boundary of the Ulleung Basin. This significant crustal thickness also implies that Anyongbok Seamount might not be related to ridge volcanism.     

The morphodynamic modelling of tidal sand waves on the shoreface

An artificial sand wave on the Dutch shoreface of the North Sea has been studied in conditions with relatively strong tidal currents in the range of 0.5 to 1 m/s and sediments in the medium sand size range of 0.2 to 0.5 mm. The sand wave is perpendicular to the tidal current and has a maximum height and length of the order of 5 m and 1 km, respectively. The sand wave is dynamically active and shows migration rates of the order of a few metres per year. A numerical morphodynamic model (DELFT3D model) has been used to simulate the morphological behaviour of the sand wave in the North Sea. This model approach is based on the numerical solution of the three-dimensional shallow water equations in combination with a surface wave propagation model (wind waves) and the advection–diffusion equation for the sediment particles with online bed updating after each time step. The model results show that the sand wave grows in the case of dominant bed-load transport (weak tidal currents; relatively coarse sediment; small roughness height; low waves) and that the sand wave decays in the case of dominant suspended transport (strong currents, relatively fine sediment, large roughness height; storm waves).     

Data assimilation for a coastal area morphodynamic model: Morecambe Bay

This paper investigates the use of data assimilation in coastal area morphodynamic modelling using Morecambe Bay as a study site. A simple model of the bay has been enhanced with a data assimilation scheme to better predict large-scale changes in bathymetry observed in the bay over a 3-year period. The 2DH decoupled morphodynamic model developed for the work is described, as is the optimal interpolation scheme used to assimilate waterline observations into the model run. Each waterline was acquired from a SAR satellite image and is essentially a contour of the bathymetry at some level within the inter-tidal zone of the bay. For model parameters calibrated against validation observations, model performance is good, even without data assimilation. However the use of data assimilation successfully compensates for a particular failing of the model, and helps to keep the model bathymetry on track. It also improves the ability of the model to predict future bathymetry. Although the benefits of data assimilation are demonstrated using waterline observations, any observations of morphology could potentially be used. These results suggest that data assimilation should be considered for use in future coastal area morphodynamic models.     

Elucidating the dynamics and mixing agents of a shallow fjord through age tracer modelling

The mixing agents and their role in the dynamics of a shallow fjord are elucidated through an Eulerian implementation of artificial tracers in a three-dimensional hydrodynamic model. The time scales of vertical mixing in this shallow estuary are short, and the artificial tracers are utilized in order to reveal information not detectable in the temperature or salinity fields. The fjord's response to external forcing is investigated through a series of model experiments in which we quantify vertical mixing, transport time scales of fresh water runoff and estuarine circulation in relation to external forcing.Using age tracers released at surface and bottom, we quantify the time scales of downward mixing of surface water and upward mixing of bottom water. Wind is shown to be the major agent for vertical mixing at nearly all depth levels in the fjord, whereas the tide or external sea level forcing is a minor agent and only occasionally more important just close to the bottom. The time scale of vertical mixing of surface water to the bottom or ventilation time scale of bottom water is estimated to be in the range 0.7 h to 9.0 days, with an average age of 2.7 days for the year 2004.The fjord receives fresh water from two streams entering the innermost part of the fjord, and the distribution and age of this water are studied using both ageing and conservative tracers. The salinity variations outside this fjord are large, and in contrast to the salinity, the artificial tracers provide a straight forward analysis of river water content. The ageing tracer is used to estimate transport time scales of river water (i.e. the time elapsed since the water left the river mouth). In May 2004, the typical age of river water leaving the fjord mouth is 5 days. As the major vertical mixing agent is wind, it controls the estuarine circulation and export of river water. When the wind stress is set to zero, the vertical mixing is reduced and the vertical salinity stratification is increased, and the river water can be effectively exported out of the fjord.We also analyse the river tracer fields and salinity field in relation to along estuary winds in order to detect signs of wind-induced straining of the along estuary density gradient. We find that events of down estuary winds are primarily associated with a reduced along estuary salinity gradient due to increased surface salinity in the innermost part of the fjord, and with an overall decrease in vertical stratification and river water content at the surface. Thus, our results show no apparent signs of wind-induced straining in this shallow fjord but instead they indicate increased levels of vertical mixing or upwelling during down estuary wind events.     

Ensuring confidence in radionuclide-based sediment chronologies and bioturbation rates

Sedimentary records of naturally occurring and fallout-derived radionuclides are widely used as tools for estimating both the ages of recent sediments and rates of sedimentation and bioturbation. Developing these records to the point of data interpretation requires careful sample collection, processing, analysis and data modeling. In this work, we document a number of potential pitfalls that can impact sediment core records and their interpretation. This paper is not intended as an exhaustive treatment of these potential problems. Rather, the emphasis is on potential problems that are not well documented in the literature, as follows: (1) the mere sampling of sediment cores at a resolution that is too coarse can result in an apparent diffusive mixing of the sedimentary record at rates comparable to diffusive bioturbation rates observed in many locations; (2) ²¹⁰Pb profiles in slowly accumulating sediments can easily be misinterpreted to be driven by sedimentation, when in fact bioturbation is the dominant control. Multiple isotopes of different half lives and/or origin may help to distinguish between these two possible interpretations; (3) apparent mixing can occur due simply to numerical artifacts inherent in the finite difference approximations of the advection diffusion equation used to model sedimentation and bioturbation. Model users need to be aware of this potential problem. Solutions to each of these potential pitfalls are offered to ensure the best possible sediment age estimates and/or sedimentation and bioturbation rates can be obtained.