Patchy Distribution Fields: Sampling Distance Unit of a Spiral Survey and Reconstruction Adequacy

A mathematical model is used to examine the effects of choosing various units of sampling distance of a spiral acoustic survey on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches (or gaps) of different shapes and spatial orientations, and an acoustic survey by a spiral of Archimedes along which a unit of sampling distance is set. For comparison, surveys are imitated by parallel transects along which the same unit of sampling distance is set. Adequacy of the reconstructed fields to those originally generated is evaluated by calculating their correlations (r). In the case of a spiral survey, the mathematical experiments conducted show that an immovable field can be reconstructed properly (r ² > 0.70) if the ratio of the units of sampling distance to the autocorrelation radius for the field averaged in various directions d/R _av < 1.0–1.5. Regarding immovable fields, the spiral surveys ensure, practically speaking, the same adequacy of the field reconstruction as do surveys by parallel transects with the same unit of sampling distance. In regard to movable fields, a comparison of the results of spiral surveys with those of surveys by parallel transects indicates that the former may ensure even higher adequacy of the field reconstruction than do the latter, provided that the units of sampling distances in these surveys are equal to each other.     

Patchy distribution fields: A spiral survey design and reconstruction adequacy

A mathematical model was used to examine the effects of a spiral survey design on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches (or gaps) of different shapes and spatial orientations, and acoustic surveys by a spiral of Archimedes; for comparison, surveys by parallel or zigzag transects are imitated. Adequacy of the reconstructed fields to those originally generated was evaluated by calculating their correlations (r). The mathematical experiments conducted showed that spiral surveys ensure, practically speaking, the same adequacy of field reconstruction (both in cases of immovable or movable patches) as do surveys by parallel or zigzag transects with greater sampling effort (overall path). In the case of a spiral survey, a patchy field can be reconstructed properly (r² = 0.70) if the overall survey path is not less than S/R_av = 20–30, where R_av is the autocorrelation radius averaged for various directions. Thus, the results obtained allow us to conclude that a spiral survey design is expedient in cases that the minimal duration of a survey is a decisive factor for its conduction and there is a priori information that no onshore-offshore gradients of fish density exist in a region under study.     

Patchy distribution fields: sampling distance unit of a zigzag survey and reconstruction adequacy

A mathematical model was used to examine the effects of choosingvarious units of sampling distance of a zigzag survey on the adequacy of reconstructing patchy distribution fields. The modelsimulates fish or plankton patches (or gaps) of different shapesand spatial orientations, and an acoustic survey by zigzag or parallel transects along which a unit of sampling distance is set. Adequacy of the reconstructed fields to those originally generated is evaluated by calculating their correlations (r). A priori information on the autocorrelation radii for the field in the directions of the survey (R_s) and perpendiculardirection (R_p) allows optimisation of the survey design andthe algorithm of data analysis. A field can be reconstructed properly (r² > 0.70) if the distance between transects D < (1.0–1.5)R_s and the unit of sampling distance d < (1.0–1.5)R_p. A posteriori determination of patchorientation allows reconstruction of the best field attainable onthe basis of the survey data. In cases of field movement, if thedimension of patches in the direction of movement exceeds that ofa surveyed area, a survey in the opposite direction gives bestresults; in contrast, if the dimension of moving patches is smaller than that of a surveyed area, it is reasonable to carryout a survey in the same direction. The criterion remains validwhen a survey is carried out by zigzag transects and a unit ofsampling distance is set along them. The results obtained indicate that, for a fixed transect spacing and a given number of sampling points on each full transect, zigzag pattern allowsless adequate reconstruction of an original distribution field (in cases of both immovable and movable fields) than corresponding parallel pattern.     

Patchy Distribution Fields: An Interleaved Survey Design and Reconstruction Adequacy

