Soil salinity evolution and its relationship with dynamics of groundwater in the oasis of inland river basins: case study from the Fubei region of Xinjiang Province, China

Soil salinization is an important worldwide environmental problem, especially in arid and semi-arid regions. Knowledge of its temporal and spatial variability is crucial for the management of oasis agriculture. The study area has experienced dramatic change in the shallow groundwater table and soil salinization during the 20th century, especially in the past two decades. Classical statistics, geostatistics and geographic information system (GIS) were applied to estimate the spatial variability of the soil salt content in relation to the shallow groundwater table and land use from 1983 to 2005. Consumption of reservoir water for agricultural irrigation was the main cause of a rise in the shallow groundwater table under intense evapotranspiration conditions, and this led indirectly to soil salinization. The area of soil salt accumulation was greater in irrigated than in non-irrigated landscape types with an increasing of 40.04% from 1983 to 2005 in cropland at ∼0.43 t ha⁻¹ year⁻¹, and an increase at ∼0.68 t ha⁻¹ year⁻¹ in saline alkaline land. Maps of the shallow groundwater table in 1985 and 2000 were used to deduce maps for 1983 and 1999, respectively, and the registration accuracy was 99%.     

Persistent Organochlorine Pesticide Residues in Soil and Surface Water of Northern Indo-Gangetic Alluvial Plains

This study reports the concentration levels and distribution pattern of the organochlorine pesticide (OCPs) residues in the soil and surface water samples collected from the northern Indo-Gangetic alluvial plains. A total of 31 soil and 23 surface water samples were collected from the study region in Unnao district covering an area of 2150 km² and analyzed for aldrin, dieldrin, endrin, HCB, HCH isomers, DDT isomers/metabolites, endosulfan isomers (α and β), endosulfan sulfate, heptachlor and its metabolites, α-chlordane, γ-chlordane and methoxychlor. In both the soil and surface water samples β- and δ-isomers of HCH were detected most frequently, whereas, methoxychlor was the least detected pesticide. The results showed contamination of soil and surface water of the region with several persistent organic pesticides. The total OCPs level ranged from 0.36–104.50 ng g^–1 and 2.63–3.72 μg L^–1 in soil and surface water samples, respectively.     

Spatial variability of soil <Superscript>137</Superscript>Cs in the South Caspian region

In a comprehensive program of environmental radioactivity survey in South Caspian region,¹³⁷Cs inventories in soil has been measured at more than 50 sites in the Iranian northern province of Guilan. This has been the first wide-range survey of soil radionuclide inventories in the narrow band sensitive ecosystem of south Caspian shore. Radioactivity measurements were carried out using HPGe gamma-spectrometry system. The activity concentration of ¹³⁷Cs in surface soil exhibits a mean value of 17.6 ± 9.4 Bq kg⁻¹, with a range of 2.3–41.7 Bq kg⁻¹. In many sites, split-level sampling method has been applied down to a depth of 20 cm. There were found generally two profiles. Most profiles exhibit a negative exponential distribution, while others revealed a clear subsurface peak in 5–10-cm layer. Cesium deposition in the study area has been estimated to be in the range of 0.38–2.9 kBq m⁻² with a mean value of 1.7 kBq m⁻². Distribution patterns of ¹³⁷Cs concentration levels and deposition values have been estimated using Kriging interpolation method. Observed hotspots in deposition pattern coincide with areas of higher precipitation.     

Assessing environmental impacts of treated wastewater through monitoring of fecal indicator bacteria and salinity in irrigated soils

