The Effect of Land Management Practices on Soil Quality Indicators in Crete

: The effects of four main practices tillage versus no-tillage, and intensive grazing versus extensive grazing, applied in characteristic agricultural and grazing lands of Crete Island were evaluated in situ using nine soil quality indicators. The following nine representative indicators of soil quality assessment were assessed using the rapid visual assessment methodology adopted at European level in the context of the EU research project iSQAPER: susceptibility to water and wind erosion, surface ponding (under cropping), formation of tillage pan, soil color, soil porosity, soil structure, susceptibility to slaking, inﬁltration rate, and biodiversity status. These indicators were measured in 48 agricultural ﬁeld-plots to adequately represent the four above-mentioned practices and the different types of geomorphological patterns existing in the area. Additionally, 38 agricultural ﬁelds were sampled in the topsoil to assess cultivation practices (tillage, no-tillage) on soil organic carbon, cation exchange capacity, exchangeable potassium, available phosphorous, and soil aggregate stability. Based on the indicators rating methodology, the appropriate statistical tests were applied and the soils under different managements were characterized in terms of their potential quality and their general agricultural value. The obtained data showed that in agricultural areas, signiﬁcant differences were detected between tillage and no-tillage management practices for the indicators of soil structure and consistency and inﬁltration rate. In grazing land, signiﬁcant differences were found for the soil quality indicators of susceptibility to erosion and inﬁltration rate for the corresponding practices of intensive and extensive grazing. Organic carbon content, exchangeable potassium content and aggregate stability were greatly affected in tillage versus no-tillage management practices.