A mathematical model was used to compare the effects of a regular (one-pass) or interleaved (two-pass) acoustic survey on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches of different shapes and spatial orientations, and a set of parallel or zigzag transects forming a regular or interleaved acoustic survey. The efficiency of a survey is determined by the adequacy of a reconstructed field to that originally generated, which is evaluated by calculating their correlations. Regarding the immovable fields, the efficiency of a regular or interleaved acoustic survey was tested with the following two alternative assumptions: (1) the entire survey was completed; (2) the survey was interrupted for some reason at the moment when one transect remained non-accomplished. In the former case, the efficiencies of both acoustic surveys were nearly the same; in the latter case, the efficiency of an interleaved survey was superior to that of a regular one. With respect to movable fields, the efficiency of the completed interleaved surveys was even higher than that of the regular ones. Thus, the results obtained allow us to conclude that an interleaved survey is expedient in cases where there is no preference regarding the position of a vessel for further work.     

Patchy distribution fields: a zigzag survey design and reconstruction adequacy

A mathematical model was used to examine the effects of a zigzagsurvey design on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches(or gaps) of different shapes and spatial orientations, and an acoustic survey by a set of transects forming a zigzag. Adequacyof the reconstructed fields to those originally generated was evaluated by calculating their correlations (r). A prioriinformation on the autocorrelation radius for the field in thedirection of a survey (R) allows optimization of survey design and the algorithm of data analysis. A patchy field can be reconstructed properly (r > 0.80) if the distance between transects D < (1.0–1.5)R. If a priori information onthe field is not available, the autocorrelation radius should bedetermined when reconstructing the field, i.e. a posteriori. In cases of field movement, the criterion for choosing a survey direction is based on the relationship betweenthe dimension of moving patches in the direction of movement andthat of the surveyed area. The results obtained indicate that, for a fixed transect spacing, zigzag pattern allow less adequatereconstruction of an original distribution field (in cases of both immovable and movable fields) than corresponding parallel pattern.     

Assessment of river water quality during snowmelt and base flow periods in two catchment areas with different land use

River water quality was evaluated with respect to eutrophication and land use during spring snowmelt and summer base flow periods in Abashiri (mixed cropland-livestock farming) and Okoppe (grassland-based dairy cattle farming), eastern Hokkaido, Japan. Water from rivers and tributaries was sampled during snowmelt and summer base flow periods in 2005, and river flow was measured. Total N (TN), NO₃–N, and Si concentrations were determined using standard methods. Total catchment and upland areas for each sampling site were determined with ArcGIS hydrology modeling software and 1:25,000-scale digital topographic maps. Specific discharge was significantly higher during snowmelt than during base flow. In both areas, TN concentrations increased, whereas Si concentrations decreased, with increased specific discharge, and were significantly higher during snowmelt. The Si:TN mole ratio decreased to below or close to the threshold value for eutrophication (2.7) in one-third of sites during snowmelt. River NO₃–N concentrations during base flow were significantly and positively correlated with the proportion of upland fields in the catchment in both the Abashiri (r = 0.88, P < 0.001) and Okoppe (r = 0.43, P < 0.01) areas. However, the regression slope, defined as the impact factor (IF) of water quality, was much higher in Abashiri (0.025) than in Okoppe (0.0094). The correlations were also significantly positive during snowmelt in both areas, but IF was four to eight times higher during snowmelt than during base flow. Higher discharge of N from upland fields and grasslands during snowmelt and the resulting eutrophication in estuaries suggest that nutrient discharge during snowmelt should be taken into account when assessing and monitoring the annual loss of nutrients from agricultural fields.     

Implications of climate change for evaporation from bare soils in a Mediterranean environment

