Non-stationary Approaches for Mapping Terrain and Assessing Prediction Uncertainty

It is well known that terrain may vary markedly over small areas and that statistics used to characterise spatial variation in terrain may be valid only over small areas. In geostatistical terminology, a non-stationary approach may be considered more appropriate than a stationary approach. In many applications, local variation is not accounted for sufficiently. This paper assesses potential benefits in using non-stationary geostatistical approaches for interpolation and for the assessment of uncertainty in predictions with implications for sampling design. Two main non-stationary approaches are employed in this paper dealing with (1) change in the mean and (2) change in the variogram across the region of interest. The relevant approaches are (1) kriging with a trend model (KT) using the variogram of residuals from local drift and (2) locally-adaptive variogram KT, both applied to a sampled photogrammetrically derived digital terrain model (DTM). The fractal dimension estimated locally from the double-log variogram is also mapped to illustrate how spatial variation changes across the data set. It is demonstrated that estimation of the variogram of residuals from local drift is worthwhile in this case for the characterisation of spatial variation. In addition, KT is shown to be useful for the assessment of uncertainty in predictions. This is shown to be true even when the sample grid is dense as is usually the case for remotely-sensed data. In addition, both ordinary kriging (OK) and KT are shown to provide more accurate predictions than inverse distance weighted (IDW) interpolation, used for comparative purposes.     

Self-organised clustering for road extraction in classified imagery

The extraction of road networks from digital imagery is a fundamental image analysis operation. Common problems encountered in automated road extraction include high sensitivity to typical scene clutter in high-resolution imagery, and inefficiency to meaningfully exploit multispectral imagery (MSI). With a ground sample distance (GSD) of less than 2 m per pixel, roads can be broadly described as elongated regions. We propose an approach of elongated region-based analysis for 2D road extraction from high-resolution imagery, which is suitable for MSI, and is insensitive to conventional edge definition. A self-organising road map (SORM) algorithm is presented, inspired from a specialised variation of Kohonen's self-organising map (SOM) neural network algorithm. A spectrally classified high-resolution image is assumed to be the input for our analysis. Our approach proceeds by performing spatial cluster analysis as a mid-level processing technique. This allows us to improve tolerance to road clutter in high-resolution images, and to minimise the effect on road extraction of common classification errors. This approach is designed in consideration of the emerging trend towards high-resolution multispectral sensors. Preliminary results demonstrate robust road extraction ability due to the non-local approach, when presented with noisy input.     

A Simple Temporal and Spatial Analysis of Flow in Complex Terrain in the Context of Wind Energy Modelling

A simple temporal and spatial analysisis done on wind speed and direction data from a number ofmeteorological towers separated by distances between roughly 1 and 100 kilometres. The analysis is done in the context of expected model error in wind energy calculations. The study first uses single point statistics to show the evolution of mean values with time. It is shown that strong seasonal signals are present and that stable means are achieved only after averaging periods of a year or more. The study then uses discrete Fourier transforms to show that significant amounts of spectral energy reside in modes with periods of a few days to less than a day. Frequency dependent cross correlation values are then derived and used to show how correlation between towers diminishes with increasing frequency. The mechanism responsible for this diminished correlation is shown through the comparison of cross-correlation phase as a function of frequency and its relationship to distance between towers. Error in wind energy estimates are shown to be strongly related to correlation and therefore distance over which the prediction is made. In summary, much of the inaccuracy in modelling flow in the context of wind energy calculations is due to a lack of scale separation between the deterministic part of the flow, which is well modelled, and that part of the flow that is stochastic at the length and time scales modelled.     

Genetic differentiation of populations of Greenlandic Arctic fox

