|
|
||||||||
Biosystems Engineering, School of Agric., Food Sci. and Vet. Medicine, Univ. College Dublin, Earlsfort Terrace, Dublin 2, Ireland
* Corresponding author (john.connolly0{at}ucd.ie).
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: Abbreviations: CORINE, Coordination of Information on the Environment DIPM, Derived Irish Peat Map DTM, digital terrain map GIS, geographic information system GSM, General Soil Map HLM, high-level montane LLA, low-level Atlantic PAVC, peat-associated vegetation class PM, Peatland Map of Ireland PSA, peat soil association RB, raised bog.
| INTRODUCTION |
|---|
|
|
|---|
As a part of the EU strategy to reduce emissions, each member state must produce, on an annual basis, a report of anthropogenic CO2 emissions by sources and removals by sinks (Gugele et al., 2003). The main source categories are classified as energy, agriculture, industrial processes, and waste, while legitimate sinks include forests and other ecosystems that absorb C (Gugele et al., 2003; European Environment Agency, 2005). Peatlands have been recognized as making significant contributions to national C fluxes (Tolonen and Turunen, 1996; Clymo et al., 1998). They are not at present counted as terrestrial sinks, however (Roulet, 2000). The reason for this is that the sourcesink status of a peatland is both spatially and temporally variable (Griffis et al., 2000; Arneth et al., 2002; Kuhry and Vitt, 1996; Bubier et al., 1998; Smith et al., 2004). Despite these issues, peatlands are a very important, and dynamic, global C pool (Gorham, 1991) and much effort is now focusing on better understanding their response to weather (Alm et al., 1997; Blodau, 2002), their long-term relationship with climate (Thormann et al., 1998; Frolking et al., 2002; Belyea and Malmer, 2004), their spatial extent and function globally (Matthews and Fung, 1987; Sheng et al., 2004), and their function and extent at national scales (Vitt et al., 2000; Gajewski et al., 2001; Sippola and Yli-Halla, 2005; Proctor, 2005).
To evaluate the role of Irish peatlands within the national C budget, it is necessary to have a reliable estimate of their spatial extent and function. In this study, a rule-based decision tree methodology was developed, based on geographic information system (GIS) map algebra and Boolean logic, that can be used for integrating national-scale soil and land resource map data to estimate the spatial extent of peatland by type and disturbance. As new resource surveys become available they can be processed using the methodology, and the approach could be used for any soil mapping task that requires making a best spatial estimate from multiple, possibly flawed, map data.
The mapping of peatlands has always been a difficult endeavor. Vitt et al. (2000) used classification and aerial photography to determine peatland distributions. In their 1996 inventory of lowland raised bogs in Great Britain, Lindsay and Immirzi (1996) used three classes of bog: primary, secondary, and archaic. This classification system was used to calculate the extent of natural, degraded, and drained areas (Lindsay and Immirzi, 1996). Milton et al. (2005) used high-resolution multispectral data and information on bog surface microrelief to create a peatland classification map. Talkkari and Nevalainen (2003) used GIS techniques and ArcMap software (ESRI, Redlands, CA) to combine multisource data consisting of soil maps of different scales and numerous other data to create a soil database for Finland (Sippola and Yli-Halla, 2005). Peatlands in Northern Ireland were mapped using aerial photographs (Cruickshank and Tomlinson, 1990). Irish peatlands were mapped in the 1970s using field surveys and the amalgamation of older maps, resulting in the Peatland Map of Ireland (Hammond, 1979). They also appear on a national-scale soil association map (Gardiner and Radford, 1980a), but often the exact spatial distribution cannot be determined from this map. Peatlands also appear as a land cover class on the Coordination of Information on the Environment (CORINE) products (O'Sullivan, 1994; Bossard et al., 2000). All available data sources in Ireland are to some degree flawed and cannot be regarded as a reliable baseline of the extent of peatlands.
Irish peatlands are conventionally classified as fens, raised bogs, and blanket bogs (Hammond, 1979; Feehan and Donovan, 1996). A bog is a surface covered in peat that receives its nutrients only from the rain (Bellamy, 1986; Charman, 2002; Schouten, 2002). Raised bogs occur predominantly in the center of Ireland (Feehan and Donovan, 1996), and can be further subclassified into true midland, transitional, or industrial types depending on location and anthropogenic use (Hammond, 1979). Blanket bogs cover large areas of the western seaboard and upland areas and hence have been classified into low-level Atlantic (LLA) blanket bogs and high-level montane (HLM) blanket bogs by Hammond (1979). Fens are peatlands that receive nutrients from surface and groundwater as well as rainfall, and are now very rare, covering just 19607 ha in Ireland (Charman, 2002; Crushnell, 2005). Due to the very limited extent of fens, no explicit attention was addressed to fens in this study.
The objective of the study was to devise a method that would permit the best estimate of the spatial extent of peatlands in Ireland to be made using combinations of existing data sources. A rule-based decision tree was implemented in a GIS to determine locations with a high likelihood of being peat, and to classify the peat into three types: raised bog (RB), HLM blanket bog, and LLA blanket bog. The methodology was formulated so that as new resource data become available, they can be integrated to make improved spatial estimates.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
The GSM used in this survey was published in 1980 (Gardiner and Radford, 1980b). The map provides a representation of 10 soil associations, including two referenced as peat: blanket peat and basin peat (raised bog) (Coulter et al., 1998). There are also a number of soil associations that the legend indicates contain between 5 and 75% peat, but these are not spatially explicit with respect to the location of the peatlands on the map. The map is based on field surveys at various scales, aerial photographic interpretation, and geological maps (Coulter et al., 1998) and, as such, is not uniform in its precision.
Each map was georeferenced and converted from a vector map into a raster image (Zeiler, 1999) with a pixel resolution of 100 m. Georeferencing is the process of establishing a relationship between a raster's coordinate system and a real world coordinate system (Zeiler, 1999; Bonham-Carter, 1994). Initially, the CORINE map was georeferenced to the Irish National Grid (a Transverse Mercator projection) (O'Sullivan, 1994) by selecting common points of reference such as road or river junctions (Zeiler, 1999). Subsequently, both the GSM and PM were georeferenced to CORINE to allow each map to be standardized in the Irish National Grid projection. It was difficult, however, to accurately georeference the source GSM digital data; therefore >60 points of reference or tic points (Zeiler, 1999) throughout Ireland were needed. Before the PM could be rasterized, minor problems identified in its attributes were corrected. Each map was rasterized and reclassified as a binary map (Fig. 1; Arctur and Zeiler, 2004), i.e., 1 and 0 indicating peatland and nonpeatland pixels, respectively.
Step 2: Rule-Based Decision Tree and the Resulting Derived Peat Map
The concept behind the methodology presented was to use raster data (pixels) and that at least two sources of evidence should be present for a pixel to be classified as peat. A rule-based decision tree was then designed for use within a GIS (Fig. 1). The decision tree takes data from the best three digital map sources available (at the time of the study: PM, CORINE, and GSM) and examines them in combination with one another to extract the pixels with a maximum likelihood of representing peatland. As new data sources become available they can be integrated into the analysis process.
The decision tree functions as a filtering system based on Boolean logic, using sequential rules to isolate pixels that are likely to be peat. There are three main rules in the decision tree (Fig. 1); Rules 1 and 2 rely on primary data from the source maps, whereas Rule 3 relies on processing map data. Output maps depicting peat pixels are produced from each rule. Rule 1 is the most rigorous, resulting in one output map of pixels representing peat on all three source maps. Pixels containing peat that do not satisfy the criteria for Rule 1 filter down and are examined in Rule 2. This rule is less rigorous and examines the combinations of any two of the three source maps and produces three output maps (PM_GSM, PM_COR, and GSM_COR). Pixels that fail to be allocated as peat by the first two rules are examined by Rule 3. This rule was devised to prevent pixels that are peat from being erroneously discarded. According to Hammond (personal communication, 2003), there are several CORINE land cover classes predominantly found on peat. The CORINE classes that were likely to be associated with peat soil were designated as a peat-associated vegetation class (PAVC). A PAVC is a class represented by a pixel that is not classed as peat in CORINE but is classed as peat on both the PM and the GSM. In total, 26 CORINE classes were associated with peatlands. Most of these covered very small areas, but six CORINE classes were found to cover a significant area of peat (>10000 ha each). They were the CORINE land cover classes transitional woodland scrub, coniferous woodlands, natural grasslands, moors and heathland, pastures and diversity of agriculture, and natural vegetation. The PAVCs were then used as part of the third rule to identify pixels that were peat on PM or GSM and PAVC in CORINE. Likewise, the GSM contains soil associations that in some cases contain very high percentages of peat and which could be considered to be peat soil associations (PSA). The soil associations not classified as peat, but in some cases containing up to 75% peat, had not been examined during previous steps of the methodology. As a large number of pixels in the CORINE peat bogs class did not coincide with peat pixels on the GSM or the PM, if a CORINE "peat bog" pixel coincided with a PSA pixel, then it was treated as peat as opposed to being discarded as part of the third rule. Rule 3 produced three output maps (PM_PAVC, GSM_PAVC, and COR_PSA).
In total, seven output maps were produced from the first three rules of the decision tree. Each map was further evaluated against two qualifier conditions: (i) was the site a cut-away industrial peatland and therefore no longer a peat C store, and (ii) was the site on a slope of >25° where peat may not develop (Bellamy, 1986; Tallis, 1998; Holden, 2005). The PM (Hammond, 1979) was used to identify industrial peatlands and a digital terrain model (DTM; McGinnity et al., 2005) was used to identify slopes that were >25°.
The decision tree was designed as a filtering system in which pixels that do not comply with Rule 1 are not discarded but are examined by the subsequent rules. If a pixel is not peat after Rule 3, however, it is discarded. The result of the process is a map derived from multisource map data by analysis, modification, or combination. As new data become available they can be readily introduced into the analysis. The seven output maps were amalgamated in ArcGIS using map algebra (ESRI, 2005; Zeiler, 1999). Map algebra is conducted using the raster calculator function in the ArcGIS Spatial Analysis Tool. This function enabled each output map to be added together, thus creating the Derived Irish Peat Map (DIPM).
Step 3: Ground Truth Data and Testing of the Decision Tree Rules
Each output map produced from the decision methodology needed to be ground tested to determine the reliability of each of the decision tree steps and thus the DIPM. In many cases, output maps covered the same area. To test each map, it was necessary to isolate the output maps from each other. This process was completed using the raster calculator to subtract the overlying maps from one another. This isolated each output map, enabling ground testing; thus pixels allocated by Rule 1 were not evaluated as part of the testing of Rules 2 and 3. Each of the seven output maps was individually tested using field data. The field testing was based on determining the probability that a peat pixel on the derived map is peat on the ground. The criterion used for verification of the map was that peat of at least 300 mm in depth (Lindsay and Immirzi, 1996; Sippola and Yli-Halla, 1999; Charman, 2002; Tomlinson, 2005) should be the predominant soil in the area represented by a sample pixel.
Field testing was performed in three selected areas along an eastwest gradient from County Mayo in the west of Ireland (Study Site 1) through County Offaly (Study Site 2) to the Wicklow Mountains in the east (Study Site 3) (Fig. 2 ). The areas were selected to be representative of LLA blanket bogs, RB, and HLM blanket bogs, respectively. One hundred sample points were allocated to each of the three sites. Due to the relative inaccessibility of many peatlands, a sampling strategy involving access routes was developed. Within each study site, a digital road map (including roads, forest tracks, and walking tracks) was created in ArcGIS (Ordnance Survey Ireland, 1995). The roads were buffered to 300 m, creating a 600-m-wide buffer zone. This buffer zone was then used to delineate the area to be randomly sampled for each of the seven output maps. Within the buffered area, each output map was ranked in terms of the area that it covered, and the allocation of sampling points was area weighted. A number of industrial peatland (Indpeat) sites were also sampled (Table 1) to aid in testing the DIPM and original map sources.
|
|
Step 4: Peatland Classification
The output maps and the DIPM were produced as binary maps: peat and nonpeat. A map legend of HLM blanket bog, LLA blanket bog, and RB was developed using precipitation and elevation data. The precipitation data came from Fitzgerald and Forrestal's (1996) monthly averages for a 30-yr period from 1961 to 1990. These 30-yr monthly averages were summed to give a 30-yr annual average for each of the 600 rainfall stations in Ireland (Keane, 2001). These data were interpolated in ArcMap to give national coverage (ESRI, 2003). The elevation data were extracted from the Environmental Protection Agency of Ireland's DTM (McGinnity et al., 2005).
Precipitation provides an independent classifying property to distinguish RB from HLM and LLA blanket bog (Hammond, 1979), and elevation permits classification of LLA and HLM blanket bogs. A preliminary study using a random-point generator (Beyer, 2004) created 157 random points in peatland areas nationwide. While testing these points, it was noted that there was a lot of misclassification along the west coast and eastern Galway and Mayo, where LLA was classed as RB. To combat this, the random-point generator (Beyer, 2004) was used to create a total of 675 random points. These points were distributed nationwide but were weighted to peatland areas in the west that were misclassified in the preliminary study. Each point recorded elevation and precipitation data. Discriminant analysis was used on the 675 points to isolate blanket bogs (HLM and LLA) from raised bogs. In the database, Points 1 to 285 were classified as blanket bog, 286 to 517 were classified as a mix of raised bog and blanket bog, and 518 to 675 were classed as raised bog. The middle mixed class was extracted and subjected to further analysis using binary logistic regression in Minitab (Minitab, 2003). Based on the variables of elevation and precipitation, a formula was derived. The raster calculator in ArcMap was used to apply the formula to the elevation and precipitation data, resulting in the classification of ambiguous peatland areas in Tipperary, Galway, and Mayo.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Decision Rules and the Derived Irish Peat Map
The rule output maps produced by this methodology are only as good as the input maps and all the maps had faults of one kind or another. In general, with each additional rule in the decision process, the reliability of the prediction of a pixel actually being peat decreased (Table 1). Two of the output maps produced from Rule 3, the PM_PAVC and GSM_PAVC, were not included in the compilation of the DIPM because they were not reliable enough (36 and 18% ground truth agreement, respectively). The COR_PSA, the third output map produced from Rule 3, was an exception to this general trend with a reliability of 62%. The reliability of the remaining four output maps produced from Rules 1 and 2 extended from 55 to 90% (Table 1) and these were included in the DIPM (Fig. 3
). There was an overall reliability of 75% for the DIPM after exclusion of industrial peatlands and sites with steep slopes. This compares with the reliability of CORINE at 65%, the GSM at 50%, and the PM at 58%, estimated using the same field data. The DIPM is more reliable than the other available data sources at the time this study was conducted so the decision process has improved the estimate of peatland spatial extent in Ireland. The DIPM legend (Fig. 3) indicates the spatial extent of peatlands, the type of peatland in three classes (RB, HLM, and LLA) and the probability of the presence of peat in each pixel.
|
Peat Classification
The discriminant analysis and binary logistic regression used to classify the peatland map into three categories (HLM, LLA, and RB) was acceptable for most of the country. There are a small number of errors in a number of counties along the western seaboard where LLA has been misclassified as RB. The misclassification may reflect the lack of resolution in the spatial interpolation of the rainfall data used for discrimination. The spatial interpolation of rainfall data along the west coast was not good in places. This may be due to the physical landscape and spatial distribution of rainfall stations, resulting in rainfall values being interpolated across larger areas than elsewhere in the country. According to Keane (2001), this could lead to problems because rainfall accumulations between stations can vary considerably across short distances.
Peatland Disturbance
According to the field survey data, 204 of the 301 sites sampled were peat; out of these, 148, or 74%, were disturbed to some degree. A high number, 71%, of the disturbed sites had evidence of anthropogenic influences including forestry, industrial and domestic cutting (mainly historical), and grazing in HLM and LLA blanket bogs and in raised bogs. This compares with somewhat higher values reported by the Irish Peatland Conservation Council of around 82 and 92% disturbance for blanket and raised bogs, respectively (Irish Peatland Conservation Council, 1996).
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
Received for publication January 23, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N. Boylan, P. Jennings, and M. Long Peat slope failure in Ireland Quarterly Journal of Engineering Geology and Hydrogeology, February 1, 2008; 41(1): 93 - 108. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||