The purpose of this study was to predict quantitative changes in evaporation from bare soils in the Mediterranean climate region of Turkey in response to the projections of a regional climate model developed in Japan (hereafter RCM). Daily RCM data for the estimation of reference evapotranspiration (ET _r) and soil evaporation were obtained for the periods of 1994–2003 and 2070–2079. Potential evaporation (E _p) from bare soils was calculated using the Penman–Monteith equation with a surface resistance of zero. Simulation of actual soil evaporation (E _a) was carried out using Aydin model (Aydin et al., Ecological Modelling 182:91–105, 2005) combined with Aydin and Uygur (2006, A model for estimating soil water potential of bare fields. In Proceedings of the 18th International Soil Meeting (ISM) on Soils Sustaining Life on Earth, Managing Soil and Technology, Sanliurfa, 477–480pp.) model of predicting soil water potential at the top surface layer of a bare soil, after performances of Aydin model (R ² = 94.0%) and Aydin and Uygur model (R ² = 97.6) were tested. The latter model is based on the relations among potential soil evaporation, hydraulic diffusivity, and soil wetness, with some simplified assumptions. Input parameters of the model are simple and easily obtainable such as climatic parameters used to compute the potential soil evaporation, average diffusivity for the drying soil, and volumetric water content at field capacity. The combination of Aydin and Aydin and Uygur models appeared to be useful in estimating water potential of soils and E _a from bare soils, with only a few parameters. Unlike ET _r and E _p projected to increase by 92 and 69 mm (equivalent to 8.0 and 7.3% increases) due to the elevated evaporative demand of the atmosphere, respectively, E _a from bare soils is projected to reduce by 50 mm (equivalent to a 16.5% decrease) in response to a decrease in rainfall by 46% in the Mediterranean region of Turkey by the 2070s predicted by RCM, and consequently, to decreased soil wetness in the future.     

Feeding of western gray whales during a seismic survey near Sakhalin Island, Russia

Exxon Neftegas Limited, as operator of the Sakhalin-1 consortium, is developing oil and gas reserves on the continental shelf off northeast Sakhalin Island, Russia. DalMorNefteGeofizika (DMNG) on behalf of the Sakhalin-1 consortium conducted a 3-D seismic survey of the Odoptu license area during 17 August–9 September 2001. A portion of the primary feeding area of the endangered western gray whale (Eschrichtius robustus) is located in the vicinity of the seismic survey. This paper presents data to assess whether western gray whale bottom feeding activity, as indicated by visible mud plumes, was affected by seismic operations. The mitigation and monitoring program associated with the seismic survey included aerial surveys during 19 July–19 November 2001. These aerial surveys documented the local and regional distribution, abundance, and bottom feeding activity of western gray whales. Data on gray whale feeding activity before, during and after the seismic survey were collected, with the whales assumed to be feeding on the benthos if mud plumes were observed on the surface. The data were used to assess the influence of seismic survey and other factors (including environmental) on feeding activity of western gray whales. A stepwise multiple regression analysis failed to find a statistically significant effect (α = 0.05) of the seismic survey on frequency of occurrence of mud plumes of western gray whales used as a proxy to evaluate bottom feeding activity in Piltun feeding area. The regression indicated that transect number (a proxy for water depth, related to distance from shore) and swell height (a proxy for sea state) were the only variables that had a significant effect on frequency of whale mud plumes. It is concluded that the 2001 seismic survey had no measurable effect (α = 0.05) on bottom feeding activity of western gray whales off Sakhalin Island. M.W. Newcomer, deceased     

Appropriate Range of Quality Control Limit for Wet Deposition Monitoring

Precipitation chemistry programs in different regions of the world have different quality control limits for ion balance parameter IPD in wet deposition monitoring (R). The range of R values was calculated by assuming sample rainwater models in this paper. It was found that R was influenced by the inorganic ion types, total ion concentrations (IS) and the accuracy of the measurements (a _i ). R was defined and calculated as a function of a _i and ion concentration C _i . R values of different types of wet deposition (including ocean type, continent type and combined oceanic and continental deposition) were different. There were also differences between the samples of the same type if the samples had different total ion concentrations. When IS>100 μeq/L, the ranges of R of ocean type deposition, continent type deposition and combined oceanic and continental deposition were 5∼7%, 5∼9% and 5∼11%, respectively. When IS<50 μeq/L, if a _i was 100%, the range of R was 33∼71% because of the lower accuracy. It was also found that R of each criterion was in the range of R as calculated in this paper when IS≥50 μeq/L, but when IS<50 μeq/L, the criteria varied greatly in their R values.     