To assess the potential for treated wastewater irrigation to impact levels of fecal indicator bacteria (FIB) and salinity in irrigated soils, levels of Escherichia coli, Enterococcus, and environmental covariates were measured in a treated wastewater holding pond (irrigation source water), water leaving the irrigation system, and in irrigated soils over 2 years in a municipal parkland in Arizona. Higher E. coli levels were measured in the pond in winter (56 CFU 100 mL⁻¹) than in summer (17 CFU 100 mL⁻¹); however, in the irrigation system, levels of FIB decreased from summer (26 CFU 100 mL⁻¹) to winter (4 CFU 100 mL⁻¹), possibly related to low winter water use and corresponding death of residual bacteria within the system. For over 2 years, no increase in FIB was found in irrigated soils, though highest E. coli levels (700 CFU g⁻¹ soil) were measured in deeper (20–25 cm) soils during summer. Measurements of water inputs vs. potential evapotranspiration indicate that irrigation levels may have been sufficient to generate bacterial percolation to deeper soil layers during summer. No overall increase in soil salinity resulting from treated wastewater irrigation was detected, but distinct seasonal peaks as high as 4 ds m⁻¹ occurred during both summers. The peaks significantly declined in winter when surface ET abated and more favorable water balances could be maintained. Monitoring of seasonal shifts in irrigation water quality and/or factors correlated with increases and decreases in FIB will aid in identification of any public health or environmental risks that could arise from the use of treated wastewater for irrigation.     

Effect of long-term application of treated sewage water on heavy metal accumulation in vegetables grown in Northern India

Use of industrial and wastewater for irrigation is on the rise in India and other developing countries because of scarcity of good-quality irrigation water. Wastewaters contain plant nutrients that favour crop growth but leave a burden of heavy metals which can enter the food chain and is a cause of great concern. The present study was undertaken on the long-term impact of irrigation with treated sewage water for growing vegetables and the potential health risk associated with consumption of such vegetable. Treated sewage water (TSW), groundwater (GW), soil and plant samples were collected from peri urban vegetable growing areas of Northern India (Varanasi) and analysed to assess the long-term effect of irrigation with TSW on Cd, Cr, Ni and Pb build-up in soils and its subsequent transfer into commonly grown vegetable crops. Results indicate that TSW was richer in essential plant nutrients but contained Cd, Cr and Ni in amounts well above the permissible limits for its use as irrigation water. Long-term application of TSW resulted in significant build-up of total and DTPA extractable Cd, Cr, Ni and Pb over GW irrigated sites. TSW also resulted in slight lowering in pH, increase in organic carbon (1.6 g kg^− 1) and cation exchange capacity (5.2 cmol kg^− 1). The tissue metal concentration and relative efficiency of transfer of heavy metals from soil to plant (transfer factor) for various groups of vegetables were worked out. Radish, turnip and spinach were grouped as hyper accumulator of heavy metals whereas brinjal and cauliflower accumulated less heavy metals. Health risk assessment by consumption of vegetables grown with TSW indicated that all the vegetables were safe for human consumption. However, significant accumulation of these heavy metals in soil and plant needs to be monitored.     

The geostatistic-based spatial distribution variations of soil salts under long-term wastewater irrigation

The purpose of this study was to determine and evaluate the spatial changes in soil salinity by using geostatistical methods. The study focused on the suburb area of Beijing, where urban development led to water shortage and accelerated wastewater reuse to farm irrigation for more than 30 years. The data were then processed by GIS using three different interpolation techniques of ordinary kriging (OK), disjunctive kriging (DK), and universal kriging (UK). The normality test and overall trend analysis were applied for each interpolation technique to select the best fitted model for soil parameters. Results showed that OK was suitable for soil sodium adsorption ratio (SAR) and Na⁺ interpolation; UK was suitable for soil Cl^? and pH; DK was suitable for soil Ca²⁺. The nugget-to-sill ratio was applied to evaluate the effects of structural and stochastic factors. The maps showed that the areas of non-saline soil and slight salinity soil accounted for 6.39 and 93.61 %, respectively. The spatial distribution and accumulation of soil salt were significantly affected by the irrigation probabilities and drainage situation under long-term wastewater irrigation.     

