IOP Conference Series: Earth and

Environmental Science

PAPER • OPEN ACCESS You may also like

- Experimental study on infiltration of
Assessment of soil erosion and landslides unsaturated soil and solution of
permeability coefficient
susceptibility based on hydrophysic soil properties Yuting Bai and Keming Tang

- Analysis of soil physical properties and

in Karangkobar catchment, Banjarnegara, infiltration rates for various landuses at
Gunung Dahu Research Forest, Bogor
Indonesia District, West Java Province
N E Saputra, C Wibowo and Y Lisnawati

- Analysis of classification hydrologic soil

To cite this article: F A Farizi et al 2019 IOP Conf. Ser.: Earth Environ. Sci. 355 012021 group distribution based on infiltration rate
of horton method in the upper Ciliwung
watershed
Firda Aulia Sartika, Dwita Sutjiningsih and
Evi Anggraheni
View the article online for updates and enhancements.

This content was downloaded from IP address 182.2.78.155 on 25/05/2024 at 17:01

The 3rd International Symposium on Agricultural and Biosystem Engineering IOP Publishing
IOP Conf. Series: Earth and Environmental Science 355 (2019) 012021 doi:10.1088/1755-1315/355/1/012021

Assessment of soil erosion and landslides susceptibility based

on hydrophysic soil properties in Karangkobar catchment,
Banjarnegara, Indonesia

F A Farizi1, Ngadisih1*, S Susanto1, H Suryatmojo2, P K V Tando1

1
Department of Agricultural and Biosystem Engineering, Universitas Gadjah Mada,
Yogyakarta, Indonesia
2
Laboratory of Watershed Management, Universitas Gadjah Mada, Yogyakarta,
Indonesia

E-mail: [email protected].

Abstract. Landslides and erosion are highly influenced by soil properties. Soil texture, structure,
particle density, bulk density, and porosity are the properties which correspond to infiltration
rates. This study aimed to assess soil erosion rate and landslide in upland of Karangkobar
catchment and to explore its correlation to hydrophysic soil properties. Material used in this
study were topographical maps, land use maps, administrative maps on scale 1:25000, and soil
samples. An unit land map was created by overlying topographic, landuse, and administrative
maps. This study found 16 unit lands. Soil sampling and infiltration test were conducted in every
unit land. Infiltration test adopted double ring infiltrometer test, while Horton's formula was
applied to determine the value of infiltration capacity. Soil hydrophysic properties in catchment
area Karangkobar: texture loam, granular structure, bulk density 0.62 – 1.09 g cm-3, particle
density 1.83 – 2.47 g cm-3, and infiltration capacity 1.5 cm hour-1 – 67.9 cm hour-1. Erosion rate
was 2 mm year-1 with the potential landslides classified as high. Soil hydrophysic properties
which have a strong correlation to erosion was bulk density (r = 0.517 *) and potential landslides
was infiltration rate (r = 0.641 **).

Introduction
Landslides are becoming common in Indonesia. Landslides struck in the Karangkobar, Banjarnegara in
2014. There are so many factors that cause landslides. Slope, soil, rainfall, land use, vibration, and
deforestation are generated of landslides [1,2]. Landslides are a type of soil erosion which transports
soil at a short time and very large volume. While causes of soil erosion are rainfall, soil, slope,
vegetation, and humans. Landslides and erosion are being related to the soil. Texture, structure, particle
density, bulk density, and porosity are influence of soil infiltration and hydrophysic soil [3].
Soil Erosion is process destroyed the soil by rainfall. After that, a soil was being transported to
another place by surface runoff [4,5]. While the process of landslides when water adsorbs into the soil
can increase the weight of the soil. When the water to penetrate to the impermeable soil and reach slip
plane, then the soil becomes slippery and move out of the slope [6].
Among all the predictive equations developed to estimate soil erosion, the most accepted, used,
convenient and suitable technique, is the Universal Soil Loss Equation (USLE). USLE is a soil erosion
model designed to predict long-term soil erosion from sheet or rill erosion [7,8]. Hydrophysic is a branch

Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
The 3rd International Symposium on Agricultural and Biosystem Engineering IOP Publishing
IOP Conf. Series: Earth and Environmental Science 355 (2019) 012021 doi:10.1088/1755-1315/355/1/012021

