# Comparing Soil Erosion Rates on Terraced and Sloping Cultivated Land in Palestine Using FRN137Cs Trace

Authors: Orwa Houshia; Moncef Benmansour; Lionel Mabit; Emil Fulajtár; Saber Abu-Jabal; Ismail Hroub; Rafat Odeh
Journal: International Journal of Analytical Chemistry
 (2022)

---

## Abstract

Soil erosion is a serious problem in Palestine. It is enhanced mainly by poor farming practices used in upland agricultural areas occupying the Central Highland of Palestine. The objective of this study is to assess the impact of terracing on soil erosion and deposition rates in the Al-Yamoun area (the Northern West Bank) using the fallout radionuclides cesium -137 (FRN137Cs). The FRN 137Cs technique, which has proved its efficiency in estimating erosion rates over the last 50–60 years, was used for the first time in Palestine to measure rates of erosion and deposition. The activity of 137Cs was measured by gamma spectrometry using an HPGe detector. For the reference site, the 137Cs inventories ranged between 2499 and 4086 Bq/m2. The average value of the reference site is 3315 ± 410 Bq/m2, which corresponds to a coefficient of variance of 12%, suggesting that the reference site is well representative for estimating 137Cs fallout. This 137Cs amount is too high for bomb-derived fallout and indicates that a significant part of the deposition is from the Chernobyl accident. The 137Cs inventories at both studied sites (terrace site and foot slope site) are significantly lower than those of the reference site. For the terrace site, the inventories are found between 1707 and 2749 Bq/m2, while for the slope site they are between 1050 and 2617 Bq/m2. The lower 137Cs values at both studied sites than values at the reference site indicate that the entire areas of both study sites are eroded and no depositional activity occurs.


---

## Body

## 1. Introduction

Agriculture is an important cultural tradition vital to the economy of the West Bank, Palestine. Also, agriculture plays a vital role in the region’s future. Palestinians are coping in creative ways [1] every day despite the challenges of establishing a livelihood or any sustainable economic activity based on agriculture. Whether the challenge lies in water scarcity, soil degradation [2], or blockade restrictions, Palestinian families implement innovative and collaborative approaches to living off the land [3]. However, the land suitable for agriculture in Palestine is shrinking due to soil erosion, land degradation, and water erosion, which in turn affect soil properties and functions [4–7] and cause the loss of soil, which is not renewable on the human timescale. Consequently, these factors represent a major threat to the soil and water resources the country needs to ensure sustainable agricultural production.Within this context, and in order to provide a comprehensive assessment of the magnitude of these problems and to support the selection of effective soil conservation measures, quantitative data on the extent and rates of soil erosion under various agroecological conditions and land use systems were needed [8]. Integrated land and water management practices improve agricultural production and enhance soil productivity and its resilience against desertification and other impacts of climate change and variability. Radionuclide and stable isotopic techniques can be used to study soil erosion and land degradation problems [9]. Data on soil erosion acquired by conventional methods takes several years to obtain and is time consuming. Furthermore, it is labor-intensive, and these conventional methods do not provide for the spatial distribution of soil erosion. On the other hand, the 137Cs isotope tracer fills the deficits of the conventional methods and facilitates the investigation of soil erosion on a much better timescale, so there is no need for costly, labor-intensive long-term monitoring.


## 2. Methodology

### 2.1. Study Area

The study area is located in the Central Highland in the Jenin Province, which receives annual precipitation of about 300 mm (Figure1). The sampling sites were selected to acquire soil samples within the Jenin Province. The reference site at an elevation of 293 m above sea level was covered with shrub vegetation and had a slope of about 2%. The other two sites were labeled Terrace Site, which at an elevation of 246 m above sea level, cultivated land on a terrace with a slope of 2%, and foot slope site, which at an elevation of 140 m above sea level, cultivated land on a foot slope that was not terraced and had a slope of 3%. The difference in slope inclination is very small, but the important difference is the slope length (field length), which is very short on the terrace (10–15 m) and several folds longer at the foot slope field (60 m).Figure 1
Location of the study area. The inset picture is an aerial view of the actual locations.


### 2.2. Field Sampling

The sampling design is as follows:


#### 2.2.1. Reference Site

Grid samplingof 1.5 m (24 samples) (3 × 8), bulk samples, 30 m depth, and a distance between two points is of 1.5 m.


#### 2.2.2. Terrace site

