Kimberly P. Van Niel

Effect of error in the DEM on environmental variables for predictive vegetation modelling

Abstract Question: Predictive vegetation modelling relies on the use of environmental variables, which are usually derived from a base data set with some level of error, and this error is propagated to any subsequently derived environmental variables. The question for this study is: What is the level of error and uncertainty in environmental variables based on the error propagated from a Digital Elevation Model (DEM) and how does it vary for both direct and indirect variables? Location: Kioloa region, New South Wales, Australia Methods: The level of error in a DEM is assessed and used to develop an error model for analysing error propagation to derived environmental variables. We tested both indirect (elevation, slope, aspect, topographic position) and direct (average air temperature, net solar radiation, and topographic wetness index) variables for their robustness to propagated error from the DEM. Results: It is shown that the direct environmental variable net solar radiation is less affected by error in the DEM than the indirect variables aspect and slope, but that regional conditions such as slope steepness and cloudiness can influence this outcome. However, the indirect environmental variable topographic position was less affected by error in the DEM than topographic wetness index. Interestingly, the results disagreed with the current assumption that indirect variables are necessarily less sensitive to propagated error because they are less derived. Conclusions: The results indicate that variables exhibit both systematic bias and instability under uncertainty. There is a clear need to consider the sensitivity of variables to error in their base data sets in addition to the question of whether to use direct or indirect variables. Abbreviations: AML = Arc/INFO Macro Language; DEM = Digital Elevation Model; GPS = Global Positioning System; HDOP = Horizontal Dilution of Precision; TWI = Topographic Wetness Index; VDOP = Vertical Dilution of Precision.

Gambling with randomness: the use of pseudo-random number generators in GIS

Analyses within the field of GIS are increasingly applying stochastic methods and systems that make use of pseudo-random number generators (PRNGs). Examples include Monte Carlo techniques, dynamic modelling, stochastic simulation, artificial life and simulated data development. PRNGs have inherent biases, and this will in turn bias any analyses using them. Therefore, the validity of stochastic analyses is reliant on the PRNG employed. Despite this, the effect of PRNGs in spatial analyses has never been completely explored, particularly a comparison of different PRNGs. Exacerbating the problem is that GIS articles applying Monte Carlo or other stochastic methods rarely report which PRNG is employed. It thus appears likely that GIS researchers rarely, if ever, check the suitability of the PRNG employed for their analyses or simulations. This paper presents a discussion of some of the characteristics of PRNGs and specific issues from a geospatial standpoint, including a demonstration of the differences in the results of a Monte Carlo analysis obtained using two different PRNGs. It then makes recommendations for the application of PRNGs in spatial analyses, including recommending specific PRNGs that have attributes appropriate for geospatial analysis. The paper concludes with a call for more research into the application of PRNGs to spatial analyses to fully understand the impact of biases, especially before they are routinely used in the wider spatial analysis community.

PREDICTIVE VEGETATION MODELING FOR CONSERVATION: IMPACT OF ERROR PROPAGATION FROM DIGITAL ELEVATION DATA

The effect of digital elevation model (DEM) error on environmental variables, and subsequently on predictive habitat models, has not been explored. Based on an error analysis of a DEM, multiple error realizations of the DEM were created and used to develop both direct and indirect environmental variables for input to predictive habitat models. The study explores the effects of DEM error and the resultant uncertainty of results on typical steps in the modeling procedure for prediction of vegetation species presence/absence. Results indicate that all of these steps and results, including the statistical significance of environmental variables, shapes of species response curves in generalized additive models (GAMs), stepwise model selection, coefficients and standard errors for generalized linear models (GLMs), prediction accuracy (Cohens kappa and AUC), and spatial extent of predictions, were greatly affected by this type of error. Error in the DEM can affect the reliability of interpretations of model results and level of accuracy in predictions, as well as the spatial extent of the predictions. We suggest that the sensitivity of DEM-derived environmental variables to error in the DEM should be considered before including them in the modeling processes.

Spatial Structure of Seagrass Suggests That Size-Dependent Plant Traits Have a Strong Influence on the Distribution and Maintenance of Tropical Multispecies Meadows

Background Seagrass species in the tropics occur in multispecies meadows. How these meadows are maintained through species co-existence and what their ecological drivers may be has been an overarching question in seagrass biogeography. In this study, we quantify the spatial structure of four co-existing species and infer potential ecological processes from these structures. Methods and Results Species presence/absence data were collected using underwater towed and dropped video cameras in Pulau Tinggi, Malaysia. The geostatistical method, utilizing semivariograms, was used to describe the spatial structure of Halophila spp, Halodule uninervis, Syringodium isoetifolium and Cymodocea serrulata. Species had spatial patterns that were oriented in the along-shore and across-shore directions, nested with larger species in meadow interiors, and consisted of multiple structures that indicate the influence of 2–3 underlying processes. The Linear Model of Coregionalization (LMC) was used to estimate the amount of variance contributing to the presence of a species at specific spatial scales. These distances were <2.5 m (micro-scale), 2.5–50 m (fine-scale) and >50 m (broad-scale) in the along-shore; and <2.5 m (micro-scale), 2.5–140 m (fine-scale) and >140 m (broad-scale) in the across-shore. The LMC suggests that smaller species (Halophila spp and H. uninervis) were most influenced by broad-scale processes such as hydrodynamics and water depth whereas large, localised species (S. isoetifolium and C. serrulata) were more influenced by finer-scale processes such as sediment burial, seagrass colonization and growth, and physical disturbance. Conclusion In this study, we provide evidence that spatial structure is distinct even when species occur in well-mixed multispecies meadows, and we suggest that size-dependent plant traits have a strong influence on the distribution and maintenance of tropical marine plant communities. This study offers a contrast from previous spatial models of seagrasses which have largely focused on monospecific temperate meadows.

