Assessing Geomorphic Processes and their Potential Relationship with Archaeological Artifact Exposure – NE Peloponnese, Greece
by
Pamela Elizabeth Tetford
A thesis submitted in conformity with the requirements for the degree of Master of Science Department of Geography University of Toronto
© Copyright by Pamela Elizabeth Tetford 2017
Assessing Geomorphic Processes and their Potential Relationship with Archaeological Artifact Exposure – NE Peloponnese, Greece
Pamela Elizabeth Tetford
Master of Science
Department of Geography University of Toronto
2017
Abstract
Landscape change can be important when encountering historical artifacts. Interactions between topography, climate and human activity shape a landscape, making sediment deposits from surface erosion and fluvial transport, and the artifacts they contain, important archives. There is a potential relationship between the rate of geomorphic processes and surface artifact density. This study compares spatially variable estimates of soil loss and stream energy, as indicators of high geomorphic activity, to surface artifact finds of the Western Argolid Regional Project (WARP).
Processes within the Inachos River watershed in the northeast Peloponnese, Greece, are quantified using the Unit Stream Power Erosion Deposition method (USPED) and the specific stream power approach in a Geographic Information Systems (GIS) environment. A statistically significant association is identified between surface erosion and artifact density, with the lowest artifact densities associated with the highest rates of soil loss. Knowledgeable interpretation of artifact distribution enables more accurate reconstruction of human settlement history.
ii
Acknowledgements
I would first like to thank my supervisor, Dr. Joseph Desloges, for his invaluable input over the past year and a half. Ongoing discussions surrounding various approaches and techniques within the field of fluvial geomorphology have challenged me and helped me develop a much stronger understanding of the discipline. His suggestions and critiques of various draft sections of this thesis have led to a much improved final version. I am also greatly appreciative for his assistance in the field, and for navigating the mountainous roads (goat trails?) of the Inachos River watershed with unparalleled confidence. His generous funding of this project has provided me with an extraordinary learning opportunity.
I would also like to thank my committee members, Drs. Joseph Desloges, George Arhonditsis and Dimitri Nakassis for providing valuable comments and suggestions that improve upon the final version of this thesis.
I also offer thanks to Tassos Venetikidis for his assistance with the translation of Greek documents and discussion regarding the tectonic activity of the northeast Peloponnese region.
I also extend my appreciation to the WARP team, for their generous hospitality during fieldwork in Greece. I thank Drs. Sarah James, Scott Gallimore and William Caraher so much for the knowledge they shared regarding the rich history and culture of the Peloponnese region, and taking the time to explain the theory and practices regarding archaeological survey methods. And to Dr. Dimitri Nakassis, I can’t begin to thank him enough for the incredible insight he provided into Greek history and culture, and his tireless assistance in the quest to find the elusive data necessary to make the completion of this project possible! His determination to find sieves and lab space to complete soil analysis was remarkable. I thank the British School of Athens for providing those laboratory facilities.
Lastly, I wish to thank my family. The collective patience and understanding that they have all shown has allowed me to devote the time necessary to complete this project. Particularly to my husband, Murray, and my daughter, Sydney, who have picked up the slack at home while “mom” was working away at the computer, and without whose support this would not have been possible…thank you so much!
iii
Table of Contents
Acknowledgements ...... iii
Table of Contents ...... iv
List of Tables ...... vii
List of Figures ...... viii
List of Appendices ...... ix
Chapter 1: Introduction ...... 1
1.1. Conceptual Understanding ...... 1
1.2. Defining the Problem ...... 2
1.3. Research Objectives ...... 4
Chapter 2: Theoretical Background ...... 5
2.1. Variable Soil Loss as an Assessment of Hillslope Stability ...... 5 2.1.1. Evolution of the Revised Universal Soil Loss Equation ...... 5 2.1.2. The RUSLE in a GIS Environment ...... 6 2.1.2.1. The Rainfall Erosivity Factor ...... 8 2.1.2.2. The Soil Erodibility Factor...... 10 2.1.2.3. Slope Length and Slope Steepness Factors ...... 11 2.1.2.4. The Vegetation Cover Factor ...... 12 2.1.2.5. The Land Management Practices/Conservation Factor ...... 13
2.2. Stream Power as an Indicator of Geomorphic Processes ...... 13 2.2.1. Downstream Discharge ...... 14 2.2.2. Bankfull Width ...... 15 2.2.3. Channel Gradient ...... 15
2.3. The Relationship between Geomorphic Processes and Artifact Exposure ...... 16
Chapter 3: Inachos River Watershed ...... 18
3.1. Location ...... 18
3.2. Climate and Hydrology ...... 19
3.3. Physiography ...... 21 iv
3.4. Current Land Use...... 22
3.5. Settlement History...... 22
Chapter 4: Methods ...... 25
4.1. Fieldwork ...... 25
4.2. Soil Analysis ...... 26
4.3. Delineation of the Watershed ...... 26
4.4. RUSLE Calculation ...... 26 4.4.1. Computation of the Rainfall Erosivity Factor, R...... 28 4.4.2. Computation of the Soil Erodibility Factor, K ...... 28 4.4.3. Computation of the Slope Length and Steepness Factor, LS ...... 28 4.4.4. Computation of the Cover Management Factor, C...... 28 4.4.5. Computation of the Support Practice Factor, P ...... 29
4.5. Specific Stream Power Calculation ...... 29 4.5.1. Determination of Slope ...... 32 4.5.2. Determination of Discharge and Width Relations...... 32 4.5.3. Mapping Specific Stream Power ...... 32
4.6. Comparing RUSLE to Artifact Density Distribution ...... 33
Chapter 5: Results ...... 35
5.1. Estimating Potential Soil Surface Erosion and Hillslope Surface Stability ...... 35 5.1.1. The R Factor (rainfall erosivity) ...... 35 5.1.2. The K Factor (soil erodibility) ...... 37 5.1.3. The LS Factor (slope length-steepness) ...... 39 5.1.4. The C Factor (land cover index) ...... 39 5.1.5. The P Factor (land management) ...... 40 5.1.6. The RUSLE ...... 40
5.2. Estimating Fluvial Geomorphic Processes ...... 41 5.2.1. Modern Bed Sediment of the Inachos River ...... 41 5.2.2. Bank Stratigraphy of the Inachos River ...... 42 5.2.3. Estimating Specific Stream Power ...... 44 5.2.3.1. Drainage Area-Discharge and Width Proxies ...... 44
v
Table of Contents 5.2.3.2. Longitudinal Profile ...... 45 5.2.3.3. Spatial Distribution of Specific Stream Power ...... 46
5.3. Assessing the Spatial Distribution of Geomorphic Processes ...... 47 5.3.1. Comparing the USPED to Artifact Distribution ...... 47 5.3.2. Regression of Density and Soil Loss using Continuous Data ...... 49 5.3.3. Analysis of Density and Soil Loss using Categorical Data ...... 50 5.3.4. Specific Stream Power Distribution Related to Artifact Density ...... 53
Chapter 6: Discussion ...... 55
6.1. Geomorphic Processes of the Upper Watershed ...... 56
6.2. Mid-Elevation Geomorphic Processes ...... 57
6.3. Geomorphic Processes of the Argolid Floodplain ...... 59
6.4. Geomorphic Processes within the WARP Survey Polygon ...... 59
Chapter 7: Conclusions ...... 62
7.1. Summary of Findings ...... 62
7.2. Future Research ...... 64
References ...... 65
Appendices ...... 70
vi
List of Tables
Table 3.1. - The chronology of cultural periods, occupation, land use and landscape .. 23
Table 5.1. - Descriptive statistics for the USPED and five factor components ...... 40
Table 5.2. - Descriptive statistics for artifact density, RUSLE and factors...... 47
Table 5.3. - Regression results of LN USPED and LN (Density +1) ...... 49
Table 5.4. - Fisher Least Significant Difference for categorized USPED ...... 51
Table 5.5. - Fisher Least Significant Difference for categorized USPED > 1 t ha-1 a-1 .. 53
vii
List of Figures
Figure 3.1. - Location of the Inachos River watershed ...... 19
Figure 3.2. - Mean monthly precipitation values for hydrograph stations ...... 20
Figure 4.1. - Location of study sites within the Inachos River watershed ...... 25
Figure 4.2. - Flowchart of RUSLE modelling process ...... 27
Figure 4.3. - Flowchart of the specific stream power modelling process ...... 31
Figure 4.4. - Artifact density per survey tract for the archaeological study polygon...... 33
Figure 5.1. - The RUSLE and individual factors as raster layers ...... 36
Figure 5.2. - Particle size distribution of small grain sediment (D < 2 mm) ...... 38
Figure 5.3. - River bank exposure (Site R5, south bank) ...... 42
Figure 5.4. - River bank exposure (Site R8, south bank) ...... 43
Figure 5.5. - Scour gullies into red-bed formations (Site R9) ...... 43
Figure 5.6. - Bankfull discharge and bankfull width proxies ...... 44
Figure 5.7. - Drainage area proxy models for bankfull width and bankfull discharge .... 45
Figure 5.8. - Spatial variability of specific stream power ...... 46
Figure 5.9. - Mean values of USPED and factors per survey tract ...... 48
Figure 5.10. - Scatterplot for LN (Density + 1) and LN USPED...... 50
Figure 5.11. - Results of one-way ANOVA for categorical USPED...... 51
Figure 5.12. - Results of one-way ANOVA for categorical USPED (A > 1.0 t ha-1 a-1) .. 52
Figure 5.13. - Spatial variability of specific stream power within the survey area ...... 54
viii
List of Appendices
Appendix A. Rainfall Erosivity, R factor ...... 70
Appendix B. Soil Erodibility, K factor ...... 70
Appendix C. Vegetative Cover, C factor ...... 72
Appendix D. Rapid Field Surveys ...... 73
Appendix E. Summary Statistics for Specific Stream Power ...... 74
Appendix F. Regression of USPED against Artifact Density...... 75
ix
Chapter 1: Introduction 1.1. Conceptual Understanding
Understanding a changing landscape, and the processes driving that change, can be important when considering the distribution of artifact finds (Fuchs, Lang & Wagner, 2004). Sediment erosion and transport are critical factors in the preservation and discovery of surface artifacts (Bevan & Conolly, 2009). The assessment of surface soil stability in survey areas assists understanding of the complex interface between human activity, surface processes, and artifact location (Gouma, van Wijnigaarden & Soetens, 2011). Better interpretation of the archaeological record and artifact dispersion enables greater accuracy in the reconstruction of human settlement history (Zananiri, Hademenos & Piteros, 2010).
Landscape form and composition are the consequence of complex interactions between topography, climate and human activity (van Andel, Zangger & Demitrack, 2010). These interactions influence processes that detach sediment, which is transported primarily by water and deposited downslope, often with large particles deposited first, then aggregates and fine particles further downslope (Garcia Rodriguez & Gimenez Suarez, 2012). The resulting sediment deposits are important archives of past events (van Andel et al., 2010). Regions of climate variability and/or extensive anthropogenic occupation are particularly vulnerable to soil erosion and subsequent sediment redistribution (Fuchs, 2007). In the Mediterranean, accelerated soil erosion and alluviation associated with early human exploitation has left ancient towns buried and has been well documented (van Andel et al., 2010). In this region, geomorphic processes often condition archaeological deposits by removing entire deposits, burying sites, or selectively transporting exposed artifacts (Byers, Hargiss & Byrd Finley, 2015). Geomorphic processes play a fundamental role in the distribution and exposure of archaeological finds (Gouma et al., 2011). Consequently, it has been suggested that there is a potentially strong relationship between the rate and intensity of geomorphic processes and site stability where relics from human settlement are exposed (Gouma et al., 2011).
1
1.2. Defining the Problem
Geomorphic investigation can offer information regarding a landscape’s morphology, and the interaction of surface processes, over a period of time (Mexia, 2015). To elucidate the mechanisms associated with long-term landscape evolution, and the behavior of earth’s materials, often requires the use of process studies from which conclusions regarding geomorphic change can be inferred (Trenhaile, 2010, p. 11). Soil erosion has been a dominant factor in the redistribution of sediment in the Mediterranean region since the inception of agriculture during the Neolithic period, with the impact of human activity intensifying the main natural physical forces such as rainfall impact, flowing water, wind, ice, temperature change, and gravity (Fuchs et al., 2004). When the rate of soil loss exceeds the rate at which soil formation may occur, land degradation will alter the physical landscape (Panagos et al., 2015b). To this end, the quantification of soil loss provides an indication of the spatial variability of potential landscape alteration occurring, while accounting for dynamic processes that entrain and transport materials across earth’s surface (Panagos et al., 2015b; Demirci & Karaburum, 2012).
Research quantifying the relationship between geomorphic work and archaeological artifact exposure is limited (Gouma et al., 2011). Bevan & Connolly (2009) argue that archaeology has been relatively slow to develop quantitative methods to examine spatially variable information regarding human deposits. However, spatial modeling can increase understanding of field survey data when developing hypotheses regarding human activities based on artifacts left behind (Bevan & Conolly, 2009). This study quantifies the relationship between soil loss and archaeological artifact exposure by examining the intensity and spatial variability of geomorphic activity in the artifact rich Argolid region of southern Greece. This semi-arid Mediterranean region is vulnerable to erosion (Garcia Rodriguez & Gimenez Suarez, 2012). Hot, dry summers and shoulder seasons, followed by intense winter rainfall, often produce rill erosion and gullying on steep, sparsely vegetated slopes (Ferreira & Panagopolous, 2014). An extensive history of human activity in the region has accelerated soil erosion beyond which would have been expected naturally (Elhag, 2015). This historically significant area contains a concentrated archive of surface artifacts dating to the Early Neolithic period (Fuchs, 2007).
2
This study identifies areas of more intense geomorphic activity (soil erosion) and compares them to the distribution of surface artifact finds in the Inachos River watershed of southern Greece. Surface soil stability is modelled, using Geographic Information Systems (GIS) and the Unit Stream Power Erosion and Deposition (USPED) approach to the Revised Universal Soil Loss Equation (RUSLE), to represent the 3-dimensional complexity of topographic processes occurring from the watershed boundary to the overland convergence that forms the main Inachos River channel. This RUSLE model has proven to effectively estimate erosion and soil loss in the semi-arid southern region of Greece (Kouli, Soupios & Vallianatos, 2009). Model inputs include a rainfall-runoff erosivity factor derived from historical precipitation records, a soil erodibility factor determined by soil sampling and soil maps, a slope length and slope steepness factor extracted from a 5 m Digital Elevation Model (DEM), and a vegetation cover factor and soil conservation practices factor derived from WorldView-3 remote sensing imagery. The RUSLE model presented here is coupled with a stream power model to address processes occurring within the main river channel. This coupling is a way of exploring the potential connections between hillslope derived material and sediment movement in the river valley. Specific stream power is modelled to assess the intensity and variability of erosion occurring within the Inachos River channel from fluvial energy. Specific stream power has been shown to effectively indicate fluvial geomorphic work (Phillips & Desloges, 2014). Model inputs include channel slope extracted from the 5 m DEM, and channel discharge and width estimated from rapid geomorphic field surveys, remote imagery and generally accepted discharge-drainage area relationships.
Modelled hillslope stability and stream channel energy are compared to artifact density records obtained from a spatially intensive, systematic fieldwalk survey conducted by the Western Argolid Regional Project (WARP). The WARP survey identifies and records surface artifact yield over a 30 km2 area of the central Inachos River watershed (WARP, 2015).
A GIS environment is used to calculate and efficiently display each element of this research on a 5 m x 5 m basis. The approach is to effectively illustrate the intensity and spatial variability of potentially high geomorphic activity and artifact yield within the Inachos River basin. The goal is to identify the spatially variability in these linkages with the aim of helping understand archaeological artifact exposure histories.
3
1.3. Research Objectives
The objective of this research project is to use generally accepted modelling techniques to demonstrate the spatial variation in soil erosion and hillslope stability, and stream energy. By doing so, this study aims to address the following questions:
1. What watershed attributes contribute the most to soil surface erosion and variations of hillslope stability in the Inachos River watershed, a mountainous Mediterranean region vulnerable to intense geomorphic processes? What are the potential soil erosion rates since occupation and how are they distributed?
2. Do the fluvial geomorphic processes of the Inachos River contribute to the erosion and transport of sediment in the watershed and how are they connected to adjacent hillslope inputs of soil and sediment?
3. Is there a relationship between the rate and intensity of geomorphic processes and site stability, and artifact distribution when artifacts from past human settlement are exposed in an artifact rich, Mediterranean region?
4
Chapter 2: Theoretical Background 2.1. Variable Soil Loss as an Assessment of Hillslope Stability
The erodibility of soil is described as the soil’s ability to resist mechanical separation and breakdown in its aggregate form (Hill & Schutt, 2000). As soil erodes, sediment is entrained, transported and deposited across a surface (Kouli et al., 2009). This is a natural process by physical forces (i.e. rain, wind, gravity, etc.) and can be accelerated by anthropogenic causes (Efthimiou, Lykouydi & Kavaritis, 2014). Regions with a Mediterranean climate are particularly vulnerable to soil loss due to erosion, with hot dry summer and shoulder seasons limiting vegetation growth followed by heavy erosive winter rainfall on steep, only partially vegetated slopes (Kouli et al., 2009). Erosion and land degradation in arid, semi-arid and sub-humid areas of the Mediterranean are consequences of climatic variation and human activities such as land clearance and agriculture (Poesen & Hooke, 1997). Landscapes of the Mediterranean have arguably suffered the most extensive degradation due to anthropogenic impact (Syvitski, 2003). It has been estimated that approximately 75% of the average sediment yield of a Mediterranean basin may be the result of human activity (Syvitski, 2003). Consequently, in an area with extensive human presence and a long agricultural history, such as the semi-arid northeast Peloponnese region of Greece, erosion is accelerated (van Andel et al., 2010), making examination of the downslope flux of sediment an appropriate diagnostic tool to predict land surface stability and geomorphic change (Tsara, Kosmas, Kirkby, Kosmas & Yassoglou, 2005).
2.1.1. Evolution of the Revised Universal Soil Loss Equation
Erosion prediction equations have been used since the mid 1940’s as a tool to assess soil erosion risk and predict soil erosion rates for various soils, land use, landscape characteristics and conservation practices (Spaeth Jr., Pierson Jr., Weltz, & Blackburn, 2003; Tsara et al., 2005). To predict long-term average soil loss (~20 years) due to sheet and rill erosion, Wischmeier and Smith (1965, 1978) developed the Universal Soil Loss Equation (USLE) to estimate long term average soil erosion rates as a conservation planning guide. Although scientifically sound, predictions using the USLE consistently produced higher values than observed values for ungrazed land (Spaeth et al., 2003). According to Renard, Foster, Weesies, and Porter (1991), the USLE did not explicitly represent the interaction of erosion processes. The USLE, as applied, did
5 not address soil deposition and gully erosion within a watershed, lacking the ability to account for deposition along hillslopes and in depressions and channels. Instead, it was designed for straight slope sections (Garcia Rodriguez & Gimenez Suarez, 2012).
Hydrology and erosion research advancements led to a revision of the USLE, the Revised Universal Soil Loss Equation (RUSLE), with indexed factors representing climate, soil, topography and land use, such that
A R K L S C P (1) where A is the average annual soil loss (t ha-1 a-1), R is the rainfall erosivity factor, K is the soil erodibility factor, L is the slope length factor, S is the slope steepness factor, C is the vegetation/crop cover management factor, and P is the soil conservation support practice factor (Kinnell, 2015). The RUSLE retains the quantification of the USLE but incorporates a number of “adjustments” in determining the weighting of individual factors in the model (Kinnell, 2015).
Improvements incorporate a 10-year frequency storm erosivity index (EI10) in the R factor, modifications to the K factor to reflect moisture extraction by growing crops, new equations for the L and S factors to reflect rill and interrill erosion, subfactors for evaluating the C factor, and new conservation practice values for the P factor (Renard & Freimund, 1994). Changes to the computation of the LS factor account for concavity and convexity using segments of irregular slopes (Mitasova, Hofierka, Zlocha & Iverson, 1996). While the RUSLE does not explicitly account for processes of detachment, transport and deposition, it does reflect the net effects of these processes in the quantification of soil loss (Spaeth, Pierson, Weltz & Blackburn, 2003). However, the methodology of subdividing landscapes to account for concavity and convexity still does not fully consider flow convergence and divergence (Oliveira, 2013). Although better guidelines are imposed for variation of terrain by the incorporation of a two dimensional profile shape, manual calculation of the RUSLE does not represent the three dimensional planform effects on erosion at the catchment scale (Garcia Rodriguez & Gimenez Suarez, 2012).
2.1.2. The RUSLE in a GIS Environment
The concept of incorporating the effect of flow discharge and sediment concentration, through calibration based on observed data, has been explored using algorithms to represent convergence
6 and divergence (Oliveira et al., 2013). In the RUSLE 3D, the upstream contributing drainage area is substituted for slope length to incorporate the impact of flow convergence (Garcia Rodriquez & Gimenez Suarez, 2012). The adoption of a GIS environment has enabled a more explicit representation of the spatial variability of erosion and deposition processes with overland flow (Garcia Rodriquez & Gimenez Suarez, 2012). The RUSLE 3D takes into account the upstream contributing area, whereas the Unit Stream Power Erosion and Deposition (USPED) model is derived using the contribution area and flow accumulation as a representative for water flow in calculation of the LS factor (Oliveira et al., 2013). The USPED assumes that the capacity of water flow determines the amount of sediment transported, and divergence and convergence are expressed by the change in sediment flow (Garcia Rodriguez & Gimenez Suarez, 2012). This model has been proven to better replicate surface processes (Gouma et al., 2011). Using GIS allows integration of topographic factors to provide effective assessment of erosion and deposition at a landscape scale (Mitasova et al., 1996). The RUSLE, therefore, is the most widely accepted and frequently used empirical soil erosion model to assess water erosion, predicting long-term average annual soil loss for specified land management classifications (Gouma et al., 2011; Garcia Rodriguez & Gimenez Suarez, 2012). When examining potential soil loss due to erosion, a watershed provides a defined hydrologic and topographic unit over which the RUSLE can be applied (Elhag, 2015). Using the RUSLE, in conjunction with GIS (and the USPED), and remote sensing, allows quantitative representation of geomorphic processes at a regional scale (Gouma et al., 2011). Remote sensing enables detection of land use change with relative ease and accuracy, and the assignment of attribute values to unsampled areas (Ferreira & Panagopoulos, 2014). Cell by cell processing allows better detailed analysis of individual contributors to soil erosion, such as soil type, slope and land use and enables quantification of a single factor’s influence on the total soil loss (Demirci & Karaburun, 2012). The RUSLE, using the USPED approach, therefore, permits better systematic examination of the spatial variability and role of model components, with an LS factor that more closely resembles expected surface runoff (Garcia Rodriguez & Gimenez Suarez, 2012; Oliveira et al., 2013).
Alternative models of several varieties have emerged in recent decades to assess sediment yield (physical models, conceptual models and empirical models) that differ primarily in complexity, data requirements, watershed characteristics and end use; however, the USPED approach to the RUSLE is useful, conceptually simple, and has fewer data requirements (Ferreira & 7
Panagopoulos, 2014). The USPED has been successfully applied in Crete (Kouli et al., 2009), a similar Mediterranean region to this study.
2.1.2.1. The Rainfall Erosivity Factor
The R factor represents the capability of rainfall to result in erosion based on rainfall intensity and depth for a given period (Kouli et al., 2009). Rainfall intensity, and the amount that occurs at each level of intensity, is determined from recorded rainfall data (Renard & Freimund, 1994). An increase in rainfall intensity, and/or amount, is reflected by an increase in the R factor (Demirci & Karaburun, 2012). The R factor is the sum of individual annual storm erosivity index values, EI, averaged over long periods of time (>20 years). It is defined as
1 n m R (E)k (I30)k (2) n j 1 k1 j
where E is the total storm kinetic energy, I30 is the maximum 30-minute rainfall intensity, j is the number of years used to produce the average, k is the number of storms in each year, n is the number of years averaged to obtain R, and m is the number of storms in each year (Renard & Freimund, 1994).
Lack of long term precipitation data in some areas can make the computation of R, and the subsequent RUSLE, difficult (Renard & Freimund, 1994). In the absence of detailed rainstorm data, mean annual and monthly rainfall have been shown to estimate the rainfall erosivity index (Arnoldus, 1980). This method uses a modified Fournier index such that
∑ = (3)
where F is the modified Fournier index value, pi is average monthly precipitation, and P is average annual precipitation (Arnoldus, 1980). The modified Fournier index has demonstrated strong linear correlation with mean annual rainfall for different European regions, even where precipitation varies seasonally (Kouli et al., 2009). Ferro, Giordano and Iovino (1999) replace the erosivity index with the modified Fournier index to propose an FF index such that
8
N 12 2 1 pi, j FF (4) N j 1 i1 Pj
where pij is the monthly rainfall depth of the year j and P is the rainfall total for the same year and N is the number of years over which data were collected. Arnoldus (1980) states, however, that the modified Fournier index must be calibrated between climatic regions. A theoretically derived constant is employed as a coefficient of proportionality (Ferro et al., 1999). Ferro et al. (1999) established a relationship between the FF index and rainfall erosivity for Sicily, Italy such that
1.56 R 0.6120 FF (5)
where FF incorporates the modified Fournier index. Kouli et al. (2009) adopted the mean R values of the Italian and a Moroccan relationship to estimate rainfall erosivity for Northwestern Crete, Greece. Ferro et al. (1999) have demonstrated that hydrologically homogenous regions can be characterized by the same R to FF relationship.
Interpolation, using various techniques, is a common application used to model spatial heterogeneity of natural data using discrete point values derived from data obtained at sampled locations (Bevan & Conolly, 2009). This application has been used successfully to model rainfall erosivity, R, in the Mediterranean, providing a continuous dataset, in raster form, that can be used for visual and computational analysis (Kouli et al., 2009; Terranova, Coscarelli & Iaquinta, 2009). Inverse distance weighting (IDW) uses a linearly weighted function to derive the contribution of neighbouring measured values to be assigned to unmeasured locations (Bevan & Conolly, 2009). This method of interpolation assumes that the influence of a variable decreases with distance. Alternatively, ordinary kriging uses autocorrelation to examine the spatial structure of measured point values for trends and directional influences (Ferreira & Panagopoulos, 2014). Ordinary kriging is considered to be more sensitive to input data than IDW, which uses pre- defined arbitrary distances. Consequently, ordinary kriging is used extensively to describe spatially heterogeneous natural phenomena, producing accurate results (Bevan & Conolly, 2009). Cross validation compares residuals between modelled values and actual point values enabling selection of the most appropriate prediction model for spatial variations in R (Ferreira & Panagopoulos, 2014).
9
2.1.2.2. The Soil Erodibility Factor
The soil erodibility factor, K, is an empirical measure of soil erodibility determined by soil texture, organic matter, structure and permeability (Kouli et al., 2009). The K factor reflects the ability for soil to separate during rainfall splash or overland flow (Demirci & Karaburun, 2012). Wischmeier and Smith (1978) uses a nomograph of five soil parameters: percent modified silt (0.002-0.1 mm), percent modified sand (0.1-2.0 mm), percent organic matter (OM), soil structure code, s, and profile permeability class, p, to determine soil erodibility. From the nomograph, Wischmeier and Smith (1978) develop an algebraic approximation for calculating soil erodibility when the silt fraction is less than 70% in representing the K factor such that
K 2.1 104 12 OMM1.14 3.25s 2 2.5p 3/100 (6) where M is the product of primary particle size fractions: (% modified silt)(% silt + % sand). K is expressed as ton per acre per erosion unit index. For conversion to SI units, K is divided by 7.59 (Renard, Foster, Weesies, McCool & Yoder, 1997). Wischmeier and Smith (1978) further develop a K factor from the geometric mean particle diameter of 225 global soil classes with less than 10% rock fragments (D > 2 mm) by weight. Soils with more than 10% rock fragments were excluded since mechanical soil analysis generally excludes rock fragments larger than 2 mm. This relationship calculates soil erodibility, K, such that