Sampling along 3 transects containing 8 points in each transect (24 soil samples), bulk samples, 30 cm depth, and a distance of 1.5 m between two points.


#### 2.2.3. Foot slope site

Sampling along 3 transects containing 8 points in each transect (24 soil samples), bulk samples, 30 cm depth, and a distance of 1.5 m between two points.


## 2.1. Study Area

The study area is located in the Central Highland in the Jenin Province, which receives annual precipitation of about 300 mm (Figure1). The sampling sites were selected to acquire soil samples within the Jenin Province. The reference site at an elevation of 293 m above sea level was covered with shrub vegetation and had a slope of about 2%. The other two sites were labeled Terrace Site, which at an elevation of 246 m above sea level, cultivated land on a terrace with a slope of 2%, and foot slope site, which at an elevation of 140 m above sea level, cultivated land on a foot slope that was not terraced and had a slope of 3%. The difference in slope inclination is very small, but the important difference is the slope length (field length), which is very short on the terrace (10–15 m) and several folds longer at the foot slope field (60 m).Figure 1
Location of the study area. The inset picture is an aerial view of the actual locations.


## 2.2. Field Sampling

The sampling design is as follows:


### 2.2.1. Reference Site

Grid samplingof 1.5 m (24 samples) (3 × 8), bulk samples, 30 m depth, and a distance between two points is of 1.5 m.


### 2.2.2. Terrace site

Sampling along 3 transects containing 8 points in each transect (24 soil samples), bulk samples, 30 cm depth, and a distance of 1.5 m between two points.


### 2.2.3. Foot slope site

Sampling along 3 transects containing 8 points in each transect (24 soil samples), bulk samples, 30 cm depth, and a distance of 1.5 m between two points.


## 2.2.1. Reference Site

Grid samplingof 1.5 m (24 samples) (3 × 8), bulk samples, 30 m depth, and a distance between two points is of 1.5 m.


## 2.2.2. Terrace site

Sampling along 3 transects containing 8 points in each transect (24 soil samples), bulk samples, 30 cm depth, and a distance of 1.5 m between two points.


## 2.2.3. Foot slope site

Sampling along 3 transects containing 8 points in each transect (24 soil samples), bulk samples, 30 cm depth, and a distance of 1.5 m between two points.


## 3. Sample Preparation and137Cs Analysis

Physical preparation of the soil samples was performed at the Palestine National Agricultural Center (NARC) including drying and light grinding. Then, samples were oven dried at 105 degrees Celsius, grinded, sieved (<2 mm), and homogenized. Radionuclide analysis of137Cs was performed by gamma spectrometry using an HPGe detector (45 efficiencies) by the Centre National de l’Energie des Sciences et des Technique Nucleaires (CNESTEN), Morocco. Tennelec/Nucleus HPGe (184 cc) planar-type coaxial intrinsic germanium detector was used. The 137Cs activity was measured by its gamma emission at 662 keV. Background counts were made at regular intervals to ensure the maintenance of low background characteristics. Each soil sample was counted for 20000 s. The Cs inventory A (Bq·m−2) was calculated using the following equation:(1)A=CMS,where C is the 137Cs activity concentration of the sample (in Bq·kg−1), M is the total dry mass of the collected soil core (in kg) and S is the cross-section of the sampling tube (in m2)


## 4. Results and Discussion