Accuracy assessment of a small-area method for estimating the spatial distribution of the degree of tree damage

Tree damage, gauged by the amount of defoliation, is one of the basic criteria used to determine treatments for protected and economic forests. Monitoring should include an assessment of the degree of tree damage in different spatial scales. Therefore, in addition to the commonly applied large-area methods, small-area methods should be used. The aim of the paper is to present the results of the accuracy assessment of a small-area method, proposed by Podlaski (2005) [Podlaski, R. (2005). Inventory of the degree of tree defoliation in small areas. Forest Ecology and Management, 215, 361–377], for monitoring the degree of tree damage. The degree of tree damage was shown in sub-blocks P₃ of the system of information on natural environment (SINUS). To estimate the spatial distribution of the degree of tree defoliation, survey sampling, based on simple random sampling with replacement (SRSWR), was used. The degree of damage to fir (Abies alba Mill.) and beech (Fagus sylvatica L.) was analysed in the Święty Krzyż forest section in the Świętokrzyski National Park. The maximum total estimation errors for the proportion of trees with a degree zero of damage, and with second and third degrees of damage together (for α = 0.05) were at most 30.8% for fir and 24.3% for beech trees. For standard, small-area evaluations, these are satisfactory values. In the Święty Krzyż forest section, the number of P₃ sub-blocks with 0.00–5.00% of undamaged trees and with 80.01–100.00% of moderately- or severely-damaged trees was significantly greater for fir than for beech. These results indicate that the fir population was unhealthier than the beech group in the study area. P₃ sub-blocks of the SINUS system, in which the proportion of the healthiest trees was highest, were situated at the forest margin, bordering on meadows and arable fields (in the case of fir) and forming dense patches consisting of several sub-blocks, or occurring singly in the whole study area (in the case of beech). The results show the significant differentiation of forest tree health in small areas.     

Seasonal Variation of the Particle Size Distribution of n-Alkanes and Polycyclic Aromatic Hydrocarbons (PAHs) in Urban Aerosol of Guangzhou, China

Seasonal aerosol samples have been collected by Andersen Hi-Vol pumping system equipped with a five stage cascade impactor and a backup filter (size range: 10–7.2 μ m, 7.2–3.0 μ m, 3.0–1.5 μ m, 1.5–0.95 μ m, 0.95–0.49 μ m, ≤0.49 μ m) in the Liwan district, Guangzhou. n-Alkanes were measured using gas chromatography and PAHs were measured using gas chromatography/mass spectrometry analysis. The bimodal log-normal distributions of n-alkanes and semi-volatile PAHs were found, while for non-volatile PAHs that was unimodal, so much as the mode of semi-volatile PAHs was similar with that of the particles. The n-alkanes and PAHs were preferably associated with fine particles. C _max (carbon number maximum) (C₂₂–C₂₆), CPI (carbon preference index) (1.12–1.21), U/R (unresolved to resolved components ratio) (7.42–10.7), wax% (0.9–3.12%) and the diagnostic ratios for PAHs revealed that vehicular emission was the major source of these organic compounds during the study periods, while the contribution of epicuticular waxes emitted by terrestrial plants was minor. CPI₂ (values for petrogenic hydrocarbons), CPI₃ (values for biogenic n-alkanes) and wax% revealed that the natural preferentially accumulated in the larger aerosol while the anthropogenic in the smaller. In addition, the different MMDs (mass median diameters) for n-alkanes and PAHs were observed in different seasons. The MMDs for n-alkanes and PAHs were higher in autumn/winter than those in spring/summer. The seasonal effect was related to the hydrocarbon content in the individual particulate fractions, showing a preferential association of n-alkanes and PAHs with larger particles in the autumn/winter season.     