Using of high-resolution topsoil magnetic screening for assessment of dust deposition: comparison of forest and arable soil datasets

Magnetic susceptibility (κ) is an easily detectable geophysical parameter that can be used as a proxy or semi-quantitative tracer of atmospheric industrial and urban dusts deposited in topsoil. An enhanced κ value of topsoil is in many cases also associated with high concentrations of soil pollutants (mostly heavy metals). High-resolution magnetic screening of topsoil in areas of high pollution influx is a useful tool for detection of pollution “hot spots”. General and regional screening maps with a grid density of 10 or 5 km have been performed on the basis of forest topsoil measurement only. The purpose of this study was to perform high-resolution magnetic screening with different grid densities in both forested and agricultural areas (arable land). Our large study area (ca. 200 km²) was located in a relatively more polluted region of the central part of Upper Silesia, and a second (small) one (ca. 100 m²) was located in the western part of Upper Silesia, with considerably lower influx of pollution. In the framework of this study, we applied a statistical comparison of data obtained in forested areas and on arable land. The arable soil showed statistically significantly lower κ values, the result of “physical dilution” of the arable layer caused by annual ploughing. Thus arable soils must be avoided during high-resolution field measurement. From semivariograms, it was clear that the spatial correlations in forest topsoil are much stronger than in arable soil, which suggests that a denser measurement grid is required in forested areas.     

Exposure Modeling of Benzene Exploiting Passive–Active Sampling Data

The objective of the present study is the exploitation of active sampling personal exposure data in assessing the factors that affect exposure to benzene in combination with the widely accepted scheme of passive sampling—time microenvironment–activity diaries (TMAD). The campaign included personal exposure measurements with both passive and active sampling in several microenvironments, evaluation of TMAD kept by the volunteers, and a variety of environmental data (ambient air benzene determination, traffic and meteorological observations). Due to the relatively elevated benzene traffic emissions, average personal exposure was determined to be equal to 8.9 μg/m³, ranging between 5 and 20 μg/m³, which is a value highly related to the average urban concentration (9.2 μg/m³). The information gained from TMAD was embedded (in terms of spatial and temporal distribution) into three zones respectively, in order to draw statistically significant conclusions about the exposure levels and the activity patterns. The contribution of the activities to the overall amount of exposure was further quantified and refined by active sampling measurements. These data revealed that driving in a traffic-congested road was the main activity leading to elevated exposure levels (up to 70 μg/m³), followed by walking on the roadside of a congested road (up to 35 μg/m³). Indoor exposure to benzene was in general lower than outdoor (indicating that traffic is the dominant source of benzene emissions in the wider area), and it was significantly affected by the presence of environmental tobacco smoke. The higher significance of the regression coefficients obtained by statistical analysis of the active sampling data was fundamental for the development of a regression-based prediction exposure model. The model was evaluated through comparison with the passive sampling data, which were considered as an unknown but realistic data exposure pattern. The model performed very well in terms of expressing the variance of the exposure data with an average score of R ² equal to 0.935. All of the above indicate that active sampling is a necessary albeit more laborious tool that needs to be used as a complement to passive sampling for precise quantification of the factors determining personal exposure patterns.     

Assessment of spatial variability in some soil properties as related to soil salinity and alkalinity in Bafra plain in northern Turkey

The objectives of this study were to assess the variability in soil properties affecting salinity and alkalinity, and to analyze spatial distribution patterns of salinity (EC) and alkalinity (ESP) in the plain, which was used irrigation agriculture with low quality waters. Soil samples were collected from 0–30cm, 30–60cm, 60–90cm and 90–120cm soil depths at 60 sampling sites. Soil pH had the minimum variability, and hydraulic conductivity (Ks) had the maximum variability at all depths. The mean values of pH, EC, ESP and Ks increased while the mean values of CEC decreased with soil depth. Values pH, EC and ESP were generally high in the east and northeastern sides. Soil properties indicated moderate to strong spatial dependence. ESP and pH were moderately spatially dependent for three of the four depths, EC exhibited moderate spatial dependence for one of the four depths, CEC had a moderate spatial dependence at all depths, and Ks exhibited a strong spatial dependence. EC, CEC, and ESP were considerably variable in small distances. The spatial variability in small distances of EC, CEC, pH and ESP generally increased with depth. All geostatistical range values were greater than 1230m. It was inferred that the strong spatial dependency of soil properties would be resulted in extrinsic factors such as ground water level, drainage, irrigation systems and microtopography.     

