Particle size distribution in soil is an important soil property that characterizes water holding capacity, water movement and other physical and chemical properties of soil. The importance of soil particles emerged from the fact that it can determine the fate of many organic and inorganic chemicals in soil environment    . Clay fraction in soil plays an important role in the retention of organic and inorganic cations in soil layers  . A number of investigations   , revealed the ability of fine soil particles to retain higher concentrations of heavy metal than larger size fractions because the finer particles have a higher surface area per unit weight. Furthermore, the finest soil particles (<10 μm) may be absorbed into human tissues  or re-suspended in the atmosphere by different mechanisms  , creating environmental problem in cities  . So far, soil particle size can create a huge diversity of micro-pores that govern the movement of water and air in soil profile   and have import roles in the ecosystem  . Large particle size creates larger pore-size, and consequently fast movement of air water molecules. On the hands, smaller particle size creates smaller pore-size or capillary pore size that retains water and increase water holding capacity of soil. So far, hydraulic conductivity in soil is directly connected with particles size distribution in soil profile. Several methods have been tested for determination of hydraulic conductivity; this included integrated global soil databases with land use/land cover databases and saturated hydraulic conductivity  , pedotransfer functions  , rainfall-induced seals  and influence of biochar  . These investigations provided valuable information that was applied in agricultural soils around the world. The abovementioned studies were limited to agricultural or forest lands. Particle size distribution, and hydraulic conductivity in coastal non-agricultural lands were poorly investigated or remains in the primary stages of investigation. The novelty of this study emerged from the fact that it targeted coastal not agricultural land and generate useful information to the environmental scientist around the world. The objectives of this investigation were to characterize the particle size distribution and hydraulic conductivity on coastal, non-agricultural lands along the coastal plain in Gaza Strip.
2. Materials and Methods
2.1. Selection of Experimental Sites
Six sites were selected for digging soil profile along the rainfall gradients in Gaza strip. These site are Khan Younis (South, site 1), Bait Lahia (North, site 2), Gaza-Shakhejjlen (sites 3 & 4), and Rafah sites (5 & 6).
1) Khan Younis (Site 1)
It is located in the southern part of Gaza. It has the following GPS information (N, 31˚26'37.56"; E, 34˚21'37.35"). The area is about 2500 m2. It is a sandy soil. It is 400 m far away from Mediterranean Sea. No agricultural activities were running in the site during the past 10 years, but the site contains a lot of wild type weeds, plants and some trees. The surrounding area of the site contains fruit trees such as Gova, palm trees and other desert shrubs. The average rainfall is about 290
mm, and far away from the site some vegetable cultivation. The information of this site was obtained from the Mayer of the city by personal communications.
2) Bait Lahia (Site 2)
Site II is located in Bait Lahia, the northern part of Gaza Strip. It has the following GPS information (N, 31˚33'55.14"; E, 34˚28'16.75"). The site is about 150 m a part from the sea. It is located in coastal area. The eastern part of the site has some agricultural activity with vegetable and greenhouse production. The historical information from the owner indicates that the land did not contain any agricultural waste, municipal or hazardous materials. The annual rainfall is about 433 mm. The information of this site was obtained from the Land Owner Mr M abu Halema by personal communications.
3) Shakhejjlen area (Site 3 & 4)
These sites are located in Shakhejjlen area, 650 meter far away from the sea sure. It has the following GPS information (N, 31˚29'0.89"; E, 34˚24"3.08"). The history of the area is an agricultural land specified with fruit production such as grapes, fig, and little citrus trees. Big trees such as Tamarex are cultivated as boarder line trees among farms.
The annual rain fall in the area is 394 mm. The history of the site indicates no hazardous waste dumpsite in the area in the past 10 year and the boarder trees, Tamarex are old trees. In this site two soil profiles were made, one directly under the stem of Tamarex and the second is 3 meters away from the stem. The information of this site was obtained from ministry of agriculture, PNA, un published results.
4) Rafah (Site 5 & 6)
These sites are located in Rafah area, 3 km far away from the sea, Rafah area is the southern part of Palestine. It has the following GPS information (N, 31˚8'38.74"; E, 34˚15"35.19"). The site is sand dunes with no human activity except grazing. It has variety of wild type vegetation. The annual rain fall in the area is about 236 mm. In this site two soil profiles were made, one directly under the stem and the second is 3 meters away from the stem. Moreover, the important data of all sites are summarized in Table 1 and their locations in Gaza Strip are presented in Figure 1. The information of this site was obtained from municipality of Rafah city, PNA, unpublished results.
2.2. Soil Sampling
Six soil profiles were digged along 45 km in the coastal plain of Gaza Stip. Each soil profile has the following dimension 100 × 200 × 150 cm (Figure 1). Samples were collected from soil profile at both sites at of depths of 0 - 10, 10 - 20, 20 - 30, 30 - 40, 40 - 50, 50 - 70, 70 - 90, 90 - 110, 110 - 130, and 130 - 150 cm.
Soil samples were air dried sieved through 2 mm mish size to remove the gravels and big stones and storage in plastic containers at room temperature for further analysis.
Table 1. Locations, distance, rainfall and GPS information of the sites.
Figure 1. Soil profile photos and locations in Gaza Strip map.
2.3. Particle Size Distribution
With modifications to the previously described procedures   , particle size distribution was determined. In this procedure ten g of each sample was transferred to a 250 ml Erlenmeyer flask, 100 ml of distilled water was added to each flask, 25 ml of H2O2 (30%) was added to each flask. The flasks were transferred to a water bath for heating up to 90˚C for 2 days to destroy all the organic matter presents in the samples. After cooling the flasks, 7 ml (1%) of HCl was added to each flask to destroy the CaCO3. After completing the CaCO3 degradation, the samples were transferred to a shaking bottles and 25 ml of Na-pyrophosphate (0.4 molar) was added to each bottle. The bottles were sealed with a plastic device and shacked for 2 h using a rotary shaker to have a homogenized soil suspension. All chemicals used in this section were HPLC grade.
The suspension in the bottle was decanted into a series of sieves with the following diameters (630, 200, 63 and 20 µm) arranged from the top to the bottom. The sieves were mechanically shacked for 10 min under continues water flow. Then each sieve was separated for the series and its content was transferred to a glass container its weight in dry form was previously collected.
The containers were transferred to the oven for an overnight drying at 105˚C. In the next day the containers were weighed. The difference between the dry empty container and dry full container equal the weight of the soil particle separated in the sieves. Fraction % was collected according to Equation (1)  .
2.4. Hydraulic Conductivity Determination
Each site was subdivided into four plots as follows: control plot (CP), and three disturbed plots (DP1, DP2, and DP3). The idea behind disturbing soil plots is to investigate the effect of disturbance on hydraulic conductivity. The differences among the three disturbed plots were in depth of disturbance. The disturbing depth was 3, 5 and 10 cm. These depths were manually achieved using special device of wood containing nails with the above mention lengths. The reason behind selecting these depths were to destroy the salt crust that may be formed in the top 3 and 5 cm due to sea water drift, whereas, and destroy soil large halls due to growth of earth worm.
Hydraulic conductivity in different soils was determined by Minidisk Infiltrometer (Decagon Devices, Inc. 2007-2017)  . On this method the Infiltrometer was placed on the top soil and measured the infiltrated amount of water in certain time (e.g., 30 second) then the measurements were repeated at least five times on the same location with different periods of time from those mentioned above. The time can be different of that mentioned above but it may be better to be attached to a certain time during the measurement.
With a slight modification of previous published work   , calculation of hydraulic conductivity required at least five readings in the same location. The time of the reading is constant and the volume of water may be changed due to hydraulic conductivity. Knowing the infiltrated water volume and diameter of the Infiltrometer enabled us to calculate the depth of water infiltration in soil. Plotting the relationship between the cumulative square root of time vs cumulative infiltration depth provided a polynomial mode of relationship. The slope of the polynomial mode indicates the hydraulic conductivity in the soil.
2.5. Statistical Analysis of the Data
Following to the procedure described in El-Nahhal et al.  , means and standard deviation of data were calculated. t-test was used to determine significant differences among values in different soil profiles. p-value below 0.05 indicate significant difference and tested parameters.
3.1. Soil Profiles
Soil profile photos are shown in Figure 1. Regardless to the differences in brightness of photos, the soil profiles are similar in color soil structure except site 2 (Bait Lahia) the soil seems to be compacted due to previous agricultural activity. So far, the differences in soil profile probably were due to time and weather during photos since the work was not achieved in the same day. The important data of the sites are summarized in Table 1. Nevertheless, the sites are different in rainfall, the highest one is Bait lahia (site 2) 433 mm, and the lowest one is Rafah (sites 5 & 6) 236 mm.
3.2. Particle Size Distributions
Particle size distribution in all sites is shown in Table 2(a) and Table 2(b). It can be seen that the trend of particle size distribution is similar in all depths at any site. For instance the largest particle size, 630 µm, and the smallest particle size, 20 µm, represented the lowest fractions at any soil layer in the soil profile. On the other hand, the particle size 200 µm and 63 µm were the largest fractions at any soil layer in the soil profile.
Nevertheless, the average of each particle size in all layers at all soil profiles was presented in Figure 2. The trend of particle size distribution is similar in all soil profile. Statistical analysis showed significant differences of each particle size. It can be noticed that the fractions of the size 630 µm are similar in all sites whereas the fractions of 200 µm are similar in sites 2, 3 and 4 and significantly different at sites 1 and 5. On the other hands the fraction of particle size 63 µm are similar in sites 1 and 2 and significantly different on the other sites. Furthermore the fraction of 20 µm is not significantly different in all sites except site 5. The data in Figure 2 were marked with different letters (a, b and c), where the statistical analysis indicated significant differences at p-value was less than 0.05.
Furthermore, it can be noticed that all soil profiles (Table 2(a), Table 2(b)) contained low fractions of the largest size 630 mm and fine particles 20 mm. The majority of particle size is 200 mm followed by 63 mm in all soil profiles. These results indicated that the structure of soils in all sites were sandy and nearly similar. Furthermore, statistical analysis of the data in Figure 2 did not detect significant differences on the particle size 630 nm found in all sites whereas significant differences were detected on the fraction of 200 nm from all sites. It is obvious that sites 2, 3 and site 4 are nearly similar but are different from site 1 and site 5. The interesting point here is the similarity on sites 3 and 4 suggesting no effects from the tree on particle size distribution. Moreover, differences were also detected on the fraction of 63 nm but sites 3 and 4 are nearly similar confirming no effects from the tree on particle size distribution.
Table 2. (a) Particles size distribution in sites 1 and 2. The data are presented in %. (b) Particles size distribution in sites 3 - 5. The data are presented in %.
Figure 2. Particle size distribution (mean ± standard derivation) in all soil depths of soil profiles in different sites. Error bars represent standard deviations. Columns have the same letter (a, b, or c) at any particle size are not significantly different at p-value = 0.05.
3.3. Hydraulic Conductivity
Hydraulic conductivity (Kc) of the soil profiles and surrounding areas are shown in Table 3. It can be noticed that the Kc values are different in the same testing site due to the disturbance depth of soil. So far, in site 1, Kc value is the smallest in the control plot soil (CP, soil was not disturbed), whereas the values of Kc in the disturbed plots are several times higher. On the other hand, Kc values are quite similar in site 2.
Furthermore, SH.E DP3 has the highest Kc value in site 3 whereas Rafah DP2 has the lowest value in site 4. These differences reflect the differences in pore size distribution in the tested soil. Nevertheless, the high value of regression coefficient (R2) indicates the precision of the presented results.
The studied land is a coastal non-agricultural land has a wide range of rainfall gradient (290 - 450 mm). The studied sites have no agricultural activity in the surrounding area for the past ten years. In general Gaza strip is considered as a semi-arid zone due to the lack of rainfall during the year. The data in Table 1 clearly shows important information of the tested soil profiles.
The data in Table 2 clearly shows the fractions of particle size in each soil layer of the soil profiles. It is obvious that the fractions of 630 nm and 20 nm are nearly low fractions in all sites, whereas the fraction of 200 nm is the dominant particle size in all sites followed by the fraction 63 nm. Statistical analysis did not detect significant differences in the particles size of 630 nm in all sites whereas significant differences were found in other fractions. These data suggest that the sites were exposed to similar weathering conditions along the coastal plain in Gaza Strip. This explanation agrees with Michel et al.  who revealed the influence of weathering conditions on soil particle size distribution. Additionally, the locations of the soil profiles are in the coastal plain of Gaza strip (Figure 1), where sands can be transferred from the coastal area to the urban area due to water movements and wend circulations, consequently similar particle size would have accumulated in the sites. Moreover, it can be suggested that differences in wind storm intensity may have resulted in different sand flow or creeping rates in the sites resulting in differences in the particle size distribution along the coastal plain. This suggestion is supported by experimental works of Bełdowska  who reported the influence of weather anomalies on metal cycling in the marine coastal zone of the southern Baltic.
An interesting outcome here is that the similarity in the fraction size 630 µm in all sites and similarity in fraction size 20 µm except site 5. On the other hands the fraction 200 µm is similar in sites 2, 3, and 4. This suggests similar weathering conditions in the corresponding sites.
Furthermore, the dominant particle size 200 nm in all sites suggests the formation of large pore size in soil structure consequently water movement in soil would strongly be affected by this fraction. This agrees with El-Nahhal et al.  , who found high water movement in sandy soil with large particle size. In addition, wide range of particle size in soil may have ecological benefits such as enhancing water movement in soil enriching the ground water aquifer, affect sorption and desorption behaviors of steroid hormones  , affect the distribution of heavy metals   . Furthermore, it can be suggested that various particle size fractions in soil profile may create specific porous size that suit the colonization and the activity of soil microorganisms. This suggestion is in accordance with Hemkemeyer et al.,  who demonstrated that the mineral composition and the particle size fractions in soil have specific and different selective effects on the abundance and diversity of Bacteria, Archaea and Fungi in natural soils with complex mineral compositions.
In addition, the dominancy fraction size 200 nm all sites followed by the fraction 63 nm suggests the formation of similar pore size in soil that can affect the hydraulic conductivity in soil. This suggestion is supported by the values of Kc presented in Table 3, regardless to the fact that Kc values in Table 3 are not equal.
However, the variations in Kc values in Table 3 can be explained by the differences in disturbing depth and differences in rainfall during the year (Table 1). This explanation is in accordance with Armenise et al.  who investigated influence of rainfall in the hydraulic conductivity and found two different trends in soil porosity vs. depth. Moreover, high rainfall may result in leaching and washing of salts from the top soil to deeper depths. This increases the free movement of water molecules in soil profile. This is in accordance with El-Nahhal and Safi   who revealed the influence of salinity in the free water molecules to dissolve organic cations.
Furthermore, land use and land cover may affect the hydraulic conductivity in soils as seen in (Table 2) in the control plot soils and disturbed plot soils. This assumption agrees with Wong et al.  who reported that biochar application increased the saturated hydraulic conductivity of compacted clay, and the increment was positively correlated to the biochar percentage. Additionally, the similarity in Kc values in site 2, probably due to the compaction of soil as seen in photos in Figure 1 and similarity in particle size distribution (Figure 2).
5. Comparisons with Previous Published Work
A comparison between the present study and those from the literature is presented in Table 4. It is obvious that the present study deals with a coastal
Table 3. Hydraulic conductivity and regression parameters in different soils.
CP, DP1-3, are control plot soil and disturbed plot soil. Y and X are cumulative infiltration depth (cm) and square root of time respectively.
Table 4. Comparing parameters with other works.
non-agricultural land and determined important soil parameters that were useful of agricultural and non-agricultural lands. Published work     dealt with agricultural, forest and wetlands and determined few parameters of soils such salinity. The previous published works were localized to certain small study site whereas the present study digged several soil profiles along a distance of 45 km. Moreover, the present study provided data on particles size distribution, hydraulic conductivity of different soil plots whereas other studies   evaluated carbonate concentrations in top soil layers. Additionally, the previous studies partially determined the salinity in soil, whereas the present study focused here on particle size and hydraulic conductivity. It can be concluded that the present study integrated the picture by determining the particle size and hydraulic conductivity in non-agricultural land. So that this study is opining a new window of investigation toward land reclamation and land use.
6. Environmental Relevance
Environmental relevance of this work emerged from the fact that particles size in soil fraction has specific surface area that plays an important role in the chemical fate of organic and inorganic molecules. For instance, the smallest particle size the largest surface area that can absorb organic cations    and create a hydrophobic zone  that can absorb organic contaminants from water. Moreover, diversity of pore size may be created and can be a suitable habitat for soil cyanobacteria that can biodegrade organic molecules and remediate soil pollution   . Recent development  in wastewater treatment used different particle size of sand for optimizing efficacy of sand filter.
7. Significance of This Work
Significance of this work emerged from the fact that it is the first time that coastal nonagricultural land being investigated for particle size distribution and hydraulic conductivity measurements in Palestine and surrounding countries. The unique point of this work is that the data will be available for the scientific community in the world for comparison or following technical procedure. Moreover, the information of this study can be of great importance to the researchers in desert lands, lands under reclamation agricultural land and/or wet lands. Furthermore, the new theory that will be arrived from this work is that soils are different and sensitive to different environmental factors that may affect the soil properties.
The rational of this work emerged from the fact that it is the first time along the coastal plain in Gaza Strip that soil profiles were digged in non-agricultural land to determine particle size distribution and hydraulic conductivity. Soil particles size analysis demonstrated that the dominancy of fraction 200 µm in all sites followed by the fraction 63 nm. The hydraulic conductivity was affected by disturbance depth of the tested soil.
JS designed the experimental work and revised the manuscript. YE did the laboratory experiments and wrote the manuscript. MS did the field work and data collection and statistically analyzed the experimental data.
The work was funded by Dousche Forshoung Gemanshaft (DFG) Germany, grant DFG: Fe 218 13/1 Special thanks go to Prof Dr Peter Felxhenengssen, at Gessen University.