of geophysics that studies the natural state of water and physical processes in the hydrosphere on the
earth. Hydrophysic soil can mean the process, nature, and condition of water in the soil. The hydrophysic
soil properties are the interaction between fluid and soil components [6,7]. Based on the results of the
study of [8] that the hydrophysic soil properties depended on infiltration. The purpose of this study is to
identify the hydrophysic soil properties, calculate the soil erosion rate and mapping potential landslides,
and analyze the correlation of potential soil erosion and landslides based on the hydrophysic soil
properties.

Methodology
The study conducted at Karangkobar catchment, Banjarnegara from Maret until November 2018. The
material used in this study were topographical maps, landuse maps, administrative maps on scale
1:25000, soil sample, and infiltration test. There were some steps: survey, secondary data, and analyze
laboratory. An unit land map created by overlying topographic, landuse, and administrative maps. This
study found 16-unit lands. Soil sampling and infiltration test were being conducted in every unit land.
Infiltration test adopted double-ring infiltrometer test, then analyzed using Horton formula to determine
the value of infiltration capacity. The methods to analyzed soil physical properties of bulk density are
ring sampler, particle density is pycnometer, texture with pipe method, structure with USB digital
microscope, and porosity were being calculated based on the formula of Hanafiah [9].

bulk density
Porosity = 1 − × 100% (1)
particle density

The predictive of soil erosion used the USLE method with ArcGIS. The soil erosion unit converted
ton/ha/year to mm/year.
𝐀 = 𝐑 × 𝐊 × 𝐋𝐒 × 𝐂 × 𝐏 (2)
Where:
A = the average annual surface erosion rate (ton/ha/year)
R = the erosive rainfalls (KJ/ha)
K = soil erodibility (ton/KJ)
LS = slope length coefficient and slope degree (%)
CP = cover-management and conservation measures (value)

The classification of landslide areas used an assessment (scoring) of trigger parameters of landslides.
The assessment based on expert judgment through consultation with related experts. Thematic maps
overlay by considering scores to get a class of distribution of landslide areas. The physical properties of
the area used as trigger parameters for landslides are being presented in Table 1.

Table 1. The score for trigger of landslides

Variable Description Value Data Source Citation
< 500 1 Automatic
BP-DAS
Rainfall 500 – 999 2 Rainfall
Jeneberang-
(D) 1000 – 1999 3 Record
Walanae, 2010
(mm/year) 2000 – 2999 4 Karangkobar
[10]
> 3000 5 catchment
0–8% 1 SK Menteri
8 – 15 % 2 Pertanian Nomor
Slope Topographic
15 -25 % 3 837/KPTS/UM/19
(A) maps
25 – 45 % 4 80
> 45 % 5 [11]
Land use Forest land 1 Land use maps

2
The 3rd International Symposium on Agricultural and Biosystem Engineering IOP Publishing
IOP Conf. Series: Earth and Environmental Science 355 (2019) 012021 doi:10.1088/1755-1315/355/1/012021

(B) Dryland 2 BP-DAS

Built-up 3 Jeneberang-
Mixing garden 4 Walanae, 2010
Cropland 5 [10]
< 0.1 1
Soil 0.1 – 0.5 2
(Infiltration 0.5 – 2 3
Arsyad (2000)
rate) 2 – 6.5 4 Infiltration test
[12]
(C) 65 – 12.5 5
(cm/hour) 12.5 – 25 6
>25 7
Calculated using the formula proposed by Dibyosaputro (1999) [13]:
C−B
I= (3)
A
Where:
I = class of interval
A = the desired number of classes (5 Class)
B = the lowest value
C = the highest value
The criteria for determining the tolerance level for classes interval is 3.6. After that, the overlay of
four maps results in areas that have potential landslides. The data obtained from observations in the area,
then analyzed by qualitative descriptive statistic methods. Finally, the hydrophysic soil properties
correlated with potential soil erosion and landslides to find out of linkages used SPSS.