Introduction
The concept of soil quality was first introduced by Mausel [1] as totally connected with the yields of the most profitable crops of the local region under high level of land management. Subsequently, Doran and Parkin [2] included biomass production as a soil quality function and linked the quality of soils with their inherent and difficult to change properties, inferred from characteristics related to those suggested by soil taxonomy. Larson and Pierce [3], considered the connection of soil quality with productivity a very narrow form that cannot embrace soil's contribution to environmental quality. Several studies carried out in early nineties [4][5][6], have supported that environmental quality as well as the human and animal health are the main three broader components that synthesize the concept of soil quality. The National Research Council of United States [7] regarded soil quality as the soil's capacity to enhance plant growth without aggravating water and air pollution. Warkentin [8], combining the previous approaches, considered soil quality as the basic concept for sustainable agriculture while Karlen et al. [9] focused on the importance of the interaction between several soil attributes and land management garded as the most representative physical indicators to soil quality under Mediterranean conditions. Rabot et al. [54] evaluated the potential of observable soil structural attributes in soil functioning assessment.
In accordance with the official definition of Sustainable Soil Management (SSM) [55] adopted by EU [56], a recommended set of indicators for the assessment of SSM has been proposed by the Food and Agriculture Organization of the United Nations [57]. Recognizing this need, the present work aims at providing a set of soil indicators for an easy and rapid assessment of the quality of agricultural soils, evaluating at the same time two different pairs of land management practices. In the framework of the EU research project iSQAPER (European Union's Horizon 2020 Programme for research & innovation under grant agreement no 635750) the following nine representative soil quality indicators were visually assessed on agricultural soils and grazing land: susceptibility to water and wind erosion, surface ponding-under cropping, formation of tillage pan, soil color, soil porosity, soil structure, susceptibility to slaking, infiltration rate, and biodiversity status. The study was carried out under four different land management practices, namely tillage versus no-tillage, and intensive grazing versus extensive grazing which are dominant in the island of Crete (Greece). Based on the indicator ratings, the soils were characterized in terms of their potential quality to estimate their general agricultural value. Additionally, a set of analytical soil data of representative agricultural soils under tillage versus no-tillage treatments were used to test the sensitivity of soil indicators in evaluating soil quality.

Study Area
The broad study area of Crete Island is bordered to the north by the Sea of Crete and to the south by the Libyan Sea. It covers an area of 8336 km 2 . Cultivated land covers 42.0% of the total island, while dryland, used mainly as pasture, is the next important ecosystem, covering 39.3% of the total area. The island is characterized as a sloping land with slopes greater than 12% in 79.5% of the area, and only 6.9% comprises lowlands with slope less than 6%. The annual rainfall ranges from 262 mm yr −1 to 978 mm yr −1 in the low areas along the coast. In the mountainous areas rainfall ranges from 551 mm yr −1 to 1836 mm yr −1 . Crete lies between the isotherms 18.5 • C to 19 • C with an annual amplitude of 14 • C to 15 • C.
The soils of the study fields present a variety of chemical and physical properties. They have been formed in a variety of parent materials such as limestone, shale, marl, conglomerates, flysch, and alluvial deposits. The soils are mainly characterized as moderately fine-textured (clay loam, silty clay loam), followed by fine (clay, silty clay), medium (loam), and coarse (sandy loam). The study fields are mainly characterized from slightly sloping (slope 2-6%) to steep (slope greater than 60%). The soils are characterized as shallow (depth 15-30 cm) to very deep (depth > 150 cm). Organic matter content in the surface A-horizon ranges from 0.5% to 8.6% depending on the land use, parent material, and land management practice. The soils present various degrees of development. Cambisols are widely extended in the study fields followed by Calcicols, Leptosols, Luvisols and Regosols. Cambisols are mainly formed on a variety of parent materials such as marl, conglomerates, flysch, and shale. Regosols are formed mainly on conglomerates parent material. Leptosols are located in mountainous areas, used mainly as grazing land. Fluvisols and Luvisols are mainly located in plain areas covering few cases of the study fields. Climatic and soil characteristics accompanied by the EU policies on subsidizing crops in the last decades have greatly favored the extensive growth of olive and vine plantations in the area providing higher farmer's income. Orange plantations and vegetables growing in greenhouses have been mainly expanded in the lowlands of the island. The intensification of agriculture resulted in accelerated rates of soil erosion in the hilly areas of the island. The following land management practices are mainly applied in the area of Crete: (a) Intensive cultivation of the land accompanied by disk harrowing and application of fertilizers and pesticides; (b) Integrated land management in olive groves by applying measures for environmental protection; (c) Organic farming in vegetables and olive groves; (d) Overgrazing in pasture land.

Collection of Quantitative and Qualitative Soil Data
A representative part of Crete island in terms of applied agricultural practices and main land uses was chosen for testing, sampling, and evaluating the soil quality indicators affected by the applied land management practices. Specifically, 48 agricultural fields (plots) (Figure 1) subjected to 4 different agricultural management practices were selected for rapid assessment of soil quality indicators. The selected study fields are located in the municipalities of Phaistos, Gortyna, Archanon, Asterousia, Heraklion, and Hersonissos of the Heraclion region. Furthermore, other 38 representative agricultural fields under tillage versus no-tillage management practices were sampled in the topsoil and used to correlate soil quality with cultivation practices. These fields are located in the western part of Crete in the municipalities of Kissamos and Platanias of the Chania region ( Figure 1).  The selected 48 field plots were chosen to be representative of the two main agricultural land uses, namely cropland and pasture land, subjected to the following four principal management practices (a) tillage, (b) no-tillage or minimum tillage in agricultural land, (c) intensive grazing, and (d) extensive grazing applied in pasture land. Furthermore, the 48 field plots have been selected to be taxonomically identical by pair, including the same topographic and geomorphic conditions, and soil taxonomical properties. Fifteen pairs of the 48 plots (30 fields) correspond to cropland under tillage and no-tillage management practices in each pair of field plots. The remaining 9 pairs of the 48 plots (18 fields) correspond to grazing land including intensive and extensive grazing management practices in each pair of field plots.
Soil quality in each field plot was evaluated using the nine representative soil quality indicators described and visually scored in the field following the methodology developed by the EU research project iSQAPER (https://www.isqaper-is.eu/soil-quality/visual-soilassessment, accessed on 3 November 2020). The scores of each indicator ranged among three evaluation levels: "0 for bad condition", "1 for moderate condition", and "2 for good condition" through a set of representative and characteristic photos to each soil indicator state. The scores of the nine indicators assessed in each pair of field plots were compared and a composite soil quality index was calculated by summing all positive and negative effects of each soil variable [58].
The selected soil quality indicators were the following: susceptibility to water and wind erosion, surface ponding (under cropping), tillage pan, soil color, soil porosity, soil structure and consistency, slaking test (aggregate stability), infiltration rate, and biodiversity (earthworms density) The scoring of the susceptibility to wind erosion was based on visual observations and whether the soil material was stable in the field or it was dispersed by the wind in the wider area considering the size of the dust plumes emanated from the cultivators on windy days. Water erosion was estimated by observing the differences in topsoil's depths between the crest and the bottom of the slope and by recording the amount of sheet and rill erosion, as well as the sedimentation into surrounding drains and streams. The period elapsed (in days) until rain water logging spots disappeared from the soil surface determined the degree of surface ponding. Tillage pan was estimated visually by the presence and degree of soil structure in combination with the compaction state of the subsoil. The visual scoring of soil color was based on the differentiation of the value of soil color, subjecting the higher score to the darker colored topsoil. Soil porosity was estimated visually by recording the presence of macropores using a set of representatives to soil structure development (photographs). The scoring of soil structure and consistency was achieved according to the distribution and the size of aggregates as well as to the degree of clodding. Soil susceptibility to slaking was tested and scored by a slaking test, assessing the ability of soil to maintain its structure when it is subjected rapid wetting placed into a container of water. The indicator infiltration rate was measured by using the experimental method developed in Bern University [59]. The earthworm density was determined by hand-counting of the population of earthworms that appeared in a soil volume of 20 cm 3 . In each study field spot, measurement of each indicator was carried out by two or three independent experts recording two or three scores for each indicator. Finally, the recorded scores for each indicator were averaged in one value.
For further characterization of the study site, laboratory analytical data, representative for tillage versus no-tillage cultivation practices, were used supplementarily, corresponding to soil samples that had been taken from the A-horizon of the above-mentioned, non taxonomically related, 38 sites ( Figure 1). The selected soil samples had been air-dried and sieved to 2 mm mesh for the subsequent laboratory analyses. The following chemical and physical soil properties have been determined: total organic carbon (TOC) determined by the Walkley-Black wet oxidation method [60]; exchangeable potassium (Ex.k) extracted by NH 4 OAc (1M, NH 4 OAc, pH 7) [61]; cation exchange capacity (CEC) determined at pH 7 by ammonium acetate [62]; aggregate stability determined by the Mean Weight Diameter (MWD) throughout wetting and sieving processes [63][64][65][66][67][68] available inorganic ortho-phosphate (Av.P) (PO 4 -P) by extraction with 0.5 N sodium bicarbonate solution adjusted to pH 8.5 [69].

Statistical Analysis
According to the used methodology, the obtained data on soil quality indicators belong to the hierarchical scale (qualitative data) of data measurement levels, given that in this subjective scale of scoring there is no objective (quantitative) distance between any two of its points. According to the nature of the data and the purpose of the present work, which was to investigate the impact of the different land management practices on soil quality indicators, the most appropriate statistical test for analyzing the qualitative data is the Wilcoxon Matched Pairs Signed Ranks Test [70]. This test is a non-parametric equivalent of the paired t-test and is used for ordinal scale data (qualitative) that they cannot support the assumptions of normal distribution. The Wilcoxon Matched Pairs Signed Ranks Test treats qualitative data as if they had arisen from an experimental design which was developed to study specific properties of a taxonomically identical couple (pair) of biological material before and after the application of a particular treatment and evaluates group differences of the paired observations with attention to medians. The null hypothesis of the Wilcoxon Matched Pairs Signed Ranks Test is that there are no statistically significant differences in the examined properties between the component units of a pair undergoing different treatments [70]. The obtained data were appropriately ordered per agricultural management practice and per soil quality index to apply the Wilcoxon Matched Pairs Signed Ranks statistical test using the open-source integrated development environment for statistical computing and graphics RStudio v. 1.2.5042 (RStudio Team (2020). RStudio: Integrated Development for R. RStudio, PBC, Boston, MA URL http://www.rstudio.com/).
In contrast to the qualitative data of the rapid assessment methodology for soil quality indicators, the laboratory analytical data of the five soil parameters (MWD, SOM, Ex.K, CEC, Av.P) are quantitative and belong to the interval scale of data measurement levels.
Thirty-eight random fields were selected with unique criterion the application of specific cultivation practices (tillage versus no-tillage) in the various but representative geomorphological and land use types of the agricultural soils in Crete. Therefore, by the nature of the specific plots' selection, the fields are independent of each other or they are not taxonomically related. To assess the effect of the aforementioned cultivation practices on soil quality, through laboratory analyses measurements of parameters, the two samples unpaired t-test was used. The specific statistical test compares the means of two independent or unrelated groups to determine significant differences [71]. The null hypothesis of the two samples unpaired t-test is that there are no statistically significant differences in the examined soil parameters between the two groups undergoing different treatments [72]. The obtained data were appropriately ordered per agricultural management practice and per soil parameter to apply the two samples unpaired t-test using the open-source integrated development environment for statistical computing and graphics RStudio v. 1.2.5042 (RStudio Team (2020). RStudio: Integrated Development for R. RStudio, PBC, Boston, MA URL http://www.rstudio.com/).

Susceptibility to Wind and Water Erosion
In agricultural land, independently of the applied land management practices, the indicator susceptibility to wind and water erosion was rated mainly as moderate (15 cases) followed by low or no susceptibility (9 cases) and high susceptibility (6 cases) ( Table 1). If the data are distinguished in terms of the applied land management, no tillage management presented mainly low or no susceptibility to wind and water erosion (8 cases), followed by moderate susceptibility class (6 cases) and high susceptibility class (1 case). However, if the tillage management practice is considered, the class of moderate susceptibility to wind and water erosion was recorded in nine cases, followed by high susceptibility in five cases, and only one case characterized with low or no susceptibility to erosion. Table 1. Soil quality indicator and soil quality scores of the nine indicators per pair of management practices in the examined fields.  SWWE  SP1  TP  SC  SP2  SSC  ST  B  IR  CSQI   25-NT17M  1  2  1  2  1  1  2  1  In grazing land, independently of management practice, the main class of wind and the water erosion was moderate (12 cases), followed by both high and low or no susceptibility with the same scores (3 cases). If the obtained data are classified according to land management practices, extensive grazing ranged from moderate susceptibility (6 cases) to low or no susceptibility (3 cases), while intensive grazing ranged from moderate susceptibility (6 cases) to high susceptibility (3 cases).

Surface Water Ponding
Surface water ponding was rated as good (cases 29 out of 30) (no ponding) in agricultural land due to the inclination of the land in the majority of the study cases, independent of tillage or no tillage management practice. Rain water in sloping land is temporarily ponded in small spots until surface water runoff is initiated. Similarly, in grazing land, no water ponding was recorded in all study field plots (18 of 18 cases) regardless of extensive or intensive grazing.

Tillage Pan
The majority of the agricultural study field plots (24 cases) have formed a tillage pan rating with a moderate degree of formation, while in some field plots (6 cases), no tillage pan was identified. The rating of the index in relation to the applied management practice showed that no tillage management mainly led to a moderate formation of tillage pan (11 cases), while in some fields (4 cases), no tillage pan formation was observed. A slightly worse rating of the index was recorded under the tillage management practice, with thirteen cases classified as having moderate tillage pan formation and only in two cases of no tillage management a pan was observed.
In grazing land, independently of management practice, no tillage pan was observed in the majority of the cases (13), or in some cases (5), a moderate degree of tillage pan formation was recorded. If the obtained data are distinguished on the basis of management practice, in both cases extensive or intensive grazing, the absence of tillage pan formation was observed in seven and six cases, respectively. In a few cases, the formation of a tillage pan (2 cases for extensive grazing and 3 cases for intensive grazing) was observed.

Soil Color
In the majority of the study field sites in agricultural areas (9 pairs), soil colour was rated by one degree higher in fields under no tillage practices compared to fields subjected to tillage operations. Under the no tillage practice, a moderate dark colour was mainly recorded (9 cases), while the rest of the study fields were classified as having a darker colour (6 cases). Fields under the tillage practice presented mainly a moderate dark soil colour (9 cases) with only two fields having a darker soil colour, while the remaining fields (4 cases) had significantly lighter soil colour.
Similarly, in most of the study fields in grazing land, the colour of the soils under extensive grazing have been characterized as being darker by one degree in six pairs of the field plots. In eight fields of the extensively grazed area, the topsoil colour was dark, and only in one field was the soil colour lighter, belonging to the moderate class. In the opposite, fields under intensive grazing showed only three cases of dark soil colour and six cases with a lighter soil colour rated in the moderate class.

Soil Porosity
The obtained data on soil porosity showed that the dominant classes in agricultural land were moderate and good recorded in thirteen and eleven cases, respectively, independent of the applied management practices. Concerning the management practice of no tillage, in nine fields, the good class was identified, characterized by many macropores between and within aggregates. In addition, under the same management practice, the number of macropores have been declined significantly in six field plots, classified in the moderate class. In contrast, under the tillage management practice, the situation of soil porosity was worse, with only two fields being classified as good, while in seven fields soil porosity was characterized as moderate and as poor in six cases (absence of soil macropores).
In general, the same trend was observed in grazing areas with moderate and good ratings in ten and eight cases, respectively. Specifically, in extensively grazed areas, seven sites were classified as having good class of soil porosity, while only two fields showed a moderate class. In contrast, intensive grazing presented eight fields of moderate soil porosity conditions and only one field was under a good condition of soil porosity.

Soil Structure and Consistency
In agricultural land, thirteen fields were characterized by good soil structural state, eleven presented a moderate class, and six fields presented with poor soil structure. The soil structure in soils under no tillage practice was rated one class higher in the majority of the study field sites (cases 10 out of 15) compared to fields under tillage practice. The poor soil structure class was identified in six cases under tillage management practice.
The soil structure condition in grazing land was characterized in nine fields as good, in eight fields as moderate, and only in one case as poor. In all cases except one, the soil structure was rated one degree higher in extensive grazing compared to intensive grazing.

Susceptibility to Slaking
The stability of soil aggregates assessed by slaking in agricultural land was characterized mainly as good in twenty-three cases, followed by the moderate class in five cases, and as poor in two cases. Soils under no tillage practice were characterized mainly by good soil aggregate stability in twelve cases and only in three cases as moderate. Almost the same pattern was observed in fields under tillage management practices, with the good class identified in eleven cases, while two cases were moderate, and two cases with poor structure stability were assessed. Only in three pairs of fields, the tillage practice had a lower rating by one degree compared to no tillage. In grazing land, the soil susceptibility to slaking was rated as good in most of the cases (8 out of 9, both in extensively and intensively grazing).

Infiltration Rate
The infiltration rate in agricultural land was mainly classified as moderate (17 cases) with several field plots (12 cases) of low ratings and only one case of a high rating. Regarding the indicator status according to the cultivation practices, no tillage management rated mostly in the moderate class (11 cases), followed by the two other classes of low (3 cases) and high (1 case). Soils under tillage management practice were mainly classified with a low score for infiltration rate (9 cases), while the class of the moderate was observed in 6 cases.
The infiltration rate in the soils of grazing land was mainly rated as moderate (12 cases), while in some fields, the infiltration rate was characterized as high (4 cases) or as low (2 cases) status. The studied grazing management practices showed different responses of the indicator. Extensive grazing ranged from moderate (5 cases) to high infiltration rate (4 cases), while under intensive grazing, the moderate infiltration rate was the dominant class (7 cases), followed by low infiltration rate recorded only in two cases.

Biodiversity
Earthworms have been recorded only in cropland in twenty-one cases. Under no tillage practice, the biodiversity was characterized as good in two cases, while the moderate class was identified in ten cases. Under tillage management practices, the moderate class of biodiversity was characterized as dominant in nine cases. In the opposite, the biodiversity in grazing land was characterized mainly as poor (seven cases in extensive grazing and eight cases in intensive grazing), with only three fields of moderate condition (2 cases in extensive grazing and 1 case in intensive grazing).

Statistical Analysis Results
The Wilcoxon matched pairs signed ranks test was applied for the statistical analysis of the obtained forty-eight (48) qualitative data of the nine (9) representative soil indicators. Since the Wilcoxon test is mainly used to cross check the differences in medians of paired observations, the distribution of the differences of the observations must be symmetrical around the median. According to the visual symmetry checking of distributions, as far as tillage and no-tillage treatments of the agricultural areas were concerned, it was inferred that in a strict manner, paired differences had come from a continuous symmetrical distribution only for the indicators soil porosity, soil structure, and infiltration rate. Significant differences were detected between the control fields (tillage practice) and the fields of no tillage application for the indicators soil porosity, soil structure and consistency, and infiltration rate ( Table 2). The application of Wilcoxon test in grazing land for intensive and extensive grazing practices detected, in a strict manner, paired differences for a continuous symmetrical distribution only for the indicators susceptibility to erosion, surface ponding, biodiversity, and infiltration rate indices. Furthermore, the indicators of susceptibility to wind and water erosion, and infiltration rate presented significant differences between the paired field plots with intensive and extensive grazing practices ( Table 3). The indicators of surface ponding and biodiversity were not statistically processed due to low variation of their values (Table 3). Based on the nine soil quality indicators, the composite soil quality index was determined. The obtained data showed that land management practice had a great impact on the general soil quality index (Table 1). In the majority of the studied field plots, no tillage has positively affected soil quality. The composite soil quality indices were greater in thirteen pairs out of fifteen pair plots in cropland. Furthermore, extensive grazing had positively affected soil quality in the majority of the study field plots (8 cases out of 9) compared to intensive grazing. For the needs of the present publication, the obtained values of the composite soil quality index were grouped in the following three classes: poor < 0.8, moderate 0.8-1.4, and good > 1.4. According to this classification, the soil quality of the study field plots in cropland can be mainly characterized as moderate in eighteen cases, followed by good in eight cases, and poor in four cases. Similarly, in the grazing land, soil quality is mainly evaluated as moderate in twelve cases, followed by good quality in six cases.
For the application of the unpaired two samples t-test, it was prerequisite to check the normality of distributions and the equality of variances [71] of the five analyzed (Table 4) soil parameters (Mean Weight Diameter-MWD, Soil Organic Matter-SOM, Cation Exchange Capacity-CEC, Exchangeable Potassium-Ex.K, and Available Phosphorus-Av.P). For that purpose, the Shapiro-Wilk test was applied for checking the normality of the data [72,73] and the F-test for assessing the equality of variances [74]. The analysis of the data by applying Shapiro-Wilk and F tests showed that only the parameter of available phosphorus was rejected (Tables 5 and 6). Therefore, the data populations for the following four parameters MWD, SOM, Ex K, and CEC are normally distributed, while the variances between the land practices of tillage and no tillage are equal for the same parameters. Therefore, it was plausible to apply the two samples unpaired t-test only for the previously mentioned four soil parameters. Finally, the application of the two samples unpaired t-test showed that the means of two land management practices (tillage versus no tillage) were statistically significantly different for the parameters MWD, Ex.K, and SOM, while the difference was insignificant for the parameter CEC (Table 7).

Discussion
The obtained data have shown that the applied land management practices in agricultural and grazing land had in various degrees affected the analyzed nine indicators. In the agricultural land, the change of agricultural practice from tillage to no tillage had mainly positively affected the studied indicators and the soil quality index. The analysis of the scores showed slight differences for the indicators formation of tillage pan, susceptibility to slaking, and biodiversity or medium for the indicators susceptibility to water and wind erosion, infiltration rate, soil color, soil porosity, soil structure, and consistency.
Tillage pan is mainly related to cultivation. The use of heavy machineries for cultivation under relatively wet soil conditions favored the formation of a compacted subsurface soil layer underlying the plough soil layer. The recorded small differences in the rating of tillage pan index between the tillage and no-tillage practices may be attributed to the tillage management practice applied almost in the whole land in the previous decades. The study area is mainly cultivated with olive trees. For various reasons, farmers have realized that it is more profitable (economically and environmentally) to keep olive plantations under a no tillage practice. However, tillage pans were already formed but not further aggravated. The same explanation could justify the small differences observed in the rating of the soil susceptibility to slaking (slaking test) between no-tillage and tillage management in agricultural areas.
The medium differences in the rating of susceptibility to water and wind erosion, between tillage and no tillage practices, can be mainly attributed to the topographic conditions and the existing land use type. The study area is mainly sloping with various degrees of inclination, cultivated mostly with olive trees. Long term soil erosion experiments carried out in the study area have shown that water erosion is greatly limited, especially under no tillage management practice in olive groves, supporting the findings of the present study [75]. Farmers used to till the soil during February or March, creating favorable conditions for water erosion, especially in the intermediate unprotected soil zone by the stem of successive trees.
Land management practices have affected the indicator infiltration rate. Medium differences of rating between tillage and no-tillage management are mainly attributed to poorer soil porosity and soil surface crusting under tillage, formed after intensive rainfall events. Furthermore, soil porosity index has been related to the soil structure development index. Soils with good macroporosity usually had good structural development. In addition, soil porosity and, therefore, soil structure and consistency were related to land management practices. No tilled soils had better porosity and soil structure than tilled soils. Statistically, only the soil structure and consistency and the infiltration rate indices were significantly different under tillage versus no-tillage management practices of agricultural areas.
Soil colour, and particularly soil darkness, is mainly related to the amount of organic matter content. Organic matter content is affected by soil characteristics (soil texture drainage conditions, etc.), climate (rainfall, air temperature), land use (annual, perennial), and land management practices. Most of the studied agricultural field sites are located on the same parent material (marls, conglomerates) under similar climatic conditions and the same land use type (olives). Therefore, land management practices had a big impact on soil colour. No tillage management practice favors both nil water erosion and the growth of understory annual vegetation, resulting in increased soil organic matter content.
Soil biodiversity in agricultural land, defined only by the number of earthworms, has been affected in the study area by various factors such as land use and management practice, soil characteristics, and climate. The climatic conditions in the study area of Crete are largely adverse for long periods of the year, particularly considering the distribution of rainfall. Due to the lack of significant amounts of rains from April until October, the soil in the upper layers remains dry, creating unfavorable conditions for earthworm growth. Therefore, the use of this indicator for assessing soil quality imposes limitations of measurement. However, statistically significant differences were found only for the indicators soil structure and consistency and infiltration rate between the management practices of tillage versus no-tillage.
The land management practices in cropland had a great impact on soil chemical and physical properties, especially in the soils' mean weight diameters and organic matter contents ( Table 4). As Figure 2 shows, a linear relation was found between these two parameters, including no-tillage versus tillage management. No tillage land management has favored higher amounts of organic matter content and better aggregate stability compared to soils in which tillage operations were carried out. However, statistically significant differences were found for the parameters organic matter content, mean weight diameter and exchangeable potassium content between the management practices of tillage versus no-tillage. Visually assessed differences in indicators ratings of grazing land were ranged from slight, for the indicators tillage pan, susceptibility to slaking, biodiversity, to moderate for the indicators susceptibility to water and wind erosion, infiltration rate, soil color, soil porosity, soil structure and consistency. Extensive grazing has mainly positively affected soil quality compared to intensive grazing, except for the indicator susceptibility to slaking for scores in moderate and low classes. The studied grazing land is mainly a sloping land covered by shrubs or cultivated with annual plants used as animal fodder. The small differences of scores in tillage pan formation index and susceptibility to slaking between extensive and intensive grazing management can be attributed to the periodic and rotational cultivation of the grazing land in combination with the pressures exerted by animal trampling. Regarding the indicator of biodiversity, the justification of the recorded small differences between extensive and intensive grazing is the same as those mentioned for agricultural land. The moderate differences in the indicator susceptibility to water and wind erosion can be attributed to the favorable conditions for soil erosion created first for few weeks in November during sowing the soil and second in early spring after intensive grazing of the corresponding fields leaving the soil partially covered by adequate vegetation.
The following indicators had big differences between extensive versus intensive grazing: infiltration rate, soil colour, soil porosity, and soil aggregate stability. Specifically, the existing higher amount of organic matter content under extensive grazing created favorable conditions for the development of a good soil structure and of a darker surface soil colour. Consequently, good soil structure in the extensively grazed areas enhanced infiltration rate and soil porosity. However, statistically significant differences were found only for the indicators susceptibility to wind and water erosion and infiltration rate between the management practices of extensive versus intensive grazing.

Conclusions
Soil quality can be considered as a composite indicator affected by several single soil indicators. In the present work, soil quality was assessed through nine qualitative and visually assessed indicators in cropland and pastures and supplementary through five quantitative indicators, only for agricultural areas. Among the most important indicators affecting soil quality were soil erosion susceptibility, soil structure and consistency, infiltration rate, aggregate stability (mean weight diameter), and soil organic matter content. However, the analyzed indicators have been affected by the land management practices of no tillage versus tillage in cropland and extensive grazing versus intensive grazing in grazing land. No tillage and extensive grazing presented a great positive effect on soil quality indicators. Tillage operations and intensive grazing tend to promote accelerated soil erosion, to enhance loss in soil organic matter content, to reduce soil aggregate stability and infiltration rate, and to deteriorate the biodiversity of soils. Qualitative soil quality indicators can be a valuable tool at the disposal of farmers and scientists for improving soil quality. The methodology of visual assessed indicators in the field is proposed to be used for specific purposes of monitoring the effects of the applied agricultural land management practices. An important advantage of using the suggested rapid methodology, of soil's general agricultural value assessment through the visually scoring of soil quality indicators, is that the required data are immediately and easily obtained in the field in contrast to indicators requiring expensive and time-consuming laboratory analyses or processing.