1. Introduction
Over the past decade, a number of high-resolution airborne geophysical surveys have been conducted across onshore UK [
1,
2]. These High Resolution Airborne Resource and Environmental (HiRES) surveys have typically acquired radiometric (gamma-ray spectroscopy), magnetic and electromagnetic (conductivity) measurements at 200 m line spacings and at low altitude (<60 m).
The UK, and recent Irish, surveys flown between 1998 and 2013 are shown in
Figure 1a. The original North Midlands survey of 1998 was largely acquired at lower spatial resolution (400 m line spacing) and at a higher elevation (90 m) than later surveys. The South West (SW) England survey is currently acquiring data. The surveys were not targeted at specific science themes (e.g., soil survey), but as a multi-faceted approach to acquiring uniform-quality, modern resource and environmental information.
Airborne radiometric data, in particular, have long been used in soil-related studies [
3,
4] due to signal interactions within the soil profile (usually < 60 cm). Across the soils of eastern England, the study by Rawlins, Webster and Lark [
5] indicated that the UK radiometric data are useful for making thematic maps of soil and confirmed that a large proportion of the thorium and potassium signals can be accounted for by parent material (bedrock). Additionally, the airborne data from Northern Ireland (NI) were used to study methods for improving current estimates of soil organic carbon (SOC) [
6]. Within the data sets acquired, it has become evident that the gamma ray attenuation considered here, displays a very evident association with organic soils, and with peat, in particular. Beamish [
7] noted that that attenuation across the soil profile is controlled by moisture content in conjunction with the density and porosity of the materials. The theory of attenuation is reprised here with special regard to peat and organic soils.
Peat assessments have become increasingly significant in relation to evaluation of sites for statutory protection, ecology and hydrogeology [
8], greenhouse gas emissions [
9] and carbon store assessments connected to climate change, land management and erosion [
10]. Conventional (non-geophysical) peat mapping within the UK may be considered relatively advanced, particularly with regard to recent assessments of soil carbon stocks [
11]. It may also be noted that many recent compilations are now regionally, and therefore spatially, devolved although UK-wide compilations exist [
12]. In order to provide a coherent assessment of the peat associations of the airborne data across the whole UK, we use a national scale mapping of deep (> 1 m) peat, generated at a 1:10 k scale and available digitally at a 1:50 k scale. These data are shown alongside the airborne survey areas in
Figure 1b.
Four survey areas at various scales, across England, Scotland, Wales and Northern Ireland were selected to demonstrate the ability of the airborne data to map peat zones. The geophysical response variations are examined across lowland areas (raised bogs, meres and fens, afforested peat) and upland areas of blanket bog, together with associated wetland zones. The four areas considered are arrowed in
Figure 1a.
2. Theory
When considering the amplitude of the radiometric response observed above the surface the simplest conceptual vertical model comprises two layers above bedrock. In the general case, the two upper layers would be defined by soil and superficial (Quaternary) deposits whose radiogenic content is assumed to be derived from the parent bedrock material. Whatever the near-surface material the observed response will be derived from a given radiometric source concentration (assumed vertically uniform in the first instance) that is primarily obtained from a shallow subsurface zone (often < 0.5 m).
A common approach to the modelling of gamma-ray fields considers the exponential absorption that characterizes the passage of electromagnetic radiation through a homogenous material:
where I
0 is the initial radiation intensity. The linear attenuation coefficient (μ) of the material is an intrinsic property of each material and would usually be associated with a specific element of given atomic number. The mass attenuation coefficient (μ
m) of the material is given by:
where ρ is density. Løvberg [
13] discusses the general principles of gamma-ray attenuation and the effect of the moisture content of materials. It is noted that Compton scattering (incoherent scattering) is the main significant attenuation interaction at the energies discussed here. In such circumstances the linear attenuation coefficient μ is proportional to the total number of electrons per unit volume of the material. Løvberg [
13] indicates that all elements with an atomic number less than 30 will have comparable mass attenuation coefficients. In the absence of water, soil, superficial and bedrock materials will have comparable attenuation coefficients at a given source concentration. Attenuation in dry materials is controlled by density alone. Any hydrogen supplied to the material as absorbed or free water (e.g., pore water) then generates an additional attenuation provided by the additional electron content. Methods of making routine airborne soil moisture measurements are described by Carroll [
14].
Soils are characterised as a three-phase system comprising solid, water and air. In the case of such mixtures it is necessary to replace the attenuation coefficient (μ) in
Equation (1) by a summation over the three phases. In the case of a uniform surface layer (a half-space), each phase becomes a fractional component and the material may be described in terms of its porosity and degree of saturation [
7]. Following Beamish [
7],
Figure 2a shows the attenuation curves obtained for two mineral soils (dry bulk densities of 1.1 to 1.6 g/cm
3). A saturation of 20% and a porosity of 20% are used for both. It is observed that 90% of the signal (a common attenuation reference level) is observed between depths of 33 and 47 cm. The attenuation curve observed for fresh water is located between those of the mineral soils and 90% attenuation is observed at a depth of 40 cm.
A standard reference vertical peat profile is provided by Clymo [
15] and this is used here when assigning dry bulk density values. The reference model displays dry density values down to 0.03 g/cm
3 in the uppermost layer (the acrotelm) which then increase to 0.12 g/cm
3 in the deeper catotelm layer (>12 cm). Here we employ a uniform value of 0.1 g/cm
3 for simplicity. A peat porosity value of 80% is taken as typical for the depth range from 0 to 100 cm following the study of Kechavarzi, Dawson and Leeds-Harrison [
16].
Figure 2a shows the attenuation observed in “dry” (20% saturation), “wet” (80% saturation) and fully saturated peat. The attenuation curves for wet peats indicate that the 90% threshold value is obtained at depths of between 50 and 60 cm. The “dry peat” attenuation curve is distinct from all other soils. The very low density and low water content provide an attenuation curve that would extend the 90% threshold value to several meters. The “wet peat” curves indicate that the at-surface radiometric response will have limited attenuation sensitivity to depths greater than 60 cm (for a uniformly radiogenic peat). The material properties, down to ∼60 cm, that control the degree of attenuation are dry density, porosity and degree of saturation. It is acknowledged that the peat profile may contain horizons such as the density variation noted in the Clymo model. Variations in other peat properties, such as carbon content, will only influence the radiometric response by virtue of any associated change in density and/or porosity. The major horizon in peat is that defined by the density variation of the Clymo model and the limiting values of this (0.03 and 1.20 g/cm
3) produce very similar behavior to peat attenuation curves shown in
Figure 2b.
The results discussed imply that the at-surface radiogenic response will have limited sensitivity to “wet peat” at depths in excess of 60 cm since a radiometric peat attenuation of >85% will exist under normal conditions of 10% at-surface soil moisture. Theory therefore suggests that the attenuation characteristics observed across wet peatlands cannot generally be used to map variations in peat depth. In the case of at-surface “wet-thin” peats (e.g., <20 cm in thickness), a localised increase in count rates may be detectable.
Again following Beamish [
7],
Figure 2b compares attenuation with degree of saturation for a range of near-surface materials. The behaviour of an intact bedrock material with bulk density of 2.7 g/cm
3 (e.g., a limestone) with a porosity of 10% is shown; the curve displays the lowest sensitivity to degree of saturation and would represent relatively intact bedrock at outcrop. The increase in attenuation when the bedrock is weathered or fractured is indicated by the 60% porosity limestone curve.
Figure 2b also shows the saturation dependence of two reference mineral soils having bulk densities of 1.1 and 1.6 g/cm
3 to a full range of saturation conditions. A porosity of 20% has been used. At full saturation, attenuation is in the range 15 to 20% and a quasi-linear decay behaviour is observed. A further curve for a higher porosity mineral soil (e.g., a clay with 60% porosity) with a bulk density of 1.2 g/cm
3 again shows a rapid attenuation and the reduction in count rate approaches 60% at full saturation. The peat behaviour (80% porosity, density 0.1 g/cm
3) is also shown. The sensitivity of the attenuation curve for peat is clearly distinct from all the other materials with a reduction in count rate of 80% at a water saturation of only 10%. Given that saturation, or soil moisture levels, of <10% might be regarded as atypical (in the UK), it is evident that peat is distinguishable in relation to the other materials studied.
The type behaviour of soils (and bedrock) summarised in
Figure 2 is intended to provide a framework for the assessment of radiometric attenuation in different types of soils. It is anticipated that soil mixtures such as organo-mineral soils will display intermediate behaviour.
Figure 2a is used to demonstrate the principal thickness of the material contributing to the soil surface gamma ray flux. Using a 90% contribution rule, the thicknesses range from 20 cm (limestone) to 60 cm (wet peat). In the unusual case of dry peat, the principal thickness may extend to several metres. Wet peat, when less than ∼20 cm thick, may begin to show significantly different (larger) fluxes than thicker zones. The same is true for the other non-organic materials but the behaviour would be confined to smaller thickness (e.g., <40 cm) as indicated in
Figure 2a.
Given sufficiently thick materials, the attenuation of gamma rays in response to the degree of saturation would be as shown in
Figure 2b. Mineral soils, with low porosity, would potentially display quite limited attenuation characteristics while any higher porosity materials (clays, weathered bedrock) would provide greater attenuation towards higher saturation levels. The clear outliers in behaviour are the high porosity, low density organic soils. The attenuation response is non-linear and appears highly sensitive to the amount of water in the interval from zero to 10% saturation. Beyond this interval (say >20% saturation) the organic soils would display essentially the same flux at all saturations. In comparative terms, the organic soils appear to be distinguishable by their gamma ray attenuation characteristics.
5. Conclusions
This study has considered gamma ray attenuation characteristics using simple theory and modern airborne radiometric data across a range of peat environments across the UK and into Ireland. The attenuation theory, reviewed here, predicts that the behaviour of wet peat is distinct from most other soil types. The case studies undertaken have largely confirmed this. The behaviour of peat is made distinct due primarily to the high water content and low density of these soils. Theory also predicts that the attenuation levels observed across wet peatlands cannot, in general, be used to map variations in peat thickness. Research in this regard [
29] continues to be further developed.
In order to demonstrate the peat mapping capability of the airborne data, a national-scale 1:50 k mapping of thick (>1 m) peat was employed. The studies broadly indicate that this mapping may be optimum in relation to the assessment of the airborne data. The variability in available control (peat-edge mapping) from alternative soil and land-use databases has also been demonstrated. A simple correlation of the geophysical information with peat zone edges is therefore not a precise procedure.
A modern UK standard for HiRES airborne acquisition uses 200 m line spacing. These data have a limited resolution in relation to small bodies of peat. Inherent sensor footprint resolution also limits lateral resolution at body edges. If finer detail is required, then 50 m line spacing can be adopted to achieve higher resolution; this resolution could potentially identify smaller peat bodies that lie outside of designated sites or that have been overlooked during geological mapping due to their small size. This type of systematic approach could help improve knowledge on the total coverage of peat across the UK.
At larger scales, both the soil and the bedrock variations (responsible for the radiometric amplitude level) require to be jointly assessed. The continuity of the airborne data particularly across wetland areas, taking in both lowland and upland contexts, is particularly useful. Given the clear association between radiometric attenuation responses and peat and wetland site boundaries, the data have the potential to identify additional areas where similar degrees of saturation exist. The methodology could therefore be used to identify new areas of peat, or rather areas of peat (e.g., small isolated bodies) that have not previously been mapped. It is also worth noting that many modern airborne geophysical data sets now exist, often over areas far larger than the UK, and across which peat and wetland mapping are probably less advanced than in UK.
It has been noted that the radiometric data do not map soils per se. Theory predicts that the attenuation characteristics observed are primarily associated with water content. Although the data contain information on the water content of all soils, the effect is most pronounced and therefore most detectable in peat soils. The intra-peat variations observed in this study are interpreted as variations in moisture content, with the lowest amplitude zones corresponding to 100% saturation. Density variations in the upper peat profile may also contribute and cannot be ruled-out. The conversion of the geophysical responses to degree of saturation would only be possible with ground calibrations performed at a scale appropriate to the airborne measurements.