Result and Discussion

Table 2 shows that the texture Karangkobar catchment dominated by loam and silt loam textures. Based
on the results of the study of Apriyono [14] that clay is often found in Indonesia comes from volcanic
eruptions. Soil structure dominated by granular dominant. Based on the results of the study of Dai et al.,
[15], that granular soil structure decreases soil erosion rate. But, according to Karnawati [16] that is
effective for trigger landslides. The highest bulk density is unit land P-1 and the lowest is unit land P-3.
This is by the opinion of Hadi et al., [17] that on land was being subjected to continuous pressure
resulting in compaction so that the contents of bulk density will be higher.
Meanwhile, the particle density in all observation points ranged from 1.82 to 2.41 grams/cm3. This
is under the opinion of Surya et al. [18] that particle density is not changed in the long term because it
is being associated with a stable of the solid composition. Table 2, the highest porosity is unit land P-3
by 68.82% and bulk density by 0.61 gram/cm3. Bulk density has negative and significant correlation
with porosity, that is (r = -0.881**). This means that the increase in bulk density, porosity is low. This
is under the opinion of Kurnia et al., [19] that bulk density affects soil porosity.
Based on Table 2, the highest infiltration rate in unit land K-1 by 43.1 cm/hour. Soil with a low bulk
density causes high porosity and will increase the infiltration rate [20]. Porosity has positive and
significant correlation with infiltration rate (r = 0.522). This means that the increase of porosity, the
infiltration rate is high.

Tabel 2. Hydrophysic soil properties at Karangkobar catchment

ρb ρs Porositsy Infiltration
No Unit land Texture Structure
(gram/cm3) (gram/cm3) (%) (cm/jam)
1 K-1 0.80 2.04 60.77 43.1 loam
2 K-2-A 0.81 1.97 58.91 8.5 loam
Granular
3 K-2-B 1.03 2.18 52.90 12.1 loam
4 K-2-C 0.78 2.02 61.27 6.2 silt loam

3
The 3rd International Symposium on Agricultural and Biosystem Engineering IOP Publishing
IOP Conf. Series: Earth and Environmental Science 355 (2019) 012021 doi:10.1088/1755-1315/355/1/012021

5 K-3 0.93 1.96 52.68 6.9 silt loam

6 P-1 1.09 2.12 48.37 1.5 loam
7 P-2-A 0.83 1.81 54.42 17.4 loam
8 P-2-B 0.89 2.47 63.99 14.3 silt loam
9 P-2-C 0.71 1.82 60.89 5.2 loam
10 P-3 0.61 1.96 68.82 22.8 loam
11 S-1 0.77 2.20 64.99 27.1 sandy loam
12 TL-1 0.66 1.86 64.44 19.2 loam
13 TL-2-A 0.63 1.94 67.54 26.3 loam
14 TL-2-B 0.80 2.03 60.76 25.7 loam
15 TL-2-C 0.62 1.94 67.91 21.3 loam
16 TL-3 0.73 1.82 59.93 22.3 loam

