The study was conducted in April 2010 on Gumara-Maksegnit watershed which has an area of 53.7 km² and located in the northwestern Amhara region, Ethiopia between 12˚24' and 12˚31' North and between 37˚33' and 37˚37' East (Figure 1). In the study watershed, the problem of soil fertility, water shortage, costly agricultural inputs and erosion are among a number of serious challenges, which makes small-scale rainfed agriculture a risky enterprise.

The mean annual rainfall in the watershed was about 1170 mm at which more than 90% of the rainfall occurs during the three month period, June to August. Meanwhile, the average daily maximum and minimum temperatures were 28.5˚C and 13.6˚C, respectively. Diverse topography is a typical characteristic of the study watershed. The elevation of the watershed varies from 1923 m at the outlet to 2851 m in the north and the mean elevation is 2237.54 m above sea level, which exhibit vast topographical variations over the research site (Figure 2). Meanwhile, the majority of the study watershed (90% of the study area) is mountainous and consists of dissected terrain with steep slopes, and the remaining 10% has an undulated topography with gentle slopes. The

slope of the study watershed ranges from nearly flat (less than 2%) to exceptionally steep (greater than 70%) in the northern part of the watershed and the mean watershed slope is 22.06%. The study watershed was mainly covered by agricultural land (63.5%) followed by forest (24.3%) and grassland 12.2%. The major crops grown in the agricultural land included sorghum, teff (EragrostisTeff), fava bean, lentil, wheat, chickpea, linseed, fenugreek, and barley. Eragrostis Teff and sorghum were the main staple crops, whereas chickpea was grown at residual moisture in the lower regions of the watershed where clay soil textural classes were dominant and this crop cannot grow at higher altitudes.

The study area, which is part of the northern Highlands of Ethiopia, belongs to the Trapp series of Tertiary volcanic eruptions [20] . The physiographic structure of the study area is a result of its geological origin. The Highlands and mountain ridges are created by uplifting of landmasses and the mountain peaks are the results of lava emissions [20] . The geomorphic processes which took place for a long period of time had significantly reshaped the area. The soil types are predominately Cambisol and Leptosol which are found in the upper and central part of the watershed, whereas Vertisol is found in the lower catchment near the main outlet. Denkez mountain ridge is the peak part of the watershed and it is a line of demarcation that separates the two main Ethiopian basin systems (Nile and Tekezi basins).

The study watershed was divided into a 500 m by 500 m square grid and approximately, at the center of each grid soil samples were collected within 184 locations over a 53.7 km² watershed using auger and core cylinder equipment. Garmin explorer GPS accuracy: (±3 m) was used for locating soil sampling points and measuring elevation in the field. However, some sampling points were shifted a bit from the center of the grid due to difficult topographic conditions. Soil samples in the study watershed were collected during the latter part of the dry season.

At each grid center, soil samples of about 2 kg were removed from the surface soil horizon, which might range from 0 - 25 cm based on the surface soil depth, with the best available tool (bucket auger) for chemical and physical analysis. Undisturbed core samples from 10 - 15 cm depth were taken to determine bulk density (ρd) following the procedure of [21] . The pH value of the soil was measured with a pH meter in the supernatant suspension of 1:2.5 ratio (sample to water mixture)measured using glass electrode pH meter as described by Peech [22] . Soil organic carbon was determined by wet oxidation method as described by De Vos et al. [23] . Available Phosphorus (AP) was extracted using sodium bicarbonate solution at pH 8.5 following the procedure described by Olsen [24] . Total nitrogen (TN) was determined using the Kjedahl method as described by Jackson and Barak [25] .

The observed data for the parameters SOC, TN, AP, pH, bulk density, and elevation were subjected to classical statistical analyses. The descriptive statistical summary for the measured soil physical and chemical properties of the study watershed is presented in Table 1. Over the sampled watershed, mean values of SOC, pH, and TN were 1.49% ± 0.07%, 6.75% ± 0.03, and 0.21% ± 0.01%, respectively. The soil analyses also showed that the value of bulk density ranged mainly from 0.83 (g・cm⁻³) to 1.70 (g・cm⁻³) with the mean value of 1.28 ± 0.01 (g・cm⁻³), while the mean value of AP was 14.20 ± 1.09 (ppm) (Table 1). SOC, AP, TN and elevation had skewness coefficient greater than 1, hence the natural logarithm was applies to stabilize the variance.