Pollution by organochlorine pesticides in Navachiste-Macapule,Sinaloa, Mexico

The lagoon system of Navachiste-Macapule is located in northern Mexico, in the state of Sinaloa, with an area of 24,000 ha. The main economic activity in the area is agriculture, and the lagoon lies next to the irrigation district ID-063 which covers 116,615 ha. The purpose of this study is to evaluate the levels of pollution generated by organochlorine pesticides (OC) in the surface sediments of the lagoon and in the agricultural drains of the ID-063 that are supposedly transported into the system as a result of agricultural activities and runoff from adjacent land. For this, between 2006 and 2007, 45 surface sediment samples were collected (warm dry, rainy and cold dry) from 15 sampling sites, during the three climatic seasons. Of these, eight were located inside the lagoon in marine conditions (salinity >31 PSU) and seven in the agricultural drains of the ID-063 in freshwater conditions (salinity <5 PSU). The average concentration of the OC in the sediments was 44.75 ng g^− 1, among which the group of the alicyclic compounds presented the greatest concentrations. The average value of the total organic carbon (TOC) in the sediments of the system was 0.90%. The sediments collected inside the lagoon had an average OC concentration of 18.97 ng g^− 1, and the predominant type of sediment was fine to very fine sand. The average OC concentration in the sediments collected in the agricultural drains was 75.69 ng g^− 1, where fine sediments (silt) were predominant. The presence of methoxychlor, endrin and heptachlor suggested that these compounds were continuously used in the system, even though their use is forbidden in Mexico.     

Fractionation and mobility of cadmium and zinc in urban vegetable gardens of Kano, Northern Nigeria

Metal fractionation provides information on mobility and stability of various metal species which can be used to evaluate the movement of such metals in soils. The effect of wastewater irrigation on the fractions, spatial distribution, and mobility of cadmium (Cd) and zinc (Zn) was investigated in five urban gardens in Kano, Nigeria. Concentration of total Zn in the surface soils (0–20 cm) ranged from 121 to 207 mg kg^− 1 while Cd concentration was 0.3–2.0 mg kg^− 1. Speciation of both heavy metals into seven operationally defined fractions indicated that the most reactive forms extracted with ammonium nitrate and ammonium acetate, also considered as the bioavailable fractions, accounted for 29–42% of total Cd and 22–54% of total Zn, respectively. The weakly bound fractions of Cd and Zn reached up to 50% of the total Cd and Zn concentrations in the soils. Such high proportions of labile Cd and Zn fractions are indicative of anthropogenic origins, arising from the application of wastewater for irrigation and municipal biosolids for soil fertility improvement. Thus, given the predominance of sandy soil textures, high concentrations of labile Cd and Zn in these garden soils represent a potential hazard for the redistribution and translocation of these metals into the food chain and aquifer.     

Spatial assessment of soil salinity in the Harran Plain using multiple kriging techniques