Karangkobar catchment has a predictive soil erosion rate was 42.8 tons/ha/year or 2 mm/year. The
results of soil erosion prediction of Karangkobar catchment shown in Figure 1. The analysis shows that
soil erosion in Karangkobar catchment is being dominated by very low 58.6%, soil erosion low 33.8%,
soil erosion medium 3.9%, and finally soil erosion high 3.7%. Figure 2 shows a comparison of the rate
of predicted soil erosion in Karangkobar catchment in unit land. The highest soil erosion prediction rate
occurred in unit land K-2-B with erosion rate of 194.6 tons/ha/year or 8.6 mm/year. The unit land K-2-
B, the soil erosion rate is being dominated by the high at 84% and medium at 16%. While the smallest
soil erosion prediction rate occurred in unit land K-2-C with erosion rate of 1,1 ton/ha/year or 0.1
mm/year. The unit land K-2-C, the soil erosion rate is being dominated by very low at 100%.
The hydrophysic soil properties cause the rate of high soil erosion in the Karangkobar catchment by
low bulk density and high porosity. Bulk density has positive and significant correlation with erosion
rate (r = 0.517*). This means that the increase of bulk density, the soil erosion is high. While porosity
has a negative and significant correlation with erosion rates of (r = -0.507*), this means that increase
porosity, the soil erosion rate is low. This is in accordance with the opinion of Hakim [21] that porosity
is the ability of the soil absorb water. The soil is harder to absorb water, the porosity is smaller.

10
60
9
Erosion rate (mm/year)

50 8
7
40
Area (%)

6
30 5
20 4
3
10 2
0 1
0
P-1

P-3
K-2-A

S-1
K-2-B
K-2-C

TL-1

TL-3
K-1

K-3

P-2-A
P-2-B
P-2-C

TL-2-A
TL-2-B
TL-2-C

Clasification soil erosion

Unit land
Figure 1. The soil erosion area at Figure 2. Soil Erosion rate at unit land of Karangkobar
Karangkobar catchment catchment

Figure 3 shows the prediction of landslides in Karangkobar catchment with a high of 67%, medium
by 18%, a very high landslide of 11%, and potential low of 4%. Based on PUPR Regulation No. 22 of
2007 [22], the landslide occurred when rainfall rate is high (> 2,500 mm/year), with steep slopes (>
40%), and earthquake-prone areas. This is in accordance with the prediction of landslides which occur
in the Karangkobar catchment on steep slopes > 40%. While the slope of less than 15% the potential is
low.

4
The 3rd International Symposium on Agricultural and Biosystem Engineering IOP Publishing
IOP Conf. Series: Earth and Environmental Science 355 (2019) 012021 doi:10.1088/1755-1315/355/1/012021

Based on Figure 4, unit land that has very high potential landslide at unit land K-1 and S-1. The unit
land K-1 and S-1 have fast infiltration rate. This is in accordance with the opinion of Karnawati [16]
that soil is easy to absorb water will effective to triggers of landslides on the slope. The infiltration rate
has positive and significant correlation with landslides (r = 0,641**). This is in accordance with the
statement of Arrozi et al. [23] that the increase of infiltration from rainfall cause saturated soil layers
and increase load on the slope.

70
100
60
50 80
Area (%)

Area (%)
40 60
30
40
20
10 20
0
0

P-II-A
K-II-A
K-II-B

P-II-B

TL-I
P-II-C
K-II-C

P-I

TL-III
k-III
K-I

S-I

TL-II-A
TL-II-B
TL-II-C
P-III
Low Medium
Clasification landslides Unit land

Figure 3. The landslides area at Figure 4. Potensial landslide at unit land of

Karangkobar catchment Karangkobar catchment

Conclusion
Karangkobar catchment is being dominated by clay texture with bulk density values (0.61 – 1.09 g.cm-
3
), particle density (1.81 – 2.47 g.cm-3), soil porosity (48 – 68%), and infiltration rate (1.5 – 43.1 cm.hour-
1
). The erosion rate in Karangkobar catchment is 42.8 tons.ha-1.year-1 or 2 mm.year-1 with a high potential
of landslides of 67%, medium by 18%, very high by 11%, and low by 4% of the total area Karangkobar
catchment. The hydrophysic soil properties have a significant correlation with soil erosion are bulk
density (r = 0.517*). Whereas the hydrophysic soil properties have significant correlation with potential
of landslides are infiltration rate (r = 0.641**).

Acknowledgement
We owe a lot to RTA grants Universitas Gadjah Mada. We also thank to Herlina and Viyata, master
student of Forestry UGM, for a good team work during field work.

