The aim of this study was to determine certain basic properties of
soils in the Batumi delta (southwestern Georgia) to determine the
relationships of studied properties and to identify differences with
regards to these properties between different sampling sites in the delta
that were selected based on the delta morphology. In this context, a total
of 125 soil samples were collected from five different sampling sites, and
the clay, silt and sand content of the samples were determined along with
their mean weight diameter (MWD) values, aggregate stability (AS) values,
amount of water retained under
Soil is a key part of the Earth system as a control of the hydrological, erosional, biochemical and biological cycles (Brevik et al., 2015; Decock et al., 2015; Smith et al., 2015). Soil contributes to basic human needs like food, clean water and clean air, and it is a major carrier for biodiversity (Keesstra et al., 2016), a consequence of soils forming at the intersection of the atmosphere, biosphere, and lithosphere (Brevik et al., 2015). Thus, soils covering the Earth vary even on small scales in terms of their properties.
Deltas, which are one of the best examples of alluvial soils, are geomorphological structures that result from the deposition and accumulation of alluvial material at estuaries (Erinç, 2001). Owing to the wetlands and high biodiversity they harbor and to their high potential for agricultural production, deltas are areas of great ecological and agricultural significance (Imentai et al., 2015; Khai and Yabe, 2015; Gillison et al., 2016). Depending on the time of accumulation of the alluvial materials, the delta morphology consists of three distinct areas, which, starting from the coast, are the active delta area, the lower deltaic plain and the upper deltaic plain (Erinç, 2001). These parts of deltas show differences in terms of soil properties (Søvik and Aagaard, 2003; Unverricht et al., 2013).
In deltas, the fact that different sections have formed at different times also leads to differences in the way these areas are utilized. The active delta section, which is the youngest section of deltas, generally lacks any notable agricultural production activity, while the upper deltaic plains, which are older in terms of their time of formation, are areas of extensive agricultural activity. The different agricultural use of different delta areas also results in the differences between their soil properties. Numerous researchers have reported that differences in the utilization of delta areas contribute to considerable differences in soil properties (Huang et al., 2012; Deng et al., 2016; Hernández et al., 2016; Madhavan et al., 2016). Conducting agricultural production in delta soils first requires reclamation activities (Krishnamoorthy et al., 2014). It is known that reclamation activities such as plowing, irrigation and fertilization cause changes in the physical (Jiao et al., 2014; Li et al., 2014) and chemical properties (Cui et al., 2012; Han et al., 2014; Jiao et al., 2014) of the soil.
The Batumi delta, which is fed with water and materials carried by the Çoruh and Adjaristsqali rivers, is one of Georgia's most important agricultural production areas. The subtropical climate of the region has promoted a diversity of agricultural production activities. Although the Batumi delta is a very important area of agricultural production, there are currently no studies that have investigated the properties of its soils. The aim of this study was (i) to determine the general properties of soils in the Batumi delta, (ii) to determine the relationships of studied properties and (iii) to compare these properties between different delta areas that are considered to be morphologically different.
Morphological position, land use type and location of sampling sites.
Location of the Batumi delta and digital elevation model (DEM) of study area.
Location of the sampling points in the Batumi delta (L1, L2, L3, L4 and L5 represent the sampling sites. L1 is located in the active delta, L2 and L3 in lower deltaic plain, and L4 and L5 in upper deltaic plain).
The presented study was conducted in the Batumi delta, located in the
south-west part of Georgia (Fig. 1). The delta covers an area
of approximately 3900 ha. The delta is flat with an elevation ranging from 2
to 50 m above the sea level (Fig. 1). The main form of land use in the
delta is agriculture, which is especially dominant in the lower deltaic
plain and upper deltaic plain, while the active deltaic plain is covered by
pastures. The most commonly grown crops are vegetables and tangerines. The
investigated area is characterized by a humid subtropical climate, with an
annual average temperature of 14.4
The Batumi delta was divided into five different sites based on morphological differences. The first site consisted of the active delta area (L1), while the second and third sites comprised the lower deltaic plain (L2 and L3), and the fourth and fifth sites comprised the upper deltaic plain (L4 and L5) (Fig. 2). Information on these sites is provided in Table 1. In each one of these sites, soils samples were collected at 25 different points from the surface layer (0–20 cm).
Soil samples were air-dried and passed through a 2 mm sieve. Particle size
distribution was determined by the hydrometer method (Gee and
Bauder, 1986), while aggregate stability (AS) for 0.25–0.5, 0.50–1.00 and
1.00–2.00 mm size aggregates was determined with the Yoder wet-sieving
method (Kemper and Rosenau, 1986). Soil organic matter (OM) content was
determined by the Walkley–Black method (Schnitzer, 1982). Mean
weight diameter (MWD) was calculated following Eq. (1) (Van Bavel,
1950):
Analysis was performed in three replications for each soil sample. Descriptive statistics, including averages, standard deviation, minimum and maximum values, and the coefficient of variation were determined for all the studied soil properties. The correlation analysis was performed to assess the relationships between variables; the analysis of variance was applied to determine the differences between the sampling sites. JMP 5.0 package software was used to conduct statistical analyses (JMP, 2007).
Summary statistics of soil properties.
MWD, mean weight diameter; OM, organic matter content; AS
Table 2 shows the descriptive statistic results regarding the soil properties examined in study. The clay content of the soils within the study area varied between 1.47 and 24.21 %, with a mean level of 8.91 %; in this respect, the clay content constituted the lowest fraction in the soils. The silt content of the soils varied between 4.90 and 53.41 %, with a mean level of 29.13 %. The sand content of the soils varied between 36.04 and 92.11 %, with a mean level of 92.11 %. Sand content thus constituted the highest fraction in the soils within the study area. The MWD of the aggregates within the study area varied between 0.26 and 1.02 mm, while the average MWD was calculated as 0.64 mm.
The correlation coefficients between field capacity and other soil
properties in sampling areas (L1, L2, L3, L4 and L5 are the sampling sites;
see Fig. 2).
The correlation coefficients between wilting point and other
properties in sampling areas (L1, L2, L3, L4 and L5 are the sampling sites;
see Fig. 2).
The AS of the studied soils was separately calculated
for each class of aggregate size. For the largest aggregate size class
(1.00–2.00 mm, AS
The moisture content retained by the soils at
Correlation matrix of investigated soil properties.
Values in bold are different from 0 with a significance level alpha
The OM content of the soils varied between 0.37 and 5.63 %, with a mean level of 2.53 %. It was determined that this level of OM content is relatively high for soils belonging to the sandy loam texture class (Marchetti et al., 2012).
Among the different properties that were examined, the lowest variation
coefficient was calculated for the AS
The results for the correlation analysis between the examined soil characteristics are shown in Table 3. Based on the analysis results, positive correlation was identified between MWD and the clay and OM of soils, while negative correlation was identified between MWD and sand content. A positive correlation was identified as well between OM and AS. It was also observed that the effect of OM on AS values decreased in parallel to the decrease in aggregate size.
Correlation analysis results evaluating the relationship between the soil water constants and the soil characteristics in question indicated statistically significant relationships between all the soil characteristics on the one hand and the FC and WP on the other.
These relationships vary among sampling areas both in direction and strength. It was determined that the clay, silt and sand contents were more effective in the water retention in L1, where the OM was lower. However, in the areas with higher OM and aggregation, the effect of the grain size distribution decreased and the effect of the OM and the MWD, albeit significantly, increased on FC and WP. In addition, evaluation of the correlation coefficients showed that the greatest effect of MWD was on FC and WP.
Correlation coefficients between soil properties and factors.
MWD, mean weight diameter; OM, organic matter content; AS
According to the factor analysis results, factors 1 and 2 explain 45.56
and 24.63 % of the variance in the soil properties, respectively (Fig. 5). Among the soil properties, clay content (
Summary of ANOVA for soil properties.
L1, L2, L3, L4 and L5 indicate sampling sites. Means with capital
letters indicate significant difference of 0.01. MWD,
mean weight diameter; OM, organic matter content; AS
The difference between the sampling sites with respect to their clay
contents was found to be statistically significant (
The MWD of the soils was also found to vary significantly depending on the
deposit area (
AS values for the different aggregate classes (AS
The level of variation between FC of the sampling sites was found to be
statistically significant (
Factor analysis of studied soil properties in the Batumi delta (MWD,
mean weight diameter; OM, organic matter content; AS
Variance analysis results showed that the variation in the OM of the
sampling sites was statistically significant (
It is observed that the information about the grain size distribution of the delta soil varied considerably in previous studies. The main cause of these differences was due to the difference in the parent material forming the delta soils, the formation conditions and the formation time. However, it was reported that delta soils predominantly had a higher sand content than silt and clay content (Elbasiouny et al., 2014; Li et al., 2015; He et al., 2017), which is similar to our study results. Although sampling areas had the same topography and the same climatic conditions, differences in land management and formation times resulted in high coefficient of variation of the studied properties. Similar to our findings, researchers reported that alluvial soil properties varied at significant levels resulting in higher standard deviation and coefficient of variation (Iqbal et al., 2005). In addition Saldana et al. (1998) reported that coefficient of variation for clay content was higher than silt and sand content, and for pH lower than other properties.
Assessments performed with the aim of determining soil aggregation have
shown that soil fractions and OM content have varying effects on
the MWD and AS. It was
determined that clay, sand and OM content have a
significant effect on the MWD of soils, while only OM has a significant
effect on AS. It is known that aggregation occurs as a result of the
reorganization, flocculation and binding soil particles
(Duiker et al., 2003) and that this process
is also mediated by the soil organic carbon, biota, ionic bridges, clay and
carbonates (Bronick and Lal, 2005). According to Tisdall and
Oades (1982), aggregates form sequentially based on a certain hierarchy,
with different binding mechanisms being involved at every stage. In this
hierarchy, clay, multivalent ions and organic materials congregate to form
micro-aggregates (< 250
Soil particle size distribution, soil OM content, and bulk
density are generally recognized as the most important factors affecting
soil water retention (Liao et al., 2011; Botula et al., 2012; Yang et al.,
2014). In this context, researchers reported that an increase in the clay
and OM content also leads to an increase in water retention
because of adsorptive effects of clay content, and affinity to water and
influence on soil structure of OM (Hillel, 1971; Kirkham, 2004;
Yang et al., 2014). In addition researchers stated that soil structure is
also related to soil water retention (Pachepsky et al., 2006; Nguyen et al.,
2015b; Hong et al., 2013; Khlosi et al., 2013). Our study findings are in
agreement with these general observations, with the increase in clay and OM
resulting in an increase in the amount of water retention at pressures of
It was observed that the clay content of the soil gradually increased moving from the active delta area towards the upper deltaic plain, while the sand content gradually decreased in the same direction. It is believed that there are two main underlying reasons for this variation: the first is the usage of the relevant areas, while the second is the morphological position of the sampling sites. Over the course of a delta's morphological development, new material is deposited in the active delta area, while the more inland areas will be converted into areas of agricultural production. It is known that the physical properties of coastal soils – such as soil particle size, aggregate and soil moisture – will eventually change as a result of reclamation actions such as plowing, irrigation and fertilization (Li et al., 2014). Similarly to our study results, other researchers also report that reclamation activities lead to a decrease in particle size (Li et al., 2012; Sun et al., 2011). It is believed that seasonal variations in the flow regime of rivers that feed the delta especially cause clay-sized particles in the active delta section to be carried away. Studies on alluvial plains also report that areas closer to the riverbed have soils with higher sand content, while the soil sand content gradually decreases at increasing distances from the riverbed (Scott, 2000; Turgut and Öztaş, 2012). It is also believed that areas of deltas closer to the sea are more exposed to wave movement, which cause clay particles to be carried away from the delta by waves. In parallel with our study results, Yu et al. (2015) determined in their study on delta soils that the sand content was higher in soils closer to the sea and gradually decreased at increasing distances from the sea.
A prerequisite for aggregation is that the clay must be flocculated (Hillel, 2003); in the field, adjacent aggregates often tend to adhere to one another, though certainly not as strongly as do the particles within each aggregate. For this reason, it was expected that the lowest clay content would be observed in soils from L1 and that the MWD values for these soils would be the lowest. It was also expected that an increase in soil clay content would be associated with a higher MWD value. Similarly to our study results, other researchers have also reported that an increase in clay content for sandy soils led to an increase in MWD values (Chrenková et al., 2014; Wang et al., 2016).
AS is defined as the resistance exhibited by aggregate against mechanical shearing forces and the dispersing effect of water (Scott, 2000). Factors which affect the formation of stable aggregates within the soil include the wetting and drying cycle, the freezing and thawing cycle, and the presence of clay, Fe and/or Al oxides and OM in the soil (Abid and Lal, 2008; Karaman et al., 2007). AS values tend to increase in parallel with increasing OM content (Bravo-Garza et al., 2010; Joseph Oyedele et al., 1999; Plante and McGill, 2002a, b; Soinne et al., 2016; Tisdall and Oades, 1982). Soil cultivation activities, on the other hand, are known to reduce AS values (Jury and Horton, 2004; Plante and McGill, 2002a; Scott, 2000). The results of our study were generally in agreement with the literature, in that higher OM was found to be associated with higher AS, while areas subject to cultivation activities generally had lower AS. Other researchers similarly report a relationship between higher OM content and higher AS (Bravo-Garza et al., 2010; Soinne et al., 2016; Turgut et al., 2015), and the relationship between higher sand content and lower AS (Chrenková et al., 2014). Various researchers also report that soil cultivation has the effect of reducing AS (Plante and McGill, 2002b; Shu et al., 2015; Soinne et al., 2016).
The soil texture and structure affects the pore size and geometry of soils, which in turn affects their moisture constants. Owing to their high water retention capacity, and their ability to improve the structural properties of soils, OM also has a positive effect on the soil moisture values (Karaman et al., 2007). For this reason, areas with high MWD, clay content and OM content values are also expected to have high field capacity and wilting point moisture contents. In parallel with our study's findings, other researchers have also reported higher field capacity and wilting point moisture content with increasing OM (Bauer and Black, 1992; Hudson, 1994; Peake et al., 2014; Nguyen et al., 2015a) and clay content (Hudson, 1994; Rawls et al., 2003).
There are two main reasons for low OM content in coastal areas: the first is high salt content (Yuan et al., 2007) and the second is wave erosion (Wong et al., 2008). It is also known that vegetation has a direct effect on the soil OM content. Soil organic content is constituted of plant and animal remains (Baldock and Nelson, 2000; Karaman et al., 2007) and is directly associated with the plant biomass (Sollins et al., 1996; Jaiarree et al., 2011; Novara et al., 2013). It is believed in the present study that the lowest OM in L1, which is close to sea, is caused by coastal erosion, high salinity, and the sparsely distributed vegetation in this area. Vegetation in L2 consisted predominantly of bushes, which led to a higher OM in the area. In the other areas (L3, L4 and L5), it was observed that due to agricultural activities that accelerate the mineralization of OM, OM was lower compared to L2. Similarly to our study results, other researchers have also reported that the soil OM content of natural pastures and prairies was higher than that of areas used for cultivation (Cates et al., 2016; Gajić, 2013; Kodešová et al., 2011).
The results of this study clearly indicated that the main factors that caused high variations in soil properties in the Batumi delta (Georgia) were the differences in parent material due to physiographical positions and land use type. Most of the measured soil properties in this study were highly correlated with soil clay content and OM content. While soils in the active delta section that were mostly sandy textured (sandy loam, loamy sand and sand) took place at the west bottom corner of the textural triangle, soils in the upper delta plain were relatively medium textured (sandy loam and loam). Clay content of soil increased from the active delta section through the upper delta plain as a function of physiographical position, and the upper plain soil consisted of, on average, about 1.5 times more clay than that of the active delta section soil. On the other hand, soil OM content showed great variability because of land use type, producing the highest OM in grassland areas fallowed by cultivated land. Therefore it is concluded that the physiographical position and land use type should be considered in evaluation of the differentiations in delta plain soil characteristics. It is expected that the results of this study provide scientific data for land use management in the Batumi delta, which had never been studied before.
The data resulting from this research were obtained by analyzing the soil samples taken by the authors of this article. You can see the detailed data in the Supplement.