The Harran Plain is located in the southeastern part of Turkey and has recently been developed for irrigation agriculture. It already faces soil salinity problems causing major yield losses. Management of the problem is hindered by the lack of information on the extent and geography of the salinization problem. A survey was carried out to delineate the spatial distribution of salt-affected areas by randomly selecting 140 locations that were sampled at two depths (0 to 30 and 30 to 60 cm) and analyzed for soil salinity variables: soil electrical conductivity (EC), soluble cations (Ca^2+, Mg²⁺, Na⁺, and K⁺), soluble anions (SO ₄ ^2? , Cl^?), exchangeable Na⁺ (me 100 g^?1) and exchangeable sodium percentage. Terrain attributes (slope, topographical wetness index) were extracted from the digital elevation model of the study area. Variogram analyses after log transformation and ordinary kriging (OK) were applied to map spatial patterns of soil salinity variables. Multivariate geostatistical methods—regression kriging (RK) and kriging with external drift (KED)—were used using elevation and soil electrical conductivity data as covariates. Performances of the three estimation methods (OK, RK, and KED) were compared using independent validation samples randomly selected from the main dataset. Soils were categorized into salinity classes using disjunctive kriging (DK) and ArcGIS, and classification accuracy was tested using the kappa statistic. Results showed that soil salinity variables all have skewed distribution and are poorly correlated with terrain indices but have strong correlations among each other. Up to 65 % improvement was obtained in the estimations of soil salinity variables using hybrid methods over OK with the best estimations obtained with RK using EC_0–30 as covariate. DK–ArcGIS successfully classified soil samples into different salinity groups with overall accuracy of 75 % and kappa of 0.55 (p?<?0.001).     

Water quality and dissolved inorganic fluxes of N,P, SO<Subscript>4</Subscript>, and K of a small catchment river in the Southwestern Coast of India

The southwestern coast of India is drained by many small rivers with lengths less than 250 km and catchment areas less than 6,500 km². These rivers are perennial and are also the major drinking water sources in the region. But, the fast pace of urbanization, industrialization, fertilizer intensive agricultural activities and rise in pilgrim tourism in the past four to five decades have imposed marked changes in water quality and solute fluxes of many of these rivers. The problems have aggravated further due to leaching of ionic constituents from the organic-rich (peaty) impervious sub-surface layers that are exposed due to channel incision resulting from indiscriminate instream mining for construction-grade sand and gravel. In this context, an attempt has been made here to evaluate the water quality and the net nutrient flux of one of the important rivers in the southwestern coast of India, the Manimala river which has a length of about 90 km and catchment area of 847 km². The river exhibits seasonal variation in most of the water quality parameters (pH, electrical conductivity, dissolved oxygen, total dissolved solids, Ca, Mg, Na, K, Fe, HCO₃, NO₂-N, NO₃-N, P _\text-inorg_{\rm \text{-}inorg}, P _\text-tot_{\rm \text{-}tot}, chloride, SO₄, and SiO₂). Except for NO₃-N and SiO₂, all the other parameters are generally enriched in non-monsoon (December–May) samples than that of monsoon (June–November). The flux estimation reveals that the Manimala river transports an amount of 2,308 t y^− 1 of dissolved inorganic nitrogen, 87 t y^− 1 dissolved inorganic phosphorus, and 9246 t y^− 1 of SO₄, and 1984 t y^− 1 K into the receiving coastal waters. These together constitute about 23% of the total dissolved fluxes transported by the Manimala river. Based on the study, a set of mitigation measures are also suggested to improve the overall water quality of small catchment rivers of the densely populated tropics in general and the south western coast in particular.     

Dynamics and quantification of dissolved heavy metals in the Mahanadi river estuarine system,India