As a standardized measure of the variance, the coefficient of variation (CV) was used to describe the shape of a frequency distribution of the observation [28] . Based on the CV values of the observed soil attributes, AP, SOC, TN were the most variable soil properties, with CV greater than 37%. Observed bulk density (ρd) was moderately variable, with CV between 14% - 27%, whereas pH was the least variable (CV less than 7%). The descriptive analyses also revealed that AP had the largest variation (CV equals to 104%). Similarly, [2] documented that AP showed the highest, whereas pH the least CV. In addition, previous researches also depicted low coefficients of variation for pH compared to the other soil properties [8] [9] [11] [12] .

Considering the soil attributes/elevation relationships, one might expect that the spatial pattern and magnitude of observed soil attributes should be strongly related to elevation of the watershed, since it dictates the flow direction and accumulation, drainage network, and slope gradient which eventually contributes to nutrient leaching, soil erosion and sediment deposition. Investigating the relationship of SOC with elevation is crucial as SOC is a major part of the terrestrial carbon reservoir which plays an important role in the global carbon (C) cycle [31] . In the study watershed, there was a positive linkage between elevation and SOC (Table 2). These findings were in line with that of Kidanemariam et al. [32] who reported a positive correlation between SOC and altitude. Soils of the upper catchments had the highest SOC contents and various factors may affect the link between SOC and elevation. Lower temperature and increased soil moisture probably reduce decomposition rates and thus the accumulation of SOC at upper catchment of the study watershed. Similarly, research conducted by Kirschbaum [33] showed that lower soil temperature usually decreased litter decay rates and soil organic matter decomposition rates [34] . Bulk density (ρd) and pH indicated an inverse (p < 0.001) and TN a direct relationship (p < 0.001)

ρd―Bulk density, SOC―Soil Organic Carbon, TN―Total Nitrogen, AP―Available Phosphorous, sd―Standard deviation, se (mean)―Standard error of mean, CV―Coefficient of Variation.

with elevation. Meanwhile, SOC was significantly correlated with TN (p < 0.001) and negatively correlated with bulk density and pH (p < 0.001) (Table 2).

The resulting semivariograms exhibits the existence of diverse spatial dependence for the study watershed (Table 3). The value of relative nugget effect or nugget to sillratio (C₀/C₀ + C) displays the spatial autocorrelation [35] . According to [16] the ratio of nugget to sill (C₀/C₀ + C) of less than 0.25, 0.25 to 0.75, and greater than 0.75 can reflect strong, moderate and weak spatial autocorrelation, respectively. The value for relative nugget effect of bulk density, pH and elevation are less than 0.25, and depicts that it has strong spatial autocorrelation. The values for relative nugget effect of SOC and TN are between 0.25 and 0.75, and indicate that they have moderate spatial autocorrelation. With regard to AP the value for relative nugget effect is more than 0.75, and depicts that it has weak spatial autocorrelation. Generally, all measured soil properties in the study site varied in their spatial dependence and revealed a positive nugget effect with low sampling error (Table 3).

The resulting semivariograms on different directions including 0, 45, 90, 135 degrees for all soil attributes showed that effective range and sill of semivariograms were uniform and there was no clear anisotropy, therefore the spatial variability for the selected soil attributes were assumed to be identical in different directions (isotropic). The semivariograms of bulk density (ρd), pH, and elevation were best fitted to Gaussian model while that of SOC and TN were best fitted to exponential function (Table 3). Although the semivariogram of AP seemed that it was best described by the exponential function, the resulting spatial structure clearly showed a pure nugget effect. Meanwhile, the lag distances for the observed soil attributes were much longer than the sampling interval (Table 3). Thus, the sampling design used in the study watershed for the selected soil attributes was adequately revealed spatial distribution. As an illustration, the experimental semivariograms of bulk density, elevation, SOC, TN, AP and pH for the best fitted model are displayed in Figure 3(a) to Figure 4(d), respectively.

Asterisks indicate correlation coefficient is significant at ^*p < 0.05; ^**p < 0.01; ^***p < 0.001 (two-sided).

RSS―Residual Sum of Squares, r²―Coefficient of determination, RNE―Relative nugget effect.