Atmosphic pollutants study of particles ionic species during high wind speed (>7 m s-1) near Taiwan strait in central Taiwan from 2004 to 2005

This study monitored atmospheric pollutants during high wind speed (> 7 m s⁻¹) at two sampling sites: Taichung Harbor (TH) and Wuci traffic (WT) during March 2004 to January 2005 in central Taiwan. The correlation coefficient (R ²) between TSP, PM_2.5, PM_2.5−10 particle concentration vs. wind speed at the TH and WT sampling site during high wind speed (< 7 m s⁻¹) were also displayed in this study. In addition, the correlation coefficients between TSP, PM_2.5 and PM_2.5−10 of ionic species vs. high wind speed were also observed. The results indicated that the correlation coefficient order was TSP > PM_2.5−10 > PM_2.5 for particle at both sampling sites near Taiwan strait. In addition, the concentration of Cl⁻, NO₃ ⁻, SO₄ ²⁻, NH₄ ⁺, Mg²⁺, Ca²⁺ and Na⁺ were also analyzed in this study.     

Interactions of Atmospheric Deposition with Coniferous Canopies in Estonia

Throughfall and open field bulk precipitation were measured at three coniferous sites during 1995–2002 in the framework of ICP Integrated Monitoring and at five coniferous sites during 1996–2002 in the framework of ICP Forests (Level II). The coniferous canopies acted as a sink for nitrate and ammonium and as a source for base cations: Ca²⁺, Mg²⁺ and K⁺. The estimated share of SO₄–S dry deposition from total deposition was 1.5–4 times higher for dormant period compared to growing period. During the study period average annual throughfall and bulk deposition of SO₄–S decreased significantly, 2.8 and 2.3 times, respectively. Throughfall enrichment with base cations increased in the order Mg < Na < Ca < K. Using Na as a tracer ion, average dry deposition and canopy leaching were calculated. Leaching was the dominant process for TF enrichment by potassium. Leaching of base cations occurred during growing as well as dormant period. The calculated internal flux of Ca²⁺ and Mg²⁺ varied in the range of 0.6–2.0 and 0.6–1.2 kg ha⁻¹ per year in spruce and pine stands, respectively. The internal circulation of K⁺ was significantly higher (8.9–10.9 kg ha⁻¹ per year) in spruce stands than in pine stands (2.7–4.4 kg ha⁻¹ per year).     

Risk Assessment of Organic Pesticides Pollution in Surface Water of Hangzhou

Advances in research on pollution of organic pesticides (OPs) in surface water, pollution survey and risk assessments of organochlorine pesticides (OCPs) and organophosphorus pesticides (OPPs) of surface water in Hangzhou are conducted. Total concentrations of dichloro-diphenyl-trichloroethane (DDT) and hexachloride-benzene (HCH) in surface water were observed to be 0–0.270 μg/L and 0–0.00625 μg/L respectively. DDE, as a metabolite of DDT and many species of OPP_S were determined in some samples of surface water. Parathion, the main pollutant among OPPs in surface water of Hangzhou, was observed to be 0–0.445 μg/L. Based on these experimental results, health risk assessments on the organic pollution are developed. It is observed that the total risk “R ^T” at present time of surface water in Hangzhou is mainly contributed by organophosphorus pesticides, especially Parathion; HCH and DDT are not the main contaminants; on the contrary, organophosphorous pesticides, especially Parathion, must be of concern at the present time.     

Spatial Scale of Autocorrelation of Assemblages of Benthic Invertebrates in Two Upland Rivers in South-Eastern Australia and Its Implications for Biomonitoring and Impact Assessment in Streams