Most microsatellites are very polymorphic. This makes them powerful markers for observing genetic differentiation between closely related populations. The population structure of the Greenlandic Arctic fox ( Alopex lagopus ) was studied genetically by analysing six polymorphic microsatellite loci of 75 foxes from four populations in different parts of Greenland. Genotypes were determined at the six loci for most of the individuals. Population differentiation was quantified in three different ways both within the total population and pairwise between all populations. The tests were Fisher's exact test, Rho estimates and Fst estimates, all of which supported a highly significant subdivision of the total population, and they showed significant differentiation in allele frequencies between all pairs of localities. It is concluded that the known long-distance migration of the Greenlandic Arctic fox has not resulted in complete genetic mixing of the populations. Fisher's exact test was also used to estimate levels of genetic differentiation between the two colour morphs: white and blue. No difference was found between allele frequencies of the two color morphs in any of the locations, and it was concluded that the white and blue morphs of the Greenlandic Arctic fox share the same habitat, at least during the mating season.     

Policy and practice in karst landscape protection: Bohol, the Philippines

The karst landscape in the interior of the Philippines' Bohol Province represents one of the world's premier kegelkarst (cone karst) environments. Government efforts to protect some of this karst, exemplified by the establishment of the Rajah Sikatuna National Park and the Chocolate Hills Natural Monument, have proven to be significant catalysts of social conflict. In Bohol there is a long history of traditional land tenure, which has recently been supplanted by a Westernized model. Protected area establishment is a response to deforestation, agricultural exploitation and uncontrolled quarrying. However, the imposition of protective legislation to prevent further degradation has disenfranchized and marginalized many local farmers and residents. The conflict between the obligation of the State to ensure environmental protection and the perceived property rights of landowners and farmers has provoked an escalation in civil unrest and armed conflict.     

The geomorphology of planetary calderas

Satellite-derived observations of the geomorphology of calderas on Earth, Mars and Venus can be used to learn more about shield volcanoes. Examples of terrestrial basaltic volcanoes from the Galapagos Islands, Hawaii, and the Comoro Islands show how these volcanoes contrast with examples found on Mars and Venus. Caldera structure, degree of infilling, and the location of vents on the flanks are used to interpret each volcano's recent history. The geometry of the caldera floor can be used to infer some of the characteristics of the magma storage system, and the orientation of the deep magma conduits. The formation of benches within the caldera and the effects of the caldera on the distribution of flank eruptions are considered, and it is evident that most calderas on the planets are/were dynamic features. Presently, deep calderas, with evidence of overflowing lavas and ponded lavas high in the caldera wall, show that these calderas were once shallow. Similarly, shallow calderas filled with ponded lavas are evidence that they were once deeper. It is probably a mistake, therefore, to place great significance on caldera depth with regard to the position, shape, or size of subsurface plumbing.     

Lunar Laser Ranging: Glorious Past And A Bright Future

Lunar Laser Ranging (LLR), a part of the NASA Apollo program, has beenon-going for more than 30 years. It provides the grist for a multi-disciplinarydata analysis mill. Results exist for solid Earth sciences, geodesy and geodynamics,solar system ephemerides, terrestrial and celestial reference frames, lunar physics,general relativity and gravitational theory. Combined with other data, it treatsprecession of the Earth's spin axis, lunar induced nutation, polar motion/Earthrotation, Earth orbit obliquity to the ecliptic, intersection of the celestial equatorwith the ecliptic, luni-solar solid body tides, lunar tidal deceleration, lunar physicaland free librations, structure of the moon and energy dissipation in the lunar interior.LLR provides input to lunar surface cartography and surveying, Earth station and lunar retroreflector location and motion, mass of the Earth-moon system, lunar and terrestrial gravity harmonics and Love numbers, relativistic geodesic precession, and the equivalence principle of general relativity. With the passive nature of the reflectors and steady improvement in observing equipment and data analysis, LLR continues to provide state-of-the-art results. Gains are steady as the data-base expands. After more than 30 years, LLR remains the only active Apollo experiment. It is important to recognize examples of efficient and cost effective progress of research. LLR is just such an example.