References
[1] Ngadisih, Yatabe R, Bhandary N P and Dahal R K 2014 Integration of statistical and heuristic
approaches for landslide risk analysis: A case of volcanic mountains in West Java Province,
Indonesia Georisk 8 29–47
[2] Ngadisih, Netra P Bhandary, Ryuichi Yatabe, Ranjan K Dahal Logistic regression and artificial
neural network models for mapping of regional-scale landslide susceptibility in volcanic
mountains of West Java (Indonesia) AIP Con. Proc.
[3] A Chay 2010 Hidrologi dan Pengelolaan Daerah Aliran Sungai (Yogyakarta: Gadjah Mada
University Press)
[4] W H. Wischmeier and J. V. Mannering 1969 Relation of Soil Properties to its Erodibility1 Soil
Sci. Soc. Am. J., 33 131–7
[5] S Hardjowigeno 2007 Ilmu Tanah (Jakarta: Akademika Pressindo)
[6] Badan Geologi 2015 Gerakan Tanah (Direktorat Vulkanologi dan Mitigasi Bencana Geologi)
[7] W H Wischmeier and S L D 1978 Predicting Rainfall-Erosion Losses: A Guide to Conservation

5
The 3rd International Symposium on Agricultural and Biosystem Engineering IOP Publishing
IOP Conf. Series: Earth and Environmental Science 355 (2019) 012021 doi:10.1088/1755-1315/355/1/012021

Planning, Agricultur.
[8] I W Sutapa 2010 Analisis potensi erosi pada daerah aliran sungai (das) di sulawesi tengah J.
SMARTek 8 169–81
[9] K A Hanafiah 2005 Dasar-dasar Ilmu Tanah (Jakarta: PT Raja Gr)
[10] BPDAS Jeneberang-Walanae 2010 Laporan Karakteristik DAS Jeneberang-Walanae 2010
[11] Kementian Pertanian 1980 Kriteria dan Tata Cara Penetapan Hutan Lindung (Indonesia)
[12] S Arsyad 2000 Konservasi Tanah dan Air (Bogor: IPB Press)
[13] S Dibyosaputro 1999 Longsor Lahan di Daerah Kecamatan Samigaluh, Kabupaten Kulon Progo,
Daerah Istimewa Yogyakarta Maj. Geogr. Indones. 16 13–34
[14] A Apriyono 2009 Analisis Penyebab Tanah Longsor di Kalitlaga Banjarnegara Din. Rekayasa 5
14–8
[15] W Dai and Y Huang 2006 Relation of soil organic matter concentration to climate and altitude in
zonal soils of China Catena 65 87–94
[16] D Karnawati 2005 Bencana Alam Gerakan Massa Tanah di Indonesia dan Upaya
Penanggulangannya (Yogyakarta: Universitas Gandjah Mada)
[17] B Hadi, Al Y Yunus, and M Idkham 2012 Analisis Sifat Fisika Tanah Akibat Lintasan Dan Bajak
Traktor Roda Empat Manaj. Sumberd. Lahan 1 43–53
[18] J A Surya and Y Nuraini 2017 Kajian porositas tanah pada pemberian beberapa jenis bahan
organik di perkebunan kopi robusta 4 463–71
[19] U Kurnia, F Agus, A Adimihadja, and A Dariah 2006 Sifat Fisik Tanah dan Metode Analisisnya
(Bogor: Badan Penelitian dan Pengembangan Pertanian)
[20] S Arsyad 2010 Konservasi Tanah dan Air (Bogor: IPB Press)
[21] N Hakim et al 1986 Dasar-dasar Ilmu Tanah (Lampung: Universitas Lampung)
[22] Kementrian PUPR 2007 Pedoman Penataan Ruang Kawasan Rawan Bencana Longsor
(Indonesia)
[23] M M F Arrozi, N S Surjandari, and N Djarwanti 2015 Menggunakan Metode Bishop
Disederhanakan ( Studi Kasus Di Dusun Pagah , Hargantoro , Tirtomoyo , Wonogiri ) e-
Journal Matriks Tek. Sipil 542–7