Spatial autocorrelation in ecological systems is a critical issue for monitoring (and a general understanding of ecological dynamics) yet there are very few data available, especially for riverine systems. Here, we report here on assemblage-level autocorrelation in the benthic-invertebrate assemblages of riffles in two adjacent, relatively pristine rivers in south-eastern Victoria, Australia (40-km reaches of the Wellington [surveys in summers of 1996 and 1997] and Wonnangatta Rivers [survey in summer of 1996 only], with 16 sites in each river). We found that analyses were similar if the data were resolved to family or to species level. Spatial autocorrelation was assessed by using Mantel-tests for the data partitioned into different sets of spatial separations of survey sites (e.g. 0–6 km, 6–12 km, etc.). We found strong small-scale (≤6 km) autocorrelation in the Wellington River, which is consistent with known dispersal abilities of many aquatic invertebrates. Surprisingly, there were strong negative correlations at longer distance classes for the Wellington River in one of the two summers (20–40 km) and the Wonnangatta River (12–20 km). That two largely unimpacted, adjacent rivers should have such different autocorrelation patterns suggests that impact assessment cannot assume dependence or independence of sites a priori. We discuss the implications of these results for use of “reference” sites to assess impacts at nominally affected sites.     

Long-Term Biomonitoring with Lichens: Comparing Data from Different Sampling Procedures

The results of 2 biomonitoring surveys carried out over a span of 7 years (1996–2003) in the Province of Ancona (central Italy) using the diversity of epiphytic lichens are reported. According to a recent standardisation process, in 2003 a different protocol was adopted, raising the problem of data comparability with the previous study. For this reason, a ‘pilot study' was carried out in which the two sampling tactics were both used. The high correlation between biodiversity values proved by the results of this study, let to compare the two lichen diversity data set and to identify the impacts on lichen communities of temporal variations of atmospheric quality in the study area, showing a reduction of stations classified as ‘altered' and an increase in stations classified as ‘semi-altered'. Furthermore, possible reduction of the sampling effort was detected performing a study on the sampling density of the survey area, indicating in a 6 × 6 km grid the best density to adopt in order to save financial resources for the future.     

Assessment of Temporal and Spatial Variation of Nitrate Removal in Riparian Zones

