Use of an expert system to map forest soils from a geographical information system

Abstract

Mapping forest soils using conventional methods is time consuming and expensive. An expert system is described and applied to the mapping of five forest soil-landscape units formed on a single granitoid parent material. Three thematic maps were considered important in influencing the distribution of soils. The first showed the distribution of nine classes of native eucalypt forests, and the second and third were derived from a digital elevation model and represented slope gradient and a soil wetness index combined with topographical position. These layers were input to a raster based geographical information system (GIS) and then geometrically co-registered to a regular 30 m grid. From a knowledge of soil distributions, the relationships between the soil-landscape units and the three data layers were quantified by an experienced soil scientist and used as rules in a rule based expert system. The thematic layers accessed from the GIS provided data for the expert system to infer the forest soil-landscape unit most likely to occur at any given pixel. The soil-landscape map output by the expert system compared favourably with a conventional soil-landscape map generated using interpretation of aerial photographs.     

Road network detection from SPOT imagery for updating geographical information systems in the rural–urban fringe

Abstract

Visual interpretation of high-resolution satellite data has been useful for mapping linear features, such as roads and updating land-use changes. However, it would be beneficial to map new road networks digitally from satellite data to update digital databases using semi-automated techniques. In this paper, an algorithm called Gradient Direction Profile Analysis (GDPA) is used to extract road networks digitally from SPOT High Resolution Visible (HRV) panchromatic data. The roads generated are compared with a visual interpretation of the SPOT HRV multispectral and panchromatic data. The technique is most effective in areas where road development is relatively recent. This is due to the spectral consistency of new road networks. As new road networks are those of most interest to the land manager, this is a useful technique for updating digital road network files within a geographical information system of urban areas.     

Evolution and sedimentology of a channel fill in the sandy braided South Saskatchewan River and its comparison to the deposits of an adjacent compound bar

The depositional stratigraphy of within‐channel deposits in sandy braided rivers is dominated by a variety of barforms (both singular ‘unit’ bars and complex ‘compound’ bars), as well as the infill of individual channels (herein termed ‘channel fills’). The deposits of bars and channel fills define the key components of facies models for braided rivers and their within‐channel heterogeneity, knowledge of which is important for reservoir characterization. However, few studies have sought to address the question of whether the deposits of bars and channel fills can be readily differentiated from each other. This paper presents the first quantitative study to achieve this aim, using aerial images of an evolving modern sandy braided river and geophysical imaging of its subsurface deposits. Aerial photographs taken between 2000 and 2004 document the abandonment and fill of a 1·3 km long, 80 m wide anabranch channel in the sandy braided South Saskatchewan River, Canada. Upstream river regulation traps the majority of very fine sediment and there is little clay (< 1%) in the bed sediments. Channel abandonment was initiated by a series of unit bars that stalled and progressively blocked the anabranch entrance, together with dune deposition and stacking at the anabranch entrance and exit. Complete channel abandonment and subsequent fill of up to 3 m of sediment took approximately two years. Thirteen kilometres of ground‐penetrating radar surveys, coupled with 18 cores, were obtained over the channel fill and an adjacent 750 m long, 400 m wide, compound bar, enabling a quantitative analysis of the channel and bar deposits. Results show that, in terms of grain‐size trends, facies proportions and scale of deposits, there are only subtle differences between the channel fill and bar deposits which, therefore, renders them indistinguishable. Thus, it may be inappropriate to assign different geometric and sedimentological attributes to channel fill and bar facies in object‐based models of sandy braided river alluvial architecture.     

ENHANCED MIGRATION OF SEISMIC DATA1

The so-called ‘enhanced migration’ which uses diffraction tomography as the ‘repair tool’ for correction of amplitudes (reflection coefficients) of migrated sections is discussed. As with any linearized procedure, diffraction tomography requires knowledge of the initial model. It is suggested that the initial model is taken as the migrated image. It will be demonstrated that diffraction tomography applied to the data residuals improves the amplitudes of the migrated images. Migration is redefined as the reconstruction of the wavefront sets of distributions (reflection interfaces), and the inversion process as tomographic correction of migrated images.     

Pressure--Temperature--Time Paths of Regional Metamorphism I. Heat Transfer during the Evolution of Regions of Thickened Continental Crust