Numerous models were established to analyze soil erosion [10]. The main postulates and necessities of the 137Cs approaches have been reported in several papers [11]. The calculation of soil erosion and deposition rates at the undisturbed site is based on 137Cs inventories and the 137Cs depth distribution. Two methods are used. The simpler method considers the fixed 137Cs fallout input and stable 137Cs profile distribution. The depth distribution over the soil profile is mathematically described [12]. If this distribution is considered by a simple mathematical function, then the soil loss can be estimated by the proportion of the 137Cs reference value removed at the examined site [13]. This method is used for the Profile Distribution Model (PDM). A more widespread method takes into consideration the time-variant processes of: (1) the 137Cs fallout and (2) the gradual postdepositional redistribution of 137Cs within the soil profile, which is caused predominantly by bioturbation and several other processes. This approach is used by the Diffusion and Migration Model (DMM) [14]. The calculation of soil loss in cultivated land can be based on two theoretical concepts. The first one called, the proportional concept, is very simple and it presumes that the removal of soil and 137Cs are directly proportional. Models based on this concept are called the Proportional Model (PM). More complex approaches involve a mass balance concept, which considers the temporal dynamics of 137Cs inputs and outputs resulting in the time-variant concentration of 137Cs in soil [12]. These changes in concentration affect considerably the relation between the 137Cs loss and soil loss caused by erosion. The 137Cs concentration in soil is controlled by several processes [13, 15–17]. Most important are (1) time-variable 137Cs fallout, (2) radioactive decay of 137Cs (i.e., 30.17 years), (3) removal of freshly deposited 137Cs by erosion prior to its incorporation into the plough horizon by tillage, and (4) incorporation of subsoil material free of 137Cs or having low 137Cs content into the eroded ploughed horizon by tillage [18].The Mass Balance Models (MBM) were developed in the mid-1980s. Different models use different approaches to handle the137Cs time-variant concentration in soil and consider some but not all processes and factors controlling it [12, 19, 20]. Moreover, three Mass Balance Models that were initially developed by Walling and He [10] are used, e.g., Mass Balance Model 1 (MBM1), Mass Balance Model 2 (MBM2), and Mass Balance Model 3 (MBM3).In order to convert inventories or areal activities (in Bq/m2) into soil erosion or deposition rates (in t/ha/yr), conversion models were used. First, the Proportional Model (PM) was used to have preliminary results about soil erosion, in accordance with equations (2), (3), and (4) obtained from references [12, 21].The basic equation (2) of the Proportional Model [10] can, therefore, be represented as follows:(2)Y=10BdX100TP,where X is the percentage reduction in total 137Cs inventory (defined as (Aref -A)/A ∗ 100), d is the depth of the plough or cultivation layer (m), B is the bulk density of the soil (kg·m−3), T the time elapsed since the start of 137 Cs accumulation (yr), Aref, represents the local reference inventory (Bq·m−2), and A is the total inventory measured at the sampling point (Bq m−2). P represents the particle size correction factor for erosion. This model does not require many parameters, but it is not realistic as he neglected some processes or phenomena related to 137Cs behavior and erosion processes in cultivated sites.In addition, the Chernobyl contribution was not included when applying this model. In this case, Mass Balance Models were used (equation (3)) and, more specifically, the Mass Balance Model 2, which can describe suitably the real situation and therefore was also applied (equation (4)). It has to be noted that the file for the annual deposition of 137Cs was modified to take the Chernobyl contribution into account. The results given by both models for both sites 1 and 2 are reported in Table1 and Table 2.(3)Y=10dB1−1−x1001/t−1963,where t represents the time since 1963, B is the bulk density of the soil (kg·m−3), and d is the depth of the plough or cultivation layer (m). X is the percentage reduction in total 137 Cs inventory (defined as (Aref -A)/A ∗ 100) (4)dAtdt=1−ΓIt−λ+PRdAt,where A (t) = cumulative 137Cs activity per unit area (Bq/m2); R = erosion rate (kg/m2. yr) d = cumulative mass depth, representing the average plough depth (kg/m−2) λ = decay constant for 137Cs (yr−1); I (t) = annual 137Cs deposition flux (Bq/m2 yr); Γ = percentage of the freshly deposited 137Cs fallout removed by erosion before being mixed into the plough layer; P = particle size correction factor.Table 1
Soil erosion rates in t/ha/yr associated with site 1 for different points of the transects using the Proportional Model (PM) and Mass Balance Model 2 (MBM2).
Foot slope site 1Distance from the top (m)PMMBM2T1T2T3T1T2T309.618.88.710.926.79.61.511.17.19.513.07.710.737.415.67.78.120.38.44.511.58.17.313.68.97.9617.48.19.323.98.910.57.516.75.48.222.45.79.095.910.06.66.211.47.110.511.17.29.512.97.810.7Table 2
Soil erosion rates t/ha/yr associated with site 2 for different points of the transects using the Proportional Model (PM) and Mass Balance Model 2 (MBM2).
Terraced site 2Distance from the top (m)PMMBM2T1T2T3T1T2T3022.923.510.040.042.412.01.522.316.59.438.223.311.2321.810.912.036.613.315.24.523.39.425.341.511.149.4622.712.924.939.616.647.87.522.611.425.439.114.249.8923.410.024.541.812.146.110.523.010.023.640.612.042.7The inventories of the samples were calculated from the137Cs concentration (in Bq/kg), the bulk density of the sample, and the depth of the core. An average bulk density was calculated for each type of site. There are about 1160 kg/m3, 1240 kg/m3, and 1190 kg/m3 for a reference site, site 1, and site 2, respectively.The results of all inventories (Bq/m2) are given in Table 3. For the reference site, the inventories ranged between 2927 and 4086 Bq/m2. The average value of the reference site is 3315 ± 410 Bq/m2, which corresponds to a coefficient of variance of about 12%. This coefficient of variance can be considered as low, confirming the low variability and low disturbance of the selected reference site. The high inventory associated with the reference site suggests that the fallout 137Cs deposited in this region is the result of both old nuclear weapon tests conducted in the 1950s–1960s and the Chernobyl disaster that occurred in 1986. Taking into account the average rainfall in this region (Jenin Province) of about 300 mm, the inventory estimated by prediction using conversion models is about 897 Bq/m2. It means that the contribution of Chernobyl is about 2418 Bq/m2, which represents around 73% of the total reference inventory. Concerning study site 1 and site 2, the inventories are significantly lower than those of the reference site (see Table 1), especially for site 2, indicating that these sites are eroded. Figure 2 shows the distribution of 137Cs inventory along 3 transects for each site and the mean reference inventory value. For site 1, inventories are found between 1707 and 2749 Bq/m2, while for the site 2, inventories are found between 1050 and 26117 Bq/m2. In comparison, all the inventories obtained are lower than the reference value (Figure 1), indicating that there is no depositional area at both sites.Table 3
Inventories (B/m2) associated with all points collected at reference and study sites.
Reference siteFoot slope site 1Terrace site 2Sample codeInventorySample codeInventorySample codeInventoryR1013804PE1012492SD1011278R1024086PE1022366SD1021328R1033543PE1032678SD1031374R1043511PE1042332SD1041239R1053564PE1051823SD1051289R1063675PE1061886SD1061303R1073790PE1072812SD1071232R1082965PE1082370SD1081264R1093800PE1091707SD1092617R1102864PE1102704SD1101217R1112791PE1111983SD1111846R1123003PE1122623SD1122345R1133250PE1132623SD1132481R1143069PE1142850SD1142167R1152979PE1152463SD1152296R1162944PE1162701SD1162420R1172499PE1172574SD1172428R1182840PE1182504SD1182478R1193108PE1192656SD1192242R1203135PE1202693SD1201057R1213616PE1212522SD1211092R1223014PE1222615SD1221050R1233369PE1232749SD1231128R1242927PE1242504SD1241210Figure 2
137Cs inventories (Bq/m2) along 3 transects associated with site 1 (a) and site 2 (b).
(a)(b)As previously reported by Nouira et al. [19], the Proportional Model does not take into account the dilution of 137 Cs concentrations in the soil within the cultivated layer, resulting from the assimilation of soil from below the originally cultivated depth due to surface lowering by erosion; thus, the results obtained by this model in the present study likely underestimate the actual rates of soil loss. The Mass Balance approach that takes into account the effect of mixing the subsoil containing no 137Cs from below the plough depth into the plough layer gave higher values for soil loss rate, as is observed in some other reports [19, 22].


## 5. Conclusions

137Cs distribution at the sites confirmed that the 137Cs tracking procedure can be effectively used to assess soil erosion and deposition rates. The application of the 137Cs tracking system enabled us to quantify the extent of soil erosion in the Northern Palestine region. For the reference site, the inventories ranged between 2927 and 4086 Bq/m2. The average value of the reference site is 3315 ± 410 Bq/m2, which corresponds to a coefficient of variance of about 12%, suggesting that the fallout 137Cs deposited in this region is a result of both old nuclear weapon tests conducted in 1950s–1960s and the Chernobyl disaster that occurred in 1986.Finally, the137 Cs tracking technique is a useful erosion model for the Northern Palestine region. This data can be used to refine erosion control guidelines influencing the selection of a vast array of management practices by an agricultural producer, such as tillage, fertilizer application, crop rotation, plant population, and structures. Thus, the future work, which will be carried out in phase–II of the project, highly recommends the selection of more locations in Gaza and the West Bank in order to generate a soil erosion map of Palestine. This result will also be presented to legislators and decision makers within the government to implement appropriate soil erosion solutions.



---

*Source: 2933661-2022-09-30.xml*