Modelling Deep Water Habitats to Develop a Spatially Explicit, Fine Scale Understanding of the Distribution of the Western Rock Lobster, Panulirus cygnus

Background The western rock lobster, Panulirus cygnus, is endemic to Western Australia and supports substantial commercial and recreational fisheries. Due to and its wide distribution and the commercial and recreational importance of the species a key component of managing western rock lobster is understanding the ecological processes and interactions that may influence lobster abundance and distribution. Using terrain analyses and distribution models of substrate and benthic biota, we assess the physical drivers that influence the distribution of lobsters at a key fishery site. Methods and Findings Using data collected from hydroacoustic and towed video surveys, 20 variables (including geophysical, substrate and biota variables) were developed to predict the distributions of substrate type (three classes of reef, rhodoliths and sand) and dominant biota (kelp, sessile invertebrates and macroalgae) within a 40 km2 area about 30 km off the west Australian coast. Lobster presence/absence data were collected within this area using georeferenced pots. These datasets were used to develop a classification tree model for predicting the distribution of the western rock lobster. Interestingly, kelp and reef were not selected as predictors. Instead, the model selected geophysical and geomorphic scalar variables, which emphasise a mix of terrain within limited distances. The model of lobster presence had an adjusted D2 of 64 and an 80% correct classification. Conclusions Species distribution models indicate that juxtaposition in fine scale terrain is most important to the western rock lobster. While key features like kelp and reef may be important to lobster distribution at a broad scale, it is the fine scale features in terrain that are likely to define its ecological niche. Determining the most appropriate landscape configuration and scale will be essential to refining niche habitats and will aid in selecting appropriate sites for protecting critical lobster habitats.

A low-altitude mountain range as an important refugium for two narrow endemics in the Southwest Australian Floristic Region biodiversity hotspot

Background and Aims Low-altitude mountains constitute important centres of diversity in landscapes with little topographic variation, such as the Southwest Australian Floristic Region (SWAFR). They also provide unique climatic and edaphic conditions that may allow them to function as refugia. We investigate whether the Porongurups (altitude 655 m) in the SWAFR will provide a refugium for the endemic Ornduffia calthifolia and O. marchantii under forecast climate change. Methods We used species distribution modelling based on WorldClim climatic data, 30-m elevation data and a 2-m-resolution LiDAR-derived digital elevation model (DEM) to predict current and future distributions of the Ornduffia species at local and regional scales based on 605 field-based abundance estimates. Future distributions were forecast using RCP2.6 and RCP4.5 projections. To determine whether local edaphic and biotic factors impact these forecasts, we tested whether soil depth and vegetation height were significant predictors of abundance using generalized additive models (GAMs). Key Results Species distribution modelling revealed the importance of elevation and topographic variables at the local scale for determining distributions of both species, which also preferred shadier locations and higher slopes. However, O. calthifolia occurred at higher (cooler) elevations with rugged, concave topography, while O. marchantii occurred in disturbed sites at lower locations with less rugged, convex topography. Under future climates both species are likely to severely contract under the milder RCP2.6 projection (approx. 2 °C of global warming), but are unlikely to persist if warming is more severe (RCP4.5). GAMs showed that soil depth and vegetation height are important predictors of O. calthifolia and O. marchantii distributions, respectively. Conclusions The Porongurups constitute an important refugium for O. calthifolia and O. marchantii, but limits to this capacity may be reached if global warming exceeds 2 °C. This capacity is moderated at local scales by biotic and edaphic factors.

Is climate change shifting the poleward limit of mangroves

Ecological (poleward) regime shifts are a predicted response to climate change and have been well documented in terrestrial and more recently ocean species. Coastal zones are amongst the most susceptible ecosystems to the impacts of climate change, yet studies particularly focused on mangroves are lacking. Recent studies have highlighted the critical ecosystem services mangroves provide, yet there is a lack of data on temporal global population response. This study tests the notion that mangroves are migrating poleward at their biogeographical limits across the globe in line with climate change. A coupled systematic approach utilising literature and land surface and air temperature data was used to determine and validate the global poleward extent of the mangrove population. Our findings indicate that whilst temperature (land and air) have both increased across the analysed time periods, the data we located showed that mangroves were not consistently extending their latitudinal range across the globe. Mangroves, unlike other marine and terrestrial taxa, do not appear to be experiencing a poleward range expansion despite warming occurring at the present distributional limits. Understanding failure for mangroves to realise the global expansion facilitated by climate warming may require a focus on local constraints, including local anthropogenic pressures and impacts, oceanographic, hydrological, and topographical conditions.

There is no good excuse for a bad random number generator: a reply to Barry

We respond to Barrys comments on Van Niel and Laffan (2003) in this journal. Barry is correct on one point, but the remainder of his critique is based on opinion and is not supported by the literature. We conclude with an update on the current application of pseudorandom number generators (PRNGs) in geospatial analysis and its implementation in current geographical information system (GIS) software systems.