The aim of this study was to monitor long-term temporal and spatial groundwater NO₃^- removal efficiencies in different riparian zones via a limited number of sampling wells. Groundwater NO₃^- concentrations were measured fortnightly or monthly over a period of two years using transects of ground water sampling wells. Depending on the level of the NO₃^- load (up to 120mgNL^-1 at the input side of the riparian zone a distance of 10 to 30m was needed to remove NO₃^- from the groundwater below 11.3mgNL^-1. Considering all seasons, the mixed vegetation and grass riparian site succeeded to remove groundwater NO₃^- efficiently (92—100% within a distance of 30m. The forested riparian zone removed 72—90% of the total NO₃^- input within a distance of 30m. Evidence emerged that NO₃^- could also be removed actively at depths up to 2m, due to the presence of organically enriched layers of alluvial deposits or roots. Our four dimensional approach (three dimensional space and time), in combination with a limited number of sampling wells, was shown to be a useful monitoring tool to assess the variability of NO₃^- removal in riparian zones.     

Assessing Respiratory Morbidity Through Pollution Status and Meteorological Conditions for Delhi

The study focuses on assessing the status of respiratory morbidity in Delhi over a four years period from 2000–2003. An attempt was made to investigate the role of important pollutants (SO₂, NO₂, SPM and RSPM) and various meteorological factors (temperature minimum & maximum, relative humidity at 0830 and 1730 hrs. and wind speed) in being responsible for respiratory admissions on account of COPD, asthma and emphysema. The study showed that winter months had greater exposure risk as pollutants often get trapped in the lower layers of atmosphere resulting in high concentrations. Statistical analysis revealed that two pollutants have significant positive correlation with the number of COPD cases viz., SPM (r = 0.474; p < 0.01) and RSPM (r = 0.353; p < 0.05), while a meteorological factor temperature (minimum) has a significant negative correlation (r = −0.318; p < 0.05) with COPD. Stepwise multiple regression analysis was performed for COPD as dependent variable and R² value of 0.33 was obtained indicating that SPM and RH(1730) were able to explain 33 percent variability in COPD. The partial correlation of SPM and RH(1730) on COPD was higher than any other combination and therefore they can be regarded as important contributing variables on COPD.     

Water quality assessment: surface water sources used for drinking and irrigation in Zaria,Nigeria are a public health hazard

We assessed the quality and pollution status of source surface waters in Zaria, Nigeria by monitoring the nature, cause and extent of pollution in Samaru stream, Kubanni River and Kubanni dam over a period of 10 months, between March and December 2002. A total of 228 water samples was collected from 12 sites and analysed for a total of ten physicochemical and one bacteriological quality indicators, using standard methods. Aesthetic water quality impairment parameters were also observed. The mean values of most water quality parameters were significantly higher (P < 0.05) in both the stream and river than in the dam. There was no significant correlation between faecal coliform counts (FCC) and water temperature (in the range 15–33°C); pH (5.77–7.32); and turbidity (1.4–567 NTU). The high FCC ranged from 2.0 × 10¹ to 1.6 × 10⁶ MPN/100 ml and exceeded the WHO standards for drinking water and water used for fresh-produce irrigation, and correlated positively (P < 0.05) with conductivity (in the range 68–1,029 μS/cm); TDS (10.0–70.0 mg/l); TSS (10.0–70.0 mg/l); Cl (7.5–181 mg/l); PO₄⁻P (0.01–0.41 mg/l); NO₃⁻N (0.6–3.8 mg/l) and BOD₅ (0.1–14.9 mg/l). The main pollution sources were municipal wastewater, stormwater runoffs, the ABU sewage treatment plant, abattoir effluents and irrigation farms treated with chemical fertilisers. We conclude that these water bodies are potentially hazardous to public health and that proper sewage treatment and river quality monitoring are needed to warn against hazards to public health.     

Assessing the effects of land use changes on soil sensitivity to erosion in a highland ecosystem of semi-arid Turkey

There has been increasing concern in highlands of semiarid Turkey that conversion of these systems results in excessive soil erosion, ecosystem degradation, and loss of sustainable resources. An increasing rate of land use/cover changes especially in semiarid mountainous areas has resulted in important effects on physical and ecological processes, causing many regions to undergo accelerated environmental degradation in terms of soil erosion, mass movement and reservoir sedimentation. This paper, therefore, explores the impact of land use changes on land degradation in a linkage to the soil erodibility, RUSLE-K, in Cankiri–Indagi Mountain Pass, Turkey. The characterization of soil erodibility in this ecosystem is important from the standpoint of conserving fragile ecosystems and planning management practices. Five adjacent land uses (cropland, grassland, woodland, plantation, and recreational land) were selected for this research. Analysis of variance showed that soil properties and RUSLE-K statistically changed with land use changes and soils of the recreational land and cropland were more sensitive to water erosion than those of the woodland, grassland, and plantation. This was mainly due to the significant decreases in soil organic matter (SOM) and hydraulic conductivity (HC) in those lands. Additionally, soil samples randomly collected from the depths of 0–10 cm (D ₁) and 10–20 cm (D ₂) with irregular intervals in an area of 1,200 by 4,200 m sufficiently characterized not only the spatial distribution of soil organic matter (SOM), hydraulic conductivity (HC), clay (C), silt (Si), sand (S) and silt plus very fine sand (Si + VFS) but also the spatial distribution of RUSLE-K as an algebraically estimate of these parameters together with field assessment of soil structure to assess the dynamic relationships between soil properties and land use types. In this study, in order to perform the spatial analyses, the mean sampling intervals were 43, 50, 64, 78, 85 m for woodland, plantation, grassland, recreation, and cropland with the sample numbers of 56, 79, 72, 13, and 69, respectively, resulting in an average interval of 64 m for whole study area. Although nugget effect and nugget effect–sill ratio gave an idea about the sampling design adequacy, the better results are undoubtedly likely by both equi-probable spatial sampling and random sampling representative of all land uses.