Dynamics of heavy metals such as Fe, Mn, Zn, Cr, Cu, Co, Ni, Pb, and Cd in surface water of Mahanadi River estuarine systems were studied taking 31 different stations and three different seasons. This study demonstrates that the elemental concentrations are extremely variable and most of them are higher than the World river average. Among the heavy metals, iron is present at highest concentration while cadmium is at the least. The spatial pattern of heavy metals suggests that their anthropogenic sources are possibly from two major fertilizer plants and municipal sewage from three major towns as well as agricultural runoff. The temporal variations for metals like Fe, Cu, and Pb exhibit higher values during the monsoon season, which are related to agricultural runoff. Concentrations of Ni, Pb, and Cd exceed the maximum permissible limits of surface water quality in some polluted stations and pose health risks. Dissolved heavy metals like Fe, Mn, Cr, Ni, and Pb exhibit a non-conservative behavior during estuarine mixing, while Zn, Cu, and Co distribution is conservative. Distribution of cadmium in the estuarine region indicates some mobilization which may be due to desorption. The enrichment ratio data suggest that various industrial wastes and municipal wastes contribute most of the dissolved metals in the Mahanadi River. The Mahanadi River transports 18.216 × 10³ t of total heavy metals into the Bay of Bengal and the calculated rate of erosion in the basin is 128.645 kg km^− 2 year^− 1.     

Identification of mangrove water quality by multivariate statistical analysis methods in Pondicherry coast,India

Different multivariate statistical analysis such as, cluster analysis, principal component analysis, and multidimensional scale plot were employed to evaluate the trophic status of water quality for four monitoring stations. The present study was carried out to determine the physicochemical parameters of water and sediment characteristics of Pondicherry mangroves—southeast coast of India, during September 2008–December 2010. Seasonal variations of different parameters investigated were as follows: salinity (10.26–35.20 psu), dissolved oxygen (3.71–5.33 mg/L), pH (7.05–8.36), electrical conductivity (26.41–41.33 ms⁻¹), sulfide (1.98–40.43 mg/L), sediment texture sand (39.54–87.31%), silt (9.89–32.97%), clay (3.06–31.20%), and organic matter (0.94–4.64%). pH, temperature, salinity, sand, silt, clay, and organic matter indicated a correlation at P < 0.01. CA grouped the four seasons in to four groups (pre-monsoon, monsoon, post-monsoon, summer) and the sampling sites in to three groups. PCA identified the spatial and temporal characteristics of trophic stations and showed that the water quality was worse in stations 3 and 4 in the Pondicherry mangroves.     

A GIS-based Study for Optimizing the Total Emission Control Strategy in Lanzhou City

A geographic information system (GIS)-based approach with spatial analysis advantages was developed to optimize the total emission control scheme on air pollution in Lanzhou, combined with Models-3 modeling system and a linear programming model. GIS system is mainly used for establishing a multifactor assessment model to quantitatively divide environmental functional zone, and also used for selecting control sites in linear programming model. The results show that most of the urban area belongs to the second type of function, which is 122.8 km² occupying about 76.5% and the third type of about 32.9 km² occupying 20.5%. To reach the air quality up to the national standard in the entire control area, some large-emission factories need to cut emissions from 19% to 27%. For small but low-emission height sources, the cut may be as high 40%. The improvement of the ecological environment is urgently needed for controlling the background particle pollution in Lanzhou city.     

Uncertainty in watershed response predictions induced by spatial variability of precipitation

Negligence to consider the spatial variability of rainfall could result in serious errors in model outputs. The objective of this study was to examine the uncertainty of both runoff and pollutant transport predictions due to the input errors of rainfall. This study used synthetic data to represent the “true” rainfall pattern, instead of interpolated precipitation. It was conducted on a synthetic case area having a total area of 20 km² with ten subbasins. Each subbasin has one rainfall gauge with synthetic precipitation records. Six rainfall storms with varied spatial distribution were generated. The average rainfall was obtained from all of the ten gauges by the arithmetic average method. The input errors of rainfall were induced by the difference between the actual rainfall pattern and estimated average rainfall. The results show that spatial variability of rainfall can cause uncertainty in modeling outputs of hydrologic, which would be transport to pollutant export predictions, when uniformity of rainfall is assumed. Since rainfall is essential information for predicting watershed responses, it is important to consider the properties of rainfall, particularly spatial rainfall variability, in the application of hydrologic and water quality models.     