The development of regional metamorphism in areas of thickenedcontinental crust is investigated in terms of the major controlson regional-scale thermal regimes. These are: the total radiogenicheat supply within the thickened crust, the supply of heat fromthe mantle, the thermal conductivity of the medium and the lengthand time scales of erosion of the continental crust. The orogenicepisode is regarded as consisting of a relatively rapid phaseof crustal thickening, during which little temperature changeoccurs in individual rocks, followed by a lengthier phase oferosion, at the end of which the crust is at its original thickness.The principal features of pressure—temperature—time(PTt) paths followed by rocks in this environment are a periodof thermal relaxation, during which the temperature rises towardsthe higher geotherm that would be supported by the thickenedcrust, followed by a period of cooling as the rock approachesthe cold land surface. The temperature increase that occursis governed by the degree of thickening of the crust, its conductivityand the time that elapses before the rock is exhumed sufficientlyto be affected by the proximity of the cold upper boundary.For much of the parameter range considered, the heating phaseencompasses a considerable portion of the exhumation (decompression)part of the PTt path. In addition to the detailed calculationof PTt paths we present an idealized model of the thickeningand exhumation process, which may be used to make simple calculationsof the amount of heating to be expected during a given thickeningand exhumation episode and of the depth at which a rock willstart to cool on its ascent path. An important feature of thesePTt paths is that most of them lie within 50 °C of the maximumtemperature attained for one third or more of the total durationof their burial and uplift, and for a geologically plausiblerange of erosion rates the rocks do not begin to cool untilthey have completed 20 to 40 per cent of the total uplift theyexperience. Considerable melting of the continental crust isa likely consequence of thickening of crust with an averagecontinental geotherm. A companion paper discusses these resultsin the context of attempts to use metamorphic petrology datato give information on tectonic processes.     

Pressure--Temperature--Time Paths of Regional Metamorphism II. Their Inference and Interpretation using Mineral Assemblages in Metamorphic Rocks

A companion paper (England & Thompson, 1984a) investigatesthe pressure-temperature-time (PTt) paths followed by rocksundergoing burial metamorphism in continental thickening events.This paper discusses problems involved in inferring such pathsfrom the petrological data available in metamorphic rocks and—oncesuch paths are determined—how they may be interpretedin terms of the thermal budgets of metamorphism. Each of theprincipal facies series (glaucophane-jadeite, andalusite- sillimaniteand kyanite-sillimanite) may be encountered by rocks involvedin the thickening and erosion of continental crust in a regimeof average continental heat flow. The inference of a minimumthermal budget required for a given metamorphism depends stronglyon a knowledge of the PTt paths followed by rocks during themetamorphism. Discrimination between possible thermal regimesis greatly enhanced if portions of PTt paths, rather than singlePT points, are available, and additional constraint is possibleif these paths are supplemented by geochronological, structuraland heat flow data.     

INVERSION OF AN UNMIGRATED STACKED SECTION TO DETERMINE AN INTERVAL VELOCITY MODEL*

A seismic inversion procedure is developed that inverts data available from an unmigrated stacked section to produce an interval velocity model. It attempts to overcome some of the limitations of existing methods by using a generalized linear inversion technique. The inversion process incorporates several features: (i) Lateral interval velocity variations are permitted, (ii) A fast accurate forward model was developed, (iii) Input data is weighted according to the accuracy with which it has been acquired. The procedure is applied to seismic data from the Gippsland Basin, an area offshore South-East Australia.     

A NEW RECEPTOR MODEL:A DIRECT TRILINEAR DECOMPOSITION FOLLOWED BY A MATRIX RECONSTRUCTION

In many cases,monitoring data for ambient airborne particles can be organized in the form of a three-way data table with one way for chemical species,one for sampling periods and one for sites.A directtrilinear decomposition followed by a matrix reconstruction(DTDMR)is developed to analyze such adata table as a whole.The three-way data set is composed into three two-way matrices by a direct trilineardecomposition(DTD).The column vectors of each of the matrices are called‘source profiles’,‘emissionpatterns’and‘site coefficients’respectively.Particulate sources are identified by examining both theirsource profiles and emission patterns.After the sources have been identified,emission patterns and sitecoefficients are used to produce a three-way matrix that gives estimates of mass contributions of sourcesto the samples collected at every site in every period.By simulation study,not only has the method beenverified,but a good indicator has been found that shows the number of factors(i.e.sources)in thesystem.Unlike other receptor models,DTDMR does not require source profile data and does not involvetrial-and-error procedures.Since DTDMR identifies sources based on variations in two dimensions,it hasa higher potential to distinguish two sources that have similar chemical compositions.The DTDMRmodel has provided excellent results with simulated data and has been applied in a real world three-waydata set.     

Trace fossils from proglacial lake sediments

Trace fossils have been discovered in Middle Pleistocene proglacial lake sediments from St. Albans, Hertfordshire, England. They indicate that this environment was capable of supporting a variety of benthonic animals. The fossils are described and tentatively assigned to several invertebrate groups.