Estimating Field Volatility of Soil Fumigants Using CHAIN_2D: Mitigation Methods and Comparison Against Chloropicrin and 1,3-Dichloropropene Field Observations

Academic, government, and industrial field researchers have generated a significant database of field studies of the volatility of soil applied fumigants. However, limited work exists in validating physical models against field volatility data sets and fully exploring the volatility parametric response surface. Field studies quantifying atmospheric flux for soil fumigants 1,3-dichloropropene and chloropicrin are validated against the United States Department of Agriculture (USDA Salinity Laboratory) soil physics model CHAIN_2D that was modified specifically for agronomic uses of soil fumigants. Comparison between model predictions and field observations for six unique field trials in five different states indicate that CHAIN_2D effectively captures the magnitude and duration of fumigant emission from soil observed experimentally with r ²  ∼ 0.14–0.96 (avg. 0.66) for peak emission, and r ²  ∼ 0.76–1.0 (avg. 0.91) for cumulative emissions (r equals the Pearson moment correlation coefficient). Correct prediction magnitudes suggest that CHAIN_2D is a useful tool for extrapolating flux predictions to diverse scenarios not addressed by field trials. Examples of mitigation strategies such as the use of agricultural films (tarps), increased soil injection depth, and management of soil water content, under near semi-infinite parameter combinations of soil, meteorological, and agronomic properties are discussed.     

Spatial variability of fallout-137Cs in the soil of a cultivated field

To obtain information on the spatial variability of fallout ¹³⁷Cs in the soil of a small area, the activity concentration of this radionuclide was determined in 100 soil samples, taken along the two diagonals of a cultivated field (150×100 m). The results show that the spatial distribution of ¹³⁷Cs in this field is at random. The frequency distribution of the values is skewed to the right, but not log-normal. The median activity concentration observed was 7.45 Bq kg^-1 dry soil, the values ranging from 4.8–15 Bq kg^-1. The spatial variability, as characterized by the relative decile deviation was 26%. If one tolerates for the establishment of baselines of global fallout ¹³⁷Cs in the soil an error in the mean activity concentration of 10% (20%) at the 95% confidence level, the minimum number of soil samples to be taken can be estimated as 14 (4). The total deposition of ¹³⁷Cs on the soil surface by fallout was determined as 3.3±0.4 kBq m^-2.     

Spatial Variability and Monitoring of Pb Contamination of Farming Soils Affected by Industry

In this study, the relationship between some physico-chemical properties of soils and lead contamination in soil due to emission from industrial operations in Samsun province of Turkey was investigated. The extent of timely contamination was studied by comparing the obtained results with the results of the study conducted in the same region in 1998. An area of 225 km² (15 km × 15 km), which was divided into 1000 × 1000 m grid squares (16 lines in the east and south directions), was selected within the industrial area. The total of 256 grid points was obtained and soil samples were collected from three depths (0–5, 5–15, and 15–30 cm) of each grid center in 2004. The total Pb concentrations of soil samples were determined as 65.84–527.04 μg g⁻¹ at 0–5 cm in depth, 58.50 – 399.54 μg g⁻¹ at 5–15 cm in depth, and 44.65–330.07 μg g⁻¹ at 15–30 cm in depth. DTPA-extractable Pb concentrations of soils were found to be in the range of 1.52–9.03 μg g⁻¹, 0.54–7.09 μg g⁻¹, 0.19–6.13 μg g⁻¹ at 0–5, 5–15, and 15–30 cm depths, respectively. There were significant relationships between both total or DTPA-extractable Pb concentrations and selected physico-chemical properties of soil. According to enrichment factor (EF) values calculated from the total Pb concentrations, 11.3% of the study area (225 km²) was enriched with Pb in high level, but 77% of the area was in significant enrichment level with Pb. The average total and DTPA-extractable Pb concentrations increased as 11 and 13%, respectively in comparison with the results of 1998.