geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 5 no. 3, december 2020 available online since 30 december 2020 at : https://jurnal.unej.ac.id/index.php/geosi/issue/view/1064 geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 5 no. 3, december 2020 available online since 30 december 2020 at : https://jurnal.unej.ac.id/index.php/geosi/issue/view/1064 geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 5 no. 3, december 2020 available online since 30 december 2020 at : https://jurnal.unej.ac.id/index.php/geosi/issue/view/1064 geosfera indonesia, vol. 5 no. 3 (2020) accredited by the ministry of research , technology , and higher education of the republic of indonesia, no. 30/e/kpt/2019 . editorial team editor in chief fahmi arif kurnianto (scopus id: 57208473928) department of geography education, university of jember, indonesia advisory international editorial boards mihai ciprian margarint (scopus id : 36698019400) department of geography, alexandru ioan cuza university of iasi, romania franck lavigne (scopus id : 15738234900) physical geography laboratory, université paris 1 panthéon-sorbonne, france fahrudi ahwan ikhsan (scopus id : 57208469257) department of geography education, university of jember, indonesia mustafa ustuner (scopus id : 56246446800) department of geomatics engineering, yildiz technical university, turkey bashkim idrizi (scopus id : 55937683800) department of geodesy, university "mother teresa" -skopje, macedonia guillermo hector re (scopus id : 7102894803) department of geology, universidad de buenos aires, buenos aires, argentina laras tursilowati (scopus id : 55317967300) indonesian national institute of aeronautics and space (lapan), indonesia kuppanagounder kumaraswamy (scopus id : 6602935596) department of geography, bharathidasan university, tiruchirappalli, india copyright (c) 2020 geosfera indonesia journal and department of geography education, university of jember focus and scope geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education : collaborative learning; comparative learning; curriculum; e-learning ; instructional technology; learning community; life skills ; remedial teaching; taxonomy of educational objectives (bloom's taxonomy); new technology; industry and education : a continous collaboration; blended learning; character; constructivist learning; disrupting innovation; expeditionary learning; flexible learning; flipped classroom; flipped learning; gamification; global view; ground up diversity; high-quality teachers; hip-hop education (hiphoped); lesson study; mobile education; problem based learning; process oriented guided inquiry lessons (pogil); project based learning (pbl); start-up; student centred learning; autodidacticism (self-teaching); informal learning; unschooling or homeschooling; pisa task, (2) physical geography : tectonics and regional structure; glacial processes and landforms; fluvial sequences; fluvial processes and landforms; mass movement; hillslopes and soil erosion; slopes processes; karst processes and landforms; aeolian processes and landforms; coastal dunes and arid landforms; coastal and marine processes; theoretical and quantitative geomorphology; soil geomorphology; soil geography; lithology; hydrogeography, (3) human geography : cultural geography; political geography; social geography; population geography; urban geography, (4) geographic information system (gis) : data collection and acquisition; data structures and algorithms; spatio-temporal databases; spatial analysis, data mining, and decision support systems; cartography; location based services; uncertainty handling in spatial data; topology; geo-computation; geo-telematics; spatial information infrastructures; interoperability and open systems; applications of geoinformation technology (all possible domains), (5) remote sensing: multi-spectral and hyperspectral remote sensing; active and passive microwave remote sensing; lidar and laser scanning; geometric reconstruction; physical modeling and signatures; change detection; image processing and pattern recognition; data fusion and data assimilation; dedicated satellite missions; operational processing facilities; spaceborne, airborne and terrestrial platforms; remote sensing applications, (6) environmental science : environmental geography; environmental education; climate change; land use and cover change; pollution; natural resources management; conservation; management and valorisation of waste; development of methods for environmental quality management; environmental system modelling and optimization; environmental analysis and assessment; social, economic and policy aspects of environmental management, (7) disaster risk reduction : risk awareness and assessment including hazard analysis and vulnerability/capacity analysis for natural disaster risk reduction; knowledge development including education, training, research and information for natural disaster risk reduction; public commitment and institutional frameworks, including organisational, policy, legislation and community action for natural disaster risk reduction. publication information secretariat of geosfera indonesia department of geography education, university of jember, fkip building, jl. kalimantan 37, jember, east java, 68121, indonesia. telp. (0331) 334988 / 330738 email : geografi.fkip@unej.ac.id website : https://jurnal.unej.ac.id/index.php/geosi copyright (c) 2020 geosfera indonesia journal and department of geography education, university of jember geosfera indonesia (geos. ind.) : | issn: 2598-9723 (print)| issn: 2614-8528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. its published three times a year in april, august, and december. geosfera indonesia is accredited by the ministry of research, technology and higher education of the republic of indonesia (ristekdikti), no. 30/e/kpt/2019. this journal has been covered by following indexing and abstracting services: (1) cabi : cab abstracts (web of science); (2) directory of open access journal (doaj); (3) ebsco; (4) google scholar (5) one search (national library of republic of indonesia); (6) sinta 2. table of contents the assessment of deforestation impact towards microclimate and environment in ilorin, nigeria 301-317toluwalope mubo agaja, elisha ademola adeleke, enekole esther adeniyi, precious temilade afolayan development of multimedia learning geography based on adobe flash to increase students’ curiosity 318-334 wahid yuda rejeki* and m. mukminan geospatial approach for the analysis of forest cover change detection using machine learning 335-351 r. sanjeeva reddy, g. anjan babu, a. rama mohan reddy emerging geospatial technologies in environmental research, education, and outreach 352-363sergio bernardes, margueritte madden, ashurst walker, andrew knight, nicholas neel, akshay mendki, dhaval bhanderi, andrew guest, shannon healy, thomas jordan a review paper on monitoring environmental consequences of land cover dynamics with the help of geo-informatics technologies copyright (c) 2020 geosfera indonesia journal and department of geography education, university of jember 364-389 ziyad ahmed abdo, satya prakash cover geosfera 5,3.pdf (p.1) geosfera indonesia akred.pdf (p.2) focus and scope.pdf (p.3-4) acknowledgement the editorial board would like to thank and express appreciation to all peer reviewers who have reviewed the manuscripts for publication of geosfera indonesia vol. 5 no. 2, august 2020. marinela istrate (scopus id : 56672996500) department of geography, universitatea alexandru ioan cuza, iasi, romania lucian rosu (scopus id : 56958158800) department of geography, universitatea alexandru ioan cuza, iasi, romania mihai niculita (scopus id : 55022909500) department of geography, alexandru ioan cuza university of iasi, romania tin lukic (scopus id : 54795557700) department of geography, university of novi sad, serbia fatih adiguzel (scopus id : 57204062132) , department of geography, nevşehir haci bektaş veli üniversitesi, nevsehir, turkey elan artono nurdin (scopus id : 57192649049) department of geography education, university of jember, indonesia bejo apriyanto (scopus id : 57208468203) department of geography education, university of jember, indonesia prima widayani (scopus id : 57188879871) faculty of geography, gadjah mada university, indonesia ivan taslim, (scopus id : 57203022349) department of geography, universitas muhammadiyah gorontalo, indonesia efdal kaya, (scopus id : 57202135756) department of architecture and urban planning, iskenderun technical university, iskenderun, turkey batuhan kilic, (scopus id : 57204675825) geomatic engineering department, yildiz technical university, istanbul, turkey nailul insani, (scopus id : 57208471498) department of geography, universitas negeri malang, indonesia rendra zainal maliki department of geography education, universitas tadulako, indonesia https://www.scopus.com/authid/detail.uri?authorid=56672996500 https://www.scopus.com/authid/detail.uri?authorid=56958158800 https://www.scopus.com/authid/detail.uri?authorid=55022909500 https://www.scopus.com/authid/detail.uri?authorid=54795557700 https://www.scopus.com/authid/detail.uri?authorid=57204062132 https://www.scopus.com/authid/detail.uri?authorid=57192649049 https://www.scopus.com/authid/detail.uri?authorid=57208468203 https://www.scopus.com/authid/detail.uri?authorid=57188879871 https://www.scopus.com/authid/detail.uri?authorid=57203022349 https://www.scopus.com/authid/detail.uri?authorid=57202135756 https://www.scopus.com/authid/detail.uri?authorid=57204675825 https://www.scopus.com/authid/detail.uri?authorid=57208471498 https://scholar.google.co.id/citations?user=qivrcmkaaaaj&hl=id geos geosfera indonesia geosfera indonesia (journal initial : geosi, journal abbreviation : geos. ind.) : issn: 2598-9723 (print) ; issn: 26148528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education, (2) geography (physical geography and human geography), (3) geographic information system (gis), (4) remote sensing, (5) environmental science, and (6) disaster risk reduction. it is published three times a year in april, august, and december. geosfera indonesia is accredited by the ministry of research, technology and higher education of the republic of indonesia (ristekdikti), no. 30/e/kpt/2019 (20172022). indexed by : geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 5 no. 2, august 2020 available online since 1 june 2020 at : https://jurnal.unej.ac.id/index.php/geosi/issue/view/952 geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 5 no. 2, august 2020 available online since 1 june 2020 at : https://jurnal.unej.ac.id/index.php/geosi/issue/view/952 geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 5 no. 2, august 2020 available online since 1 june 2020 at : https://jurnal.unej.ac.id/index.php/geosi/issue/view/952 geosfera indonesia, vol. 5 no. 2 (2020) accredited by the ministry of research , technology , and higher education of the republic of indonesia, no. 30/e/kpt/2019 . editorial team editor in chief fahmi arif kurnianto (scopus id: 57208473928) department of geography education, university of jember, indonesia advisory international editorial boards mihai ciprian margarint (scopus id : 36698019400) department of geography, alexandru ioan cuza university of iasi, romania franck lavigne (scopus id : 15738234900) physical geography laboratory, université paris 1 panthéon-sorbonne, france fahrudi ahwan ikhsan (scopus id : 57208469257) department of geography education, university of jember, indonesia mustafa ustuner (scopus id : 56246446800) department of geomatics engineering, yildiz technical university, turkey bashkim idrizi (scopus id : 55937683800) department of geodesy, university "mother teresa" -skopje, macedonia guillermo hector re (scopus id : 7102894803) department of geology, universidad de buenos aires, buenos aires, argentina laras tursilowati (scopus id : 55317967300) indonesian national institute of aeronautics and space (lapan), indonesia kuppanagounder kumaraswamy (scopus id : 6602935596) department of geography, bharathidasan university, tiruchirappalli, india copyright (c) 2020 geosfera indonesia journal and department of geography education, university of jember focus and scope geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education : collaborative learning; comparative learning; curriculum; e-learning ; instructional technology; learning community; life skills ; remedial teaching; taxonomy of educational objectives (bloom's taxonomy); new technology; industry and education : a continous collaboration; blended learning; character; constructivist learning; disrupting innovation; expeditionary learning; flexible learning; flipped classroom; flipped learning; gamification; global view; ground up diversity; high-quality teachers; hip-hop education (hiphoped); lesson study; mobile education; problem based learning; process oriented guided inquiry lessons (pogil); project based learning (pbl); start-up; student centred learning; autodidacticism (self-teaching); informal learning; unschooling or homeschooling; pisa task, (2) physical geography : tectonics and regional structure; glacial processes and landforms; fluvial sequences; fluvial processes and landforms; mass movement; hillslopes and soil erosion; slopes processes; karst processes and landforms; aeolian processes and landforms; coastal dunes and arid landforms; coastal and marine processes; theoretical and quantitative geomorphology; soil geomorphology; soil geography; lithology; hydrogeography, (3) human geography : cultural geography; political geography; social geography; population geography; urban geography, (4) geographic information system (gis) : data collection and acquisition; data structures and algorithms; spatio-temporal databases; spatial analysis, data mining, and decision support systems; cartography; location based services; uncertainty handling in spatial data; topology; geo-computation; geo-telematics; spatial information infrastructures; interoperability and open systems; applications of geoinformation technology (all possible domains), (5) remote sensing: multi-spectral and hyperspectral remote sensing; active and passive microwave remote sensing; lidar and laser scanning; geometric reconstruction; physical modeling and signatures; change detection; image processing and pattern recognition; data fusion and data assimilation; dedicated satellite missions; operational processing facilities; spaceborne, airborne and terrestrial platforms; remote sensing applications, (6) environmental science : environmental geography; environmental education; climate change; land use and cover change; pollution; natural resources management; conservation; management and valorisation of waste; development of methods for environmental quality management; environmental system modelling and optimization; environmental analysis and assessment; social, economic and policy aspects of environmental management, (7) disaster risk reduction : risk awareness and assessment including hazard analysis and vulnerability/capacity analysis for natural disaster risk reduction; knowledge development including education, training, research and information for natural disaster risk reduction; public commitment and institutional frameworks, including organisational, policy, legislation and community action for natural disaster risk reduction. publication information secretariat of geosfera indonesia department of geography education, university of jember, fkip building, jl. kalimantan 37, jember, east java, 68121, indonesia. telp. (0331) 334988 / 330738 email : geografi.fkip@unej.ac.id website : https://jurnal.unej.ac.id/index.php/geosi copyright (c) 2020 geosfera indonesia journal and department of geography education, university of jember geosfera indonesia (geos. ind.) : | issn: 2598-9723 (print)| issn: 2614-8528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. its published three times a year in april, august, and december. geosfera indonesia is accredited by the ministry of research, technology and higher education of the republic of indonesia (ristekdikti), no. 30/e/kpt/2019. this journal has been covered by following indexing and abstracting services: (1) cabi : cab abstracts (web of science); (2) directory of open access journal (doaj); (3) ebsco; (4) google scholar (5) one search (national library of republic of indonesia); (6) sinta 2. table of contents a new algorithm for the grid cell-based runoff routing model based on travel time concept 160-185 baina afkril, m. pramono hadi, slamet suprayogi a preliminary study on tsunami disaster in yogyakarta: identification of vulnerability order and components 186-195 lina wahyuni, muh. aris marfai, m. pramono hadi soil zonation and the shaking table test of the embankment on clayey soil 196-209 ripon hore, sudipta chakraborty, md. fayjul bari, ayaz mahmud shuvon, mehedi ahmed ansary modeling land use and land cover dynamic using geographic information system and markov-ca 210-225 millary agung widiawaty, arif ismail, moh. dede, n. nurhanifah spatio-statistical analysis of rainfall and temperature distribution, anomaly and trend in nigeria 226-249 elisha ademola adeleke and eniola aminat orebayo the facies and metamorphism types determination of metamorphic rock in the part of mekongga complex 268-287 m. musnajam, ahmad tarmizi abd karim , n. nurfadillah, fahrudi ahwan ikhsan, andri estining sejati building density level of urban slum area in jakarta tenty melvianti legarias , renny nurhasana, edy irwansyah land value potential zonation : implication towards urban planning 288-300 revi mainaki, anita eka putri, dwiyono hari utomo copyright (c) 2020 geosfera indonesia journal and department of geography education, university of jember 250-267 cover geosfera 5,2 malam.pdf (p.1) geosfera indonesia akredcover.pdf (p.2) focus and scope2.pdf (p.3-4) acknowledgement the editorial board would like to thank and express appreciation to all peer reviewers who have reviewed the manuscripts for publication of geosfera indonesia. marinela istrate (scopus id : 56672996500) department of geography, universitatea alexandru ioan cuza, iasi, romania lucian rosu (scopus id : 56958158800) department of geography, universitatea alexandru ioan cuza, iasi, romania mihai niculita (scopus id : 55022909500) department of geography, alexandru ioan cuza university of iasi, romania tin lukic (scopus id : 54795557700) department of geography, university of novi sad, serbia fatih adiguzel (scopus id : 57204062132) , department of geography, nevşehir haci bektaş veli üniversitesi, nevsehir, turkey elan artono nurdin (scopus id : 57192649049) department of geography education, university of jember, indonesia bejo apriyanto (scopus id : 57208468203) department of geography education, university of jember, indonesia prima widayani (scopus id : 57188879871) faculty of geography, gadjah mada university, indonesia ivan taslim (scopus id : 57203022349) department of geography, universitas muhammadiyah gorontalo, indonesia efdal kaya (scopus id : 57202135756) department of architecture and urban planning, iskenderun technical university, iskenderun, turkey batuhan kilic (scopus id : 57204675825) geomatic engineering department, yildiz technical university, istanbul, turkey nailul insani (scopus id : 57208471498) department of geography, universitas negeri malang, indonesia rendra zainal maliki department of geography education, universitas tadulako, indonesia faisal arif setiawan department of geography education, universitas lambung mangkurat, indonesia aulia ulfa farahdiba (scopus id : 57208130620) , department of environmental engineering, universitas pembangunan nasional "veteran" jawa timur, indonesia guruh samodra (scopus id : 55053491600) faculty of geography, universitas gadjah mada, yogyakarta, indonesia andri estining sejati (scopus id : 57211280452) geography education study program, universitas sembilanbelas november, kolaka, indonesia wahid akhsin budi nur sidiq (scopus id : 57201676728), department of geography, universitas negeri semarang, semarang indonesia pertiwi andarani (scopus id : 55959123900), department of environmental engineering, universitas diponegoro, semarang, indonesia prama ardha aryaguna (scopus id : 56236908900), department of survey and mapping, faculty of engineering, universitas esa unggul, jakarta, indonesia ionuț minea (scopus id : 56951000900), faculty of geography and geology, department of geography, alexandru ioan cuza university, 20 a, carol i bd., 700505 iasi, romania rosmadi fauzi (scopus id : 6508044388), department of geography, university of malaya, kuala lumpur, malaysia geos geosfera indonesia geosfera indonesia (journal initial : geosi, journal abbreviation : geos. ind.) : issn: 2598-9723 (print) ; issn: 26148528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education, (2) geography (physical geography and human geography), (3) geographic information system (gis), (4) remote sensing, (5) environmental science, and (6) disaster risk reduction. it is published three times a year in april, august, and december. geosfera indonesia is accredited (sinta 2) by the ministry of research and technology of the republic of indonesia (ristekbrin), no. 200/m/kpt/2020 (2020-2025). indexed by : geos geosfera indonesia geosfera indonesia (journal initial : geosi, journal abbreviation : geos. ind.) : issn: 2598-9723 (print) ; issn: 26148528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education, (2) geography (physical geography and human geography), (3) geographic information system (gis), (4) remote sensing, (5) environmental science, and (6) disaster risk reduction. it is published three times a year in april, august, and december. geosfera indonesia is accredited (sinta 2) by the ministry of research and technology of the republic of indonesia (ristekbrin), no. 200/m/kpt/2020 (2020-2025). indexed by : geos geosfera indonesia geosfera indonesia (journal initial : geosi, journal abbreviation : geos. ind.) : issn: 2598-9723 (print) ; issn: 26148528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education, (2) geography (physical geography and human geography), (3) geographic information system (gis), (4) remote sensing, (5) environmental science, and (6) disaster risk reduction. it is published three times a year in april, august, and december. geosfera indonesia is accredited (sinta 2) by the ministry of research and technology of the republic of indonesia (ristekbrin), no. 200/m/kpt/2020 (2020-2025). indexed by : acknowledgement the editorial board would like to thank and express appreciation to all peer reviewers who have reviewed the manuscripts for publication of geosfera indonesia. marinela istrate (scopus id : 56672996500) department of geography, universitatea alexandru ioan cuza, iasi, romania lucian rosu (scopus id : 56958158800) department of geography, universitatea alexandru ioan cuza, iasi, romania mihai niculita (scopus id : 55022909500) department of geography, alexandru ioan cuza university of iasi, romania tin lukic (scopus id : 54795557700) department of geography, university of novi sad, serbia fatih adiguzel (scopus id : 57204062132) , department of geography, nevşehir haci bektaş veli üniversitesi, nevsehir, turkey elan artono nurdin (scopus id : 57192649049) department of geography education, university of jember, indonesia bejo apriyanto (scopus id : 57208468203) department of geography education, university of jember, indonesia prima widayani (scopus id : 57188879871) faculty of geography, gadjah mada university, indonesia ivan taslim, (scopus id : 57203022349) department of geography, universitas muhammadiyah gorontalo, indonesia efdal kaya, (scopus id : 57202135756) department of architecture and urban planning, iskenderun technical university, iskenderun, turkey batuhan kilic, (scopus id : 57204675825) geomatic engineering department, yildiz technical university, istanbul, turkey nailul insani, (scopus id : 57208471498) department of geography, universitas negeri malang, indonesia rendra zainal maliki department of geography education, universitas tadulako, indonesia faisal arif setiawan department of geography education, universitas lambung mangkurat, indonesia geos geosfera indonesia geosfera indonesia (journal initial : geosi, journal abbreviation : geos. ind.) : issn: 2598-9723 (print) ; issn: 26148528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education, (2) geography (physical geography and human geography), (3) geographic information system (gis), (4) remote sensing, (5) environmental science, and (6) disaster risk reduction. it is published three times a year in april, august, and december. geosfera indonesia is accredited by the ministry of research, technology and higher education of the republic of indonesia (ristekdikti), no. 30/e/kpt/2019 (20172022). indexed by : geos geosfera indonesia geosfera indonesia (journal initial : geosi, journal abbreviation : geos. ind.) : issn: 2598-9723 (print) ; issn: 26148528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education, (2) geography (physical geography and human geography), (3) geographic information system (gis), (4) remote sensing, (5) environmental science, and (6) disaster risk reduction. it is published three times a year in april, august, and december. geosfera indonesia is accredited by the ministry of research, technology and higher education of the republic of indonesia (ristekdikti), no. 30/e/kpt/2019 (20172022). indexed by : geos geosfera indonesia geosfera indonesia (journal initial : geosi, journal abbreviation : geos. ind.) : issn: 2598-9723 (print) ; issn: 26148528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education, (2) geography (physical geography and human geography), (3) geographic information system (gis), (4) remote sensing, (5) environmental science, and (6) disaster risk reduction. it is published three times a year in april, august, and december. geosfera indonesia is accredited by the ministry of research, technology and higher education of the republic of indonesia (ristekdikti), no. 30/e/kpt/2019 (20172022). indexed by : geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 6 no. 1, april 2021 available online since 25 april 2021 at : https://jurnal.unej.ac.id/index.php/geosi geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 6 no. 1, april 2021 available online since 25 april 2021 at : https://jurnal.unej.ac.id/index.php/geosi geosfera indonesia published by: department of geography education university of jember e-issn 2614-8528 p-issn 2598-9723 geos. ind. vol. 6 no. 1, april 2021 available online since 25 april 2021 at : https://jurnal.unej.ac.id/index.php/geosi geosfera indonesia, vol. 6 no. 1 (2021) accredited (sinta 2) by the ministry of research and technology of the republic of indonesia (ristekbrin), no. 200/m/kpt/2020 . editorial team editor in chief fahmi arif kurnianto (scopus id: 57208473928) department of geography education, university of jember, indonesia advisory international editorial boards mihai ciprian margarint (scopus id : 36698019400) department of geography, alexandru ioan cuza university of iasi, romania fahrudi ahwan ikhsan (scopus id : 57208469257) department of geography education, university of jember, indonesia mustafa ustuner (scopus id : 56246446800) department of geomatics engineering, yildiz technical university, turkey bashkim idrizi (scopus id : 55937683800) department of geodesy, university "mother teresa" -skopje, macedonia guillermo hector re (scopus id : 7102894803) department of geology, universidad de buenos aires, buenos aires, argentina laras tursilowati (scopus id : 55317967300) indonesian national institute of aeronautics and space (lapan), indonesia m. asyroful mujib (scopus id : 57217104177) department of geography education, university of jember, indonesia lyubka pashova (scopus id : 6506546577) national institute of geophysics, geodesy and geography, bulgarian academy of sciencesdisabled, sofia, bulgaria layout editors yucha risdarani (department of geography education, university of jember, indonesia) m. mitasari (department of geography education, university of jember, indonesia) giofani ginolla ardiyanto (department of geography education, university of jember, indonesia) copyright (c) 2021 geosfera indonesia and department of geography education, university of jember https://drive.google.com/file/d/1qarud7btkxdsqlrwobcji-eacl76rta8/view?usp=sharing https://www.scopus.com/authid/detail.uri?authorid=57208473928 https://www.scopus.com/authid/detail.uri?authorid=36698019400 https://www.scopus.com/authid/detail.uri?authorid=57208469257 https://www.scopus.com/authid/detail.uri?authorid=56246446800 https://www.scopus.com/authid/detail.uri?authorid=55937683800 https://www.scopus.com/authid/detail.uri?authorid=7102894803 https://www.scopus.com/authid/detail.uri?authorid=55317967300 https://www.scopus.com/authid/detail.uri?authorid=57217104177 https://www.scopus.com/authid/detail.uri?authorid=6506546577 focus and scope geosfera indonesia welcomes high quality original research articles, short communications, and review articles written by researchers, academicians, professional, and practitioners from all over the world about : (1) geography education : collaborative learning; comparative learning; curriculum; e-learning ; instructional technology; learning community; life skills ; remedial teaching; taxonomy of educational objectives (bloom's taxonomy); new technology; industry and education : a continous collaboration; blended learning; character; constructivist learning; disrupting innovation; expeditionary learning; flexible learning; flipped classroom; flipped learning; gamification; global view; ground up diversity; high-quality teachers; hip-hop education (hiphoped); lesson study; mobile education; problem based learning; process oriented guided inquiry lessons (pogil); project based learning (pbl); start-up; student centred learning; autodidacticism (self-teaching); informal learning; unschooling or homeschooling; pisa task, (2) physical geography : tectonics and regional structure; glacial processes and landforms; fluvial sequences; fluvial processes and landforms; mass movement; hillslopes and soil erosion; slopes processes; karst processes and landforms; aeolian processes and landforms; coastal dunes and arid landforms; coastal and marine processes; theoretical and quantitative geomorphology; soil geomorphology; soil geography; lithology; hydrogeography, (3) human geography : cultural geography; political geography; social geography; population geography; urban geography, (4) geographic information system (gis) : data collection and acquisition; data structures and algorithms; spatio-temporal databases; spatial analysis, data mining, and decision support systems; cartography; location based services; uncertainty handling in spatial data; topology; geo-computation; geo-telematics; spatial information infrastructures; interoperability and open systems; applications of geoinformation technology (all possible domains), (5) remote sensing: multi-spectral and hyperspectral remote sensing; active and passive microwave remote sensing; lidar and laser scanning; geometric reconstruction; physical modeling and signatures; change detection; image processing and pattern recognition; data fusion and data assimilation; dedicated satellite missions; operational processing facilities; spaceborne, airborne and terrestrial platforms; remote sensing applications, (6) environmental science : environmental geography; environmental education; climate change; land use and cover change; pollution; natural resources management; conservation; management and valorisation of waste; development of methods for environmental quality management; environmental system modelling and optimization; environmental analysis and assessment; social, economic and policy aspects of environmental management, (7) disaster risk reduction : risk awareness and assessment including hazard analysis and vulnerability/capacity analysis for natural disaster risk reduction; knowledge development including education, training, research and information for natural disaster risk reduction; public commitment and institutional frameworks, including organisational, policy, legislation and community action for natural disaster risk reduction. publication information geosfera indonesia (geos. ind.) : | issn: 2598-9723 (print)| issn: 2614-8528 (online) is an international open access and peer-reviewed journal, published by department of geography education, university of jember, indonesia. its published three times a year in april, august, and december. geosfera indonesia is accredited (sinta 2) by the ministry of research and technology of the republic of indonesia (ristekbrin), no. 200/m/kpt/2020. this journal has been covered by following indexing and abstracting services: (1) cabi : cab abstracts ; (2) directory of open access journal (doaj); (3) ebsco; (4) google scholar (5) one search (national library of republic of indonesia), (6) web of science. secretariat of geosfera indonesia department of geography education, university of jember, fkip building, jl. kalimantan 37, jember, east java, 68121, indonesia. telp. (0331) 334988 / 330738 email : geografi.fkip@unej.ac.id website : https://jurnal.unej.ac.id/index.php/geosi copyright (c) 2021 geosfera indonesia and department of geography education, university of jember table of contents rethinking urbanization: a transit-information-communication – technology-oriented development path for the developing countries and post-industrial towns 1-19 schuman lam, heng li, ann t.w. yu landslide hazard analysis using a multilayered approach based on various input data configurations 20-39 ilyas ahmad huqqani, tay lea tien, junita mohamad-saleh quantifying the significance of distance to temporal dynamics of covid-19 cases in nigeria using a geographic information system 40-54 ifeyinwa sarah obuekwe, umar saleh anka, sodiq opeyemi ibrahim, usman ahmad adam assessment of water balance at mayang watershed, east java 55-76ariska mia christiwarda sihombing, indarto indarto, sri wahyuningsih development of web-based gis alert system for informing environmental risk of dengue infections in major cities of pakistan 77-95 naureen zainab, aqil tariq, saima siddiqui 96-126 assessing the impacts of climate variability on rural households in agricultural land through the application of livelihood vulnerability index ginjo gitima, abiyot legesse, dereje biru copyright (c) 2021 geosfera indonesia and department of geography education, university of jember cover geosfera 61rev.pdf (p.1) 1.pdf (p.2) focus and scoperev.pdf (p.3-4) sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 291 the effect of google earth utilization on students' spatial thinking ability sri rahayu, murjainah, m. idris department of geography education, university of pgri palembang jln. jend. a. yani, lr gotong royong 9/10 ulu, 30116, palembang, indonesia email : srirahayu9796@gmail.com received 1 august 2019/ revised 2 december 2019/ accepted 5 december 2019/ published 12 december 2019 abstract the ability to think spatially in geography learning is essential, so it requires technology-based learning resources in the form of google earth, which can facilitate students in imagining or visualizing images in mind. in this regard, this study aims to determine the effect of the use of google earth on the spatial thinking abilities of students in the class x geography of sma pgri 2 palembang. this study used an experimental research method (posttest-only control design), because this design is suitable to use if the pre-test is not possible or pre-test can influence the experimental. the sample data collection technique uses purposive sampling, which is based on considerations or criteria that must be met by the sample used in the study. the sample in this study is class x ips 1 as the experimental class and x ips 2 as the control class. data collection techniques used documentation and tests. for data analysis techniques, normality test, homogeneity test, and hypothesis testing using the ibm spss statistics 20 formula for windows. based on the results of the study, the average value of the experimental class's superior post-test was 82.92, and the results of the posttest control class were 66.39. it shows that there are differences in the spatial thinking ability of the experimental group students who were treated using google earth during the learning process. the significance of the results of the posttest t-test from the two experimental and control groups was 0.000, and then the null hypothesis ho was declared rejected because based on the t-test criteria, the significance value was <0.05 or the sig (2-tailed) value of 0,000 was obtained <0.05. so it can be concluded that there is a significant influence between the use of google earth on the spatial thinking ability of students in the class x geography subject of sma pgri 2 palembang. keywords: google earth, spatial thinking ability, geography. 1. introduction spatial thinking is recognized as a collection of three cognitive skills about the nature and concepts of space (such as distance, closeness, and distribution), about the representation of spatial information (such as maps and graphs), and the process of spatial reasoning (such as decision making) (support committee for thinking thinking spatial; liu, et al. 2019). spatial geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.4 no. 3 (2019), 291-301, december, 2019 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v4i3.13350 accredited by the ministry of research, technology and higher education of the republic of indonesia, no. 30/e/kpt/2019. https://jurnal.unej.ac.id/index.php/geosi https://drive.google.com/file/d/1rsnvas6cuhowhl5bj87cl2l6k5dqz7s6/view sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 292 ability according to albert and golledge; setiawan (2015) consists of spatial visualization, spatial orientation and spatial relationships. in fact, according to golledge & stimson; halpern; aliman, mutia & yustesia (2018) spatial thinking is the ability of human reasoning to recognize spaces that can develop due to input, processing and output processes. according to the national research council, 2006; setiawan (2015) spatial thinking is one form of thought among other types of view, such as verbal, logical, statistical, hypothetical and so on. furthermore, spatial thinking is an important character in geography learning activities. according to setiawan (2015), the study of geographic phenomena not only explains the existence of a phenomenon and the process of occurrence of this phenomenon on the surface of the earth but also the shape, size, direction, pattern of phenomena and their relationship with other phenomena. according to hidayat et al (2017), spatial thinking is a basic skill that can be accessed by everyone to different degrees in different contexts to solve problems in various contexts. meanwhile according to lee jongwon & bednars s robert; hidayat et al (2017) spatial thinking requires three related components, namely: the concept of space, the method used to represent spatial information, and the process of spatial reasoning. therefore, these three components are interrelated, mutually supportive and inseparable. the importance of spatial thinking can and must be taught at all levels in the education system. the goal is that each individual has the good spatial ability. gersmehl & gersmehl; oktavianto et al (2017) define spatial thinking as an ability that can be used by a geographer to analyze spatial relationships on earth. this ability will be very useful for students when deciding or making decisions from things that are very simple to complications related to space or location. when someone travels, he must know about distance and direction, so he can predict the time of arrival and not get lost. therefore, in learning geography, it is very important to emphasize spatial thinking, not only information about geographic phenomena, but students must have the ability to analyze spatial aspects, because the ability of spatial thinking in geography can affect students' ability to imagine or visualize images in the mind. based on data obtained from geography subject teachers in class x sma pgri 2 palembang, said that there were 4 class x social studies, where 60% of students in grade 10 had reached kkm with a score above 70, while 40% of students had not yet reached kkm. core competencies of geography subject in class x are managing, reasoning, and presenting in the abstract realm related to the development of what is learned in schools independently and can sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 293 use methods that are following curriculum rules. the basic competency used is to present the results of an analysis of the relationship between humans and their environment as the influence of atmospheric dynamics in the form of narratives, tables, graphs, graphs, illustrated images, and concept maps. then, indicators identify the type of inland waters, identify the use of inland waters, analyze conservation of inland waters, and watersheds. a watershed is a part of the earth's surface where water flows into certain rivers. in other words, a watershed is a rainwater reservoir that enters the watershed. watershed consists of 3 types, namely upstream, downstream and middle watersheds. in this connection, the teacher's role becomes important in learning to improve students' spatial abilities. the teacher is expected to be able to provide stimulation to students and interesting innovations and learning strategies are needed so that students can understand the geography concept about watersheds in southern sumatera, so that it can be understood more easily. recognizing that all technological facilities are needed, that can improve spatial thinking skills. so, it is needs technology that can support learning, using google earth. according to yousman (2008) google earth is an interactive mapping application released by google. google earth displays globe maps, topography, satellite photos, terrain that can be overlaid with roads, buildings, locations, or other geographical information. with google earth, we can plan trips, find tourist attractions, motorbikes, restaurants, hotels, hospitals, schools and more where we can get latitude and longitude coordinates. google earth can display low-resolution satellite photos that depict mountains, seas, forests, to high-resolution satellite photos that can depict objects such as roads, office homes. for certain areas that are already equipped with 3d building views. google earth provides an application for educators to display images of the earth visually. google earth also provides opportunities for students to see every side of the world. google earth helps students see distance and other geographical features. google earth is a free software package available for anyone who has a computer and an internet connection. this is an online resource available in the classroom and can be used by students at home. studies have shown that the use of online resources has helped increase students' understanding of key concepts and skills while also helping students gain confidence in their knowledge of geographical problems (solem and gersmehl; cuviello, 2010). according to bodzin et al. (2009); oktavianto et al. (2017) geospatial-based sites such as google earth can accelerate the improvement of spatial thinking skills in a variety of students. this is in line with the main material in this study that sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 294 requires digital technology facilities with easy access to information that is relatively faster without having to be present directly on the object to be addressed. the novelty in this study is analysis regarding google earth for conveying material taught specifically in geography. bearing in mind that many students find it difficult to understand an object or map if they learn to use ordinary maps, so the teacher is required to more attention for help students develop their own spatial thinking skills so as to achieve the expected competencies in learning geography in school. based on these problems, the reference for researcher to conducting research on the effect of using google earth to students' spatial thinking abilities. the purpose of this study was to determine the effect of using google earth to students' spatial thinking abilities. 2. methods the method used in this study is the experimental research method in the form of posttest only control design (sugiyono, 2010). data collection of this study uses test techniques. the test is given at the end of the meeting, used to obtain data on the ability of students to solve the questions given after using google earth. the test used is a multiple-choice form test.validity test using the ibm spss statistics 20 formula for windows. testing the validity of the instrument in this study uses the validity of construction, namely the pearson productmoment correlation formula as follows: rxy = 𝑛𝛴𝑋𝑌‐(𝛴𝑋)(𝛴𝑌) √{(𝛴𝑋2)‐(𝛴𝑋)2}{𝑁𝛴𝑌2‐(𝛴𝑌)2} (1) (arikunto,2010) reliability testing is done using ibm spss statistics 20 for windows with the cronbach's alpha model, which is measured on the cronbach's alpha scale 0 to 1. to find the reliability value using alpha formula: (arikunto,2010) the data analysis technique used in this study is the statistical test parameter t (t test). the t statistical test is used to test the rejection or acceptance of the null hypothesis, provided that the sample is homogeneous and normally distributed. the value of the geography learning test (2) http://1.bp.blogspot.com/-1itfstw-fbw/udeqtuw1ayi/aaaaaaaaaio/47lrmv2ymcw/s1600/reliabilitas.jpg sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 295 results obtained in the experimental class and the control class. according to sudjana test t-test (2005) as follows: (sudjana, 2005) normality test to analyze data by testing whether the data obtained is normal or not. data normality test needs to be done to find out whether the data analyzed is normal or not. data is said to be normal if the km price is located between 1 (-1≤ km ≤ 1). the data created in the frequency distribution table is tested for normality using the curve slope normality test formula as follows: km = 𝑋 − 𝑀𝑜 𝑆 homogeneity test data is done to prove the similarity of group variance, where samples taken are from the same population. to test the sample using the bertlett test with the chi-square equation. at the beginning of the meeting in the experimental class, students look confused in the following learning because the delivery of learning has never used google earth so that the beginning of learning researchers are active to guide students during the learning process by using google earth. after researchers explain how to use google earth, students become more active. the use of google earth helps researchers in delivering material well so that students become more focused on visualizing images. the water inland described by researchers in the form of rivers, lakes and swamps makes students interested in listening carefully. this is proven by the good results when students work on the posttest. while in the control class applying learning using conventional methods, in this class, students are asked to pay attention to the explanations of researchers with conventional methods, without using google earth, so that makes students bored and tired. students only get information from the teacher. some students ask questions, but many students are passive during the learning process. so it is not enough to pay more attention to assist students in developing sgab = √ (𝑛₁‐1)𝑆1 2+ (𝑛21)𝑆1 2 𝑛1+𝑛2−2 (3) (4) sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 296 spatial thinking skills and cause understanding of the material in the control class can not be maximized, this is seen from the results of students' posttest. this study uses two data collection techniques, namely, documentation and test techniques. the test questions that were given in the sample class had previously been tested for validity and reliability. test questions are given as many as 20 multiple-choice questions. the question is made by adjusting the indicators of spatial thinking ability consisting of location, condition, connection, comparison, aura, region, hierarchy, transition, analogy, pattern, spatial association. 3. results and discussion tests are given at the end of learning (posttest). meanwhile, documentation is used to support research data in the form of test result data and learning documentation. the results of the experimental and control class posttest can be seen in figure 1 below. figure 1. post test results students in the experimental and control classes based on figure 1. shows that in the experimental class and the control class gets a variety of values. however, in the control class, many students received grades at intervals of 5565, and none of the students received grades at intervals of 88-98. students in the experimental class get the highest value of 95, the lowest value of 65, and the average value obtained by the 0 5 10 15 20 25 55-65 66-76 77-87 88-98 experiment class control class sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 297 experimental class is 82.92. meanwhile, the control class test results with the highest score of 85, the lowest score of 55 and the average value obtained by the control class is 66.39. this is proof that the results of experimental class tests using google earth can affect student learning outcomes. in contrast, verma and estaville; jo and hong (2018) revealed that currently, there is no evidence of empirical research that shows that learning geography helps students develop spatial thinking skills. to find out whether the data is normally distributed or not, then the normality test is done using shapiro wilk, with the help of a computer with the statistical package for social science (spss) version 20. the data criteria are said to be normal if significance> 0.05. the results of data processing using the shapiro wilk technique can be seen in table 1 below. table 1. distribution of normality test results with the shapiro wilk test class kolmogorovsmirnova shapiro-wilk statistic df sig. statistic df sig. spatial thinking ability posttest experimen .131 36 .123 .945 36 .073 posttest control .150 36 .039 .928 36 .021 a. lilliefors significance correction based on table 1. the normality of the test results above is known that the results of the experimental group posttest significance value (sig) on the posttest score of the experimental class 0.73> 0.05 while the posttest score of the control class 0.21> 0.05. this shows that the data is normally distributed because the significance is 0.21> 0.05. therefore it can be concluded that both are normally distributed. next, to find out whether or not some of the variants of the research data were tested for homogeneity. criteria for decision making 0.05. in this homogeneity test, the researchers used spss 20. the results of the homogeneity test can be seen in the following table. sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 298 table 2. distribution of homogeneity test results levene statistic df1 df2 sig. spatial thinking ability based on mean .000 1 70 .983 based on median .013 1 70 .908 based on median and with adjusted df .013 1 69.59 1 .908 based on trimmed mean .005 1 70 .942 based on table 2, it is known that the results of the posttest of the experimental and control groups, obtained a significance value of 0.983> 0.05, it can be concluded that the variants of the experimental and control groups are homogeneous. after observing the characteristics of the variables that have been studied and the requirements of the analysis, then testing the hypothesis. for the purposes of the hypothesis, inferential statistics are used with the help of spss version 20, namely t-test statistics. the decision-making criteria are as follows: if sig> 0.05, ho is accepted if sig <0.05 then ho is rejected hypothesis testing, the step taken is to analyze the results of the t-test. the results of the t test analysis can be seen in table 3. table 3. experiment and control posttest t-test results mean std. deviation df thitung ttabel sig (2tailed) conclusion experimen 82,92 7.962 70 8.651 1.666 .000 h0rejected control 66,39 8.247 ardyodyantoro's research results (2014) on "utilization of google earth in geography learning to improve student learning outcomes of class x high school widya kutoarjo" shows that learning using google earth media is effective in improving learning outcomes in geography. learning outcomes with google earth media are higher than learning outcomes with lectures. the mean value of learning outcomes with google earth media 83,397, while with lectures 78,348. improved learning outcomes indicated by the achievement score of 0.68 sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 299 improvements in the experimental class and 0.58 in the control class. p value of learning outcomes 0.01 <0.05, then ha is accepted and h0 is rejected. this proves that the use of google earth media in learning geography is effective in improving the learning outcomes of class x high school students widya kutoarjo. isnaini's research results (2015) on "comparative use of google earth media with digital maps on class xi ips fauna distribution materials in state high school 1 semarang" shows that (right-side t-test) shows a tcount of 2,433 and a table of 1.67 with a significance level of 5 % and dk = 31 + 31 2 = 60, because tcount> ttable, it was concluded that the experimental group taught using the google earth media the learning outcomes were better than the control group taught using digital map media. this means that the research hypothesis was accepted. meanwhile, nofirman (2018) conducted a study on the geographical spatial ability of class xii students of sma negeri 6 bengkulu city showing that the results of data processing it was found that the spatial abilities of class xii students at sman 6 bengkulu there is in the largest group (43.55%). the potential geographic spatial ability of class xii students in city 6 of bengkulu city on the group with the largest number of 38.71%. the average position is in the highest score group. furthermore, the results of ervina, asyik and mizwar (2012) research on "the influence of the use of google earth and maps media in the improvement of geography learning outcomes in material of southeast asian regional at sma negeri 14 bandar lampung in the academic year 2011/2012" shows that there are differences in the increase in results student learning and the value of student learning outcomes on the use of google earth media is higher than the use of media maps. meanwhile, the results of oktavianto, et al (2017) research on "the effect of google earth assist project based learning to spatial thinking ability" shows that google earth-aided project-based learning has a significant effect on students' spatial thinking skills. in addition, also found several advantages of google earth-based project-based learning, including: (1) encouraging students to be solve real problems through project activities, (2) students are more active in learning, (3) student performance in completing projects are more organized, (4) students have more freedom to complete projects, (5) students are motivated to compete to produce the best products, and (6) students experience increased spatial thinking skills. sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 300 the use of google earth in learning geography has an effect to students on spatial thinking ability indicated by an increase in student learning outcomes after using google earth while learning. in addition, student activities and attention increase and students' difficulty in understanding spatial decreases. 4. conclusion based on the results of research, it can be concluded that there is a significant influence between the use of google earth on the spatial thinking abilities of students in class x sma pgri 2 palembang. judging from the average value of the experimental class's posttest is 82.92 and the posttest of the control class is 66.39. this has been proven by examiners that the t-count ha is accepted, that's indicating that there are differences in the spatial thinking ability of the experimental group students who are treated using google earth during the learning process. the significance of the results of the posttest t-test from the two experimental and control groups was 0.000, and then the null hypothesis ho was declared rejected because, based on the t-test criteria, the significance value was <0.05 or the sig (2-tailed) value of 0,000 was obtained <0.05. conflict of interest the authors declare that there is no conflict of interest with any financial organization regarding the material discussed in the article. references aliman, mutia, & yustesia. (2018). integrasi kebangsaan dalam tes berpikir spasial. jurnal geografi fkip ump , 82-89. arikunto, s. (2010). prosedur penelitian suatu pendekatan praktik. jakarta: rineka cipta. ardyodyantoro, gatty. (2014). pemanfaatan google earth dalam pembelajaran geografi untuk meningkatkan hasil belajar siswa kelas x sma widya kutoarjo. skripsi. program studi pendidikan geografi fakultas ilmu sosial universitas negeri yogyakarta. cuviello, matthew p. (2010). evaluating google earth in the classroom. new york : center for teaching excellence sri rahayu et al / geosi vol 4 no 3 (2019) 291-301 301 ervina, e., asyik, b., & miswar, d. (2012). pengaruh penggunaan media google earth dan peta terhadap peningkatan hasil belajar geografi. jpg (jurnal penelitian geografi), 1(1). hidayat, k. n., & fiantika, f. r. (2017). analisis proses berfikir spasial siswa pada materi geometri. prosiding si manis (seminar nasional integrasi matematika dan nilai islami) , 385-394. isnaini, n. (2018). komparasi penggunaan media google earth dengan peta digital pada materi persebaran fauna kelas xi ips di sma negeri 1 semarang. jurnal geografi: media informasi pengembangan dan profesi kegeografian, 12(1), 52-61. jo, i., & hong, j. e. (2018). geography education, spatial thinking, and geospatial technologies: introduction to the special issue. international journal of geospatial and environmental research, 5(3), 1. liu, r., greene, r., li, x., wang, t., lu, m., & xu, y. (2019). comparing geoinformation and geography students’ spatial thinking skills with a human-geography pedagogical approach in a chinese context. sustainability, 11(20), 5573. doi:10.3390/su11205573 nofirman, n. (2019). studi kemampuan spasial geografi siswa kelas xii sma negeri 6 kota bengkulu. jurnal georafflesia: artikel ilmiah pendidikan geografi, 3(2), 11-24. oktavianto, d. a. (2017). pengaruh pembelajaran berbasis proyek berbantuan google earth terhadap keterampilan berpikir spasial. jurnal teknodik, 21(1), 059. patterson, t. c. (2007). google earth as a (not just) geography education tool. journal of geography, 106(4), 145–152. doi:10.1080/00221340701678032 setiawan, i. (2016). peran sistem informasi geografis (sig) dalam meningkatkan kemampuan berpikir spasial (spatial thinking). jurnal geografi gea, 15(1). doi:10.17509/gea.v15i1.4187 sugiyono. (2010). metode penelitian pendidikan pendekatan kuantitatif, kualitatif, dan r&d. bandung: alfabeta. sudjana. (2005). metode statistika. bandung: tarsito. yousman, y. (2008). google earth. yogyakarta: c.v andi. la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 280 mapping of subsurface geological structure and land cover using microgravity techniques for geography and geophysic surveys: a case study of maluri park, malaysia la ode nursalam1, a arisona2, ramli1, la harudu1, eko harianto3, sitti kasmiati4, fahrudi ahwan ikhsan5, andri estining sejati6 1department of geography education, halu oleo university, kendari, sulawesi tenggara, 93232, indonesia 2department of geology engineering, halu oleo university, kendari, sulawesi tenggara, 93232, indonesia 3department of agrobusiness, kendari terbuka university, kendari, sulawesi tenggara, 93232, indonesia 4department of archeology, halu oleo university, kendari, sulawesi tenggara 93232, indonesia 5department of geography education, jember university, east java, 68121, indonesia 6department of geography education, sembilanbelas november kolaka university, kolaka, sulawesi tenggara 93561, indonesia email: laodenursalam77@gmail.com received 13 september 2019/ revised 9 november 2019/ accepted 22 november 2019/ published 1 december 2019/ available online 25 november 2019 abstract a microgravity investigation on bedrock topography was conducted at maluri park reference level in kuala lumpur, malaysia. the study aim to mapping the near-surface structure and soil and land cover distribution for geography and geophysics surveys. two types of crosssection modeling of the residual anomaly generated the maluribouguer anomaly model for site-1 and site-2 at maluri park. the 2d microgravity models produced the contour map, displaying the characterization due to density contrast in rock types while mapping the subsurface geological structure at different depths. moreover, a synthetic model was initiated with the assumption of lateral distance on the left and right sides taken at 50 m and a depth of 60 m. the results of modeling confirmed that the soil and rock type composition on both models site tests are topsoil (1.1 to 1.92 g/cm3), soil (1.8 g/cm3), clay (1.63 g/cm3), gravel (1.7 g/cm3), sand (2.0 g/cm3), shale (2.4 g/cm3), sandstone (2.76 g/cm3), and limestone (2.9 g/cm3). the 2d gravity modeling using two model site tests obtained a correspondence with the observed microgravity data. keywords: bouguer anomaly, limestone, microgravity, soil structure, topography. 1. introduction the survey area is underlain by limestones. the survey area in maluri park in eastern kuala lumpur, malaysia, near cheras. the survey location can be seen in the fig.1. it is commonly understood that limestones pose a certain threat to soil structures. as limestones are generally easily dissolved by water, voids and cavities are common internal structures of this type of rock (arisona et al., 2018). identification of such structures can help to ensure sustainability, especially for land cover and strategic structures. the area under study is geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.4 no. 3 (2019), 280-290, december, 2019 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v4i3.13738 accredited by the ministry of research, technology and higher education of the republic of indonesia, no. 30/e/kpt/2019 la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 281 planned for both establishment of transportation facilities and an economic and urbanization projects. the motivation of the present work is to help reducing the threats to soil structures in the mentioned projects. a significant application of geography and geophysical methods in mapping practice is to determine the bedrock and the characteristics of soil structure (grandjean, 2009; hiltunen, 2012). the gravity method is widely used for geography and geophysical surveys, especially in the detection of subsurface geological features and land covers. gravity' anomalities controlled by the lateral variation of densities or other words, the lateral density contrast. generally speaking, modeling of gravity data in 2d sections is useful for the determination of the depth of various features and can be done by either forward or inverse algorithms. in similarity with other geography and geophysical methods, the interpretation of gravity data is non-unique because many possible models could result in the same gravity anomaly. constraints from borehole data can help to reduce the uncertainties greatly. the kind of useful information for this purpose in the depths of rock boundaries and rock densities. the irregular structure of limestones and land cover produces gravity anomalies, which are our target for gravity analysis and modeling. in other words, gravity anomalies are modeled to determine the subsurface structure and land cover problems. figure1. survey location: maluri park, malaysia la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 282 as mentioned, the problem targeted by the present work is the presence of voids or cavities underneath the sites under the survey area. these cavities are either empty (i.e. filled with air or water) or filled with loose sediments. in both cases, the lateral density contrast is large, which favors the application of microgravity techniques. as the density of the cavities or voids is low compared to the surrounding host limestone bed, the anomalies will be negative. hence, the analysis focusses on these negative anomalies for modeling their dimension and depth(zabidi et al., 2011). in this study, we perform microgravity modeling using physical parameters from real data that are robust. the bedrock within the survey area is a limestone formation. limestone formations are well known for their highly unusual karstic features (tan, 2005). as a consequence, the depth of the limestone bedrock is highly irregular. the overburden soils above limestone formation are mainly silty sands with significantly variable thickness due to the irregular topography of the limestone bedrock (zabidi et al., 2011). the soil thickness in the survey is that trends show variations of overburden soil between 3 m and 5 m. the soil is comprised mainly of sandy silt with embedded layers of soft clay. moreover, the residual soil above the limestone bedrock is mainly loose fine-grained materials. the soil color in the survey area is light greyish-brown and mostly sandy. they are described and named according to the grain size classification as silty clay, clayey silt, silty sand, and sandy clay. some are identified as fill and slime materials. they are essentially loose and soft material and are probably transported materials (yusoff et al., 2016). the purpose of the research is mapping the near-surface structure and soil and land cover distribution for geography and geophysical surveys. 2. methods microgravity data from the survey area were processed using surfer ® 13 software , which reduced bouguer anomaly values at each station of the microgravity survey. fig.2 shows that the modeling has been well constrained because the parameters required to obtained the bedrock topography were well defined from the borehole data (figure 3). according to samsuddin (2003), modeling enables the determination of the presence of ridge and valley features. additionally, the modeling indicates that there can be smaller features such as cracks or steep, narrow valleys within the ridges. this result explains the differences in the limestone depth in nearby boreholes. the post-processing procedure checked the microgravity instrument corrections for latitude and longitude, diurnal variations, and instrument drift using base station polynomial drift values and relative elevation. this procedure is merged with the respective gravity station la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 283 topographic survey data, which is being modeled. the same density value of 1.8 g∕cm3 was used to calculate the bouguer correction for all survey data sets. table 1 presents the density of rock types that were reviewed in this study. table 1. rock types density values rock type density range (gr/cm3) average (gr/cm3) source sediment rock overburden (topsoil) 1.92 soil 1.20 to 2.40 1.92 clay 1.63 to 2.60 2.21 gravel 1.70 to 2.40 2.00 sand 1.70 to 2.30 2.00 telford et al.,1990 sandstone 1.61 to 2.76 2.35 shale 1.77 to 3.20 2.40 limestone 1.93 to 2.90 2.55 dolomite 2.28 to 2.90 2.70 subsequent processing was the elevation correction to address the variation in data points due to the topography. this evaluation is necessary for the 2d geotechnical modeling (pringle,2012). furthermore, the removal of regional values (low frequency and high amplitude) to express the residual anomalies (high frequency and low amplitude) was performed. the modeling of the residual anomaly was generated in cooper ™ grav2dc v.2.10 software. however, qualitative interpretation using geological and geography maps are used only as additional information. according to (amaluddin et al., 2019) geological formation or structure for instance can use the maps which poured in the data field or 2d. the qualitative interpretation explains the anomaly by geological and geophysical information. on this basis, the geological structure and distribution of masses of different densities may be delineated. the gravity anomaly in the study area might be generated in the field due to the following factors; variations in the thickness, density differences of the subcrustal matter (crustal thickness), density variation within the basement rocks, thickness variations of the sedimentary rocks, and density variations within the sedimentary rocks. the difference in density values can be correlated to divergent material types such as soil, rock, and voids (samsuddin,2003). la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 284 figure 2.a 2d microgravity contour map at maluri station compiled from topography and bouguer anomalies at site test-1 and site test-2. the microgravity countur color information above is on the right and adjusted the rock conditions in the study site. figure 3. borehole (bh) logs showing bedrock rl at maluri station la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 285 3. results and discussion the bedrock topography for the maluri station was compiled from bh results, secant bored piles (sbp) and kingpost, and the gravity model. the density values were not uniform, which showed that some lines in the study area lack good control. the bedrock topography map is shown in fig. 4was correlated with the 31boreholes, where the average topography height of about 39.91mwas determined (figure 3). two sinkholes were detected in the study area in the southeast of the survey area. the sinkholes occurred near the lrt bridge, which might occur under the devotional of clay construction site tests. the cause of the sinkhole is attributed to the dissolution of the limestone bedrock and the subsequent raveling of the overburden soil cover. figure 4.maluri station bedrock topography model. la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 286 the gravity pattern correlates subsurface topography of bedrock (tajuddin and lat, 2004). the results of this investigation confirmed earlier borehole results which indicated thepresence of the cavities (samsudin, 2003). however, the uneven distribution and clustering of the data necessitate the use of an interpolation algorithm to create a uniformly spaced grid. all data processing in the contour map was generated with cooper ™ surfer ®13software. the survey results were represented in contour maps for delineating anomalies varying from negative to a positive value. according to kamal et al. (2010), negative values are interpreted as low density subsurface layers and for the possibility of the existence of cavities. no density measurements were made during this work. the density values used are based on the investigations of soils and foundations sdn. bhd., employing a superficial deposit density of1.8 g/cm3. table 1presents the density of rock types that were reviewed in this study. fig 5 shows the results of the microgravity data of maluri site for the bouguer, regional, and residual anomalies. fig. 5a displays the tendency of response towards positive anomalies, and yet it was not significant to influence the gravity anomaly around model site test-1. the small variation in gravity responses may be due to the effect of the distribution of rainwater around the area and the inhomogeneity in the soil types. fig.5b shows the gravity anomaly at model site test-2, characterized by negative values, probably due to the inhomogeneous geo-materials consisting of a mixture of clay and silt with grain sizes, which is from fine to medium. additional geological and geography mapping showed that the complexity of the subsurface profile regarding the geomaterials, geological structure, and water seepage influences the contrast in the zones, thus resulting in some of the inconsistency of the gravity values. figure 5. results of gravity field measurements at maluri site for (a) profiles site test-1 and (b) profiles site test-2 using extracted techniques from anomaly contour maps. la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 287 figs. 5a and 5b show results of the 2d geotechnical model generated with cooper ™ grav2dc v.2.10 software. the two curves for the model at the sites confirmed the product of a misfit of4.44 % and 3.18 %between calculated curves with observed curves. furthermore, modelcalculated gravity confirmed density contrast at both model sites as shown in table 2. table 2. estimated density contrast from gravity field measurements at maluri site for (a) profiles site test-1 and (b) profiles site test-2 model site test-1 model site test-2 density (g/cm3) rock type density (g/cm3) rock type 1.80 soil 1.10 dry density 1.63 clay 1.80 clayey, sandy silt 2.00 gravel 1.20 clay 1.70 sand 1.92 overburden (topsoil) 2.40 shale 2.76 sandstone 2.90 limestone rms error = 4.44 % rms error = 3.18 % fig. 5 shows the zones associated with negative anomalies. an important to observation is the decrease of soil compactness due to the cavity and existence of groundwater around the area. in other words, the model site test encountered a unsewed soil at an adepth of about 5 to 10 m filled with water. additionally, different factors include changes in gravity values caused by the dynamics around the observation points, such as variations in the depth of the groundwater level and land subsidence. the occurrence of strain and shear failure due to loading at the top layer of soil could be another reason. in general, loading in the ground layers produces strain in the sediments (fig.5), a decline in groundwater level at a depth of 15 m to 35 m). strain observed in this area could be due to variation in the composition of the soil and water in the cavities. the loading effect is often referred to as consolidation. as a result, regional anomalies when compared with bouguer anomalies in a horizontal plane, the obtained residual anomalies are still significant. this difference is probably caused by the material homogeneity and similarity of the geological structure. the residual anomaly values are not affected by topography. anomalies appear in the northeast and northwest. that result in line with the research (arisona et al, 2018), difference assumed that this contrast is affected by density variations of the host rock and the possibility of rock density. the difference location implies that the anomaly is a result of different depth to the bedrock, different thickness of the overlying material la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 288 the bouguer anomaly map result characterized the density contrast due to rock type and mapped subsurface geological structure and land cover at different depth. the research (samsuddin, 2003) describes that the modeling enables the determination of the presence of ridge and valley features. additionally, the modeling indicates that there can be smaller features such as cracks or steep, narrow valleys within the ridges. the result is eight rock types from the profiling from different and depth. it is shown that the mapping for a geophysics survey can obtain information matching with the mapping aim. the research (wanjohi, 2014) geophysical field mapping is the process in selecting the interest area, then identifying all geophysical aspect matching with the mapping aim. the mapping aim in the research (wanjohi, 2014) is to understand all physical parameters of a geothermal field. the survey results were represented in contour maps for delineating anomalies varying from negative to a positive value. mapping in a geophysical method in the research (georgsson, 2009) uses contouring to see the object interest phenomena. the contour or line data combined with other data like a polygon in the form of coloring. microgravity method solves problems targeted by the present work is the presence of voids or cavities underneath the sites under the survey area. according to the previous research (tuckwell, grossey, owen, & stearns, 2008) that micrografity establihed as technique for detection natural or man-made cavilities. in the case mounchel parkman on behalf of herfordshire country council, a doline had opened up within a school playground. that example in the natural voids in the limestone bedrock. 4. conclusion the result obtained from the grav2dc v.2.10 software correlates the model-calculated gravity, and the corrected gravity data in site tests have minimal percentage errors. the results of modeling showed that there are eight rock types from the gravity profiles; topsoil (1.1 to 1.92 g/cm3), soil (1.8 g/cm3), clay (1.63 g/cm3), gravel (1.7 g/cm3), sand (2.0 g/cm3), shale (2.4 g/cm3), sandstone (2.76 g/cm3) and limestone (2.9 g/cm3). the utilization of extracted technique characterized the density contrast due to rock type and mapped subsurface geological structure and land cover for geography at different and depth. models of microgravity distribution in the ground could be useful for the mapping of variations in soil composition. the changes in gravity anomaly observed throughout the sections were due to the heterogeneities in the composition of the subsurface materials and density contrasts in the study area. the information on the background lithology is paramount for an acceptable interpretation in a microgravity modeling of soil structure. la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 289 conflict of interest the authors declare that there is no conflict of interest with any financial organization regarding the material discussed in the article. acknowledgment this study is a collaboration with a graduate student of the physics and geography education study program, at halu oleo university. we thank to a colleague at halu oleo university, which has been giving us for their guidance, comment, and correction on this paper. references amaluddin, l. o., rahmat, r., surdin, s., ramadhan, m. i., hidayat, d. n., purwana, i. g., & fayanto, s. (2019). the effectiveness of outdoor learning in improving spatial intelligence. journal for the education of gifted young scientists, 7(3), 667–680. https://doi.org/10.17478/jegys.613987 arisona,a., mohd n., amin e.k., &abdullahi, a.(2018).assessment of microgravity anomalies of soil structure for geotechnical 2d models.journal of geoscience, engineering, environment, and technology (jgeet)3(3), 151-154. georgsson, l.s. (2009). geophysical methotds used in geothermal exploration. presented at exploration for geothermal resources, 1-22 november 2009, 1-16. grandjean, g. (2009). from geophysical parameters to soil characteristics.florida: report n°brgm/fp7-digisoil project deliverable 2.1, final reportdepartment of civil and coastal engineeringuniversity of florida. hiltunen, d.r., hudyma,n.,tran,k.t.,&sarno,a.i. (2012).geophysical testing of rock and its relationthipsto physical properties.florida:final reportdepartment ofcivil and coastal engineeringuniversity offlorida. kirsch,r. (2006).groundwatergeophysics, atool for hydrogeology.new york: springer. kamal,h.,taha,m.,&al-sanad,s. (2010). geoenvironmental engineering and geotechnics, geoshanghai 2010 international conference. (accessed 02.03.17) lilie, r.j. (1999).whole earth geophysics: an introductory textbook for geologists and geophysicists. new jersey:prentice-hallinc. pringle, j.k., styles, p., howell, c.p.,branston, m.w., furner, r., &toon,s.m. (2012). longterm time-lapse microgravity and geotechnical monitoring of relict salt mines, marston, cheshire, uk. geophysic77(6), 165-171. samsudin, h.t.(2003).a microgravity survey over deep limestone bedrock.bulletin of geological society of malaysia4(6), 201-208. la ode nursalam et al / geosi vol 4 no 3 (2019) 280-290 290 tan, s.m. (2005). karsticfeatures of kualalumpur limestone. bulletin of the institution of enginnermalaysia 4(7), 6-11. tajuddin, a.&lat, c.n. (2004).detecting subsurfacevoids using the microgravity method, a case study from kualalipis, pahang.bulletin of geological society of malaysia 3(48), 31-35. tuckwell, g., grossey, t., owen, s., & stearns, p. (2008). the use of microgravity to detect small distributed voids and low-density ground. quarterly journal of engineering geology and hydrogeology, 41(3), 371–380. https://doi.org/10.1144/1470-9236/07-224 wanjohi, a.w. (2014). geophysical field mapping. presented at exploration for geothermal resources, 2-23 november 2014, 1-9. yusoff , z.m., raju,g. &nahazanan, h.(2016).static and dynamic behaviour of kualalumpur limestone. malaysian journal of civil engineering special issue vol.28 (1), p.:18-25. zabidi, h. & de freitas, m.h. (2011).re-evaluation of rock core logging for the prediction of preferred orientations of karst in the kualalumpur limestone formation. engineering geology, 117(3-4), p.: 159–169. 196 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 soil zonation and the shaking table test of the embankment on clayey soil ripon hore1*, sudipta chakraborty2, md. fayjul bari1, ayaz mahmud shuvon2, and mehedi ahmed ansary1 1department of civil engineering, bangladesh university of engineering and technology (buet), dhaka, 1000, bangladesh 2bangladesh network office for urban safety (bnus), bangladesh university of engineering and technology, dhaka, 1000, bangladesh *corresponding author : riponhore@gmail.com received 22 march 2020/ revised 5 june 2020 / accepted 21 june 2020/ available online 28 june 2020 abstract the main objective of this research was to model the zonation of wrap faced embankment on soft clay foundation, by applying a shake table test. also, to investigate the dynamic behaviors of clay soil, such as acceleration amplification, displacement and pore water pressure of wrap faced embankment. this was done with respect to changes in frequencies of 1 hz, 3 hz, 5 hz, 10 hz, 12 hz and 15 hz respectively. constant acceleration (0.1 g) and surcharge (19 kg) were applied by using a laminar box, placed on a shake table testing machine. the main elements of this research were the laboratory test, which was used for preparing reconstitute soil samples, and wrap faced embankment with frequency arrangement. after applying all test parameters, dynamic parameters were increased by rise in elevation with respect to frequency. the result shows that the maximum dynamic parameters were found at the frequency of 10 hz. it is beneficial to the relative performances of the wrap faced embankment, which is the updated design parameter. keywords: seismic; clay soil; frequency; shake table test; wrap faced; soil zonation 1. introduction the soil-foundation formed from soft clay becomes the focus of seismic engineering. in some cases, the foundation on soft clay is creating a problem for the design and construction of any type of structure. in bangladesh, the southern part of the country, an geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.5 no. 2 (2020), 196-209, august, 2020 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v5i2.17873 accredited by the ministry of research , technology , and higher education of the republic of indonesia, no. 30/e/kpt/2019. https://jurnal.unej.ac.id/index.php/geosi https://drive.google.com/file/d/1rsnvas6cuhowhl5bj87cl2l6k5dqz7s6/view 197 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 excessive amount of soft clay is found in the khulna and bagerhat districts and also in the surroundings of dhaka city (hore et al., 2019). bangladesh has the largest delta in the world. therefore, very large alluviums are deposited on its surface. the oldest deposits are the barend, madhapur and lamaicregion clay. the sediments deposited are not evenly distributed throughout the country. at the northern part, it is about 128 m thick and this is where granite is extracted for construction purposes. conversely, the thickness was gradually increased towards the south. at the centre, the capital city, dhaka of bangladesh has a sediment covering of over 22 km (al zaman & jahan monira, 2017; alam & islam, 2009; bazlar rashid et al., 2018; haque et al., 2013; hore et al., 2019). earthquakes or any seismic effects may cause a major damage on this type of soil. the earthquakes or any sinusoidal waves (railway vibration) have often created significant problems during the design and construction of the embankment. moreover, it may be damaged due to the softening of soilfoundation (ering & sivakumar babu, 2020; krishna & latha, 2007; kumar et al., 2020; latha & krishna, 2006; zhou et al., 2020). the soil and structure interaction (ssi) system was considered in this experimental study, to simulate actual soil-foundation (bullock et al., 2019; çelebi et al., 2019; he & jiang, 2019; hore et al., 2020; srilatha, et al., 2013; srinivasan et al., 2016). soft soils in dhaka, bangladesh were used to build the soilfoundation. krishna & latha (2007) determined the result of shake table tests on geotextilereinforced wrap faced soil-retaining walls. a total number of 9 model tests were described, such as development, experiment methodology, and outcomes. srilatha et al., (2013) illustrated the effect of frequency of base shakes on the dynamic response of unreinforced and reinforced soil slopes. xiao et al., (2014) investigated the earthquake response of a slurry wall and presented a minimized scale shake table test. in this research, soil-cement-bentonite (scb) was evaluated. this is also a common type of slurry wall. reinforced soil is only slightly damaged during the recent seismic disasters in japan (suzuki et al., 2015). the seismic behaviour of this soil wall is formed from clay, cement-treated clay, sand based on shake table tests and results of the pull-out. fleming et al., (2016) conducted a contemporary phenomenon on soft clay soil, to determine the consequences of soil improvement on the seismic resistance of piles. yazdandoust (2017) conducted a recent laboratory test to assess the character of 1-g shake table tests. many researchers also performed a shake table test to study the seismic responses of different soil structures on soft soil (beskhyroun et al., 2011; 198 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 edinçliler & toksoy, 2017; hamayoon et al., 2016; hassan & pal, 2018; helwany et al., 2017; madhavi latha & manju, 2016). few researchers focused on shake table tests on soft clay soil in south asian regions. in bangladesh, there have been limited studies related to shake table tests of wrap faced retaining walls, on soft clay. therefore, the purpose of this research was to study the effects of frequency on reinforced soil wall models. the objectives of this study were to: a) measure the response of wrap faced and reinforced soil retaining wall, which was subjected to dynamic loading through shake table; b) investigate the acceleration amplification, deformation, and pore water pressure with respect to frequency response; c) to zone the wrap faced embankment on the soft clay foundation, by applying a shake table test, and d) draw the contour maps using plaxis 3d software. 2. methods 2.1 area of study the clay soil sample was collected at a depth of 1.5 m below the existing ground level from the dhaka city of bangladesh as shown in figure 1. these homogeneous, stiff, reddish brown samples were at first oven-dried. subsequently, the dry lumps were then powdered gently by using a wooden hammer. it was finally sieved through a 200 standard sieve to obtain clean clay-like soil powder. this type of clay is very dominant in catchment areas (areas around ariver) with a thickness from 1 to 20 m. a total number of 1000 drill holes were carried out for the standard penetration test (spt) in and around bangladesh. furthermore, 470 drill holes were chosen for the formation of soft clay layers. in addition, a map identified with zoning for spt values were from 1 to 5. the total area was divided into five subsections according to the thickness of soft clay soil. the lower range of soft clay thickness is 0 to 1 and the higher range is between 10 to 20. from the thickness map and bore log, spt n zone was realized. the target of this study was to obtain the accuracy of the seismic design of roadway and rail sub structure, based on the aforementioned necessities. on this model embankment, 100 shake table tests were experimented. the analysis implemented a repeated loading and unloading process. 199 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 figure 1. study area and the thickness map for soft soil in bangladesh 2.2 using equipment and materials the shake table facility, which is computer-controlled, was used to simulate the horizontal shake action associated with dynamics. the platform of testing was a square, with a dimension of 2.5×2.5 m² and an approximate payload capacity of 1100 kg, made up of steel plates. the range of acceleration is 0.05 g to 2 g. the frequency range is 0.05 to 50 hz with a maximum amplitude of ±200 mm. the highest velocity is 0.03 m/s. the shake table test machine is shown in figure 2. the laminar box is positioned on the shake table as presented in figure 3. figure 2. shake table test apparatus figure 3. laminar box mounted on shake table a laminar box was constructed on the shake table apparatus to reduce boundary effects wherever possible. the laminar box neither resists nor promotes soil displacement to 200 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 accommodate the movement of soil. the laminar shear box has 24 hollow aluminum layers of frames developed at bangladesh university of engineering and technology (buet). every layer consists of an internal structure with internal dimensions of 915 mm × 1220 mm × 1220 mm. a laminar box was constructed on the shake table apparatus to reduce boundary effects wherever possible. the laminar box neither resists nor promotes soil displacement to accommodate the movement of soil. the laminar shear box has 24 hollow aluminum layers of frames developed at buet. every layer consists of an internal structure with internal dimensions of 915 mm × 1220 mm × 1220 mm. sylhet sand was used as the backfill material that was available locally. the unified soil classification system classified the sand as poorly graded sand (s.p). the maximum and minimum dry densities were 18 kn/m³ and 16 kn/m³m3, respectively. the specific gravity of the sand particles was 2.34. the relative density of sand was 60%. the soft clay soil in dhaka was used in this research work.the liquid limit and water content of this soil sample were found at 40% and 23% respectively. the soil sample used was prepared, by using 50% of water content (1.25 times of liquid limit). cohesion was obtained by 14.8 kn/m² and friction was obtained by 1.0 from the direct shear test. the water content of the soil sample was 14%, and the unconfined compressive strength (qu) was 19 kpa, after the loading had been done. the thickness of the clay layer in the soil sample was 6 m. reconstituted soil sample preparation was displayed in figure 4. figure 5 illustrates the schematic diagram of the test configuration. figure 4. preparation of clay layer figure 5. schematic illustration of shake table test 2.3 testing procedure the dimension of the laminar box was 1220 mm deep, and the size of the plan area was 915 mm × 1220 mm. in lifts of equal height, the model was constructed while 201 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 reinforcing each lift with a layer of woven geotextile. to measure pore water pressure, two pore water sensors (p1 and p2) were placed as shown in figure 5. besides this, 6 acceleration sensors (a1, a2, a3, a4, a5, and a6) were placed on different points of the model, to measure acceleration. a1 and a2 sensors were placed in clay soil layers. three displacement sensors (lvdt1, lvdt2, and lvdt3), were placed on different locations of the model to measure the displacement of the embankment. table 1 shows the test sequence. table 1. test sequence test name acceleration amax(g) frequency hz relative density (%) surcharge (kg) ft1 0.1 1 60 19 ft2 0.1 3 60 19 ft3 0.1 5 60 19 ft4 0.1 10 60 19 ft5 0.1 12 60 19 ft6 0.1 15 60 19 3. results and discussion the total number of 100 shaking table tests were conducted for this research. in this section, 6 shake table tests were described among the 100 shaking tests to evaluate the seismic response of the model retaining wall. the chosen base accelerations for this research were 0.1 g, 0.2 g, 0.3 g, 0.4 g, and 0.5 g respectively. the natural frequency in this shake table test was determined to be 16 hz. therefore, input frequency should be less than the natural frequency of the model. the input frequencies were 1 hz, 3 hz, 5 hz, 10 hz, 12 hz and 15 hz respectively. the surcharge pressures selected for this study were 19 kg, 34 kg, and 49 kg. a total number of 470 layers of soft clay soil were also presented in the gis interface map. the standard penetration test (spt) zonation map was prepared based on spt n value. the spt zone map has four sub sections. the lower range of spt zonation map is 1 to 2 and the higher range is greater than 5. the n value is less than five, showing the existence of soft clay content of that area. this influences the dynamic behaviours of the embankment, (hore et al., 2019). from the map (figure 6), the green colour shows the spt value between 1 to 2. the capital of bangladesh, dhaka, is occupied with soft clay zone containing spt value 1 to 2. the output result was the dynamic behaviour like acceleration amplification, displacement and pore water pressure. this was shown in the contour map and the graph. the research 202 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 showed that dynamic behaviour of the soil depends on the soil standard penetration test result. the zonation map using standard penetration test (n) value was shown in figure 6. figure 6. zonation map using standard penetration test (spt) 3.1 acceleration response the dynamic parameters for acceleration amplification response were ft1, ft2, ft 3, ft4, ft5, ft6 with frequencies of 1 hz , 3 hz , 5 hz , 10 hz , 12 hz and 15 hz respectively. these were conducted at 0.1 g base acceleration and 19 kg surcharge pressure, as shown in figure 7. figure 8 shows the acceleration amplification variation with respect to frequency without clay layers. figure 7 displays the two sensors in the clay soil sample layer for the different frequencies of 1 hz, 3 hz, 5 hz, 10 hz, 12 hz and 15 hz from ft1, ft2, ft3, ft4, ft5, and ft6 model tests respectively. from the figure, it can be seen that acceleration amplitude decreases, with increasing normalized elevation,for the frequency of 1 hz. conversely, acceleration amplitude increases with rise innormalized elevation for other frequencies. the previous figure shows that acceleration amplification and frequency are not directly proportional. in fact, within the range of tests conducted, accelerations were amplified less for 1 hz, 3 hz, and 5 hz and more for 10 hz and 12 hz compared to that of 15 hz at all elevations. moreover, accelerations at normalized elevations of 0.25, 0.5, 075, and 1 203 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 were amplified closer or slightly more than 1 for the frequency of 1 hz. the differences in acceleration amplification for various frequencies were increased with increase in wall height. at a normalized height of 1, for 1 hz, 3 hz, 5 hz, 10 hz, 12 hz, and 15 hz, the values for acceleration amplification were 1.04, 1.24, 1.46, 3.26, 2.47, and 1.80 respectively. these test results of acceleration are partially similar to the (cai et al., 2019; fleming et al., 2016; hore et al., 2020). 3.2 face displacement response figure 9 shows the displacement profiles observed for tests ft1, ft2, ft3, ft4, ft5, and ft6 with frequencies 1 hz, 3 hz, 5 hz, 10 hz, 12 hz and 15 hz respectively. it was observed that at the highest elevation (z/h=0.875), the displacement was maximum. the highest normalized displacement of 2.01% was observed for 12 hz frequency. the corresponding values for the frequencies of magnitudes 1 hz, 3 hz, 5 hz, 10 hz and 15 hz were 0.02 %, 1.99 %, 2.05 %, 2.01 % and 1.93 % respectively. these test results of face displacement response are partially similar to the (krishna & latha, 2007; srinivasan et al., 2016; suzuki et al., 2015). figure 7. acceleration amplification (clay layer) response on frequency figure 8. acceleration amplification response on effect of frequency frequ whole.grf acceleration amplification n o rm a li se d e le v a ti o n , z/ h 0 0.5 1 1.5 2 2.5 3 3.5 0 0.2 0.4 0.6 0.8 1 test no. freq. ft1 1hz ft2 3hz ft3 5hz ft4 10hz ft5 12hz ft6 15hz frequu.grf acceleration amplification n or rm al is ed e le va ti on , z /h 0 0.5 1 1.5 2 2.5 3 3.5 0 0.2 0.4 0.6 0.8 1 test no. freq. ft1 1hz ft2 3hz ft3 5hz ft4 10hz ft5 12hz ft6 15hz 204 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 3.3 variations of pore water pressure figure 10 shows the effect of frequency for fixed base acceleration (0.1 g) and surcharge (19 kg) on acceleration amplification, strain and pore water pressures. the variations were for tests ft1, ft2, ft3, ft4, ft5 and ft6 with the frequency of 1 hz, 3 hz, 5 hz, 10 hz, 12 hz and 15 hz respectively for 0.1 g base accelerations and 19 kg surcharge. the pore water pressures increase with increasing elevation. the pore water pressure was 0.07 kpa at frequency of 10 hz. the maximum pore water pressure for model tests ft1, ft2, ft3, ft4, ft5, and ft6 were 0.02 kpa, 0.03 kpa, 0.03 kpa, 0.11 kpa, 0.09 kpa, and 0.07 kpa respectively. figure 11 shows the contour map (plaxis 3d output results) of acceleration, displacement, and pore water pressure response. figure 9. displacement profile response on frequency figure 10. pore water pressure variations 0.1g_19kg new.grf normalised face displacement, h/h (%) n o r m a li s e d e le v a t io n , z /h 1.7 1.75 1.8 1.85 1.9 1.95 2 2.05 2.1 2.15 0.3 0.36 0.42 0.48 0.54 0.6 0.66 0.72 0.78 0.84 0.9 test no. freq. ft1 1 hz ft2 3 hz ft3 5 hz ft4 10 hz ft5 12 hz ft6 15 hz 0.1g_19 kg.grf pore water pressure, kpa n o r m a li s e d e le v a t io n , z /h 0 0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 0 0.08 0.16 0.24 0.32 0.4 0.48 0.56 0.64 0.72 0.8 test no. freq. ft1 1 hz ft2 3 hz ft3 5 hz ft4 10 hz ft5 12 hz ft6 15 hz 205 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 figure 11. contour map of acceleration response figure 12. contour map of displacement response figure 13. contour map of pore water pressure response this research discusses the results found from 6 different shake table tests from 100 combinations, on the embankment with soft clay models. the acceleration of thistest varied from time to time. the acceleration of the test were as follows: 0.1 g, 0.2 g, 0.3 g, 0.4 g and 0.5 g. during the shake table test, the natural frequency was the first calculated and its test result was 16 hz. therefore, the other frequencies should be less than the natural frequency. the other frequencies were 1 hz, 3 hz, 5 hz, 10 hz, 12 hz and 15 hz. this means that the 206 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 frequencies ranged from 1 hz to 15 hz. during the shake table test, surcharge pressure varied for each test. changing the surcharge pressure and wrap faced embankment was newly created. the surcharge pressures were 19 kg, 34 kg and 49 kg. after conducting the unconfined compressive test, the result was 20 kpa. the wrap faced embankment height was 4 m and therefore, the model was developed. the length (l) of the geotextile reinforcement was 3.75 m. also, 20 cycles of sinusoidal shakes were subjected to the model wall. all of the present model walls were later constructed with sand placed on the same medium to loosen density. the sand embankment properties were as follows: i) average unit weight ii) the relative density which were 18 kn/m3 and 60% respectively. these results obtained were similar to that of other researcher swho only used sand embankment without clay soil (eric et al., 2013; goktepe et al., 2019; srilatha et al., 2013). the cyclic behaviour of the soft clay has been analysed based on the different input frequency. the face deformations were high for low-frequency shaking, low for surcharge pressures and high for base accelerations. the pore water pressures were observed to increase with a rise in base motion frequency. three contour maps were drawn to investigate the acceleration amplification, deformation, and pore water pressure with respect to frequency response. these maps also showed dynamic parameters such as acceleration amplification, deformation, and pore water pressure representations of the soft soil. these results were very important for observing the dynamic behaviour of wrap faced soil retaining walls on the soft clay layer. this can be applied to improve incorporating dynamic loading, considering the design specification of this type of retaining wall (railway and road embankment). 4. conclusion the shake table test on wrap faced embankment on the clay soil foundation, is a new form of test. the lower range of spt zonation map is 1 to 2 and the maximum range is greater than 5. in bangladesh, embankment on soft clay soil plays a very vital role in seismic perspective. therefore, the research opportunity of this area will bring about an upgrade to seismic design specifications. from the test results, it was discovered that face deformation increases with rise in elevation. at high elevations and frequencies, pore water pressure is also high. these test results are beneficial in understanding the relative performance of the wrap faced embankment, which is the indicative direction to the design process of embankment. 207 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 conflict of interest the authors declare that there is no conflict of interest with any financial, personal, or other relationships with other people or organizations related to the material discussed in the article. references al zaman, m. a., & jahan monira, n. (2017). a study of earthquakes in bangladesh and the data analysis of the earthquakes that were generated in bangladesh and its’ very clo se regions for the last forty years (1976-2016). journal of geology & geophysics, 06(04), 6–10. https://doi.org/10.4172/2381-8719.1000300. alam, m. j., & islam, m. s. (2009). geological aspects of soil formation of bangladesh. conference: bangladesh geotechnical conference 2009 (bgc-2009), dhaka, bongladesh. bazlar rashid, m., rashid, b., & islam, b. (2018). structure and lineaments of the northwestern part of bangladesh and evolution of the barind tract. american journal of earth sciences, 5(3), 26–36. beskhyroun, s., wegner, l. d., & sparling, b. f. (2011). integral resonant control scheme for cancelling human-induced vibrations in light-weight pedestrian structures. structural control and health monitoring, may 2011, n/a-n/a. https://doi.org/10.1002/stc. bullock, z., dashti, s., liel, a., porter, k., & karimi, z. (2019). this preprint is a pdf of a manuscript that has been accepted for. earthquake spectra (preprint), 1–33. https://doi.org/10.1016/j.ssresearch.2016.09.015. cai, w., yu, b., & kaewunruen, s. (2019). shaking table tests of suspended structures equipped with viscous dampers. applied sciences (switzerland), 9(13). https://doi.org/10.3390/app9132616. çelebi, e., göktepe, f., & omid, a. j. (2019). seismic soil response of scaled geotechnical test model on small shaking table. arabian journal of geosciences, 12(2). https://doi.org/10.1007/s12517-018-4197-8. edinçliler, a., & toksoy, y. s. (2017). shake table tests to measure the dynamic performance of geotextile-reinforced embankment. periodica polytechnica civil engineering, 61(4), 803–814. https://doi.org/10.3311/ppci.10540. eric, r., moss, s., & crosariol, v. a. (2013). scale model shake table testing of an underground tunnel cross section in soft clay. 29(4), 1413–1440. https://doi.org/10.1193/070611eqs162m. ering, p., & sivakumar babu, g. l. (2020). effect of spatial variability of earthquake ground motions on the reliability of road system. soil dynamics and earthquake engineering, 136(december 2019), 106207. https://doi.org/10.1016/j.soildyn.2020.106207. fleming, b. j., eeri, m., sritharan, s., & miller, g. a. (2016). full-scale seismic testing of piles in improved and unimproved soft clay. 32(1), 239–265. 208 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 https://doi.org/10.1193/012714eqs018m. goktepe, f., celebi, e., & omid, a. j. (2019). numerical and experimental study on scaled soil-structure model for small shaking table tests. soil dynamics and earthquake engineering, 119(august 2018), 308–319. https://doi.org/10.1016/j.soildyn.2019.01.016 hamayoon, k., morikawa, y., oka, r., & zhang, f. (2016). 3d dynamic finite element analyses and 1 g shaking table tests on seismic performance of existing group-pile foundation in partially improved grounds under dry condition. soil dynamics and earthquake engineering, 90, 196–210. https://doi.org/10.1016/j.soildyn.2016.08.032. haque, m. e., sayem, h. m., & kamal, m. (2013). assessment of some engineering geological aspects of the sub-soil of ganakbari, dhaka, bangladesh. jahangirnagar university environmental bulletin, 2(september), 61–70. https://doi.org/10.3329/jueb.v2i0.16331. hassan, a., & pal, s. (2018). effect of soil condition on seismic response of isolated base buildings. international journal of advanced structural engineering, 10(3), 249–261. https://doi.org/10.1007/s40091-018-0195-z. he, t., & jiang, n. (2019). substructure shake table test for equipment-adjacent structure– soil interaction based on the branch mode method. structural design of tall and special buildings, 28(4), 1–19. https://doi.org/10.1002/tal.1573. helwany, s., wu, j. t. h., meinholz, p., alizadeh, v., & ghaderi, r. (2017). seismic behavior of grs bridge abutments with concrete block facing: an experimental study. transportation infrastructure geotechnology, 4(4), 85–105. https://doi.org/10.1007/s40515-017-0040-z. hore, r., arefin, m. r., & ansary, m. a. (2019). development of zonation map based on soft clay for bangladesh. journal of engineering, 10(1), 13–18. hore, r., chakraborty, s., shuvon, a. m., & ansary, m. a. (2020). effect of acceleration on wrap faced reinforced soil retaining wall on soft clay by performing shaking table test. 15, 1–11. krishna, a. m., & latha, g. m. (2007). seismic response of wrap-faced reinforced soilretaining wall models using shaking table tests. geosynthetics international, 14(6), 355– 364. https://doi.org/10.1680/gein.2007.14.6.355. kumar, v., chopra, s., choudhury, p., & kumar, d. (2020). estimation of near surface attenuation parameter kappa (κ) in northwest and northeast himalaya region. soil dynamics and earthquake engineering, 136(september 2019), 106237. https://doi.org/10.1016/j.soildyn.2020.106237. latha, g., & krishna, a. (2006). shaking table studies on reinforced soil retaining walls. indian geotechnical journal, 36(4), 321–333. madhavi latha, g., & manju, g. s. (2016). seismic response of geocell retaining walls through shaking table tests. international journal of geosynthetics and ground engineering, 2(1), 1–15. https://doi.org/10.1007/s40891-016-0048-4. srilatha, n., latha, g. m., & puttappa, c. g. (2013). effect of frequency on seismic response of reinforced soil slopes in shaking table tests. geotextiles and geomembranes, 209 ripon hore et al / geosi vol 5 no 2 (2020) 196-209 36(september 2014), 27–32. https://doi.org/10.1016/j.geotexmem.2012.10.004. srinivasan, v., srivastava, s., & ghosh, p. (2016). optimization and parametrical investigation to assess the reconstitution of different types of indian sand using portable travelling pluviator. geotechnical and geological engineering, 34(1), 59–73. https://doi.org/10.1007/s10706-015-9928-y. suzuki, m., shimura, n., fukumura, t., yoneda, o., & tasaka, y. (2015). seismic performance of reinforced soil wall with untreated and cement-treated soils as backfill using a 1-g shaking table. soils and foundations, 55(3), 626–636. https://doi.org/10.1016/j.sandf.2015.04.013. xiao, m., ledezma, m., & hartman, c. (2014). shake table test to investigate seismic response of a slurry wall. geotechnical special publication, 234 gsp, 1234–1243. https://doi.org/10.1061/9780784413272.120. yazdandoust, m. (2017). investigation on the seismic performance of steel-strip reinforcedsoil retaining walls using shaking table test. soil dynamics and earthquake engineering, 97(november 2016), 216–232. https://doi.org/10.1016/j.soildyn.2017.03.011. ye, j., & jiang, l. (2018). simplified analytical model and shaking table test validation for seismic analysis of mid-rise cold-formed steel composite shear wall building. sustainability (switzerland), 10(9). https://doi.org/10.3390/su10093188. zhou, y., wang, h., wen, r., ren, y., & ji, k. (2020). insights on nonlinear soil behavior and its variation with time at strong-motion stations during the mw7.8 kaikōura, new zealand earthquake. soil dynamics and earthquake engineering, 136(29), 106215. https://doi.org/10.1016/j.soildyn.2020.106215. 288 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 land value potential zonation : implication towards urban planning revi mainaki1, anita eka putri1*, dwiyono hari utomo2 1department of geography education, faculty of teacher training and education, universitas siliwangi, jl. siliwangi no. 24, kota tasikmalaya 46115, indonesia 2department of geography education, faculty of social science, universitas negeri malang, jl. semarang 5, malang, 65145, indonesia *corresponding author : anita.eka@unsil.ac.id received 15 april 2020/ revised 7 august 2020 / accepted 15 august 2020/ available online 22 august 2020 abstract potential land prices are strongly influenced by various factors, cimahi city has three subdistricts displaying unique characteristics, since it is located between two districts and one large city, which affects the potential price of land. the potential price of land is crucial to identify, especially in determining the policies of related agencies, the purpose of this research was to zoning potential land prices in cimahi city. this study engaged a quantitative approach utilizing data collection technique in the form of observation, literature, documentation, and interviews, then analysis was conducted using a gis which composed of assessment, weighting, coating, and buffering. the study population was all sub-districts in cimahi. the research samples were taken from several sub-districts which were influenced by districts and cities characteristics indicators such as accessibility, land usage, and land ownership status. the results showed land prices potential was classified as low, medium, and high which was derived by several indicators. it could be concluded that the potential land price is strongly influenced by the surrounding area features, especially urban and district infrastructure. . keywords: land value potential; zonation; urban planing 1. introduction land is needed by humans for all forms of activities, especially for living and farming to meet their daily needs. land ownership in indonesia have been regulated in government policy related to the start of laws, legislation, government regulations, presidential regulations, presidential decrees and other technical rules (virgo et al., 2011). cimahi city presents unique characteristics to investigate how the land price is predicted, since cimahi city is one of the bearers and directly adjacent to the city of bandung with proximity of 12 km as the capital of west java province. cimahi city is bordered by bandung regency and west bandung regency in the other side. ownership along with rapid development are determining the value of a particular currency, the more strategic or the closer residence to meet human needs, certainly the higher it gets to be paid. land prices are mainly drived by several factors (anggraini et al., 2015; geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.5 no. 2 (2020), 288-300, august, 2020 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v5i2.17442 accredited by the ministry of research , technology , and higher education of the republic of indonesia, no. 30/e/kpt/2019. mailto:anita.eka@unsil.ac.id https://jurnal.unej.ac.id/index.php/geosi https://drive.google.com/file/d/1rsnvas6cuhowhl5bj87cl2l6k5dqz7s6/view 289 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 glaeser et al., 2005) namely (1) physical factors such as soil type, slope, elevation, land area and land utilization; (2) economic factorssuch as purchasing power, interests, amount of benefits offered, interest rates and land beneficial; (3) social factors namely population, security level and community lifestyle; (4) government factors namely taxes and related policies; (5) location and accessibility such asaccessibility towards education, health facilities, main roads, markets, transportation, mobility pace and etc; (6) availability of facilities for electricity networks, clean water, telephone networks, religious amenities, education, health and the tranquil spots. a study result indicated that land prices are drawn by several factors (fakhirah, 2010; hidayanti, 2013) as follows (1) distance aspect upon road accessibility; (2) elevation aspect which correlates to temperature and water; (3) distance towards donwtown or main activities spots; (4) zoning; (5) mobility pace; (6) land topography; (7) land area, and (8) land usage types. the land classification scheme to determine the price of land in a certain parameter including land purposed for (1) settlements; (2) trade and services; (3) industry; (4) transportation, communication and facilities; (5) a relatively complex trading industry; 6) urban areas within downtown and (7) urban areas within urban areas (rusdi, 2013; prabowo et al., 2016; sasono & susetyo, 2018; prihandoko, 2018). land prices would possibly increase in accordance to the occurrence of economic activity that is correlated towards infrastructure such as roads, government financial and other activities in centrals (rynjani & ragil, 2015). the significant of zoning as a form of reconstruction is to estimate the potential price and to determine certain area performing a spatial approach. a research showed how zoning could be utilized as basic decision-making theories of the region's level of urgency (hermawan & mainaki, 2019). zoning could also be conducted using basic relics and other object indicators with reconstruction principle to determine the potential or forecasts of a particular are (mainaki & hermawan, 2019). once the zoning results are identified, the results are not only beneficial for decision making, but also could be used as a learning resource, as the map becomes a basis in identifying an area as a whole, to compherend the thoughts that would foster spatial intelligence. the results of this study could be used as one of the additional indicators in determining the purchase of shelter for the community in cimahi as one of the growing centers (hilman & mainaki, 2020; rahayu & haryatiningsih, 2013). the other side of cimahi is bordered by west bandung regency, which is a new administrative area that has expanded, only a few years since the district was formed. this 290 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 feature establishes the city of cimahi occupies urban areas both within the middle of urban and in the suburbs. it has a relatively diverse area with social characteristics and physical potential, therefore land utilization determines the price of the land is relatively varied and has not been zoned in a wider scope area. hence, land price predictions in cimahi city which are essential information for urban planning and management (koka et al., 2013). a few researchers focused on development, since the potential price of land turns high or low. there have been limited studies concerned on a smaller scale and rarely utilized a regional approach, which was a draw back if a study attempted to analyze particular zone or region. therefore, this current study aimed to describe the potential price of land in the form of reports and zoning maps engaging scoring and weighting on gis, the purpose of this research was to zoning potential land prices in cimahi city. 2. methods in order to analyze the potential of land prices in the city of cimahi, this study employed a quantitative research approach using procedural steps based on the geographic information system (gis). the basic analysis was performed for potential prices which later used as a feasibility study, as the indicators used in the study were based on relevant agency documentation research (putri et al., 2020) namely (1) preliminary study to determine the topic, issues formulation, objectives and results to be obtained; (2) determine the research methodology and instruments development in accordance to collected data; (3 ) data collection was gained by documentation method were taken from (a) national land agency (bpn) of cimahi city to obtain data regarding land ownership status, land use and land area, (b) cimahi city transportation agency to obtain data relating to the existing road network in the city; supported by literature review and field observations based on documentation data to identify and strengthen research results referring to various scientific references, hence the basis of zoning grows vivid as it is founded upon theory (mainaki & putri, 2020) data processing obtained relation upon the geographic information system (gis) includes (a) scoring; (b) weighting; (c) determining the maximum score and weight to identify the range of scores (kusumo & nursari, 2016; gunawan et al., 2014; siagiaan et al., 2015). the analysis unit in this study was urban villages, thus the population was all urban villages in cimahi city and the sample was several urban villages that displayed similar characteristics of one to two district. the potential land prices were determined using scoring 291 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 and weighting tables (table 1, table 2, table 3, table 4, table 5), the higher score means the higher potential land price and vice versa as the basis of map overlay is presented below: table 1. the process of scoring data based on the area of the land or building no. land use large large (>10000 m²) is being (1000-10000 m²) narrow (<1000m 2 ) 1. rice fields 3 2 1 2. plantation 3 2 1 3. offices/schools/hospitals/health centers 3 2 1 4. factory/industry 3 2 1 5. shops/supermarkets/markets 3 2 1 6. settlements 3 2 1 7. empty land 1 1 1 (source: 2019, analysis results) table 2. the process of scoring data based on land utilization no. land use score 1. rice fields, plantations and vacant land 1 2. offices/schools/hospitals/health centers and settlements 2 3. shops/supermarkets/markets and factories/industries 3 (source: 2019, analysis results) table 3. the process of scoring (scoring) data based on accessibility no. road type score close (<50 m) is being (50-100 m) far (> 100 m) 1. distance to frontage road 3 2 1 2. distance to boulevard 3 2 1 3. distance to provincial road 3 2 1 (source: 2019, analysis results) table 4. process of scoring data based on mastery of land ownership no. criteria score 1. right to use (tni hp) 1 2. right of use (mobile) 2 3. right to build (hgb) 2 4. customary property/rights (hm/hma) 3 (source: 2019,analysis results) table 5. the process of weighting (weight) based on map effect no map weight information 1 land use 3 the type of use and area of land greatly influences the value of land and impacts on prices and determination of tax rates. 2 road network 3 affordability of land affects the value of the land and the calculation of the tax rate. 3 land tenure 2 to identify land tenure (property rights or state property) in determining taxes. (source: 2019, analysis results) 292 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 the overlay calculation for zonation is performed as follows: minimum total score = (minimum score of land area x minimum weight of land use map) + (minimum land use score + minimum weight of land use map) + (minimal score of adjacency to road x minimum weight of road network map) + (minimum land mastery score x minimum weight of status map land) minimum total score = (1x3) + (1x3) + (1x3) + (1x2) = 3 + 3 + 3 + 2 = 11 maximum total score = (maximum land area score x maximum weight of land use map) + (maximum land use score + maximum weight of land use map) + (maximum score of adjacency to road x maximum weight of road network map) + (maximum score of land mastery x maximum weight of status map land) maximum total score = (3x3) + (3x3) + (3x3) + (3x2) = 9 + 9 + 9 + 6 = 33 zonation range : = ∑𝑀𝑎𝑥𝑖𝑚𝑎𝑙 𝑆𝑐𝑜𝑟𝑒 − ∑ 𝑀𝑖𝑛𝑖𝑚𝑢𝑚 𝑆𝑐𝑜𝑟𝑒 𝐶𝑎𝑡𝑒𝑔𝑜𝑟𝑦 (1) = 33−11 3 = 7.33 rounded 7 hence, the zonation range is: 11 to 18 = low land price potential; 19 to 26 = potential medium land prices; 27 to 33 = potential high land prices. 3. results and discussion land utilization in the city of cimahi (figure 1) presented practically similar characteristics of land use in each village that was dominated with settlements, rice fields and plantations. however, in the south side were relatively large industrial and shopping areas which indicated on how the southern part gained the potential for higher land prices compared to central and northern parts. the northern part is bordered by west bandung regency while the southern part is bordered by bandung city. hence, considering both the land utilization and location, the south side seemed promising higher price potential and vice versa occurred upon the north that obtained lower price potential. the road network in cimahi (figure 1) has been relatively wide-spread, as each place was connected with local avenue and highway. the central part was cut across over provincial road in a single direction, that way the south was remarkably more strategic than north. the south was through fare towards provincial road one-way and two-way, making it more strategic than the central and northern parts of cimahi. this phenomenon explains why 293 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 the south has higher potential land prices according to the road network adjacency and accessibility upon low to high, namely the north, center and south. according to land ownership status (figure 1) in the north and south, most of them earned the land ownership/customary ownership rights status. in the southern part obtained land utilization rights, where the central part obtained of land use rights status (military),as cimahi as a military basis in the sense of the community implying that the middle part was inventory or land of country that was only function for military. the northern and southern displayed interests for the benefit upon other human activities, means the northern and southern parts will likely draw a higher potential land price compared to the central part of cimahi city. potential land prices (figure 2) in the north or cimahi utara sub-district used pasirkaliki sub-district as sample possed regular or irregular residential or irregular land use conditions and rice fields, the land use status of all customary property rights. this urban village presented practically similar characteristics as other urban villages in north cimahi sub-district. it endowed potential high land prices of customary property rights status, in term of residential or regular land utilization and relatively adjacent towards the road based on distance factor, while areas with potential taxes were customary land/customary land status with regular and relatively moderate or far-off residential land use by road. based on the analysis conducted to the middle section with setiamanah village as a sample that displayed identical characteristics towards other urban villages within central cimahi sub-district, it could be concluded that setiamanah and other urban villages within the central part of cimahi city enjoyed high and medium land prices, most of which were potentially due to land use of military status of usufructuary rights (government land). the other minority gained high potential price, since those earned customary land/customary rights status. the sample of the southern part or the southern district of cimahi in the current study was the cibeber village, with land status, land use and various road conditions as appeared in other villages in the southern part of cimahi city. the administrative map, land status, land use and road network of cimahi city is presented be on figure 1. 294 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 figure 1. administrative map, land status, land use and road network of cimahi city revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 figure 2. map of land use, land status, potential land prices in pasirkaliki, setiamanah and cibeber villages in cimahi city 295 296 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 cibeber sub district demonstrated a relatively diverse land prices potential (figure 2). the potential high price was expected for land with customary property rights status, strategically adjacent towards one or two-way provincial roads, supported with residential and industrial/market land use. land use is the land with right to establish building or military base, while the potential for low land prices was presented with the condition upon usufructuary land status of employment for paddy fields or plantations and relatively remote towards the main road in accordance with the direction of proximity to the road. considering the influence of land status, land use, adjacency, and land accessibility, cimahi city endowed potential of land prices across north to south. in the north side was potentially possess low prices, while in the south showed higher value, and in middle had more diverse potential land prices. the map of landuse, land status, and potential land prices in pasirkaliki, setiamanah and cibeber village in cimahi city is presented on figure 2. the results showed that potential for the low, medium and high land potential prices within cimahi city were strongly influenced whether by road network or land utilization in the term of urban village surroundings. land status affects the land prices potential, however not by two other aspects. this result is in accordance with research conducted by fakhirah (2010) and hidayanti (2013) which land prices were drawn by (1) distance aspect towards road accessibility; (2) elevation aspect that relates to temperature and water; (3) accessibility towards center of activities or activities; (4) zoning; (5) mobility pace; (6) land topography; (7) land area and (8) land use type. the result also comformable towards rusdi (2013); prasetya & sunaryo, (2013); sasono & susetyo (2018); prihandoko (2018) which described the land classification scheme to determine the land price in a certain parameter, namely land use for (1) settlements; (2) trade and services; (3) industry; (4) transportation, communication and facilities; (5) a relatively complex trading industry; (6) urban areas within downtown and (7) urban within the urban areas. the results showed the correlation between several studies in a single result with a specific area and broader approach (regional approach). this result is in accordance with research conducted by elmanisa et al. (2017) examined land prices across jakarta, bogor, depok, tanggerang and bekasi or abbreviated as jabodetabek illustrated how land prices in these areas were relatively high due to resources and infrastructure, especially within metropolitan area of the capital. in indonesia, there are several areas with zoning potential for research land, as urban areas in cimahi city gained relatively high land prices which tended to be close to bandung city as the capital area of 297 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 west java province and otherwise, it tended to below or moderate once the areas were rather remote towards bandung city or closer to urban areas. this research is also in line with the research of hidayanti (2013) that engaged highresolution satellite image analysis to zoning land prices in jetis district, yogyakarta utilizing relatively similar indicators. the results are quite equal, which strengthened the research results of cimahi city with relatively similar results are obtained. therefore, the results of this study could be utilized as a general reference for upcoming research regarding zoning potential land prices in various regions. in contrast to research of masykuroh & rudiarto (2016) which attempted on changes in land prices, as there was an effect highway development, in this case is accessibility aspect. there was an increase of land prices towards the area nearby the road or highway gate of ungaran. the result supports the research undertaken,since accessibility aspect significantly affects the potential price of land. research of rahati et al. (2015) is also in line with this study, as the occurrence of land damage due to natural disasters had greatly reduced land prices, certainly, it is related to the physical aspect that is being affected by the damage and causes material loss. based on the research discussions according to several other studies, this study presents a parallel and general position, what distinguishes this study from other research is that the location demonstrates intriguing characteristics. in this study of cimahi city, which represented an area that obtain characteristics established upon cities and districts. 4. conclusion all over areas in cimahi city obtained a road construction network that was relatively evenly crossed by city roads, sub-district roads and provincial roads. the southern part of cimahi city possed an area that was supported with bandung city namely factories, shophouses and industrial areas, where the land status was dominated by ownership and user rights,therefore most land prices were relatively high. the central area was dominated by settlements and military areas utilization with land status dominated by military rights. hence, most areas displayed the potential for moderate land prices as those were the central controlled by municipalities or regencies, while the northern part was dominated by open land use such as plantations and land. rice fields ownership and use rights, thus most of the potential land prices were low and moderate due to excess, by the district area. 298 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 conflict of interests the authors declare that there is no conflict of interest with any financial, personal, or other relationships with other people or organizations related to the material discussed in the article. acknowledgment we would like to express our gratitude to universitas pendidikan indonesia, specifically bachelor degree programs that have provided us opportunities to perform research and formulated the journal articles. references anggraini, s., prawira, y. c. s., untung, m., & indonesia, p. (2015). estimated the reasonable land price in industrial estate : case study of an industrial estate in jakarta. fig working week. sofia, bulgaria. elmanisa, a. m., kartiva, a. a., fernando, a., arianto, r., winarso, h., & zulkaidi, d. (2017). land price mapping of jabodetabek, indonesia. geoplanning: journal of geomatics and planning, 4(1), 53-62. fakhirah. (2010). identifikasi faktor yang mempengaruhi nilai jual lahan dan bangunan pada perumahan tipe sederhana. jurnal smartek, 8(4), 251–269. glaeser, e. l., gyourko, j., & saks, r. (2005). why is manhattan so expensive? regulation and the rise in housing prices. journal of law and economics, 48(2), 331–369. https://doi.org/10.1086/429979. gunawan, d; hartanti, s & maulana, y. (2014). rancang bangun aplikasi analisis kredit menggunakan metode skoring pada bintang jaya variasi audio. jurnal sistem informasi, 3(2), 97–103. hermawan, i & mainaki, r. (2019). pemetaan jalur dan tinggalan perkeretaapian masa kolonial belanda di wilayah cirebon timur. jurnal sosioteknologi, 18(3), 560–574. hidayanti, i. (2013). analisis harga lahan berdasarkan citra penginderaan jauh resolusi tinggi. jurnal pendidikan geografi gea, 8(4), 57–71. hilman, i & mainaki, r. (2020). advantage of map as geography learning media to enhance students spatial intelligence. internasional journal of geomate, 18(68), 225– 232. koka, n; monkkonen, p & quigleyc, g. (2013). land use regulations and the value of land and housing; an intra-metropolitan analysis. journals of urban economic, 81(2), 136–148. https://doi.org/https://doi.org/10.1016/j.jue.2014.03.004. kusumo, p & nursari, e. (2016). zonasi tingkat kerawanan banjir dengan sistem informasi 299 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 geografis pada das cidurian kabupaten serang, banten. jurnal string, 1(1), 29–38. rahati, m., kahar, s., & subiyanto, s. (2015). analisis perubahan zona nilai tanah kaitannya dengan banjir di kecamatan pedurungan kota semarang. jurnal geodesi undip, 4(1), 117-128. mainaki, r & hermawan, i. (2019). perkeretaapian masa kolonial belanda di wilayah indramayu: pemetaan jalur dan bukti tinggalan arkeologis. jurnal walennae, 17(2), 125–142. mainaki, r & putri, a. (2020). paleogeografi: perkembangan keanekaragaman hayati dalam ruang dan waktu (biodiversity in time and space). jurnal geografi, edukasi dan lingkungan, 4(1), 17–24. masykuroh, d. k., & rudiarto, i. (2016). kajian perubahan penggunaan lahan dan harga lahan di wilayah sekitar pintu tol ungaran. tataloka, 18(1), 53-66. prabowo, s., hamid, d., & prasetya, a. (2016). analisis partisipasi masyarakat dalam pengembangan desa wisata (studi pada desa pujonkidul kecamatan pujon kabupaten malang). jurnal administrasi bisnis s1 universitas brawijaya, 33(2), 18–24. prasetya, n & sunaryo, b. (2013). faktor-faktor yang mempengaruhi harga lahan di kawasan banjarsari kelurahan tembalang, semarang. jurnal teknik pwk, 2(2), 223– 232. prihandoko, t. (2018). analisis nilai lahan kecamatan mergangsan kota yogyakarta menggunakan aplikasi sig dan penginderaan jauh. jurnal pendidikan geografi ums, 13(1), 1–14. putri, a. e., utomo, d. h., & mainaki, r. (2020). analisis kesesuaian lahan rawa untuk pengembangan kawasan permukiman di kecamatan gondang kabupaten tulungagung. jurnal pendidikan geografi: kajian, teori, dan praktek dalam bidang pendidikan dan ilmu geografi, 25(2), 170-178. rahayu, m & haryatiningsih, r. (2013). persepsi masyarakat dalam menentukan tempat tinggal : kota cimahi sebagai pusat pertumbuhan. jurnal kajian ekonomi dan pembangunan, 9(2), 58–71. https://doi.org/doi: https://doi.org/10.29313/.v9i2.2305. rusdi, m. (2013). faktor-faktor yang mempengaruhi harga dan penggunaan lahan di sekitar jalan lingkar salatiga. jurnal pembangunan wilayah dan kota, 9(3), 317–329. https://doi.org/https://doi.org/10.14710/pwk.v9i3.6543. rynjani, g & ragil, h. (2015). kajian harga tanah dan penggunaan lahan di kawasan perdagangan dan jasa kelurahan lamper kidul, kota semarang. jurnal teknik pwk, 4(3), 417–427. sasono, m & susetyo, c. (2018). analisis potensi perubahan pemanfaatan lahan berdasarkan model spasial harga lahan di kecamatan tembelang kabupaten jombang. jurnal teknik its, 7(1), 60–65. 300 revi mainaki et al / geosi vol 5 no 2 (2020) 288-300 siagiaan, t; sudarsono, b & wijaya, a. (2015). evaluasi kriteria kesesuaian lahan permukiman dengan analitycal hierarchy process (studi kasus :kecamatan boja dan kecamatan limbangan di kabupaten kendal). jurnal geodesi undip, 5(1), 107–115. virgo, e; mahendra, a; rusmawar, w & jp, t. (2011). land value indexation in indonesia: a pilot study in pontianak city. facing the challenges-building the capacity congress. proceeding fig congress. fix-dikonversi_3.pdf (p.1-6) 1 hal fix (1) edited.pdf (p.7) 2 hal fix edited.pdf (p.8) fix-dikonversi_2.pdf (p.9-13) r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 335 geospatial approach for the analysis of forest cover change detection using machine learning r. sanjeeva reddy1,*, g. anjan babu1, a. rama mohan reddy2 1 department of computer science, sri venkateswara university, tirupati, andhra pradesh, 517502, india 2 department of computer science & engineering, sri venkateswara university, tirupati, andhra pradesh, 517502, india *corresponding author : haisanjeevareddy@gmail.com received 13 october 2020/ revised 11 december 2020 / accepted 24 december 2020/ published 30 december 2020 abstract spatial data classification is famous over recent years in order to extract knowledge and insights into the data. it occurs because vast experimentation was used with various classifiers, and significant improvement was examined in accuracy and performance. this study aimed to analyze forest cover change detection using machine learning. supervised and unsupervised learning methods were used to analyze spatial data. a vector machine was used to support the supervised learning, and a neural network method was used to support unsupervised learning. the normalized difference vegetation index (ndvi) was used to identify the bands and extract pixel information relevant to the vegetation. the supervised method shows better results because of its robust performance and better analysis of spatial data classification using vegetation index. the proposed system experimentation was implemented by analyzing the results obtained from support vector machine (svm) and nn (neural network) methods. it is demonstrated in the results that the use of ndvi mainly enhances the performance and increases the classifier's accuracy to a greater extent. keywords: spatial data; normalized difference vegetation index; ndvi;vegetation index, support vector machine; neural network; forest cover change 1. introduction the knowledge of various social, economic, and cultural aspects is considered a correct viewpoint in land management and its planning. in various scenarios, landscape changes did not show up ultimately. therefore, to give exceptional improvement to the lands, different geographic tools are used such as geographic information system and photogrammetry. the point judges of the criteria of forest management that the population is decreasing in some regions. therefore, most of the land is covered with trees and bushes, seriously affecting the landscape (an et al., 2007). to describe the data and gather useful geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.5 no. 3 (2020), 335-351, december, 2020 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v5i3.20157 accredited by the ministry of research , technology , and higher education of the republic of indonesia, no. 30/e/kpt/2019. geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.5 no. 3 (2020), 301-317, december, 2020 https://jurnal.unej.ac.id/index.php/geosi doi : accredited by the ministry of research , technology , and higher education of the republic of indonesia, no. 30/e/kpt/2019. https://jurnal.unej.ac.id/index.php/geosi https://drive.google.com/file/d/1rsnvas6cuhowhl5bj87cl2l6k5dqz7s6/view https://jurnal.unej.ac.id/index.php/geosi https://drive.google.com/file/d/1rsnvas6cuhowhl5bj87cl2l6k5dqz7s6/view r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 336 information regarding the earth's surface, geographic information has been considered as a fruitful source. there are various applications involved, such as digital image analysis, analysis, and detection of a change in environmental conditions, science, education. however, these areas are a source of the right domain to conduct adequate research. geospatial approach for the analysis of forest cover change detection using machine learning. the spatial data is gathered from satellites that include images and define information regarding the image's pixels. the data collected seems unstructured and complex, and then this data is evaluated to get that hidden information. this process is mainly called spatial data analysis, and to get and relocate the landscape in spatial data is known as geospatial data analysis. these landscape types are identified and classified by utilizing techniques that involve deep learning and machine learning. svm is one of the adequate mechanisms used in ml, and in it, the kernel function is activated to conduct the descriptive analysis on the dataset of images. the features are extracted and based on those features, and the machine classifies the landscape types. the main landscape types involved are bare soil, urban land, waterbody, natural vegetation, and forest area. the spatial data gathered was very difficult to cater to, and it involves various critical issues regarding orientation, structure, and other atmospheric conditions (aubrecht et al., 2009; addink et al., 2007). the ongoing writing study on distant detecting information, characterization utilizing ai techniques incorporates the rich data regarding the spatial information focal points, natural biology, accuracy farming, science and building, and military use. recently, the distant detecting information characterization has been finished utilizing better ai and profound learning approaches (pozdnoukhov & kanevski, 2006). the diverse order strategies were utilized on far off detecting information. in high dimensional space, the impediment of dimensionality may yield excellent outcomes. the high dimensional information dealing with is a fundamental task in the enhancement issues (gangappa et al., 2016). subsequently, the enhancement technique, such as svm may regularly be unaware of high dimensional space (acharya & yang, 2015). there is numerous grouping calculation which is regularly used to anticipate the class objects in the spatial data. the regulated learning techniques for spatial information are neural networks, choice tree technique, irregular timberlands order techniques, k-implies bunching and arrangement strategy, harsh set based information decrease and characterization technique, and fuzzy rough set based order method (singh et al., 2016; chi et al., 2008; pal, 2005). the fluffy rationale and neural networks are utilized in spatial information characterization. the r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 337 adequate measure of research work commitment has been on fluffy and unpleasant based components (shanthini et al., 2017; foddy & mathur, 2004; al-obeidat et al., 2015; ham et al., 2005). numerous different calculations have been utilized in advancement issues. in the graphic models, the preparation information with class marks was given at preparing a calculation. in request to get prepared, these techniques (rawat & kumar, 2015) utilize some spatial information highlights, for example, spatial goal, entropy, mean eleven and mean slop, and other pertinent highlights in the input information. we firmly contend that expectation exactness depends on the noteworthy highlights utilized in that model. various machine learning techniques were used to identify different landscape images, and they are considered supervised learning. it involves training data, and the classifier learns from the training data, and other decisions are mainly based on the lear ning of the classifier. the classifier was done its training based on the features extracted. sometimes more the faster features are the classifier's efficiency; however, the feature's amount may affect the classifier performance in some cases. as we can say, the classifier takes more time to validate those features of the dataset gathered, and the performance of the system falls. this may also lead to affect the accuracy of the classifier badly. therefore, for this purpose and feature extraction, some feature reduction techniques are also being involved in the classification process. dimensionality reduction is used to make the dimension space more accurate for the spatial data. ultimately, we emphasize the evaluation of machine learning methods suitable for spatial data pixels. this study aimed to analyze forest cover change detection using machine learning. this research also aimed to find out suitable machine learning techniques that efficiently extract the features, reduce the features according to the demand, and distinguish the features based on the dataset's pixel information. 2. methods 2.1 supervised and unsupervised learning in supervised learning, labeling is involved, in which there is a specific outcome against each entity. furthermore, based on those outcomes, the algorithm accuracy is computed. while on the other hand, in an unsupervised approach, unlabeled data is provided in which the patterns are formed based on the features extracted. this study aimed to analyze forest cover change detection using machine learning. r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 338 2.1.1 supervised learning supervised learning contains a properly labeled dataset and then train the algorithm based on this labeled data (mahmon & ya'acob, 2014). the term fully labeled means that the training dataset contains answers to each question or query. a complete illustration of labeled data and its supervision through supervised learning is shown in figure 1. for example, as related to this study, the forest images are labeled according to the years and their specific characteristics, and then after classification, it is to be identified which ones belong to specific families. when the model is fully trained, it is tested on a new set of images, and the models have to predict values against each set of images fed into it. figure 1. algorithm learns labelled data with supervised learning the support vector machine (svm) is a supervised learning method, and they automatically analyze data, make classes, and put each object into a class by using some rules (gangappa et al., 2017). in svm, labeled data is manipulated for all the classes to be the data classified. the following equations (equation 1 and equation 2) for the svm are as follows: 1 1i iu w l for v      (1) 1 1i iu w l for v      (2) r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 339 supervised learning is mainly used in two scenarios, one in the classification problems and the other in the regression problems. in classification problems, the prediction of values is made by the classifier in which data is recognized based on class. while in regression problems gather continuous data, and in it, the effect of one variable on the variable is identified, such as for a particular value x, what would be the expected value of the variable y. 2.1.2 unsupervised learning in contrast to supervised learning, unsupervised learning contains a deep learning model with a set of instructions on what to do next. in the training dataset, no labeling is involved, and the dataset is without any desired outcome. the network automatically gets the useful features and then analyzes the structure based on the features extracted. the illustration of unsupervised learning is given in figure 2. figure 2. illustration of unsupervised learning the unsupervised scenario data is arranged in two ways, either by arranging clustering or association, which are deeply explained as follows: (1) clustering: in clustering, the data is divided into clusters or groups having the same properties. based on those features or properties, the images are identified and put in their respective groups, which is mathematically explained in equation 3. for example, it is somehow not possible in different scenarios to distinguish the plants or trees based on their appearance. such as all plants have left and stems branches. furthermore, in clustering, these features (equation 3) are locked in different groups. 2 1 1 ( ) (|| ||) mzz m n m n j b a b     (3) || ||m na b  distance among ix and jv , mz total data points, z total cluster centers. r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 340 (2) association: in this way, the algorithm tries to learn without the data being labeled. the algorithm takes some different decisions, such as the forest images are not labeled according to their specific characteristics and fed into the classifier. this is a case of an association, in which highlights of an information test associate with different highlights. by taking a gander at a couple of crucial characteristics of an information point, a solo learning model can foresee different properties with which they are ordinarily related (maity, 2016). the dataset used in this study was landsat 8, and it consists of earth images based on two mechanisms, which are operational land imager oli and thermal infrared sensor trs (lu & yang, 2009). data were collected near the infrared and panchromatic band. the dataset collected consists of various landscape types such as vegetation area, water area, and bare land. more details of the dataset used are illustrated in table 1. table 1. dataset used with respect to its attributes different attributes of dataset used values related to those attributes image format geo tiff orientation north up (map) pixel size 30 meters the data is gathered through remote satellites by using special remote sensors. they mainly reflect the energy on the earth's surface. different types of remote sensors are used, such as spot, landsat, and sar. these sensors contain spectral resolution, and the wavelength represented on the spectral resolution plane is known as the band. different remote sensing systems have different bands. however, in this research work, we have used landsat dataset sensing images, which consists of 11 electromagnetic spectral bands and mainly have 30 meters of spectral resolution. this section presents a framework for spatial data. it contains essential information regarding handling data, evaluating and validating the data to compute useful results. the complete procedure is illustrated in figure 3, as given below. r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 341 figure 3. the pipeline of the main classification framework r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 342 first, pre-processing is done to make the instructed data in a structured format. in it, the dataset is considered as input. at first, the images are digitized to represent the intensity of each pixel in the spectral band. before starting the primary procedure, the raw data need some extra techniques to be used in order to correct errors and avoid noise. the main techniques involved are radiometric correction, geometric correction, and noise removal. after implying these techniques, the data is aligned to the real-world coordinates. in the next step, the vegetation indices are computed to be implied in different scenarios like climate change, detection, monitoring, and modeling of vegetation studies. this procedure helps in combining the information of different bands, and it is also very fruitful to find ndvi by using the following equation (equation 4). nir red ndvi nir red        (4) the final step is to get the region of interest of the image shape file is created. the training and testing samples are also extracted for classification. then the training model is built by using machine learning techniques. 3. results and discussion when the land data was pre-processed, class labels are generated geometrically, such as water, bare land, and forest. then region of interest is extracted from the image and is saved in the shapefile. moreover, from the shapefile, training samples are collected and used in the building of the model. all the features were combined in a shapefile, and then they are used for classification purposes. the features are were represented by the pixel values having specific intensity. the figures (figure 4 – figure 7) show the innumerable spectral assets, which are extricated from the spatial dataset starting from the year 2005 to 2020. r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 343 figure 4. area of land-use in 2005 figure 5. area of land-use in 2010 figure 6. area of land-used in 2015 r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 344 figure 7. area of land-used in 2020 the spatial trials are divided into two training and testing sections. the training and testing sections are collected based on 70% and 30% phenomena. the svm and nn models were used in the classification process. the detailed performance of land cover categories of 2005 and 2020 is deeply illustrated in table 2. the model accuracy and other statistics related to ndvi are recorded and shown in figure 8. the mean of land-use (lu)/land cover (lc) and fcc are shown in figures (figure 9 figure 13). figure 8. accuracy and other statistics related to ndvi -0.015 -0.01 -0.005 0 0.005 0.01 0.015 1995 2000 2005 2010 2015 2020 2025 ndvi r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 345 table 2. land cover categories related to the area square kilometres and area in percentage from 2000 to 2020 the use of ndvi mainly enhances the performance and increases the classifier's accuracy. this result is in line with the results reported by somching et al. (2020), ramos et al. (2020), ahmad et al. (2020), nachapa et al. (2020), gumma et al. (2020), zurqani et al. 2020), liu et al. (2020) that machine learning can improve accuracy in analyzing land-use and land cover change. land-use change from forest to bare land or rocky land is the most dominant. this result is supported by the findings of mirici et al. (2020), yu et al. (2020), devi et al. (2019), that land-use change from forest to bare land and rocky land is becoming dominant due to urbanization and increasing human needs. land-use change from forest to bare land or rocky land occurred mostly in the northern region (figure 9-figure 13). these changes are caused by human activities that require wood as building material. urbanization also causes the human need for forest wood to increase, especially to build new settlements. furthermore, the accessibility of forest locations is also a factor that causes the conversion of forest land to bare land or rocky land. figure 9. mean of land-use/land cover and fcc of 2000 years 2000 2005 2010 2015 2020 landuse/land cover categories area in (sq km) area in (%) area in (sq km) area in (%) area in (sq km) area in (%) area in (sq km) area in (%) area in (sq km) area in (%) less forested 563.0274 21.4 578.387 22.0 559.8693 21.3 632.1591 24.0 504.6606 19.2 rocky 739.6398 28.1 647.5645 24.6 785.8458 29.9 700.2333 26.6 731.9502 27.8 bare land 644.6781 24.5 564.9028 21.5 634.9212 24.1 573.1749 21.8 530.2575 20.2 water body 60.1938 2.3 52.10608 2.0 51.5898 2.0 43.7445 1.7 37.2285 1.4 forested 621.9 23.7 786.5129 29.9 597.213 22.7 680.1273 25.9 825.3423 31.4 total 2629.4391 100 2629.47328 100 2629.4391 100 2629.4391 100 2629.4391 100 r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 346 figure 10. mean of land-use/land cover and fcc of 2005 figure 11. mean of land-use/land cover and fcc of 2010 figure 12. mean of land-use/land cover and fcc of 2015 r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 347 ndvi with false color composite (fcc) analyzes land changes with very high accuracy. this result is in line with the results reported by oon et al. (2019), das & pandey (2019), murad & pearse (2018), that the fcc can improve accuracy for vegetation analysis. moreover, the fcc is also able to identify bare land in green color. this result is in line with the results reported by khalil & haque (2018), rajani & varadarajan (2018), sharma et al. (2016), that that the fcc is also able to identify bare land. 4. conclusion the accuracy tends to increase when a new band formed in the slack of spatial image. the svm model performed better than nn method, so we can say that supervised learning increases the system's accuracy and performance compared to unsupervised learning. the classification is performed with the spatial vegetation index related to ndvi. it is also observed that if different vegetation indices are used and evaluated on different datasets, there are great chances that the system's accuracy improves to a large extent. figure 13. mean of land-use/land cover and fcc of 2020 r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 348 conflict of interest the authors declare that there is no conflict of interest with any financial, personal, or other relationships with other people or organizations related to the material discussed in the article. references acharya, t. d., & yang, i. (2015). exploring landsat 8. international journal of it, engineering and applied sciences research, 4(4), 4-10. addink, e. a., de jong, s. m., & pebesma, e. j. (2007). the importance of scale in object based mapping of vegetation parameters with hyperspectral imagery. photogrammetric engineering and remote sensing, 73(8), 905–912. https://doi.org/10.14358/pers.73.8.905. ahmad, i., singh, a., fahad, m., & waqas, m. m. (2020). remote sensing-based framework to predict and assess the interannual variability of maize yields in pakistan using landsat imagery. computers and electronics in agriculture, 178 doi:10.1016/j.compag.2020.105732. al-obeidat, f., al-taani, a. t., belacel, n., feltrin, l., & banerjee, n. (2015). a fuzzy decision tree for processing satellite images and landsat data. procedia computer science, 52(c), 1192–1197. https://doi.org/10.1016/j.procs.2015.05.157. an, k., zhang, j., & xiao, y. (2007). object-oriented urban dynamic monitoring a case study of haidian district of beijing. chinese geographical science, 17(3), 236–242. https://doi.org/10.1007/s11769-007-0236-1. aubrecht, c., steinnocher, k., hollaus, m., & wagner, w. (2009). integrating earth observation and giscience for high resolution spatial and functional modeling of urban land-use. computers, environment and urban systems, 33(1), 15–25. https://doi.org/10.1016/j.compenvurbsys.2008.09.007. chi, m., feng, r., & bruzzone, l. (2008). classification of hyperspectral remote-sensing data with primal svm for small-sized training dataset problem. advances in space research, 41(11), 1793-1799. doi:10.1016/j.asr.2008.02.012. das, p., & pandey, v. (2019). use of logistic regression in land-cover classification with moderate-resolution multispectral data. journal of the indian society of remote sensing, 47(8), 1443-1454. doi:10.1007/s12524-019-00986-8. devi, m.s., dewangan, s., ambashta, s. k., jaiswal, a., & kondapalli, s. (2019). recognition of forest fire spruce type tagging using machine learning classification. international journal of recent technology and engineering, 8(3), 4309-4313. doi:10.35940/ijrte.c5176.098319. foody, g. m., & mathur, a. (2004). a relative evaluation of multiclass image classification by support vector machines. ieee transactions on geoscience and remote sensing, 42(6), 1335–1343. https://doi.org/10.1109/tgrs.2004.827257. https://doi.org/10.14358/pers.73.8.905 https://doi.org/10.1007/s11769-007-0236-1 https://doi.org/10.1016/j.compenvurbsys.2008.09.007 r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 349 gangappa, m., mai, d., & sammulal, d. (2016). esdas: explorative spatial data analysis scale to predict spatial structure of landscape. international journal of scientific and engineering research, 7(12), 466–472. https://doi.org/10.14299/ijser.2016.12.001. gangappa, m., mai, d. c. k., & sammulal, d. p. (2017). analysis of advanced data mining prototypes in spatial data analysis'. journal of engineering and applied sciences, 12(12), 3213-3219. gumma, m. k., thenkabail, p. s., teluguntla, p. g., oliphant, a., xiong, j., giri, c., . . . whitbread, a. m. (2020). agricultural cropland extent and areas of south asia derived using landsat satellite 30-m time-series big-data using random forest machine learning algorithms on the google earth engine cloud. giscience and remote sensing, 57(3), 302-322. doi:10.1080/15481603.2019.1690780. ham, j. s., chen, y., crawford, m. m., & ghosh, j. (2005). investigation of the random forest framework for classification of hyperspectral data. ieee transactions on geoscience and remote sensing, 43(3), 492–501. https://doi.org/10.1109/tgrs.2004.842481. oon, a., mohd shafri, h. z., lechner, a. m., & azhar, b. (2019). discriminating between large-scale oil palm plantations and smallholdings on tropical peatlands using vegetation indices and supervised classification of landsat-8. international journal of remote sensing, 40(19), 7312-7328. doi:10.1080/01431161.2019.1579944. pal, m. (2005). random forest classifier for remote sensing classification. international journal of remote sensing, 26 (1), 217–222. https://doi.org/10.1080/01431160412331269698. pozdnoukhov, a., & kanevski, m. (2006). monitoring network optimisation for spatial data classification using support vector machines. international journal of environment and pollution, 28(3–4), 465–484. https://doi.org/10.1504/ijep.2006.011223. rajani, a., & varadarajan, s. (2020). lu/lc change detection using ndvi & mlc through remote sensing and gis for kadapa region. cognitive informatics and soft computing (pp. 215-223). springer, singapore. doi:10.1007/978-981-15-1451-7_24. ramos, a. p. m., osco, l. p., furuya, d.e.g., gonçalves, w.n., santana, d.c., teodoro, l.p.r, . . . pistori, h. (2020). a random forest ranking approach to predict yield in maize with uav-based vegetation spectral indices. computers and electronics in agriculture, 178 doi:10.1016/j.compag.2020.105791 rawat, j. s., & kumar, m. (2015). monitoring land-use/cover change using remote sensing and gis techniques: a case study of hawalbagh block, district almora, uttarakhand, india. the egyptian journal of remote sensing and space science, 18(1), 77–84. https://doi.org/10.1016/j.ejrs.2015.02.002 shanthini, d., shanthi, m., & bhuvaneswari, m. (2017). a comparative study of svm kernel functions based on polynomial. int. j. eng. comput. sci., 6(3). https://doi.org/10.14299/ijser.2016.12.001 https://doi.org/10.1109/tgrs.2004.842481 https://doi.org/10.1504/ijep.2006.011223 r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 350 sharma, r. c., tateishi, r., & hara, k. (2016). a biophysical image compositing technique for the global-scale extraction and mapping of barren lands. isprs international journal of geo-information, 5(12) doi:10.3390/ijgi5120225. singh, r. p., singh, n., singh, s., & mukherjee, s. (2016). normalized difference vegetation index (ndvi) based classification to assess the change in land-use/land cover (lulc) in lower assam, india. international journal of advanced remote sensing and gis, 5(1), 1963–1970. doi:10.23953/cloud.ijarsg.74. somching, n., wongsai, s., wongsai, n., & koedsin, w. (2020). using machine learning algorithm and landsat time series to identify establishment year of para rubber plantations: a case study in thalang district, phuket island, thailand. international journal of remote sensing, 41(23), 9075-9100. doi:10.1080/01431161.2020.1799450. khalil, r.z., & haque, s. (2018). insar coherence-based land cover classification of okara, pakistan. egyptian journal of remote sensing and space science, 21, s23-s28. doi:10.1016/j.ejrs.2017.08.005. liu, b., gao, l., li, b., marcos-martinez, r., & bryan, b. a. (2020). nonparametric machine learning for mapping forest cover and exploring influential factors. landscape ecology, 35(7), 1683-1699. doi:10.1007/s10980-020-01046-0. lu, k. c., & yang, d. l. (2009). image processing and image mining using decision trees. journal of information science and engineering, 25(4), 989–1003. https://doi.org/10.6688/jise.2009.25.4.2. maity, a. (2016). supervised classification of radarsat-2 polarimetric data for different land features. arxiv preprint arxiv:1608.00501. mahmon, n. a., & ya'acob, n. (2014). a review on classification of satellite image using artificial neural network (ann). in 2014 ieee 5th control and system graduate research colloquium (pp. 153-157). ieee. mirici, m. e., satir, o., & berberoglu, s. (2020). monitoring the mediterranean type forests and land-use/cover changes using appropriate landscape metrics and hybrid classification approach in eastern mediterranean of turkey. environmental earth sciences, 79(21) doi:10.1007/s12665-020-09239-1. murad, c. a., & pearse, j. (2018). landsat study of deforestation in the amazon region of colombia: departments of caquetá and putumayo. remote sensing applications: society and environment, 11, 161-171. doi:10.1016/j.rsase.2018.07.003. nachappa, t. g., ghorbanzadeh, o., gholamnia, k., & blaschke, t. (2020). multi-hazard exposure mapping using machine learning for the state of salzburg, austria. remote sensing, 12(17) doi:10.3390/rs12172757. yu, j., li, f., wang, y., lin, y., peng, z., & cheng, k. (2020). spatiotemporal evolution of tropical forest degradation and its impact on ecological sensitivity: a case study in jinghong, xishuangbanna, china. science of the total environment, 727 doi:10.1016/j.scitotenv.2020.138678. r. sanjeeva reddy et al /geosi vol 5 no 3 (2020) 335-351 351 zurqani, h. a., post, c. j., mikhailova, e. a., cope, m. p., allen, j. s., & lytle, b. a. (2020). evaluating the integrity of forested riparian buffers over a large area using lidar data and google earth engine. scientific reports, 10(1) doi:10.1038/s41598020-69743-z. the final step is to get the region of interest of the image shape file is created. the training and testing samples are also extracted for classification. then the training model is built by using machine learning techniques. 3. results and discussion when the land data was pre-processed, class labels are generated geometrically, such as water, bare land, and forest. then region of interest is extracted from the image and is saved in the shapefile. moreover, from the shapefile, training samples are... the accuracy tends to increase when a new band formed in the slack of spatial image. the svm model performed better than nn method, so we can say that supervised learning increases the system's accuracy and performance compared to unsupervised learnin... 186 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 a preliminary study on tsunami disaster in yogyakarta: identification of vulnerability order and components lina wahyuni1*, muh. aris marfai2, m. pramono hadi2 1doctoral program of regional development study program, faculty of geography, gadjah mada university, jl. kaliurang, sekip utara, bulaksumur, yogyakarta, 55281, indonesia 2 regional development study program, faculty of geography, gadjah mada university, jl. kaliurang, sekip utara, bulaksumur, yogyakarta, 55281, indonesia *corresponding author : linawahyuniananda@gmail.com received 10 march 2020/ revised 14 may 2020 / accepted 23 may 2020/ available online 2 june 2020 abstract a tsunami is a disaster that can be hardly estimated. it is a significant concern un since more than 60% of the world's population lives in coastal areas prone to tsunamis, including indonesia. the county community with complex and dynamic plate requires mastering of mitigation strategies as a tsunami preventive effort. understanding the vulnerable elements in risky areas is critical. however, the magnitude of potential disasters cannot be minimized. this study analyzes the tsunami vulnerability in bantul, special region of yogyakarta (diy). the analysis was based on a description of assessment parameters such as land use, the physical condition of the area, social conditions, and availability of infrastructure. the results show that social vulnerability had the most significant impact. keywords: vulnerability, tsunami, bantul, diy 1. introduction a tsunami refers to waves that are faster, taller, and stronger than wind or storm surge (chen & cheng, 2016; rangel-buitrago et al., 2020). it has frequently occurred in the last decade, damaging coastal structures (nandasena et al., 2012). its incidencet on a large scale is relatively less frequent compared to hydrometeorological disasters.the associated waves are unpredictable because they are caused by sudden significant volcanic displacements, initially triggered by earthquakes, landslides, volcanic eruptions, or meteors (al-faesly et al., 2012). on december 26, 2004, a tsunami disaster was triggered by the magnitude of an earthquake with a strength of 9.1 ritcher scale (sr) in the indian ocean. the maximum height of the waves was 30 meters, causing more than 200,000 deaths and massive geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.5 no. 2 (2020), 186-195, august, 2020 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v5i2.17006 accredited by the ministry of research , technology , and higher education of the republic of indonesia, no. 30/e/kpt/2019. mailto:linawahyuniananda@gmail.com https://jurnal.unej.ac.id/index.php/geosi https://drive.google.com/file/d/1rsnvas6cuhowhl5bj87cl2l6k5dqz7s6/view 187 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 destruction of property in more than ten countries bordering the indian ocean (grilli 2007, leonard & lucinda, 2014; iverson & prasad, 2007; roshan et al., 2016). on february 27, 2010, a tsunami disaster was triggered by an 8.8 magnitude earthquake of the coast of chile. the waves reached local run-ups of 29 meters high on coastal cliffs (fritz, 2010). on march 11, 2011, a magnitude 9.0 earthquake stroked near the coast of northeastern japan and swept along the coast, penetrating the land with a maximum height of 40 meters (yeh et al., 2013). on september 16, 2015, an 8.3 magnitude earthquake occurred off chile's central coast and triggered a tsunami with a maximum runoff height of 13 meters (contreras-lopez et al., 2016). in thailand, the dock plates at the port of khao lak and the fishing port of the ban nam kem deck were severely damaged by the uplifted pressure due to the indian ocean tsunami of 2004 (ghobarah et al., 2006). the same incident occurred in japan and damaged the sendai port in the tohoku region during the 2011 tsunami (suppasri, 2012). based on the risk analysis conducted by the national disaster management agency (bnpb) in 2012, four major areas have high risk and probability of tsunami, including mentawai, sunda strait, and southern part of java, megathrust south of bali and nusa tenggara, and northern papua region. of the four areas, the south java coast or pansela, has the largest population (bnpb, 2012). due to a large number of residents in this region, the spatial planning along the south coast of java should be based on coastal area disaster mitigation. the coastal area of the bantul regency is prone to tsunami because it is a low lying area designed as one of the national strategic tourism area (hadipour et al., 2019; mcguire, 2020). the existence of south cross road (jjls) increase the strategic value of the coastal areas. it connects the southern coast of java island to the coast of bantul regency. also, the new airport development plan in the coastal area of the kulonprogo regency is next to the regency, which is integrated with national tourism strategic area (kspn) borobudur and road along jjls. generally, airports are encouraged to support the development of kspn borobudur and surrounding areas. this is stated in the national tourism development master plan 2010-2025. figure 1 shows the diy southern coast strategic region. the implementation of spatial planning should be carried out comprehensively, holistically, well-coordinated, integrated, effective, and efficient, focusing on political, economic, social, culture, defense, security, and environmental sustainability (ibrahim & hegazy, 2013). spatial planning needs to be based on the system approach, main function, 188 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 administration, activity, and strategic value areas, taking into account the disaster factor (muta'ali, 2014). figure 1. kspn pansela yogyakarta area. source: government regulation number 50 the year 2011 concerning master plan of national tourism development year 2010 2025. the development plan of the southern coastal region, including the bantul area, has encouraged unity in the spatial planning following the risk of the tsunami disaster (balasundareshwaran et al., 2020). this study analyzes the tsunami vulnerability orderbased on the assessment parameters, including land use, the physical condition, social state,and infrastructure availability. 189 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 2. methods the identification of dangerous elements in the disaster-prone area is part of mitigation. this aspect was investigated by tanaka, (2008); tanaka et al., (2010); liyanage & lee, (2012); freire et al., (2012); & shibayama et al., (2012). in this research, the risky component was considered a significant factor in mitigation since the magnitude of the disaster cannot be reduced. the assessment of the damage was conducted by identifying and calculating the vulnerability order, including land use, the physical and social condition of the area, availability of infrastructure, and economic. the unit of analysis was the village administration. theidentification of elements at risk within the tsunami danger zone, including physical, social, and economic elements, is a significant step in determining vulnerability order. the parameters used are based on the national disaster management agency (bnpb) regulation no. 2 of 2012 and several previous studies. the research flows as follows. figure 2. research flow 3. results and discussion indonesia is a country prone to the tsunamis, especially in the coastal area that directly faces the meeting layer of eurasia plate, indo-australia, and the pacific, including the western part of sumatra island, the southern part of java island, nusa tenggara, the northern part of papua, sulawesi, maluku, and the eastern part of kalimantan island (bmkg, 190 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 2012; mcguire, 2020).the common disaster occurringis a close-range tsunami of around 200 km from the earthquake epicenter. local tsunamis can be caused by earthquakes,slide, and volcanic eruptions (bmkg, 2012). table 2.1 depicts the tsunami history occurring in indonesia. table 1. significant tsunami event in indonesia no year location magnitude total victims 1 1883 krakatau volcano 36,000 2 1833 west sumatera, bengkulu, and lampung 8.8 unreported 3 1938 kal island, bangka 8.5 unreported 4 1967 tinambung 58 5 1968 tambu, southeast sulawesi 6.0 200 6 1977 sumbawa 6.1 161 7 1992 flores 6.8 2,080 8 1994 banyuwangi 7.2 377 9 1996 toli toli 7.0 9 10 1996 biak 8.2 166 11 2000 banggal 7.3 50 12 2004 nangro aceh darussalam 9.0 250,000 13 2006 pangandaraan 7.2 >600 14 2010 mentawai 7.7 >400 source : mardiatno (2008 ) the greatest tsunami in the history of indonesia occurred in aceh on december 26, 2004. it started by the earthquake magnitude of 9.3 sr, which caused a strong shock and fault, stretching from aceh to andaman. the tsunami was attributed to the earthquake with huge losses and 250,000 deaths (mardiatno, 2008). almost all of the tsunami disasters led to material losses and claimed many lives. according to table 2.1, the most recent tsunami occurred in october 25th 2010 in mentawai island, west sumatra. it started with an earthquake of magnitude 7.7 sr, followed by tsunami waves of 3-10 meters. this caused destruction of 77 villages and more than 400 deaths (mardiatno, 2008). bantul regency has ahigh vulnerability because it directly faces the indian ocean. additionally, the coastal typology tends to be flat (trihatmoko, 2017; mcguire, 2020). when a tsunami strikes, it is likely to damage the physical and social aspects, as well as the existing infrastructure. social vulnerability should be the first concern since it relates to the number of people affected (koroglu et al., 2019; liu et al., 2020; malherbe et al., 2020). the readiness of every resident in the face of disasters significantly affects vulnerability. in case the community is ready to face any disaster, the severity can be reduced. 191 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 social vulnerability is the ability to recover from the impact of natural disasters based on age and sex group of the populations (dawyer, 2004. in zulkarnaen, 2012). it is based on the understanding of the disaster and the resulting conditions. this includes the ability to evaluate when it occurs and the recovery process. the population of women, children, and the elderly is considered the most vulnerable (subarkah, 2009). based on the law of the republic of indonesia number 24, 2007 on disaster countermeasure, the vulnerable groups include infants, toddlers, and children of pregnant and lactating mothers, the disabled and the elderly. according to subarkah. (2009), the components used in calculating social vulnerability include age and sex groups. the vulnerability is assessed based on an understanding of current and post-disaster conditions based on the evacuation capability. in this assessment, the population of women, children, and the elderly are targeted. table 2 shows the population in bantul regency, which can be potentially affected in case the tsunami strikes. this data includes the number of people living in the district directly facing the indian ocean. population density also affects the vulnerability of a region. table 2 shows the pandak district has the highest population and density. table 2. population and population density based on district which directly facing the indian ocean in bantul regency district population population density srandakan 29,130 1,590 sanden 30,114 1,300 kretek 30,111 1,125 pundong 32,321 1,365 bambang lipuro 38,206 1,684 pandak 48,950 2,014 source : statistics of bantul regency (2016) economic vulnerability is the risk of damage, negative impact, or external shock resistance due to unexpected events (koroglu et al., 2019; liu et al., 2020). the calculation losses can be a good indicator of the economy (gulllaumom, 1999, in zulkarnain 2012). physical vulnerability is the last aspect describing the extent of damage to physical infrastructures exposed to hazards, such as residential buildings (ishtiaque et al., 2019). it affects the local community's structural readiness and the condition of structures (prasstiya, 192 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 2013). the main physical vulnerability in bantul regency include jjls and the existence of airport . the most important thing to do in adjusting the development of risk management strategies is to assess vulnerability to potential tsunami damage in the proper order (lantz et al., 2020). however, vulnerability assessment has never been formulated in an appropriate legal document. building vulnerability means calculating thestructures' capacity in the horizontal pressure of tsunami flows and its susceptibility to water (dall'osso et al., 2009; koroglu et al., 2019). 4. conclusion in this study, the ranking of the element at risk involves social, economic, and physical vulnerabilities. social vulnerability is considered to be the most powerful aspect associated with detailed and varied subcomponents. the economic vulnerability comes second in terms of the economic conditions of a society. it is closely related tocommunities’specific capacity to survive. physical vulnerability is rated last since rebuilding of infrastructures is easier than social or economic development. conflict of interest the authors declare that there is no conflict of interest with any financial, personal, or other relationships with other people or organizations related to the material discussed in the article. references al-faesly, t., palermo, d., nistor, i., and a. cornett, a. (2012). experimental modeling of extreme hydrodynamic forces on structural models. int. j. protect. struct, vol. 3 (4), pp. 477–506. balasundareshwaran, a., kumaraswamy, k., balasubramani, k. (2020). multi-hazard zonation for effective management of disasters in tamil nadu. geosfera indonesia, vol. 5(1), doi: https://doi.org/10.19184/geosi.v5i1.16710. bmkg. (2012). pedoman pelayanan peringatan dini tsunami. inatews. edisi kedua. bnpb. (2012). indeks risiko bencana indonesia. badan nasional penanggulangan bencana. jakarta. https://doi.org/10.19184/geosi.v5i1.16710 193 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 chen & cheng. (2016). experimental study of uplift loads due to tsunami bore impact on a wharf model. international journal of coastal engineering, vol. 117 (2016), pp. 126–137. contreras-lopez, m., winckler, p., sepúlveda, i, andaur-álvarez, a, cortés-molina, f., guerrero, c. j., mizobe, c.e., igualt, f., breuer, w., beyá, vergara, h., figueroasterquel, r. (2016). field survey of the 2015 chile tsunami with emphasis on coastal wetland and conservation areas. pure appl. geophys. vol. 173 (2), pp. 349367. dall’osso, f.m., gabbianelli, g., withycombe, g., goff, j., dominey-howes, d.a. (2009). revised (ptva) model for assesing the vulnerability of building to tsunami damage. natural hazard earth science. 9:1557-1565. freire, s., aubrecht, c., wegscheider , s. (2012). advancing tsunami risk assessment by improving spatio-temporal population exposure and evacuation modeling. natural hazards 68:1311-1324. fritz, h.m. (2010). field survey of the 27 february 2010 chile tsunami. pure appl. geophys, vol. 168 (11), pp. 1989–2010. ghobarah, a., saatcioglu, m., nistor, i. (2006). the impact of the 26 december 2004 earthquake and tsunami on structures and infrastructure. eng. struct, vol. 28 (2), pp. 312–326. government regulation number 50. (2011) concerning master plan of national tourism development year 2010 2025. jakarta : president of the republic of indonesia grilli, s.t. (2007). “source constraints and model simulation of the december 26, 2004, indian ocean tsunami”. port coast. ocean eng., vol. 133 (6), pp. 414–428. haipour, v., vafaie, f., kerle, n. (2019). an indicator-based approach to assess social vulnerability of coastal areas to sea-level rise and flooding: a case study of bandar abbas city, iran. ocean & coastal management, vol 188., 105077. ibrahim, h.s., hegazy, i. (2013). decentralization in the egyptian coastal management. journal of coastal zone management, vol 16(2), pp. 102-113. ishtiaque, a., eakin, h., chhetri, n., myint, s.w., dewan, a., kamruzzam, m. (2019). examination of coastal vulnerability framings at multiple levels of governance using spatial mcda approach. ocean & coastal management, vol. 171, pp. 66-79. iverson, l.r., & prasad, a.m. (2007). using landscape analysis to assess and model tsunami damage in aceh province, sumatra. landscape ecol, vol. 22, pp. 323–331. koroglu, a., ranasinghe, r., jimenez, j.a., dastgheib, a. (2019). comparison of coastal vulnerability index applications for barcelona province. ocean & coastal management, vol. 178, 104799. 194 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 lantz, t.c., moffat, n.d., jones, b.m., chen, q., tweedie, c.e. (2020). mapping exposure to flooding in three coastal communities on the north slope of alaska using airborne lidar. coastal management, vol. 48(2), pp. 96-117. leonard, l. & lucinda (2014). tsunami hazard assessment of canada. nat hazards, vol. 70, pp. 237–274, 2014. liu, r., pu, l., zhu, m., huang, s., jiang, y. (2020). coastal resource-environmental carrying capacity assessment: a comprehensive and trade-off analysis of the case study in jiangsu coastal zone, eastern china. ocean & coastal management, vol. 186, 105092. liyanage, a.n. & lee, h. (2012). a schematic approach for gis application for tsunami disaster management. asea/drbc, ccis 340, pp. 395-400. marherbe, w., sauer, w., aswani, s. (2020). social capital reduces vulnerability in rural coastal communities of solomon islands. ocean & coastal management, vol. 191, 105186. mardiatno, d. (2008) risiko tsunami di pantai selatan jawa: belajar dari kejadian tsunami di banyuwangi pada tahun 1994 dan di pengandaran pada tahun 2006. jurnal kebencanaan indonesia, vol. 6 (2), pp. 15-32 mcguire, c.j. (2020). coastal resiliency through value capture and transfer: a framework proposal. coastal management, vol. 48(2), pp. 57-76. muta’ali, l. (2014). development planning based on disaster risk reduction, (in bahasa perencanaan pengembangan wilayah berbasis pengurangan risiko bencana). badan penerbit fakultas geografi, universitas gadjah mada, yogyakarta. nandasena, n., sasaki, y., tanaka, n. (2012). “modeling field observations of the 2011 great east japan tsunami: efficacy of artificial and natural structures on tsunami mitigation. coast. eng, vo. 67, pp. 1–13. prasstiya, r.t. (2013). manajemen risiko tsunami di pulau-pulau kecil (gilitrawangan) kecamatan pemenang kabupaten lombok utara. master thesis.universitas gadjah mada, 2013. rangel-buitrago, n., neal, w.j., bonetti, j., anfuso, g., de jonge, v.n. (2020). vulnerability assessments as a tool for the coastal and marine hazards management: an overview. ocean & coastal management, vol. 189, 105134. roshan, a.d., basu, p.c., jangi, r.s. (2016). tsunami hazard assessment of indian coast. nat hazards, vol. 82, pp. 733–762. shibayama, t., esteban, m., nistor, i., takagi,h., thao, n.d., matsumaru, r., mikami, t., aranguiz, r., jayaratne, r., ohira k. (2012). classification of tsunami and evacuation areas. natural hazards 67:365-386. statistic of bantul regency. (2016). bantul dalam angka. bantul: bps. 195 lina wahyuni et al / geosi vol 5 no 2 (2020) 186-195 subarkah. (2009). spatial multi criteria evaluation for tsunamis vulnerability: case study of coastal area parangtritis, yogyakarta, indonesia. master thesis, universitas gadjah mada, yogyakarta. suppasri, a. (2012). damage characteristic and field survey of the 2011 great east japan tsunami in miyagi prefecture. coast. eng. j., vol. 54 (01). tanaka, n. (2008). vegetation bioshields for tsunami mitigation: review of effectiveness, limitations, construction, and sustainable management. landscape ecol eng 5:71-79. tanaka, n., jinadasa, k. b. s. n., mowjood, m.i.m., fasly, m.s.m. (2010). coastal vegetation planting projects for tsunami disaster mitigation: effectiveness evaluation of new establishments. landscape ecol eng 7:127-135. trihatmoko, e. (2017). proses dan dampak dinamika wilayah kepesisiran jawa tengah dan daerah istimewa yogyakarta. master thesis. universitas gadjah mada, yogykarta. yeh, h., sato, s., tajima, y. (2013). the 11 march 2011 east japan earthquake and tsunami: tsunami effects on coastal infrastructure and buildings. pure appl. geophys, vol. 170 (6–8), pp. 1019–1031. zulkarnain, m. w. d.(2012). evaluasi multi-kriteria keruangan untuk penilaian risiko total tsunami di pacitan. master thesis. universitas gadjah mada, yogyakarta. luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 175 improving the learning outcomes of knowledge and inquiry skill domain of third grade students at smp negeri 14 kendari using guided inquiry learning model assisted by science kit luh sukariasih1*, i gede purwana edi saputra2, fahrudi ahwan ikhsan3, andri estining sejati4, khaerun nisa1 1department of physics education, halu oleo university, indonesia 2department of physics education, sembilanbelas november kolaka university, indonesia 3department of geography education, university of jember, indonesia 4department of geography education, sembilanbelas november kolaka university, indonesia *email: luhsukariasih@yahoo.com received 29 march 2019/ revised 13 august 2019/ accepted 18 august 2019/ published 23 august 2019 abstract the study aims to improve the learning outcomes in the field of knowledge and inquiry skill in class viii 5 smp negeri 14 kendari on the subject matter of light in atmosphere as the effect of applying the guided inquiry learning model assisted by science kit. the method of the study used a classroom action research with research design is cycle model. the research subject is the students of class viii 5 smp negeri 14 kendari in the academic year 2016/2017 which consist of 26 students. the learning data achievements of the learners' realm were obtained through the learning result test (cycle test), the skill data of the learners were obtained through the inquiry sheet, and then was analyzed used the descriptive statistics. results of data analysis are: 1) learning outcomes increased from 60,31 in cycle i to 75 in cycle ii; 2) the students group inquiry skill increased form average value 2.68 (enough category) in the cycle i to 3.15 (good category) in cycle ii; 3) the students mastery learning percentage increase from 42.31% (11 students) in cycle i to 77% (20 students) in cycle ii. it could be concluded that the implementation of guided inquiry learning model assisted by science kit could improve the learning outcomes of knowledge and inquiry skill domain on class viii 5 smp negeri 14 kendari in the subject matter of light in atmosphere. keywords: guided inquiry, inquiry skills, learning outcomes,science kit. 1. introduction education science is a branch of science that is built on observation and classification of data, and is usually compiled and verified in quantitative laws, which involve the application of mathematical reasoning and data analysis to natural phenomena. natural events are included in the study of physical geography. in essence, education science is a science of natural phenomena that is poured in the form of facts, concepts, principles, and laws that are validated and through a series of activities in the scientific method (depdiknas, 2004). according to aksa et al. (2019), geography spesification are physical geography, social geography, and technical geography. geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.4 no. 2 (2019), 175-187, august, 2019 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v4i2.10097 luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 176 one of physical geography study way through mathematics or counting. this is often feared and tends to be disliked by most students which has an impact on the low level of understanding of teaching material. according to sohibun (2014), learning with a count is still to understanding in the form of conceptsand formulas. the results of preliminary observations made on viii class students in smp negeri 14 kendari on 6th december 2016, found low learning outcomes in the realm of students' knowledge especially in the subject matter of light in atmospherein class viii 5. this can be seen from the average daily test results obtained by the value of 67.3. 13 out of 23 students (57%) scored below the minimum completion criteria (68 points) and only 10 students (43%) completed. the results of interviews with teachers in class viii 5 obtained learning outcomes in the realm of students' knowledge more on factual knowledge, while for conceptual knowledge, procedural knowledge, and metacognitive knowledge were only a few that students could understand. students rarely ask questions in the learning process, meaning that the ability of students to be involved in formulating problems or questions is also lacking. experimental activities in atmosphere light material that lack make competent skills in this case the skills to carry out experiments, collecting observational data and formulating conclusions are also lacking. according to anderson & krathwohl (2015), the knowledge dimension consists of: four types: factual knowledge, conceptual knowledge, procedural knowledge, and metacognitive knowledge. the low average learning outcomes of atmosphere light material according to observation and interviews viii 5 class smp negeri 14 kendari caused in the learning process the model used by the teacher is teacher centered. teacher is very active while the students tend to be passive in the learning process which makes learning monotonous and less student participation. this learning process makes it students difficult to understands the material taught by the teacher. teacher teach less of experiment activities on material that requires experimental activities made students think abstractly about the material. according to sanjaya (2014), in conventional learning the position of students as objects of learning with the role as recipients of information passively. the students are only faced with concepts and formulas without any learning activities that demand the full activity of students in understanding the knowledge taught. the learning process at the junior secondary level is still dominated by teachers compared to students with one-way communication especially in the viii 5 class smp negeri 14 kendari. this is thought to be one of the factors causing low understanding of subject matter for junior high luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 177 school students. according to amaliana (2017), in teacher-centeder learning, teachers play important roles in the learning process as information providers or evaluator, and students are viewed as learners who passively receive information. the solution about the lack of student learning process and outcome is using leaning model. one of learning model that can make students more active and the learning process more interesting is the guided inquiry. the guided inquiry emphasizes students to be active and find their own knowledge. this relates to conceptual knowledge and procedural knowledge is easier to understand if students can find their own knowledge through experiments. a good learning process is that students are more active than teachers, because the learning process like this is more interesting and the material is easier to understand. according to hidayati et al.(2016), in formulate problems and make hypotheses stages the student become active thinking direcly confronts the issues to be resolved. according to hardianti & kuswanto(2017), learning through inquiry gives learnenrs independence by encouraging them to have a more active and responsible role in various stages of investigation. guided inquiry can improve the inquiry skills and students learning outcome. guided inquiry can more effective if use experimental equipment in the school like science kit tools. there is science kit tools in the science laboratory in smp negeri 14 kendari. based on the analysis of the problem, the teacher obtained the appropriate model to solve the student problem in the class with the guided inquiry learning model assisted by science kit. according to greenwald & quitadamo (2014), inquiry learning model combine with clinical case namely ibcc in the health matter in the neuro anatomy cursus students learn by using process skills, attitudes, and rational thinking knowledges. according to kuhlthau et al.(2015), inquiry directly connects between matter and the real world. many people think that the light in atmosphere material is quite difficult because some of the sub-material makes students think abstractly so that there is needed for direct activity to finding concepts. the use of guided inquiry learning assisted by science kit is very suitable in learning atmosphere light material. guided inquiry learning model is learning that involves all the ability of students to search and investigate systematically and logically so that they can formulate their own findings. guided inquiry learning will be more effective with practical activities in this case utilizing science kit. application of this model is carried out in two cycles of learning in the classroom to improve learning outcomes in the realm of students' knowledge and skills. guided inquiry learning model that is a model that emphasizes more on students to actively train courage, luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 178 communicate and try to find their own knowledge to solve problems faced through an experiment with direct guidance by the teacher. the combination of inquiry and science kit models is very suitable, because this model will make students find their own knowledge through direct activities of students in this case the experiment using science kit teaching aids. according to piaget (1970), by using real experience, a person's cognitive development will be better than just using language to communicate. according to sapriya (2014), guided inquiry learning has the advantage of developing thinking skills, knowledge, attitudes, and values in students compared to the classical approach. the research by ahmadi (2015) conducted at sdn 1 telaga, gorontalo district showed that learning by using kit teaching aids in science learning can make students more easily understand the material being taught, and students are more enthusiastic in learning learning. based on this description, the purpose of this study is to improve learning outcomes in the realm of knowledge and skills in the form of students ofclassviii 5 of smp 14 kendari in learning the subject matter of light in atmosphere through a guided inquiry guided model science kit. 2. the methods this study used a type of classroom action research. the research was conducted from 10 may to 24 may 2017 in 2016/2017 academic year with the subject matter of light in atmosphere. the subjects of this study were all students on class viii 5 in the even semester of smp negeri 14 kendari with 26 students. the research design is a cycle model. before each cycle phase is carried out first a preliminary study is carried out. the data in this research were analyzed using descriptive statistics. determination of the value of the real learning outcomes of students’ knowledge of the range of scores used for the description test in this study is 0 to 100 with the formula (tim direktorat pembinaan smp, 2017). the average value of learning outcomes in the knowledge class ( ), standard deviation (sd),% complete, and categorizing the value of learning outcomes with formula (sudjana, 2014). normalized gain (n-gain) is an increase in learning outcomes in the realm of students knowledge and skills in the first cycle and second cycle of inquiry are determined using equations from (riduwan, 2015). n-gain criteria ≥ 0.7 (high), 0.3 ≤ n-gain> 0.7 (medium), n-gain <0.3 (low). calculating and classifying the average score of inquirt skills associated with students with formula (sugiyono, 2014).below is reseach flow diagram. luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 179 figure 1. classroom action research flow diagram (arikunto, 2016) 3. results and discussion 3.1. inquiry skills of student groups the summary of the results of data analysis in the inquiry group skills of students in the teaching-learning process in each cycle can be seen in table 1. table 1. results of data analysis of inquiry skills for student groups per cycle no observed aspects the average value of inquiry skills student groups per cycle n-gain cyclei categories cycle ii categories 1. formulate problems 2,91 satisfactory 3,12 good 0,19 2. formulate a hypothesis 2,25 satisfactory 3 good 0,43 3. experiment 2,67 satisfactory 3,25 good 0,44 4. collecting data 3,08 good 3,37 good 0,32 5. formulating conclusions 2,5 satisfactory 3 good 0,33 the average value of inquiry skills in groups of students 2,68 3,15 0,34 categories satisfactory good medium resources:riduwan (2015), sanjaya (2014), sugiyono (2014) based on table 1, it can be seen that in the first cycle there were four aspects of inquiry skills in groups of students who obtained an average value with satisfactory beginning reflextion planning action observation reflextion cycle 1 planning action observation reflextion cycle 2 luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 180 categories, namely aspects of formulating problems, aspects of formulating hypotheses, aspects of conducting experiments and aspects of formulating conclusions. this shows that in the first cycle, the skills of the inquiry group of students was not optimal. in the second cycle the skills of inquiry students groups in each aspect had increased. the lowest average value of inquiry group skills of students in the cicle i, that’s mean the aspect of formulating a hypothesis increases in the second cycle to 0.75. while the average value of inquiry skills in the highest students group in the first cycle, that’s mean the aspect of collecting data increased by cycle ii to 0.29. in the second cycle of the 5 aspects of inquiry skills the observed group of students had an average score of 3.15 which good categorized. the average value of inquiry group skills in the first cycle and second cycle can be seen in table 1. from the table, the increase in student inquiry skills is shown by the average n-gain which is categorized as a moderate increase and the average value of inquiry skills group of students experienced an increase of 0.47. 3.2 teacher activity an overview of teacher activities in managing learning using a guided inkuri learning model assisted by science kit, teacher activities in each cycle can be seen in table 2. table2.results of data analysis of teacher activities on each cycle no aspects observed value average cycle i cycle ii a. introduction 1 opening the lessons and check the readiness of students 4 4 2 give apperception to students 3 4 3 deliver/write down topics and learning objectives. 3 3 b. main activity 4 presenting problems so students can formulate problems 4 4 5 direct students to gather information according to what they see at the stage of presenting the problem 3 4 6 direct students to get information through experiments 2,5 3 7 direct students to formulate explanations based on the results of the experiment 2,5 3 8 direct students to analyze experimental data in the form of conclusions 3 4 c. post activity 9 conclusion on the results of the activity 2,5 4 10 the teacher informs the material for the next meeting 3 3 d. class atmosfer 11 enthusiastic students 2,5 3 12 enthusiastic teacher 3 3 13 time according to allocation 3 4 14 kbm according to rpp 2,5 3 average value 2,9 3,5 category enough good luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 181 resousces: sanjaya (2014), sugiyono (2014) based on table 2 it can be seen that in the first cycle the average value of the teacher' lowest activity was the activity of the teacher directing students to obtain information through experiments, directing students to formulate explanations based on the results of the experiment and the conclusion of the results of the activity. while the highest teacher activity in the first cycle was when opening lessons and checking the readiness of students and presenting problems so that students can formulate problems. in the second cycle, it was seen that every aspect observed was maintained and experienced an increase. one aspect of teacher activity with the lowest average value in cycle i is that the teacher directs students to obtain information through experiments increasing in the second cycle of 0.5. in addition, the lowest average value of other teacher activities in the first cycle is the conclusion of the results of increased activities in the second cycle of 1.5. table 2, shows an increase in teacher activity from cycle i to cycle ii where the average value of teacher activity has increased by 0.6. 3.3 learning outcomes of students knowledge data on the realm learning outcomes of students' knowledge is obtained by using learning outcomes tests. based on the descriptive analysis of the learning outcomes of the realm of students 'knowledge on the subject matter of atmosphere light shown in the form of cycle tests consisting of cycle i tests and cycle ii tests, a summary of the results of data analysis of the learning outcomes of students' knowledge per cycle is as follows. table 3.learning outcomes data analysis results of student knowledge sphere per cycle no. value cycle n-gain i ii 1 minimum 20 22 -0,272 2 maximum 78 95 0,785 3 average 60,31 75 0,345 4 deviation standard 17,26 14 0,254 amount of completion 11 20 amount of incomplete 15 6 % complete 42,31 77 % incomplee 57,69 23 resources: riduwan ( 2015), sudjana (2014) tim direktorat pembinaan smp (2017) from table 3 above, it can be seen that the knowledge learning outcomes of viii5 students of smp negeri 14 kendari on the subject matter of atmosphere light learning luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 182 through the application of a guided inquiry learning model assisted by science kit indicate an increase from cycle i to cycle ii. this is seen in the value of the average learning outcomes in the realm of knowledge increasing from cycle i to cycle ii, it can be seen from the average n-gain of 0.345 in the medium category. the percentage of learning completeness from cycle i increased to cycle ii to be above 75% who had achieved kkm, showing completeness classically from classroom action research has been fulfilled which means that the model of guided inquiry teaching assisted by science kit can solve learning problems of viii5 students in smp negeri 14 kendari. 3.4 cycle i the results showed that in the first cycle there were still some aspects of the skills of the inquiry group of students that needed to be improved. the skills of inquiry groups of students who are still lacking and need to be improved are in the aspects of formulating problems, formulating hypotheses, conducting experiments and formulating conclusions. this is because students are not used to formulating problems, formulating hypotheses and formulating conclusions with learning models with the help of applied science kit. the low average value of group inquiry skills is influenced also by the teacher who still lacks directing students in conducting experiments and formulating conclusions during the learning process. the average value of the problem formulation is higher than the average value of the aspects of formulating the hypothesis, this is due to students' initial knowledge of the concepts related to the material in presenting problems that are still lacking so that students have difficulty when formulating hypotheses. according to hidayati et al. (2017), teachers should be able to prepare the questions about the problem problem is close to student environment. in the first cycle based on descriptive analysis, the teacher's activity shows the average value of teacher activity in a sufficient category. where teacher activity is still quite based on. table 2 of which is the teacher is still lacking in directing students to obtain information through experiments, directing students to formulate explanations based on the results of the experiment and formulating conclusions on the results of the activities. based on the results of reflection on teacher activity, by knowing the shortcomings in the first cycle, the teacher improved the way to teach learning material in accordance with the guided inquiry learning model assisted by science kit, so that it was expected that in the next meeting an increase in teacher activity. according to sitorus et al. (2017), in guided learning teacher can manage and implement a learning to students which contains scientific steps. luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 183 based on the first problem, namely how to improve the learning outcomes of knowledge in class viii5 smp negeri 14 kendari as the effect of applying a guided inquiry learning model assisted by science kit on the subject matter of light in atmosphere, it can be explained that based on the results of descriptive analysis of learning outcomes in the knowledge domain each cycle tends to increase towards the better. according to zaini (2016), there is a significant effect that is given by using guided inquiry based learning with the cognitive process learning outcome. the percentage of completeness in this first cycle has not yet reached the research target, which is to achieve learning completeness in a classical minimum of 75%. this is because students have not been able to recall the material being taught, besides that students pay less attention and listen to the initial problems conveyed by the teacher, students are not accustomed to formulating conclusions, and during the experiment the students are less cooperative with their group friends. in addition, the teacher is not optimal in delivering direction and guidance to students in conducting experiments and formulating explanations in the form of conclusions so that students experience obstacles to conducting experiments and formulating conclusions on the results of experiments. according to wardani et al. (2015), cooperation can be increased through guided inquiry learning. after analyzing and reflecting on the first cycle, the teacher in this case the researcher and observer make improvements in learning through guided inquiry learning assisted by science kit to be applied in the second cycle, so that the realm of learning outcomes of knowledge of viii5 students in smp negeri 14 kendari can increase as expected in cycle ii. 3.5 cycle ii from the results of the descriptive analysis of the skills of inquiry groups of students in the second cycle showed an increase in the skills of inquiry groups of students from cycle i. this is as shown in table 1, the average skill scores of groups in the second cycle were 3.15 by category well. the skills aspect of inquiry groups in the second cycle on average experienced an increase from the first cycle. improving the skills of the best group inquiry inquiry was on the aspect of formulating a hypothesis, where the aspect of formulating this hypothesis increased by 0.75 higher than the increase in other aspects of inquiry skills. . this is because students are getting used to formulating hypotheses, assisted by teachers with clearer direction and presenting clearer problems also in the learning process using a guided inquiry learning model assisted by science kit. according to sarwi&prayitno (2016), there is luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 184 a significant application of guided learning model to the students inquiry 'understanding of the concept and being able to improve practice in inquiry syntax. the inquiry process is characterized by the presentation of problems that can be solved in groups through experiments that develop ideas and thoughts and activities centered on problem solving skills. in this learning model students make questions that guide the next investigation, including activities to form hypotheses, carry out experiments, collect experimental data, and formulate conclusions experiments. the teacher acts as a facilitator. the increase in the average value of the inquiry group skills indicates that the weaknesses or weaknesses found in the first cycle can be overcome so that the skills of the inquiry group of students are obtained as expected. according to putra et al. (2016) the application of guided inquiry trains students to pass the learning syntax well and is able to improve students science literary skills. in the second cycle the teaching activity of the teacher showed a significant increase, where in the second cycle the average value of the teacher's activity obtained an average value of 3.5 which was categorized as good. the results of the analysis and observations in the second cycle showed an increase in teacher activity by implementing a guided inquiry learning model assisted by science kit. increasing the average value of teacher activity indicates that the weaknesses contained in the first cycle can be overcome so that the teacher can manage learning using a guided inquiry learning model assisted by science kit. according to niana et al.(2016), the teacher's activities in applying the guided inquiry model were able to improve student learning success. from the results of descriptive analysis of the learning outcomes of the realm of students' knowledge shows an increase in learning outcomes in the realm of knowledge from cycle i to cycle ii. this can be seen with the average value obtained by students in the first cycle increased in the second cycle. improving the learning outcomes of the students 'knowledge in the second cycle showed an increase in students' mastery of learning material and the motivation of students to attend learning until the last meeting. according to mulyana et al. (2018), there is a significant influence guided inquiry learning the student learning result. the increase in learning outcomes is also because the teacher has been able to manage the learning process, besides that students are also familiar with the guided inqury learning model assisted by science kit in the learning process, because in the guided inquiry learning model assisted by science kit students can observe directly the concepts that are taught by conducting experiments and can be more creative looking for answers to the problems given. luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 185 so that students are not difficult to understand the concept of science, especially in the subject matter of atmosphere light and relate it to their experience in everyday life. according to yewang et al. (2016) guided inquiry is better than free inquiry in guiding students in investigations, observations and hypothesestests. 4. conclusion the conclusion in this study is the application of a guided inquiry learning model assisted by science kit can improve learning outcomes in the aspect of knowledge, skills and learning completeness of students of class viii 5 of smp negeri 14 kendari on the subject matter of atmosphere light. this conclusion is supported by several specific conclusions of the results of descriptive empirical data analysis as follows. a. the skills of inquiry in groups of viii 5 students in smp negeri 14 kendari through the assisted guided inquiry learning model of kit ipa have increased. this can be seen from the average value of the skills aspects of the inquiry group of students in the first cycle and second cycle. where in the first cycle obtained an average value of 2.68 sufficient categories, increased in the second cycle to 3.15 good categories. b. the realm of learning outcomes of knowledge of grade viii 5 students of smp negeri 14 kendari who learned through assisted guided inquiry learning model kit ipa increased from the first cycle to the second cycle which was indicated by the average value of the realm learning outcomes of students' knowledge increased from 60.31 to 75, the standard deviation decreased from 17.26 to 14 and the average n-gain was 0.34 medium category. c. the learning completeness of grade viii 5 students of smp negeri 14 kendari who learned through a guided inquiry learning model assisted by science kit showed that there was an increase in the percentage of learning completeness from cycle i to cycle ii, in the first cycle the completeness percentage was 42.31% or 11 students had reached minimum completeness criteriaand in cycle ii the percentage of completeness increased to 77% or 20 students had reached. references ahmadi, l. (2015). pemanfaatan alat kit pada pembelajaran.gorontalo: universitas negeri gorontalo. aksa, f.i., utaya, s., & bachri, s. (2019). geografi dalam perspektif filsafat ilmu. majalah geografi indonesia,33(1), 43-37. luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 186 amaliana, i. (2017). teacher-centered or student-centered learning approach to promote learning?. jurnal sosial humaniora, 10(2), 59-70. anderson, l.w., & krathwohl, d.r. (2015). kerangka landasan untuk pembelajaran, pengajaran, dan asesmen revisi taksonomi pendidikan bloom (a. priantoro, trans.). yogyakarta: pustaka pelajar. arikunto, s. (2016). prosedur penelitian suatu pendekatan praktik. jakarta: rineka cipta. depdiknas. (2004). kurikulum mata pelajaransains smp danmts.jakarta: depdiknas. greenwald, r.r.,&quitadamo, i.j. (2014). a mind of their own: using inquiry-based teaching to build critical thinking skills and intellectual engagement in an undergraduate neuroanatomy course. the journal of undergraduate neuroscience education, 12(2), 100-106. hardianti, t., & kuswanto, h. (2017). difference among levels of inquiry: process skills improvement at senior high school in indonesia. international journal of instruction, 10(2), 119-130. hidayati, d.n., amaluddin, l.o., & surdin. (2016). the effect guided inquiry to critical thinking ability to build student character in geography subject. social science, education and humanities research,9(1), 367-371. kuhlthau, c.c., maniotes, l.k., & caspari, a.k. (2015). guided inquiry: learning in the 21st century(2nd ed.). california: libraries unlimeted. mulyana, s., rusdi, & vivanti, d. (2018). the effect of guided inquiry learning model and scientific performance on student learning outcome. indonesian journal of science and education, 2(1), 105-109. niana, r., sarwanto, & ekawati, e.y. (2016). the application of guided inquiry model on physic learning to improve scientific attitude and students analysis ability. proceedings of the 2nd international conference on teacher training and education sebelas maret university,2(1), 605-615. piaget, j. (1970). science of education and the psychology of the child. new york: wiley. putra, m.i.s., widodo, w., & jatmiko, b. (2016). the development of guided inquiry science learning materials to improve science literacy skill of prospective mi teachers.jpii, 5(1), 83-93. riduwan. (2015). dasar-dasar statistika. bandung: tarsito. sanjaya, w. (2014). strategi pembelajaran. jakarta: kencana prenada media group. sapriya. (2014). pendidikan ips. bandung: remajarosdakarya. sarwi, s.,& prayitno, w.w. (2016). implementation of guided inquiry physics instruction to increase an understanding concept and to develop the students character conservation. jurnalpendidikanfisika indonesia,12(1), 1-7. luh sukariasih et al/ geosi vol. 4 no. 2 (2019) 175-187 187 sitorus, h.h., hasruddin, & edi, s. (2017). the influence of inquiry learning model on student’s scientific attitudes in ecosystem topic at mts. daarul hikmah sei alim (islamic junior high school) asahan. international journal of humanities social sciences and education (ijhsse), 4(11), 170-175. sohibun. (2014). penerapan strategi belajar dengar lihat kerjakan (delikan) berbasis laboratorium mini terhadap ketrampilan proses sains (kps) siswa sma kelas x mia. jurnal imliah edu research, 3(1), 53-67. sudjana, n. (2014). penelitian hasil proses belajar mengajar. bandung: remaja rosdakarya. sugiyono. (2014). metode penelitian pendidikan pendekatan kuantitatif, kualitatif, dan r&d. bandung: alfabeta. tim direktorat pembinaan smp. (2017). panduan penilaian oleh pendidik dan satuan pendidikan sekolah menengah pertama. jakarta: kementerian pendidikan dan kebudayaan direktorat jenderal pendidikan dasar dan menengah direktorat pembinaan sekolah menengah pertama. wardani, s., nurhayati, s., & safitri, a. (2015). the effectiveness of the guided inquiry learning module towards students’ character and concept understanding. international journal of science and research (ijsr), 5(6), 1589-1594. yewang, m.u.k., degeng, i.n.s., setyosari, p., & sulton. (2016). the effect of guided inquiry learning method vs free inquiry against learning outcomes. international conference on educationum, 561-568. zaini, m. (2016). guided inquiry based learning on the concept of ecosystem toward learning outcomes and critical thinking skills of high school student. iosr journal of research & method in education (iosr-jrme), 6(6), 50-55. foreword 65 multi-hazard zonation for effective management of disasters in tamil nadu a. balasundareshwaran1, k. kumaraswamy 2 and k. balasubramani 1* 1 department of geography, central university of tamil nadu, thiruvarur, tamil nadu 610005, india 2 department of geography, bharathidasan university, tiruchirappalli, tamil nadu 620024, india *corresponding author: geobalas@cutn.ac.id received 6 february 2020/ revised 16 february 2020/ accepted 1 march 2020/ published 10 april 2020 abstract natural hazards are a long existing threat to human and their surroundings which may occur throughout the world. tamil nadu is one of the indian states with a number of natural hazard incidences. the occurrence of natural hazards, such as cyclone, storm surge, flood, drought, landslide, forest fire etc., has increased manifold in the recent decades.the multi-hazard zonation is one of the preliminary studies in disaster management scenario, which is used to understand the product of all prominent natural hazards. at the state level, it is imperative for the government to know the regions affected by multiple hazards to help them prepare the management plans appropriately to protect the local communities and infrastructures. however, such systematic hazard assessment and integration in an administrative unit is largely missing in tamil nadu. further, t he utilisation of geoinformatics in the preparation of multi-hazard zonation helps to identify the most endangered areas of the state precisely and offers insights to detailed studies in highly risk zones. this paper attempts on these lines to prepare multi hazard zones (mhz) based on natural hazards viz. earthquake, landslide, cyclone, storm surge, flood, drought and forest fire of tamil nadu. the data for the study were generated from multiple sources, which were all generalised and integrated in a normalised scale. the occurrences, intensities and frequencies of hazards, namely seismic, landslide, and forest fire are the reasons for a very high multi-hazard in hilly tracts of the nilgiris and parts of shayadhri hills in coimbatore and theni districts, whereas cyclone, storm surge, and flood caused a very high risk along the coastal stretch of chennai, kancheepuram, cuddalore and ramanathapuram districts. the segregation of the results into administrative division’s which was then categorised in an order of high risk zones may provide a powerful tool to the state authorities to allocate fund and resources. the output of this study also offers zonation for immediate knowledge, policy briefs, and proper disaster management plan at state level. keywords : disaster management; natural hazards; geoinformatics; tamil nadu 1. introduction according to the disaster statistics from undp’s disaster risk reduction programme (2012), the average disasters per year have increased for more than 60 per cent in the last decade. geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.5 no. 1 (2020), 65-79, april, 2020 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v5i1.16710 accredited by the ministry of research , technology , and higher education of the republic of indonesia, no. 30/e/kpt/2019. 65 mailto:geobalas@cutn.ac.id https://jurnal.unej.ac.id/index.php/geosi https://drive.google.com/file/d/1rsnvas6cuhowhl5bj87cl2l6k5dqz7s6/view a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 asia has the highest number of victims due to natural disasters. it is not peculiar at all that the increase in the frequency of disaster, its magnitude, complexity, economic impact and number of deaths, registers in the low human development region and less prepared communities. according to the ndma (2016), hazards pose threats to people and assume serious proportions in the under developed countries with dense population. the regional study of hazards helps us understand the influencing factors, analysing them without exclusion of nearby quantifiable factors, provides applicable regional alike solutions. there is a coarse approach followed in disaster plans at state and local level in india where national level zonation considered for all kinds of disaster risk management without detailed regional or local attempts. tamil nadu, the southern-most state of india, nestles in the indian peninsula between the bay of bengal in the east, the indian ocean in the south, and the western ghats and the arabian sea on the west. in the north and west, the state adjoins karnataka, andhra pradesh, and kerala (figure 1). the eastern extremity of the state is point calimere situated at 80°20’ e longitude, while the western tip is the mudumalai sanctuary at 71°15’ e longitude. the northern and southern extremities are defined by pulicat lake (13°35’ n latitude) and cape comorin in kanyakumari (08°50’ n latitude). an effective disaster management program always requires a detailed hazard zone, which is largely omitted in the state. at the same time, natural hazards that rattle the vulnerable communities of tamil nadu are increasing manifold and often occur simultaneously. the natural hazards of tamil nadu occur mostly isolated with seasonal characteristics and seldom occured all together that poseda maximum threat in a combined form within the same region. drought is considered one of the reoccurring disasters of tamil nadu, which hampers the development of the region. even though drought is non-structural, it spreads over a larger geographical area and posesmore serious threats to population than the resultsof other natural hazards such as floods, tropical storms, and earthquakes (balasubramani, 2014). balasundareshwaran, et al., (2019a) evaluate drought within an administrative unit and suggested that, the drought restricts the growth of a region. similarly debated that, flood both gains and losses, where it either sweeps the entire regions or deposits minerals for next seasonal irrigation. likewise, cyclones bring in gusty winds to the coastal parts of tamil nadu, especially during the north-east monsoon season (october-december), which influence the sea level to swell and cause coastal inundation, storm surge, and flooding in the inland rivers where the impacts are multi-fold (ramkumar, 2009). a spatio-temporal pattern of earthquake occurrences in south india was epitomised by gangai et al., (2009) who revealed that seventy-seven occurrences of earthquakes with magnitude above 3.0 were recorded within and adjoined regions of tamil nadu in the last two centuries. ramkumar (2019) argues that the net effects of tectonic quiescent regions and 66 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 geomorphometric processes induce the earthquake occurrences in tamil nadu. natarajan, et al.,. (1988) studied the landslides of the nilgiri regionand described that the deforestation had resulted into massive soil erosion along the slopes, which increased the seepage and ,thereby, resulted into an increase in pore pressure leading to landslides in many hilly topography of the state. further, forest fire (wildfires) also increases the air pollutants and the instantaneous concentration of particulate matters in the atmosphere during the fires and poses serious health risks (sofowote & dempsey, 2015). geoinformatics aids in identification, demarcation, assessment, and management of hazards. nowadays, geoinformatics plays a key role in extracting potential information from an extensive region or inaccessible location. satellite remote sensing has successfully proven itself as a valuable information generator for various hazards studies (bhanumurthy et al., 2010). almost all of the natural hazards that can be studied are mapped with the help of geoinformatics. many studies proved that the utilisation of geoinformatics increases the value of datasets as highly reliable, accurate, and cost effective. many attempts, including behanzin et al., (2015), roy and blaschke (2015); shankar et al., (2015); udo (2015); udani & mathur (2016); prasad & narayanan (2016); hoque (2017); oluwasegun (2017); chigbu et al., (2018); sharma et al., (2018), help us to understand the hazard scenario of the state. however, the integration of multi-hazards and delineation of accurate zones using geoinformatics is still missing in tamil nadu. even though arguably the hazard hotspot regions covers the coastal and urban regions (shrinarayan, 2015; balasundareshwaran et al., 2019), mhz covering the entire state with reliable datasets is still scarce. hence, the present study attempted to identify the multi-hazard prone zones of tamil nadu and showcase the use of geoinformatics in a systematic hazard assessment. 2. methods the research primarily deals with the use of geoinformatics in mhz for tamil nadu. in the initial phase, the spatial and temporal dimensions of natural hazards of tamil nadu were examined. all the possible natural hazards were considered and ranked according to the occurrence and impacts in the past, namely drought, flood, cyclone, storm surge, earthquake, landslides, and forest fire. 67 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 figure 1. location and study area of tamil nadu the drought datasets used for the study, mainly comprised of three categories (1) meteorological records acquired from ground stations, (2) products derived from the satellite sensors, and (3) ancillary data from various sources. high temporal satellite datasets considering vegetation, temperature, and precipitation were used to compute proxies for drought hazard assessment viz. precipitation condition index (pci), temperature condition index (tci) and vegetation condition index (vci). terra modis surface reflectance mod13q1 (250m) (20002016) used for vci, modis lst of mod11a2 (1km) (2000-2016) was for tci and trmm 3b43 (0.25ο×0.25ο) precipitation estimate (2000-2016) was used for pci. finally, the scaled drought condition index (sdci) a multi-sensor drought index proposed by rhee et al. (2010) was used to calculate composed drought index as given in the equation: vcipcitcisdci 25.05.025.0  (1) the low values of sdci imply serious condition of drought. under a drought process, the sdci is close or equal to 0 and at wet conditions the sdci is close to 1. 68 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 the cyclone hazard was assessed by utilising the geographic information system (gis) through a long-term best tracks cyclone database from indian meteorological department (imd) and joint typhoon waring center (jtwc) of bay of bengal for the period of 127 years from 1891 to 2018. storm surge computed by imd based on 50-year return data was adopted for the preparation of storm surge hazard layer. both the layers were compared with the building material and technology promotion council’s (bmptc) vulnerability atlas of india (2006) for conformity and classification. the flood hazard, compared to other hazards,was very difficult to compute on a wider geographical scale, thus it was extracted from the bhuvan geoportal where tile wise flood vulnerability index was computed based on multiple origin of floods (riverine, coastal, urban, flash flood and cyclonic floods). the spatio-temporal characteristic of seismic activities was studied using the data sources of united states geological survey (usgs), coastal geodetic survey (cgs), amateur seismic centre (asc) and imd for tamil nadu and its surrounding regions from 1807 to 2018. based on these datasets and bmptc vulnerability atlas of india, the earthquake zones were delineated. reported landslides and bmptc atlas were used to classify the landslides zones. forest fire hazard map was prepared from the data extracted from tamil nadu forest department and modis satellite data. all the generated hazard layers (earthquake, landslide, cyclone, storm surge, flood, drought, and forest fire) were reclassified into four levels of hazard,namely high, moderate, low, and nonhazard with a cell size of 250 m. by utilising raster calculator tool in arcgis, all the hazard layers were intersected on a normalised scale and multi hazard zones were prepared. based on the combined intensity values, the multi hazard zones were classified into four classes, namely very high, high, moderate, and low. 3. results and discussion 3.1 earthquake the spatio-temporal characteristic of seismic activities in tamil nadu shows concentration around coimbatore, madurai, hosur, salem, shevaroy, villupuram, pondicherry, chennai and dharmapuri (figure 2). the magnitude of the past earthquakes in the study area ranged from 3 to 6. it was observed that the earthquakes of maximum intensity 5-6 had occurred in the coimbatore and surrounding districts. also,it was seen that the areas of maximum seismicity were also the ones that are densely populated. whereas, in spite of the sparse distribution of seismic locations in central and western parts, the frequency tendedto be greater along with magnitude in the western parts in comparison to the central part of the region. 3.2 drought the scaled drought condition index (sdci) is categorised into extreme drought (<0.2), severe drought (0.2-0.3), moderate drought (0.3-0.4), abnormal dry (0.4-0.5), and no drought (>0.5). during the dry year (2016), about half of the state was witnessed the extreme drought. severe drought (0.2 to 0.3) condition was observed in about 30 percent of the state. in 2016, the state was badly hit by the worst annual rainfall in 140 years as it received just 543 mm of rain against the annual average of 920 mm. a weak la-nina over the equatorial pacific that followed a year of 69 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 strong el-nino was cited as one of the reasons for this drought. by considering drought condition during june to december in dry year (2016) and normalcy in base year (2007), the sdci drought index was prepared as shown in figure 3. the interior northern district of tamil nadu, such as karur, nammakal, salem, dindigul and erode faced acute drought. figure 2. earthquake zones in tamil nadu figure 3. scaled drought condition index for tamil nadu 70 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 z 3.3 cyclone and storm surge the cyclone affecting tamil nadu are largely from bay of bengal. the very high hazard appeared in the northern coastal region of tamil nadu. the moderate zone was observed in central part of tamil nadu, extending as tongue from cauvery delta. the november cyclones generally move west or north-westward and strike northern coasts of tamil nadu (siddiki et al., 2012). the storm surge is often accompanied with the cyclone or sometimes with high wind, which affects the coastal margins, that vary up to 5 km from the coast (figure 5). the hazard map shows a very high figure 4.cyclone zones in tamil nadu figure 5. storm surge in coastal parts of tamil nadu 71 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 storm surge (> 4.5 m) in ramanathapuram district’s coast. a high storm surge (3.5 4.5 m) was observed around pulicat lake (northern tamil nadu) which is understood for its tidal activities. 3.4 landslide the western and eastern ghats of the state are prone to landslides. the severe up to very high class was observed in the steep slopes of the western ghats and very confined to the nilgiris, anamalai, kodaikanal, courtallam and agasthyamalai hills. eastern ghat sections, especially shayadri hills and kolli hills,fell under high to moderate landslide hazard (figure 6). 3.5 flood the flood hazard layer shows a combined effect of all flood events of tamil nadu that devastated settlements as well as vegetation. flood vulnerability index for the whole tamil nadu is categorised into five classes, namely very high, high, moderate, low, and very low hazard zones (figure 7). the flood vulnerability zone followed the deltaic regions of majorrivers of tamil nadu, especially between palar to cauvery deltaic region. the upper reaches of tributaries of cauvery and tamirabharani rivers were also found with very high hazard zones. 3.6 forest fire the forest fire is highly concentrated in the western and eastern ghats of the state where forest cover is intersected with settlements and transport lines (figure 8). the forest fire occurred between 2006 and 2015 are classified into five classes of frequency; very low (<4), low (5-9), moderate (10-14), high (15-19), and very high (>20). the very high forest fire occurrences were at mudumalai and gundri in the east kalvarayan, jawathu, kanamangalam, and villapakkam in the west and northern parts of tamil nadu. a cluster of forest fires was observed at valpari, anamalai, palani, kodaikanal, and agamalai hills at eastern central margins of tamil nadu and southern dense vegetation in agasthyamalai hills. 3.7 multi-hazard zonation multi-hazard zones reflect the complex nature of interaction between the hazards. it is very difficult to quantify the interaction but it can be relatively compared (kappes et al., 2012). the natural hazards occur in the same location at different interval cause more damage than the location affected by a single hazard. the multi-hazard zones (mhz) sometimes are found with a domino effect, where one hazard triggers the next possible hazard, which is extremely dangerous. hence, the assessment of mhz is indispensable for disaster risk management and governance of critical time. based on the interaction of all major natural hazards, mhzare delimited and classified into four classes, namely very high, high, moderate, and low (figure 9). the very high mhz is concentrated at hilly stretches of the nilgiris and parts of shayadhri hills in coimbatore and theni districts. further, the very high mhz is found along the coastal stretch of chennai, kancheepuram, cuddalore, and ramanathapuram districts. 72 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 figure 6. landslides in western and eastern and ghats of tamil nadu figure 7. tiled flood index figure 8. forest fire hazard in tamil nadu 73 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 figure 9. multi-hazard zonation (mhz) for tamil nadu 74 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 the high mhz is observed in the central regions of the state engulfing rest of the nilgiris, villupuram, thiruvanamalai, vellore and karur districts, and all along the central and northern coastal tracts of tamil nadu. the occurrences, intensities and frequencies of hazards, namely seismic, landslide and forest fire are the reasons for very high multi-hazard in the hilly tracts, whereas cyclone, storm surge and flood is the major reason for a very high risk along the coastal regions. the moderate mhz vastly spreads over the rest of the state except for a few parts in the northern, eastern, and southern region where hazards with high intensities are absent. the moderate mhz is the direct cause of higher intense of drought but moderate or low intense of cyclone and earthquake hazards. the very high-to-high mhz occurs in the western regions of hilly terrains, central plains, and the eastern margins of coastal regions. the result of this study coincides with case studies attempted for different parts of the state. for instance, agricultural drought assessment in palar basin (krishna et al., 2009) and drought condition over the cauvery basin (balasubramani, 2014) indicate similar prevalence of high mhz at central and northeastern part of the state. the southern part of the nilgirs falls under high and very high vulnerability as the area is under steep slope (rahaman & aruchamy, 2017). moreover, the fire risk assessment at kurangani of theni district forest debated by balaguru et al., (2018) propose that the highly vulnerable zones falls at steep slopes and high altitude regions. comparatively, a study on macrozonation of seismic hazard pursued by menon et al., (2009) exhibits moderate vulnerability in northern and western parts of tn and much lower in the southern half of the state. the study also suggest that the hazard in the northwestern parts, some central districts of the state are clearly underestimated i.e., the hazard levels classified in seismic indian standard zone iii could be marginally higher, especially for chennai, coimbatore, udhagamandalam, and krishnagiri. these discussions support the very highto-high mhz of western regions of hilly terrains and central plains of state. similarly, the study on habitat risk assessment for coastal taluks of tamil nadu (balasundareshwaran et al., 2019) reflects on the assessment of coastal vulnerability utilising the geospatial technology. the study implies the need of urgent measures to thwart negative impacts at chennai, cuddalore, nagapattinam, thoothukudi, and agatheeswaram areas, which are highly vulnerable zones. in the same way, the disaster risk management plan for karur district attempted by balasundareshwaran et al., (2019a) reveals that, during northeast monsoon, heavy water flowing in cauvery and its tributaries, has affected the low-lying areas of the riverbanks with floods. the results of this present study show all the flood occurrences attributed to larger basins of cauvery 75 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 and palar rivers. in an integrated study, classifying coastal tamil nadu (rajan et al. 2019) reveals that the very high vulnerability zone covers nagapattinam coast. further, very high to high vulnerable zones cascade over cuddalore, kancheepuram pondicherry, karaikkal, nagapattinam, and thiruvarur. the physical vulnerability analysis conducted by rajan and vengadasalam (2019) along the coastal tamil nadu observed that coastal zone is still in underdeveloped condition. the study proposed that the ramanathapuram district appears to be a better coast and kancheepuram as a poor. consequently, these deliberations correspond and represent the very high-to-high mhz senario of coastal regions. the major research gap exists in developing appropriate approaches to evolve disaster plans for the administrative units like districts. the hazards demarcated at the state level will be effective for planning at district and sub-district level planning and, thus, this study can be used as a pioneer material. the best use of geoinformatics in the digital era, especially in disaster scenarios, is very much needed, and the present work attempts to fulfil this gap. 4. conclusion this study is akin to a first-hand atlas to the planners and decision makers. the study combined different hazards and their intensities in a single map with varying intensities using the power of geoinformatics. these research highlights need for immediate and proper action plans at high multi-hazard prone zones. the very high mhz is observed at the nilgiris and parts of shayadhri hills in coimbatore and theni districts in the western region and along the coastal stretch of chennai, kancheepuram, cuddalore, and ramanathapuram districts in eastern region of tamil nadu. the following recommendations are suggested for the very high -high mhz of natural disasters of tamil nadu: (1) identification and maintenance of the location specific disaster trigger based response system need to be given a higher priority, especially in the eastern hilly terrains and central plains; (2) preparing region specific disaster management plans along with training for professionals / technical institutions to minimise the ill-effects for both administrative and natural regions, especially for the eastern coastal margins. this study offers scope for immediate knowledge, effective regionalisation, and possibilities of preparing disaster management plans. the segregation of results into administrative division’s categorised in an order of highly risk zones will provide a power tool to the state authorities to allocate fund and resources. the further improvement of accuracy of certain datasets, e.g. flood, could precise the zonation and provide robust tool to prepare long-term plans for safeguarding the vulnerable communities, mitigating the physio-socio-economic losses, and attaining sustainable development 76 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 conflict of interest the authors declare that there is no conflict of interest with any financial, personal, or other relationships with other people or organizations related to the material discussed in the article. acknowledgement the authors thank the council of scientific and industrial research university grant commission (csir-ugc) and indian council of social science research (icssr) for the financial support to complete the study. the authors also thank ugc-sap-drs-ii programme of the department of geography, bharathidasan university, tiruchirappalli for the infrastructural and technical supports. references balaguru, m., navammuniyammal, m., vidhya, r. and sathyavathi, g., (2018). assessing forest fire prone area in kurangani, tamil nadu using remote sensing and gis. international journal of applied engineering research, 13(7), pp. 264-269 balasubramani, k., (2014). “assessing vegetative drought from satellite data: a case study of cauvery basin, south india”. the konkan geographer, 8, pp.12-18. balasundareshwaran, a., abdul rahaman, s., balasubramani, k., kumaraswamy, k. and ramkumar, m., (2019). “habitat risk assessment along coastal tamil nadu”, india an integrated methodology for mitigating coastal hazards”. coastal zone management, pp. 515-542. doi.org/10.1016/b978-0-12-814350-6.00023-9 balasundareshwaran, a., kumaraswamy, k., balasubramani, k. and abdul rahaman, s. (2019a). “geoinformatics catalyst for preparation of district disaster risk management plan: a case of karur district, tamil nadu”, indian geographical journal, 94 (1), pp. 1-14. behanzin, i.d., thiel, m., szarzynski, j., and boko, m., (2015). “gis-based mapping of flood vulnerability and risk in the benin niger river valley”. international journal of geomatics and geosciences, 6 (3), pp. 16531668. bhanumurthy, v., manjusree, p., and srinivasa rao, g., (2010). “flood disaster management”. remote sensing applications, pp. 283-296. chigbu, n., mbah, e. o., and anyanwu, i., (2018). “assessment of flood risk of cross river state using geographic information system (gis)”. embracing our smart world where the continents connect: enhancing the geospatial maturity of societies, istanbul. gangai, p., selvam, k. and kumaraswamy, k., (2009). “a spatio-temporal patterns of earthquake occurrences in south india”. indian national geographer, 24 (1), pp.48-54. hoque, m.a., (2017). “cyclone disaster mapping, monitoring and management using satellite remote sensing and spatial analysis”. (doctoral dissertation) https://espace.library.uq.edu.au kappes, m. s., keiler, m., von elverfeldt, k., glade, t., (2012). “challenges of analysing multihazard risk: a review”. natural hazards, 64 (2), pp. 1925-1958. 77 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 krishna, t.m., ravikumar, g. and krishnaveni, m., (2009). “remote sensing based agricultural drought assessment in palar basin of tamil nadu state, india”. journal of the indian society of remote sensing, 37, pp.9–20. doi.org/10.1007/s12524-009-0008-8 menon, a., ornthammarath, t., corigliano, m. and lai, c.g. (2009). probabilistic seismic hazard macrozonation of tamil nadu in southern india. bulletin of the seismological society of america, 100(3), pp.1320–1341. doi.org/10.1785/0120090071 natarajan, t.k., murty, a.v.s.r., and chandra, d., (1988). “control of surfacial slides by different erosion control techniques”. proceedings 2”. international conference in case histories in geotechnical engineering, rolla, 1, pp. 361-364. ndma, (2016). national disaster management plan, goi, report. oluwasegun, a. h., (2017). “gis analysis of flood vulnerable area in benin-owena river basin nigeria”, indonesian journal of geography, 49 (1), pp. 27-33. doi.org/10.22146/ijg.12777 prasad, n.n.r., and narayanan, p., (2016). “vulnerability assessment of flood-affected locations of bangalore by using multi-criteria evaluation”, annals of gis, 22 (2), pp. 151-162, doi.org/10.1080/19475683.2016.1144649 rahaman, s.a. and aruchamy s., (2017). geoinformatics based landslide vulnerable zonation mapping using analytical hierarchy process (ahp), a study of kallar river sub watershed, kallar watershed, bhavani basin, tamil nadu. modeling earth systems and environment, 3 (1) pp.1-13. doi.org/10.1007/s40808-017-0298-8 rajan s.m.p., and vengadasalam r., (2019): swot and pci based grading of coastal zone: a model emulation from tamil nadu coast, india. human and ecological risk assessment: an international journal. doi.org/10.1080/10807039.2019.1678382 rajan s.m.p., nellayaputhenpeedika m., tiwari s.p. and vengadasalam r., (2019): mapping and analysis of the physical vulnerability of coastal tamil nadu. human and ecological risk assessment: an international journal. doi.org/10.1080/10807039.2019. 1602752 ramkumar, m., santosh, m., abdul rahaman, s., balasundareshwaran, a., balasubramani, k., mathew, m. j., sautter, b., siddiqui, n., usha, k. p., sreerhishya, k., prithiviraj, g., nagarajan, r., thirukumaran, v., menier, d. and kumaraswamy, k., (2019). “tectonomorphological evolution of the cauvery, vaigai, and thamirabarani river basins: implications on timing, stratigraphic markers, relative roles of intrinsic and extrinsic factors, and transience of southern indian landscape”, geological journal, pp.1-42. (issn: 0072-1050) doi.org/10.1002/gj.3520 ramkumar, mu., (2009). “types, causes and strategies for mitigation of geological hazards”, in: ramkumar mu., (ed.) geological hazards: causes, consequences & methods of containment. new india publishing agency, pp. 1-22. roy, d.c. and blaschke, t., (2015). “spatial vulnerability assessment of floods in the coastal regions of bangladesh”, geomatics, natural hazards and risk, 6 (1), pp. 21-44. doi.org/10.1080 /19475705.2013.816785 shankar, v.s., dharanirajan, k., and narsimulu, g., (2015). “cyclone vulnerability zonation of southern part of south andaman, india using multi-criteria weighted overlay techniques”, indian journal of geo-marine sciences, 44 (8), pp. 11811190. 78 a. balasundareshwaran et al / geosi vol 5 no 1 (2020) 65-79 sharma, s.v.s., roy, p. s., chakravarthi, v. and rao, g., (2018). “flood risk assessment using multi-criteria analysis: a case study from kopili river basin, assam, india”. geomatics, natural hazards and risk, 9 (1), pp. 79-93. doi.org/10.1080/19475705.2017.1408705 shrinarayan, j. n., (2015). “study of multi-hazards on urban habitats of india: a case study of ahmedabad city”, thesis submitted to gujarat technological university, ahmedabad. siddiki, u. r., islam, m. n., and aziz, m. n., (2012). “cyclonic track analysis using gis over the bay of bengal”, international journal of applied sciences and engineering research, 1(5), pp. 689-701. doi: 10.6088/ijaser.0020101070. sofowote, u., and dempsey, f., (2015). “impacts of forest fires on ambient near–real–time pm2 in ontario, canada: meteorological analyses and source apportionment of the july 2011–2013 episodes. atmospheric pollution research, 6(1), 1-10. udani, p.m., and mathur, d.k., (2016). “flood hazard vulnerability mapping using remote sensing and gis: a case study along the villages of anand taluka”, advances in applied science research, 7 (3), pp. 214-221. udo, e.a., ojinnaka, o.c., baywood, c.n., and gift. u.a., (2015). “flood hazard analysis and damage assessment of 2012 flood in anambra state using gis and remote sensing approach”, american journal of geographic information system, 4 (1), pp. 38-51. undp (2012). “disaster risk reduction programme”, united national development programme. (date accessed: 28/07/2018) https://www.in.undp.org/content/india/en/home/operations/projects/closed/goiundp_ disaster riskreductionproject.html 79 multi-hazard zonation for effective management of disasters in tamil nadu 1. introduction 2. methods 3. results and discussion 3.1 earthquake 3.2 drought 3.3 cyclone and storm surge 3.4 landslide 3.5 flood 3.6 forest fire 3.7 multi-hazard zonation 4. conclusion references 77 research article development of web-based gis alert system for informing environmental risk of dengue infections in major cities of pakistan naureen zainab 1, aqil tariq 2,* , saima siddiqui 3 1department of computer software engineering, military college of signals, national university of science and technology (nust), islamabad, 44000, pakistan 2state key laboratory of information engineering in surveying mapping and remote sensing (liesmars), wuhan university, wuhan, 430079, china 3department of geography, university of punjab, lahore, 54590, pakistan received 23 november 2020/revised 10 march 2021/accepted 19 march 2021/published 25 april 2021 abstract dengue is one of the emerging major public health problems, and its incidence varies with climate conditions. it affects millions of people's lives owing to unusual socioeconomic conditions and epidemiological factors. this study was designed to build a web-based gis alert system for dengue data management and analysis which would centralize information and make it accessible to all relevant stakeholders before, during, and after crises. three geographical regions were selected in this study. the user interface of the dengue alert system was developed based upon mapguide. results indicate that risk level was mainly associated with breteau index. karachi and lahore were at their highest risk, i.e., level 4. islamabad and chakwal were also at the highest risk, i.e., level 4. attock had high risk, i.e., level 3 followed by haripur with minimal level 1. the high breteau index showed a direct relationship to high potential transmission of dengue outbreaks, a more significant peak of dengue was the result of monsoons, while smaller peaks were observed due to domestic water storage. hence, it was concluded that monsoon is the best suitable season for the development of dengue. web-based gis alert system for dengue data management and analysis was developed, centralizing information and making it accessible to all relevant stakeholders before, during & after a crisis. this program creation will provide a more analytical forum for advising multiple levels of risk and an experimental method for measuring the effect of different factors on risk level distribution by adjusting the component's weighting. keywords : dengue; gis analysis; gui; alert system; breteau index; weighted overlay 1. introduction dengue is a mosquito-borne virus whose prevalence differs with temperature and weather (chang et al., 2009; gubler, 2006). dengue is rising as one of pakistan's most significant public health problems (asif et al., 2013). in october 2005, after ten years, dengue affected karachi again, and 21 deaths were reported out of 103 confirmed cases (mukhtar et al., 2011). geosfera indonesia *corresponding author. email address : aqiltariq@whu.edu.cn (aqil tariq) vol. 6 no. 1, april 2021, 77-95 p-issn 2598-9723, e-issn 2614-8528 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v6i1.20792 https://orcid.org/0000-0003-1196-1248 https://orcid.org/0000-0003-3020-0233 mailto:aqiltariq@whu.edu.cn 78 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 the year 2011 was pakistan's worst year for dengue. in 2013 the most affected area was the dengue fever (df) outbreak in pakistan. aedes-aegypti and aedes-albopictus is considered to be the main dengue vectors in south asia, including pakistan. the identification of disease outbreaks is very significant (sirisena et al., 2017). epidemics are silent in contrast to explosions. outbreaks kill or cause sickness before detection (kahn et al., 1975). disease outbreaks can easily cause such hurt, and they can spsread rapidly, too. in the worst scenario, the window of opportunity to minimize this harm could be limited to a few days (thompson et al., 2016). the united states of america spends billions of dollars a year on different forms of safety surveillance. the essential expenses include patient infection control, public health surveillance (phs), air and water inspection, preparation, improved public health, and science services in information technology (kahn et al., 1975). proper and prompt treatment of incidents of sickness can save many meaningful humans lives. the breteau index calculates the number of positive containers per hundred surveyed homes, which in turn represents the distribution of aedine mosquitoes, the dengue vector ( udayanga et al., 2018). an alert system was established in pakistan from a synthesis of geospatial data on the breteau map. the interrelationship was rendered between the breteau index and the temperature. it forms the basis for creating weighted overlays to determine risk levels (attaway et al., 2016). dengue alert system generally predicts four different risk levels for the risk of dengue infections, i.e., highest, high, medium, and minimal ( olubadewo-joshua & ugom, 2019; tran et al., 2020). hence, a web-based application was created from a synthesis of geospatial data related to breteau index and temperature effect in pakistan's major cities. using various weighting factors on the previous history, appropriate weighting factor was evaluated and used to generate the alert for informing environmental risk of dengue infection in major cities of pakistan. it forms the basis for creating weighted overlays to determine levels of risk. the weighting can be adjusted to adjust alarm device sensitivity (bowman et al., 2014). every year pakistan experiences a lot of rain during the mild summer (monsoon) season. the daytime temperature is greater than the night temperature. its good climatic conditions for the growth of dengue mosquitos. in this study, we target six major cities of pakistan because pakistan faces a lot of problems in these regions every year. in this study, we clearly describe the spatio-temporal distributions of dengue cases. the study is novel because its first time correlates the local predominant conditions with dengue prevalence and builds a web-based gis alert system. therefore, this study was designed to build a web-based gis alert system for dengue data management and analysis which would centralize 79 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 information and make it accessible to all relevant stakeholders before, during, and after crises. 2. methods 2.1. study area the study area includes rawalpindi, haripur, attock, chakwal, lahore, and karachi (figure 1). all metropolitan regions are located in different latitudes and longitude in pakistan. these cities' climate is characterized by four seasons: dry, wet, hot, and cold (abbas, 2013). rawalpindi and lahore experience a monsoonal environment with rainy, hot summers and cool dry winters (wet and dry season); rainfall is typical of pakistan's semi-arid region. karachi also receives monsoon rainfall, but the landscape is different compared to rawalpindi and pakistan (burney et al., 2018). karachi harbor is a secure and majestic natural harbor on the shores of which the town is located. a low-lying coastal area stretches along the harbor's edge. the malir river, a seasonal current, flows through the eastern part of the city, and the seasonal layari river runs through the northern section that is most densely populated. several ridges and small hills occur; the maximum elevation, mango pir, is 585 feet high. figure 1. study area of dengue epidemic-prone districts in pakistan 80 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 2.2. system design, development, and assumption the alert system had been developed to forecast the risk of infection with dengue. it was presumed that the two variables that were used in the method, namely temperature and breteau scale, were strongly associated with dengue fever infection sensitivity and are of significant geographical significance (kaya et al., 2019). 2.3 data source data on breteau index and affected people from 2006 to 2013 were collected from the national institute of health and governments health departments. the monthly temperature data from 2006 to 2013 was acquired from the pakistan meteorological department of islamabad. 2.4 weighted overlay a web-based application was developed from a synthesis of geospatial data related to the breteau index and temperature (cetin et al., 2019).using various weighting factors on the previous history, appropriate weighting factor was evaluated and was used to generate the alert for informing environmental risk of dengue infection in major cities of pakistan. the breteau index was calculated according to eq. (1) as described by gubler et al. ( 2014), while the risk level was calculated based on the formula generated as eq. (2). breteau index (bi) = number of positive containers houses inspected × 100 (1) risk level = bi (60 + temperature) 40 (2) when the breteau index value is zero, the risk level will be zero too, even in the case of high temperature. for instance, in similar temperatures in two different cities, two different risk levels can be observed due to different breteau indexes. 2.5 risk levels dengue risk levels were classified (table 1) into four different groups based on the developed formula of the risk level. the breteau index against these four risk levels is also analyzed. it was assumed that if breteau index is high, then the risk level is also high as breteau index and risk levels are directly proportional to each other (table 2). 81 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 table 1. risk levels (gubler et al., 2014) classification risk level level 1 01200 table 2. breteau index according to risk levels (gubler et al., 2014) classification risk level breteau index level 1 01200 bi >=20 2.6 prevention or management of dengue borne vector for confined risk level the adoption of integrated pest management did the management of dengue vectors at confined four risk levels. integrated pest management is the pest population's level to check the pest below the economic injury level using all possible control measures (novotny et al., 2007). recommendations regarding each risk level were given on the basis of ipsm strategies. there are specific preventive and control measures recommended for each alert level shown in table 3. table 3. preventive measures against each risk level (novotny et al., 2007) classification breteau index risk level actions level 1 0< bi<=4 0=20 risk level>1200 launch a task-oriented force to create awareness among people regarding dengue elimination programmed and give them responsibilities to restrain the potential breeding places by allotting them specific territories. private competitive pest control and management contractor must be employed to control and eradicate the arising potential mosquito problem. 82 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 2.7 conceptual design the present study's concept was to build a web-based gis alert system both online and offline for dengue data management and analysis, centralizing the information and making it accessible to all relevant stakeholders before, during, and after a crisis. the system aimed to ensure timely availability of information on health care services. the project data were mapped accurately and with geographic features displayed. maps were drawn according to pre-set parameters in the present framework, and then the web browser showed the map in a .jpeg format. the users can create and display a new map by adjusting the parameters. on the server-side, this sort creates a heavy load. the server-side had regional info, gis tools, and a system with hundreds of complex reports on the interface (cetin, 2015). the systematic conceptual diagram of the web-based gis alert system is shown in figure 2. figure 2. conceptual design of web-based gis alert system. 2.8 user interface (ui) mapguide maestro is a mapguide open-source map authoring tool. mapguide maestro strives to support the capabilities of the open-source application mapguide. both spatial and non-spatial details have been contained in sql server express, which is a free database service that fits well with any web application platform (zavlavsky, 2000). the graphical user interphase (mapguide) was systematically developed, as shown in figure 3. 83 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 figure 3. graphical user interphase (mapguide) 2.9 component of web-based gis alert system three components of the web-based gis alert system (figure 4) were used as: (1) dengue information databases (2) gis layers and maps (3) queries and reports figure 4. components of web-based gis alert system overall, the system describes the risk levels of dengue infection in major cities of pakistan, classifying them into three geographical regions. a web-based application was created from a synthesis of geospatial data related to the breteau index and temperature effect in pakistan's major cities. currently, statistics on the breteau index and temperature are 84 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 displayed in a table format that is difficult for interpretation by general users. this project allows the visualization of the breteau index and temperature information in a spatial pattern and association with four risk levels, which enhances public awareness. dengue alert system describes the risk level of dengue infections from 2006 to 2013 in major cities of pakistan. it can also predict the current risk level of dengue in any district of pakistan. it can predict four different risk levels for the risk of dengue infections, i.e., highest, high, medium, and minimal. 3. results and discussion 3.1 temporal changes in temperature and breteaux index a good association between the two variables has been identified from the temporal change in temperature and the breteau scale (bozdogan et al., 2021). higher breteau index focused primarily in the summer and autumn when higher temperatures prevailed. results shown in figure 5a indicate the temporal changes in breteau index in september 2010 in selected study sites. the breteau index of lahore and islamabad represented by green color was high. the breteau index of karachi denoted by blue color was medium. the breteau index of chakwal and haripur represented by yellow color was minimal, and the breteau index of attock represented by gray showed no risk. similarly, the temporal changes in temperature in major cities of pakistan figure 5b indicate that karachi had high temperature while the remaining cities showed low temperature. the high temperature was represented by dark gray color, and low temperature was represented by light gray. the effect of temperature on the weighted overlay product can be calculated when the temperature and the breteaux index are overlaid. the temperature was not considered to be the main factor deciding the level of risk according to the given weighting, and when the breteau scale equals zero, the risk level was also zero, irrespective of the temperature increase. the cities with a high breteau index are at a higher risk level. lahore and islamabad's temperature was low compared to karachi, but they had a high-risk level due to the high breteau index. breteau index of attock was zero, so there was no risk; however, the attock temperature was 28oc, similar to chakwal 27.6oc. 85 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 figure 5a. breteau index of september 2010 in major cities of pakistan, figure 4b. temperature of september 2010 in major cities of pakistan, figure 4c. breteau index of september 2011 in major cities of pakistan, figure 4d. temperature of september 2011 in major cities of pakistan. figure 5c and figure 5d show the breteau index and temperature of september 2011 in selected study sites. attock and haripur very high followed the breteau index of karachi, lahore, islamabad, rawalpindi, and chakwal. the temperature of karachi and lahore was high as compared to other cities. figure 5d represents that karachi and lahore were at their highest risk and the temperature of these two cities was also high. rawalpindi, islamabad, and chakwal were also at the highest risk, but they had a low temperature as compared to karachi and lahore. attock had high risk, and the temperature was low. haripur showed minimal risk and had a low temperature. rawalpindi and attock's temperature showed a similar trend but had different risk levels due to different breteau index. web-based generated map of the risk levels of dengueshows that risk level was high in those cities which had high breteau index. results indicate that karachi and lahore were at their highest risk, i.e., level 4. islamabad and chakwal were also at the highest risk, 86 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 i.e., level 4. attock had high risk, i.e., level 3 followed by haripur with minimal level 1. hence, it is concluded that risk level was mainly associated with breteau index, if breteau index was high, then the risk level was also high, so this is the main reason for which the highest weightage is assigned to the breteau index (gungor et al., 2020). temperature also plays a significant role because the risk level was high in mostly the summer months, which had a high temperature. if the temperature was high, mosquito's growth rates were also high because high temperature helps in breeding aedes aegypti and aedes albopictus mosquitoes. 3.2 peaks of dengue based on the analysis of risk levels from 20062013 in major cities of pakistan results from the analysis of dengue risk levels from 2006 to 2013 in major cities of pakistan indicated that the high risk of dengue from 2006 to 2013 was mostly observed in the month of september, october, and november (table 4). the dengue-related deaths were reported between the months from september to november. dengue case-load registered during these months, and the reports declined rapidly during and after december. the high breteau index showed a major direct relationship to the high potential transmission of dengue outbreaks in the study sites. two peaks of dengue were shown in the form of risk levels. the formation of a more significant peak of dengue resulted from monsoon rains in which water was filled at the low-lying places. smaller peaks of dengue were formed in the dry season due to domestic storage of water. it was seen that transmission of dengue during the monsoon and after monsoon was maximum than that of the dry season. finally, it was concluded that monsoon is the best suitable season for the reproduction, fecundity, survival, growth, and development of dengue. dengue fever is estimated to affect 50 to 100 million people with 1/2 million lifethreatening infections globally in a year as of 2010. it has risen 30-fold significantly in frequency between 1960 and 2010. in pakistan in the year 2010, there were more than 7000 positive confirmed cases of dengue.with rising epidemics, dengue fever has become a significant disease in pakistan. given pakistan's government's efforts, especially in punjab, the high treatment cost has hindered pakistan's ability to control epidemics. dengue-fever mortality in pakistan in the summer of 2011 was over 300 people, and epidemic incidence was over 14,000 infections. the outbreaks mainly occurred in the region of lahore, punjab, pakistan. 87 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 table 4. results from the analysis of dengue risk levels from 2006 to 2013 in major cities of pakistan city dengue risk 2006 2007 2008 2009 2010 2011 2012 2013 k a r a c h i duration jan apr jul aug jan, sep to dec aug to dec jan apr aug to dec jan apr to dec jan to dec all over the year minimum jan apr jul aug jul dec jan jan aug sep dec jan apr aug dec jan apr to jul jan to sep jan to aug dec moderate sep sep sep dec sep aug nov dec dec sep highest oct nov oct nov oct nov oct nov oct nov sep oct oct nov oct nov l a h o r e duration jan aug to dec sep dec jul dec apr dec apr aug to dec mar to dec minimum jul aug sep dec dec jan, aug dec dec jul aug dec apr jul apr aug dec mar to july sep moderate oct sep oct sep sep oct sep aug nov dec sep aug dec highest nov nov oct nov nov oct nov sep oct oct nov oct nov a tt o c k duration oct to dec oct to dec oct nov july sep to dec minimum sep sep dec dec jul dec nov moderate oct nov oct nov oct oct oct sep nov highest nov nov nov oct continued 88 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 city dengue risk 2006 2007 2008 2009 2010 2011 2012 2013 r a w a lp in d i duration sep to dec sep to dec sep dec jun to dec jun to dec minimum jul aug aug dec dec sep dec jun jul aug jul aug dec moderate oct sep oct sep sep oct nov dec sep highest nov nov oct nov nov oct nov sep oct oct nov is la m a b a d duration aug to dec sep to dec sep to nov sep to dec jun to dec aug to dec minimum dec sep dec jun jul aug dec aug to nov aug sep dec moderate sep oct oct sep sep oct nov oct highest nov nov oct nov oct nov nov sep oct nov h a r ip u r duration sep to nov oct nov sep to nov apr jul to dec minimum dec dec dec sep apr jul aug sep dec sep oct sep oct nov moderate oct sep oct oct highest nov oct, nov nov oct nov nov c h a k w a l duration no case jan nov sep to nov sep to nov minimum oct, nov sep sep nov moderate oct, nov jan oct oct nov highest nov nov sep note : = no case 89 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 3.3 application web-based gis alert system for web-based gis alert system application, a user can select url localhost/dengue. the dengue support system will display. for login purposes, the user can click on the login option. before getting detailed information, the customer will be asked to enter the id and password. if the user identification and password are wrongly entered, an error message will pop up, "your attempt to login is not successful. please try again," and the user will be asked to register until the correct information is typed in. this system would help to protect data security under study. the existing database had included eight years' data. the application was developed monthly, so a first select year and then the month, and then click/check the risk levels. database contents and collections could be expanded in the future. after selecting the month and the year, the selected year's risk level will be displayed. users can also zoom the areas concerned by clicking on the option zoom rectangle. this feature provides the consumer a general understanding of the environment impacted and theoretically affected. one of the most important functions of this system was to find the current risk level of any district based upon the breteau index and temperature. by entering the breteau index and temperature, the system can find the current risk level of dengue of any pakistan district by selecting the district in the drop-down list. results shows that the user selected district faisalabad entered breteau index 7.53 and temperature 21.4 and then clicked on select, the risk level of faisalabad was displayed. the blue color represents that it was at level 2. the dengue alert system was capable of generating the reports of risk levels of each month. the reports of the risk level of dengue in september 2011 along with "actions to be taken" at each level. the system was also capable of generating the graphical representation of the risk level of each month (figure 6). figure 6. graphical representation of risk level of september 2011 1741.9 254.15 2044.5 1456.18 1457.7 1342.7 0 500 1000 1500 2000 2500 chakwal haripur islamabad karachi lahore rawalpindi series1risk level 90 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 it has been seen from the data and observed that high dengue indicators are in august, september, october, and november, i.e., after july, which have a high temperature in pakistan. if the breteau scale is 0, the risk level will also be zero, regardless of how high the temperature is. it is concluded that risk level is mainly associated with breteau index, if breteau index is high, then risk level is also high, so this is the main reason due to which highest weightage is assigned to breteau index, but temperature also plays a significant role because risk level is high in mostly the months which have a high temperature. if the temperature is high, mosquito growth chances are also high because high temperature helps in the breeding of aedine mosquitoes. there is no chance even if the temperature is very high without the existence of aedine mosquitoes. for instance, if the temperature is similar for two cities, they have very different risk levels. the device may be combined with other potential factors. human population growth is a significant dengue consideration because the desert region containing many mosquitoes is not dangerous. if the weather around the house is suitable for dengue, then dengue infections will be high. the moisture and containers may be important influences. eventually, mosquito species often play a part. the system is not without limitations. it is a web-based application, and for this application, internet connectivity must require. if the internet speed is slow, this application will not work properly because it requires high bandwidth. according to dengue outbreaks, the likelihood of dengue is caused by many variables in real life including the virus and the climate. the application of temperature and breteau index offers one element of risk assessment, concentrating on exposure impact calculation alone. certain considerations might be applied to the system's robustness. similar results/an online alert system (gis) were developed. the interrelationship between various indices and temperature was measured by wong et al., (2007), who observed that the inference forms the rationale for the generation of weighted overlays to define risk levels. weighting can be adjusted to set the sensitivity of the alert system. they concluded that the alert system offers one objective means to defining the risk of dengue in a society, which would not be affected by the incidence of the infection itself similarly, shen et al. (2015) conducted a study to explore the associations between the monthly number of dengue fever (df) cases and possible risk factors in guangzhou, a subtropical city of china. in their study, a total of 39,697 df cases were detected in guangzhou. dengue fever incidence showed an obvious seasonal pattern, with most cases occurring from june to november. they observed that the current month's breteau index, average temperature, previous month's minimum temperature, and tave were positively 91 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 associated with dengue fever incidence. they also concluded that the mosquito density, tave, and tmin play a critical role in dengue fever transmission in guangzhou. in another study, liyanage et al., (2019) used an interrupted time-series design with a non-linear extension, and evaluated the impact of vector control interventions from june 22, 2014, to dec 29, 2016, in panadura, a high-risk moh division in western province, sri lanka. they used dengue notification and larval survey data to estimate the reduction in breteau index and dengue incidence before and after the intervention using two separate models, adjusting for time-varying confounding variables (i.e., rainfall, temperature, and the oceanic niño index). they found that breteau index is an essential index for controlling and measuring the extent of vector control in dengue transmission. bajwala et al (2020) conducted a study using mean temperature (°c), relative humidity (%), and precipitation (rainfall in mm) as climatic parameters affecting vector ecology. they recorded the peak number of cases was recorded in the post-monsoon period. they proposed that presence of some serologically positive cases even during dry months in this study could probably be reflective of the year-round activity of the dengue vector. moreover, some advanced techniques, i.e., artificial intelligence-based mathematical model and fuzzy logic, could be a suitable option to implement control measures in dengueprone areas. adak & jana (2021) developed an artificial intelligence-based mathematical model taking the three stegomyia indices, namely house index, breteau index, and container index. the project's concept was to build a web-based gis alert system for dengue data management and analysis, centralizing information and making it accessible to all relevant stakeholders before, during & after a crisis situation. the system aims to ensure timely availability of information on health care services and project data, mapped accurately and with geographic features displayed, designed an alert system for users to visualize the spatial distribution of the risk of exposure to aedine mosquito in the local setting. the system complements information on the incidence of infection and is a more objective way of defining risk levels. overall, the system describes the risk levels of dengue infection in major cities of pakistan, classifying them into three geographical regions, and can predict all districts of pakistan's risk levels. a web-based application was created from a synthesis of geospatial data related to breteau index and temperature effect in pakistan's major cities. future alignment of this program with practices of prevention and control of dengue fever will be important. this pilot project may be associated with the various government departments involved. similar systems could be built and adapted to estimate exposure to other diseases, such as bird flu, for which environmental factors are critical for propagation. 92 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 the program could become a tool for local, regional, and even national health authorities and other stakeholders in the long run. 4. conclusion the user interface of the dengue alert system was developed based upon mapguide. the inter-relationship between breteau index and the temperature was developed using a web-based gis application. four different patterns mapping were developed for the risk of dengue, i.e., highest, high, medium, and minimal. maps were created using pre-set criteria, and the map was shown as an image by the web browser.a high risk of dengue was found in august, september, october, and november, i.e., after july, in which had high temperature was noted. risk level was found mainly associated with breteau index; if breteau index was high, then risk level was also high. when the breteau index equals zero, the risk level was zero, irrespective of the temperature rise. temperature also played a significant role because the risk level was high, mainly in the months with high temperatures. there was no risk without aedine mosquitoes' presence even if the temperature was very high. for instance, if the temperature was similar for two cities, they depicted very different risk levels. this study's significance is introducing a new method for expressing the alert system on a spatial scale. the web-based alert system should no longer be expressed in statistical and textual formats, but that spatial graphics could be incorporated. in addition, the weighted overlay allows the users to weigh the contribution of factors affecting the exposure risk to dengue. conflict of interest the authors declare that there is no conflict of interest. acknowledgments authors gratefully acknowledge the support received from dr. mukhtar senior scientific officer (entomology), ex. head of department of zoonotic and vector-borne diseases and thanks to dr jan muhammad director pakistan meteorological department for monthly mean temperature data. we are thankful to dr. m. farrukh sultan, dpc, edoh, lahore. the authors would also like to express their appreciation to rana naveed mustafa (narc), naveed yaqoob, department of mathematics, qau islamabad and azmat ali sr. project manager erp (it/mis) at agribusiness support fund for their kind cooperation. this research was supported by the national natural science foundation of china (no. 41871345). 93 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 references abbas, f. (2013). analysis of a historical (1981-2010) temperature record of the punjab province of pakistan. earth interactions, 17(15), 1–23. https://doi.org/10.1175/2013ei000528.1. adak, s., & jana, s. (2021). a study on stegomyia indices in dengue control: a fuzzy approach. soft computing a fusion of foundations, methodologies and applications, 25(1), 699–709. https://doi.org/10.1007/s00500-020-05179-x. asif, m., tripathi, n. k., & ahmed, s. (2013). towards near real time public health surveillance (a decision support system for public health surveillance). international journal of computer applications, 61(21), 45–50. attaway, d. f., jacobsen, k. h., falconer, a., manca, g., & waters, n. m. (2016). risk analysis for dengue suitability in africa using the arcgis predictive analysis tools (pa tools). acta tropica, 158, 248-257. https://doi.org/10.1016/j.actatropica.2016.02.018. bajwala, v. r., john, d., rajasekar, d., eapen, a., & murhekar, m. v. (2020). burden of dengue with related entomological and climatic characteristics in surat city, gujarat, india, 2011–2016: an analysis of surveillance data. the american journal of tropical medicine and hygiene, 103(1), 142-148. https://doi.org/10.4269/ajtmh.19-0967. bowman, l. r., runge-ranzinger, s., & mccall, p. j. (2014). assessing the relationship between vector indices and dengue transmission: a systematic review of the evidence. plos neglected tropical diseases, 8(5). https://doi.org/10.1371/journal.pntd.0002848. bozdogan sert, e., kaya, e., adiguzel, f., cetin, m., gungor, s., zeren cetin, i., & dinc, y. (2021). effect of the surface temperature of surface materials on thermal comfort: a case study of iskenderun (hatay, turkey). theoretical and applied climatology, 1–11. https://doi.org/10.1007/s00704-021-03524-0. cetin, m. (2019). the effect of urban planning on urban formations determining bioclimatic comfort area’s effect using satellitia imagines on air quality: a case study of bursa city. air quality, atmosphere & health, 12(10), 1237-1249. https://doi.org/10.1007/s11869-019-00742-4. cetin, m., adiguzel, f., gungor, s., kaya, e., & sancar, m. c. (2019). evaluation of thermal climatic region areas in terms of building density in urban management and planning for burdur, turkey. air quality, atmosphere & health, 12(9), 1103-1112. https://doi.org/10.1007/s11869-019-00727-3. cetin, m. (2015). using gis analysis to assess urban green space in terms of accessibility: case study in kutahya. international journal of sustainable development & world ecology, 22(5), 420-424. https://doi.org/10.1080/13504509.2015.1061066. chang, a. y., parrales, m. e., jimenez, j., sobieszczyk, m. e., hammer, s. m., copenhaver, d. j., & kulkarni, r. p. (2009). combining google earth and gis mapping technologies in a dengue surveillance system for developing countries. international journal of https://www.ijcaonline.org/archives/volume61/number21/10216-5083 94 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 health geographics, 8(1), 1–11. https://doi.org/10.1186/1476-072x-8-49. gubler, d. j. (2006). dengue/dengue haemorrhagic fever: history and current status. in novartis foundation symposium (vol. 277, p. 3). chichester; new york; john wiley. gubler, d. j., ooi, e. e., vasudevan, s., & farrar, j. (2014). dengue and dengue hemorrhagic fever. oxfordshire : cabi. gungor, s., cetin, m., & adiguzel, f. (2020). calculation of comfortable thermal conditions for mersin urban city planning in turkey. air quality, atmosphere and health, 1–8. https://doi.org/10.1007/s11869-020-00955-y. kahn, t. c., cameron, j. t., & giffen, m. b. (1975). methods and evaluation in clinical and counseling psychology. elmsford, ny: pergamon press. kaya, e., agca, m., adiguzel, f., & cetin, m. (2019). spatial data analysis with r programming for environment. human and ecological risk assessment, 25(6), 1521– 1530. https://doi.org/10.1080/10807039.2018.1470896. liyanage, p., rocklöv, j., tissera, h., palihawadana, p., wilder-smith, a., & tozan, y. (2019). evaluation of intensified dengue control measures with interrupted time series analysis in the panadura medical officer of health division in sri lanka: a case study and cost-effectiveness analysis. the lancet. planetary health, 3(5), e211–e218. https://doi.org/10.1016/s2542-5196(19)30057-9. mukhtar, m., tahir, z., baloch, t., mansoor, f., & kamran, j. (2011). entomological investigations of dengue vectors in epidemic-prone districts of pakistan during 2006– 2010. world health organization regional publications. south east asia series. novotny, v., miller, s. e., hulcr, j., drew, r. a. i., basset, y., janda, m., … weiblen, g. d. (2007). low beta diversity of herbivorous insects in tropical forests. nature, 448(7154), 692–695. https://doi.org/10.1038/nature06021. olubadewo-joshua, o., & ugom, k. m. (2019). application of geospatial techniques in the locational planning of health care centers in minna, nigeria. geosfera indonesia, 3(3), 59-72. https://doi.org/10.19184/geosi.v3i3.8754. shen, j. c., luo, l., li, l., jing, q. l., ou, c. q., yang, z. c., & chen, x. g. (2015). the impacts of mosquito density and meteorological factors on dengue fever epidemics in guangzhou, china, 2006-2014: a time-series analysis. biomedical and environmental sciences, 28(5), 321–329. https://doi.org/10.3967/bes2015.046. sirisena, p., noordeen, f., kurukulasuriya, h., romesh, t. a., & fernando, l. k. (2017). effect of climatic factors and population density on the distribution of dengue in sri lanka: a gis based evaluation for prediction of outbreaks. plos one, 12(1), e0166806. https://doi.org/10.1371/journal.pone.0166806. burney, s. a., barakzai, m. a. k., & james, s. e. (2020). forecasting monthly maximum temperature of karachi city using time series analysis. pakistan journal of engineering, technology & science, 7(2). https://doi.org/10.22555/pjets.v7i2.2439. https://onlinelibrary.wiley.com/doi/book/10.1002/0470058005#page=13 https://books.google.co.id/books?hl=id&lr=&id=tl_ybaaaqbaj&oi=fnd&pg=pr5&dq=gubler+dengue+2014&ots=ctlavgemxi&sig=poapgum45gpxtfhkxsje09knnje&redir_esc=y#v=onepage&q=gubler%20dengue%202014&f=false https://books.google.co.id/books?hl=id&lr=&id=tl_ybaaaqbaj&oi=fnd&pg=pr5&dq=gubler+dengue+2014&ots=ctlavgemxi&sig=poapgum45gpxtfhkxsje09knnje&redir_esc=y#v=onepage&q=gubler%20dengue%202014&f=false https://www.sciencedirect.com/book/9780080178622/methods-and-evaluation-in-clinical-and-counseling-psychology https://www.sciencedirect.com/book/9780080178622/methods-and-evaluation-in-clinical-and-counseling-psychology https://apps.who.int/iris/bitstream/handle/10665/171002/db2011v35p99.pdf;jsessionid=455d9654c6994338d775c931e795e6ef?sequence=1 https://apps.who.int/iris/bitstream/handle/10665/171002/db2011v35p99.pdf;jsessionid=455d9654c6994338d775c931e795e6ef?sequence=1 https://apps.who.int/iris/bitstream/handle/10665/171002/db2011v35p99.pdf;jsessionid=455d9654c6994338d775c931e795e6ef?sequence=1 95 naureen zainab et al. / geosfera indonesia 6 (1), 2021, 77-95 thompson, r. n., gilligan, c. a., & cunniffe, n. j. (2016). detecting presymptomatic infection is necessary to forecast major epidemics in the earliest stages of infectious disease outbreaks. plos computational biology, 12(4). https://doi.org/10.1371/journal.pcbi.1004836. tran, b. l., tseng, w. c., chen, c. c., & liao, s. y. (2020). estimating the threshold effects of climate on dengue: a case study of taiwan. international journal of environmental research and public health, 17(4), 1–17. https://doi.org/10.3390/ijerph17041392. udayanga, l., gunathilaka, n., iqbal, m. c. m., najim, m. m. m., pahalagedara, k., & abeyewickreme, w. (2018). empirical optimization of risk thresholds for dengue: an approach towards entomological management of aedes mosquitoes based on larval indices in the kandy district of sri lanka. parasites and vectors, 11(1), 368. https://doi.org/10.1186/s13071-018-2961-y. wong, n. s., law, c. y., lee, m. k., lee, s. s., & lin, h. (2007). an alert system for informing environmental risk of dengue infections. in p. c. lai & a. s. h. mak (eds.), gis for health and the environment: development in the asia-pacific region with 110 figures (pp. 171–183). springer berlin heidelberg. https://doi.org/10.1007/978-3-540-71318-0_12. zavlavsky, i. (2000). a new technology for interactive online mapping with vector markup and xml. cartographic perspectives, 0(37), 65-77–77. https://doi.org/10.14714/cp37.810. 40 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 research article quantifying the significance of distance to temporal dynamics of covid-19 cases in nigeria using a geographic information system ifeyinwa sarah obuekwe1, 2* , umar saleh anka1,3 , sodiq opeyemi ibrahim1,3 , usman ahmad adam4 1nigerian environmental society, nigeria 2department of microbiology, faculty of life sciences, university of benin, benin city, edo state, pmb. 1154, nigeria 3 department of geography, faculty of earth and environmental sciences, bayero university, kano state, nigeria 4 department of geography, sa’adatu rimi college of education, kumbotso, kano, nigeria received 16 december 2020 /revised 8 april 2021 /accepted 17 april 2021 /published 25 april 2021 abstract the coronavirus disease 2019 (covid-19) is caused by a new strain of coronavirus that spreads primarily by close contact. although nigeria adopted lockdown measures, no defined strategies were used in setting the distance threshold for these lockdowns. hence, understanding the drivers of covid-19 is pivotal to an informed decision for containment measures in the absence of vaccines. spatial and temporal analyses are crucial drivers to apprehending the pattern of diseases over space and time. thus, this study aimed to quantify the significance of distance to the temporal dynamics of covid-19 cases in nigeria using the geographic information system. incremental spatial autocorrelation was used to analyze datasets of each month in arcgis. march, april, may, and june exhibited patterns with no significant peaks, while july and august exhibited patterns with two statistically significant peaks. the first and second peaks of july were 301,338.39 and 365,947.83 meters, respectively, while august was 301,338.39 and 336,128.09 meters, respectively. therefore, a significant difference in the clustering of covid19 over distances between july and august was established. this indicated that progression in the spread of the virus increased the virus's spatial coverage while the distance of risk of exposure decreased. this study's findings could be utilized to establish maximum movement restriction areas to contain the spread of covid-19. keywords: distance; incremental spatial autocorrelation; covid-19; disease; nigeria 1. introduction covid-19 is a respiratory disease that emerged in late 2019 caused by a new coronavirus that presents pneumonia-like symptoms. therefore, it has high transmission competence and spreads primarily when people are in close contact (dhama et al., 2020). geosfera indonesia *corresponding author. email address : ifeyinwa.obuekwe@uniben.edu (ifeyinwa sarah obuekwe) vol. 6 no. 1, april 2021, 40-54 p-issn 2598-9723, e-issn 2614-8528 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v6i1.21405 https://orcid.org/0000-0002-0187-7731 https://orcid.org/0000-0002-7984-0947 https://orcid.org/0000-0002-4635-8521 mailto:ifeyinwa.obuekwe@uniben.edu 41 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 covid-19 is causing an ongoing pandemic in many countries and territories, and therefore, it is a global health crisis (mckee & stuckler, 2020). infected people are accommodated and isolated in designated health facilities across the world for treatment. in the absence of a vaccine to prevent or aid speedy recovery, a significant concern is that, as the number of infected people increases, the healthcare system will be overwhelmed (murray et al., 2020). nevertheless, efforts to prevent further spreading are ongoing in many nations as vaccines, and antiviral drugs are awaited. in response to the worldwide covid-19 outbreak, many sub-saharan africa countries have implemented strict lockdown measures to emulate high-income countries to contain the spread of the virus (teachout & zipfel, 2020). however, severe restrictions placed on people's travel may impact different livelihoods, especially the world’s poorest people (woodhill, 2020). gradually, the impact of the lockdown is unfolding and remains largely guesswork. despite the likelihood of severe impact of lockdown on livelihoods, there is no adequate information to help guide the threshold with which movement restriction can be placed. details about places where cases may arise concerning identified cases could help enhance movement restriction decisions. wang et al. (2021) claimed that the temporal trends of covid-19 across space require deliberation when considering lifting lockdown. this makes the need for a timely knowledge of the temporal and spatial patterns of covid-19 transmission. however, knowledge of the geographical mean distances between the locations of cases and the direct estimation of distances between sequential cases requires both the recognition of cases and their infectors (salje et al., 2016). such contact tracing efforts are nearly impossible because of two major restrains, cost and time. thus, limited data hinders the ability to characterize the space and time scales where cases tend to occur. estimating mean transmission distances have been possible previously only in situations where there is available data for most cases in a transmission network or detailed epidemiological data on who infected whom (memish et al., 2014; marziano et al., 2017; yang et al., 2020). in the absence of data on where covid-19 cases were discovered in nigeria, there is the need to seek alternatives that can help in decision-making before the spread of the virus is beyond management capacity. hence, this study's novelty was that it demonstrated an approach to quantify the significance of temporal dynamics in the mean transmission distance of recorded covid-19 cases in nigeria. this was based on the use of geographic information system 42 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 software, data of the monthly covid-19 recorded cases, and the coordinates of the capital of the thirty-six (36) states in nigeria and the federal capital territory. this study aimed to quantify the significance of distance to the temporal dynamics of covid-19 cases in nigeria using the geographic information system. 2. methods 2.1 description of study area nigeria covers an area of 923,769 square kilometers, and it’s bordered by different countries on the north, south, east and west. on the north is niger and chad, south is gulf of guinea while the east and west are cameroon, benin and niger respectively (national bureau of statistics, 2011). 2.2 study area ecosystem the coast of nigeria is surrounded by flowering trees in marshes which are traversed by estuaries and rivers as well as the great niger basin (national bureau of statistics, 2011).eastwards however, lie successive belts of tropical rain forests and the undulating plateau, rising from 609.6 meters on the average to 1,828.8 meters. grassland vegetation interspersed with trees and shrubs is seen midway north of the country however, this stops at the north-east desert which is the grassy coast region. river niger receives its benue counterpart at lokoja, and from there flows into the atlantic ocean for about 547 kilometers. its tributaries include include sokoto, kaduna and anambra rivers (national bureau of statistics, 2011). however, the benue river is fed by two rivers namely katsina-ala and gongola. nigeria has other rivers which include cross river, benin and ogun rivers (national bureau of statistics, 2011). 2.3 study area population and geopolitical zones the population of nigeria in mid-year of 2020 is estimated at 206,139,589 people based on united nation data (worldometer, 2020). nigeria practices democratic system of government and each state and the federal capital territory has their elected officials. the elected officials are from political parties within the country. 43 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 figure 1. nigeria 36 states and the federal capital territory in their geopolitical zones (gis lab geography department, 2020) 2.4 research design the study was an explorative research. due to the newness of covid-19, explorative research has been adopted in several studies (scarpone et al., 2020) as it offers researchers the opportunity to study a problem that has not been thoroughly studied in the past. 2.5 data and sources of data the study made use of publicly available covid-19 and spatial datasets. the covid-19 dataset used include monthly recorded cases from march to august, 2020. the data were downloaded from nigeria centre for disease control (2020) website (https://covid19.ncdc.gov.ng/). this data is updated once a day (at 23:59 wat).the spatial data used for this study was shapefile of nigeria including states and coordinates of state capitals. the spatial data were sourced from the geographic information system laboratory, department of geography bayero university, kano. https://covid19.ncdc.gov.ng/ 44 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 2.6 data analysis the monthly covid-19 recorded cases for the 36 states and the federal capital of nigeria were subjected to incremental spatial autocorrelation analysis in arcgis 10.5 with the use of a conceptualized fixed distance. the incremental spatial autocorrelation analysis work in such a way that it calculates the moran’s i index and z score of a single data set at multiple distances. z-score is a standard deviation which measures the distance of a raw score from the mean. also, the rule of thumb for pattern analysis like spatial autocorrelation is that the z-score and p-values obtained determine the confidence level which can be used to accept or reject a null hypothesis. hence, the threshold for determining the confidence level of a z-score and p-value is shown in table 1. hence, analyzed fixed distances were derived from each of the possible distance between a data point (a state capital) and all other data points (other state capitals). as in the case of this study, moran’s i index was run for all the possible distances that existed between the 36 states and the federal capital of nigeria. instead of using the moran’s i index value, the distance was used to plot the z scores. this enabled a standardized distance-based comparison of statistical significance using the threshold in table 1. larger positive z scores indicate high significant clustering (table 1). thus, in each of the plots, distances of the first peak and maximum peak were identified, as the first peak indicates lower statistically significant clustering while the maximum peak indicates the distance that clustering or spatial autocorrelation was of most statistical significance in the data set. also, different colors were used to indicate the z-scores each distant point belong to as well as which distance is characterized as a peak. table 1. critical p-values and z-scores for different confident levels (moran, 1948) z-score (standard deviations) p-value (probability) confidence level (%) < -1.65 or > + 1.65 < 0.10 90 < -1.96 or > + 1.96 < 0.05 95 < -2.58 or > + 2.58 < 0.01 99 furthermore, the distance of each state capital to other state capital and the federal capital was first obtained in arcgis via point distance analysis. the obtained distance was then exported to excel where the mean was determined. in addition, bar graph was used to visualize the mean distance and distances where covid 19 clustering were significant. afterwards, 45 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 distances with significant peaks was compared with the mean distance of each state capital and the federal capital to other state capital and federal capital. the flowchart of the methodology for the study is summarized in figure 2. figure 2. flowchart of the method of the study 3. results and discussion after the index case of covid 19 in february, by the end of march, lagos has recorded 62 cases and was the highest state in nigeria with recorded cases (figure 3). this was followed by the federal capital territory, which had a recorded 9-28 victims. south western states that were closer to lagos were also observed to have recorded 4-8 victims. the cases recorded in kaduna, bauchi, ekiti, edo, benue, enugu, and rivers were within 1-3 while the other states had no recorded cases. alkhamis et al. (2020) also observed that after about a month from the first reported case of covid-19 in kuwait, the size was small despite sporadic infections. gayawan et al. (2020) also stated that covid-19 initially had a slow spread across africa, only for the situation to escalate in the last week of march. in april, although lagos still had the highest recorded cases, the spread of the virus in the north within the addition of one month is a course to worry about (figure 5). kano state that had no recorded case in march had the same range of 77-219 recorded cases with the federal capital territory. this implies that the rate of spread in the north was high in april, despite the federal government imposing a lockdown on lagos, ogun, and abuja on march 30, 2020. the high 46 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 spread indicates the lockdown was not as effective as that of china's strict lockdown (prem et al., 2020). may's spread exceeds that of march and april combined with kano having more recorded cases than the federal capital territory. despite the lockdown, border porosity, which is an issue in most african countries (okunade & ogunnubi, 2019), makes it hard to restrict mobility effectively. the ability to contain human mobility has been an issue worldwide as migrant workers and travelers have been the mostly affected and responsible for spreading the virus (kuwait government online, 2020; gayawan et al.,2020). there was an increase in the spread of covid-19 between june and july. however, the range of increase in abuja and oyo for july exceeded kano, edo, delta, and rivers state (figure 9 and figure 11). it was also observed that cross river was the last state to record the covid19 case in nigeria (figure 11). as of august, all the states in nigeria had a recorded case of more than 52 except kogi. also, states like kano, kaduna, plateau, oyo, ogun, ondo, edo delta, river lagos, and the federal capital territory each had more than 1000 recorded cases. the variability in the monthly pattern of covid-19 that can be observed from 3,5,7,9, 11 and 13 was how clustered the covid-19 cases are. in march, the first cluster was observed where states around lagos had high recorded cases (figure 3). the cluster that stood out in april, may, and june was that of kano state and the state around it. by july and august, more cases have been recorded in most states, thus resulting in a smaller number of states forming a cluster. studies of wang et al. (2021); kim & castro (2020), and xie et al. (2020) also observed spatial clustering about covid-19 distribution at a national scale. 47 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 figure 3. covid 19 recorded cases-march figure 4. spatial autocorrelation by distance-march figure 5. covid 19 recorded cases-april figure 6. spatial autocorrelation by distance-april ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 48 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 figure 7. covid 19 recorded cases-may figure 8. spatial autocorrelation by distance-may figure 9. covid 19 recorded cases-june figure 10. spatial autocorrelation by distance-june ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 49 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 figure 11. covid 19 recorded cases-july figure 12. spatial autocorrelation by distance-july figure 13. covid 19 recorded cases-august figure 14. spatial autocorrelation by distance-august figure 13. covid 19 recorded cases-august figure 14. spatial autocorrelation by distance-august ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 50 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 figure 15. mean distance of nigeria state capitals and the federal capital to other states capital and federal capital with peak category 0 100000 200000 300000 400000 500000 600000 700000 800000 900000 m e a n d is ta n c e i n m e te r state capital and fct with peak category ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 despite more visible clustering, the trend in the clustering of the covid-19 over distance shows no significant peak from march to june (figures 4, 6, 8, and 10). however, that of july and august exhibited patterns with two statistically significant peaks (figure 12 and figure 14). the first and second peak of july were 301,338.39 meter (z score =2.098; p<0.036) and 365,947.83 meter (z score = 2.299; p=0.022) respectively while that of august were 301,338.39 meter (z score =1.918; p=0.055) and 336,128.09 meter (z score =2.013; p=0.044) respectively (figure 12 and figure 14). the result shows that the mean distance of july and august's first peak was the same (301,338.39 meter). therefore, the distance where the spatial processes responsible for clustering is the most pronounced was about 301,338.39 meter. also, the second peak of august (336,128.09 meter) was less than that of july (365,947.83 meter). this implies that there was a significant drop in the mean transmission distance between july and august 2020. this was attributed to the fact that the spread has increased among neighboring state as more states recorded a similar range of covid-19 cases. blank et al. (2016) used a similar approach, observed spatial clustering of bacterial canker severity, and showed that disease severity had significant spatial autocorrelation with dominant pattern recognized as clustered. a drop in mean transmission distance throughout foot-and-mouth outbreaks was similarly reported by salje et al. (2016). kiang et al. (2020) also found a significant spatial correlation associated with covid-19 infections in china. each state capital and the federal capital had a mean distance of more than 300,000 meters to other states and federal capital (figure 15). although, only yola, maiduguri, sokoto, and damaturu have a mean distance above 600,000 meters (figure 15). this is because these states are in the far north of the country, bordering the niger republic. also, each of the state and the federal capital had a mean distance that is above the first and second peak of july and august except for lokoja that has a mean distance (361,924 meters) that was less than the second peak of july (365,947.83 meters) (figure 15). this implies that the interstate lockdown placed to curb the spread of covid-19 was justifiable as any movement within the minimum and maximum significant peak distance can fast-track the spread of covid-19. this corresponds with salje et al. (2016) 's assertion, who claimed that a fundamental property of infectious disease dispersal is the average spatial distance between transmission-linked cases. although this study does not compare covid-19 spread with population density nor socioeconomic activities. however, covid-19 cases have been clustered around two states in nigeria (kano and lagos). moreover, this cannot be farfetched from the fact that kano and lagos happen to be the most populous states in nigeria with vibrant industrial and 51 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 122 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 commercial activities. rapone (2020); rocklöv & sjödin (2020) also pointed out that high population density and socioeconomic status are essential drivers of the spread of covid-19. the disbelief regarding the presence of the virus also results in people taking preventive measures lightly. even six months after the first case of coronavirus was reported in nigeria, many nigerians still doubt its existence (csr-in-action, 2020). negligence could be detrimental to the spread to other states. kiang et al. (2020) also reported that people with a history of travel in epidemic areas fail to isolate themselves on time, thus accelerating transmission risk. thus, a study of this capacity that was able to showcase the significance of the temporal dynamics of covid-19 can help guide decision-making regarding how a movement should be restricted to reduce the clustering phenomenon observed across nigeria. 4. conclusion this study has shown that the point location of each state capital in nigeria and recorded cases for each of the 36 nigeria states and the federal capital territory, abuja displayed the significance of temporal dynamics in the mean transmission distance of recorded covid-19 cases in nigeria. this indicates that progression in the spread of the virus increased the spatial coverage of the virus while the distance of risk of exposure decreased. this information could be utilized to establish maximum movement restrictions to contain the spread of covid-19 through lockdown as implemented by the nigeria government. this method can also be applied to well-detailed data and a small proportion of adequately collected data. conflicts of interest the authors declare no conflict of interest. references alkhamis, m. a., al youha, s., khajah, m. m., ben haider, n., alhardan, s., nabeel, a., … al-sabah, s. k. (2020). spatiotemporal dynamics of the covid-19 pandemic in the state of kuwait. international journal of infectious diseases, 98(c), 153–160. https://doi.org/10.1016/j.ijid.2020.06.078. blank, l., cohen, y., borenstein, m., shulhani, r., lofthouse, m., sofer, m., & shtienberg, d. (2016). variables associated with severity of bacterial canker and wilt caused by clavibacter michiganensis subsp. michiganensis in tomato greenhouses. phytopathology, 106(3), 254–261. https://doi.org/10.1094/phyto-07-15-0159-r. 52 https://doi.org/10.1016/j.ijid.2020.06.078 https://doi.org/10.1094/phyto-07-15-0159-r 123 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 csr-in-action. (2020). the covid-19 disbelief amongst the nigerian citizenry. retrieved from https://www.csr-in-action.org/the-covid-19-disbelief-amongst-the-nigerian citizenry/. dhama, k., khan, s., tiwari, r., sircar, s., bhat, s., malik, y. s., ... & rodriguez-morales, a. j. (2020). coronavirus disease 2019–covid-19. clinical microbiology reviews, 33(4). https://doi.org/10.1128/cmr.00028-20. gayawan, e., awe, o., oseni, b., uzochukwu, i., adekunle, a., samuel, g., … adegboye, o. (2020). the spatio-temporal epidemic dynamics of covid-19 outbreak in africa. epidemiology and infection. https://doi.org/10.1101/2020.04.21.20074435. gis laboratory geography department. (2020). map of nigeria. kano : bayero university. kiang, m., chin, e., huynh, b., chapman, l., rodríguez-barraquer, i., greenhouse, b., … lo, n. (2020). routine asymptomatic testing strategies for airline travel during the covid-19 pandemic: a simulation analysis. the lancet infectious diseases. https://doi.org/10.1101/2020.12.08.20246132. kim, s., & castro, m. c. (2020). spatiotemporal pattern of covid-19 and government response in south korea (as of may 31, 2020). international journal of infectious diseases, 98(c), 328–333. https://doi.org/10.1016/j.ijid.2020.07.004. kuwait government online. (2020) covid-19 updates, state of kuwait. retrieved from https://e.gov.kw/sites/kgoarabic/pages/homepage.aspx. marziano, v., pugliese, a., merler, s., & ajelli, m. (2017). detecting a surprisingly low transmission distance in the early phase of the 2009 influenza pandemic. scientific reports, 7(1), 1–9. https://doi.org/10.1038/s41598-017-12415-2. mckee, m., & stuckler, d. (2020). if the world fails to protect the economy, covid-19 will damage health not just now but also in the future. nature medicine, 26(5), 640–642. https://doi.org/10.1038/s41591-020-0863-y. memish, z. a., cotten, m., meyer, b., watson, s. j., alsahafi, a. j., al rabeeah, a. a., … drosten, c. (2014). human infection with mers coronavirus after exposure to infected camels, saudi arabia, 2013. emerging infectious diseases, 20(6), 1012– 1015. https://doi.org/10.3201/eid2006.140402. moran, p. a. (1948). the interpretation of statistical maps. journal of the royal statistical society: series b (methodological), 10(2), 243-251. murray, c. j. (2020). forecasting covid-19 impact on hospital bed-days, icu-days, ventilator-days and deaths by us state in the next 4 months. medrxiv. https://doi.org/10.1101/2020.03.27.20043752. national bureau of statistics (2011). annual abstract of statistics. abuja : national bureau of statistics. nigeria centre for disease control (2020). covid-19 data. retrieved from https://covid19.ncdc.gov.ng/. 53 https://www.csr-in-action.org/the-covid-19-disbelief-amongst-the-nigerian%20citizenry/ https://www.csr-in-action.org/the-covid-19-disbelief-amongst-the-nigerian%20citizenry/ https://doi.org/10.1128/cmr.00028-20 https://doi.org/10.1101/2020.04.21.20074435 https://doi.org/10.1101/2020.12.08.20246132 https://doi.org/10.1016/j.ijid.2020.07.004 https://e.gov.kw/sites/kgoarabic/pages/homepage.aspx https://doi.org/10.1038/s41598-017-12415-2 https://doi.org/10.1038/s41591-020-0863-y https://doi.org/10.3201/eid2006.140402 https://rss.onlinelibrary.wiley.com/doi/abs/10.1111/j.2517-6161.1948.tb00012.x https://rss.onlinelibrary.wiley.com/doi/abs/10.1111/j.2517-6161.1948.tb00012.x https://doi.org/10.1101/2020.03.27.20043752 https://nigerianstat.gov.ng/pdfuploads/annual_abstract_of_statistics_2011.pdf https://nigerianstat.gov.ng/pdfuploads/annual_abstract_of_statistics_2011.pdf https://covid19.ncdc.gov.ng/ 124 ifeyinwa sarah obuekwe et al. / geosfera indonesia 6 (1), 2021, 40-54 okunade, s. k., & ogunnubi, o. (2019). the african union protocol on free movement: a panacea to end border porosity? journal of african union studies, 8(1), 73–91. https://doi.org/10.31920/2050-4306/2019/v8n1a4. prem, k., liu, y., russell, t. w., kucharski, a. j., eggo, r. m., davies, n., … klepac, p. (2020). the effect of control strategies to reduce social mixing on outcomes of the covid-19 epidemic in wuhan, china: a modelling study. the lancet. public health, 5(5), e261–e270. https://doi.org/10.1016/s2468-2667(20)30073-6. rapone, c. (2020); migrant workers and the covid 19 pandemic. united nations food and agriculture organization. retrieved from http://www.fao.org/documents/card/en/c/ca8559en/. rocklöv, j., & sjödin, h. (2020). high population densities catalyze the spread of covid19. journal of travel medicine. https://doi.org/10.1093/jtm/taaa038/5807719. salje, h., cummings, d. a. t., & lessler, j. (2016). estimating infectious disease transmission distances using the overall distribution of cases. epidemics, 17(c), 10–18. https://doi.org/10.1016/j.epidem.2016.10.001. scarpone, c., brinkmann, s. t., große, t., sonnenwald, d., fuchs, m., & walker, b. b. (2020). a multimethod approach for county-scale geospatial analysis of emerging infectious diseases: a cross-sectional case study of covid-19 incidence in germany. international journal of health geographics, 19(1), 1–17. https://doi.org/10.1186/s12942-020-00225-1. teachout, m., & zipel, c., (2020). the economic impact of covid-19 lockdowns in subsaharan africa. retrieved from https://www.theigc.org/wpcontent/uploads/2020/05/teachout-and-zipfel-2020-policy-brief-.pdf. wang, y., liu, y., struthers, j., & lian, m. (2021). spatiotemporal characteristics of the covid-19 epidemic in the united states. clinical infectious diseases, 72(4), 643– 651. https://doi.org/10.1093/cid/ciaa934. woodhill, j (2020). responding to the impact of covid-19 on rural people and food systems. retrieved from https://www.foresight4food.net/wp-content/uploads/2020/05/impactof-covid-19-on-rural-poverty-and-food-systems-v2.pdf. worldometer (2020). nigeria population. retrieved from https://www.worldometers.info/world-population/nigeria-population/. xie, z., qin, y., li, y., shen, w., zheng, z., & liu, s. (2020). spatial and temporal differentiation of covid-19 epidemic spread in mainland china and its influencing factors. science of the total environment, 744(c), 140929. https://doi.org/10.1016/j.scitotenv.2020.140929. yang, f., pahlavan, a. a., mendez, s., abkarian, m., & stone, h. a. (2020). towards improved social distancing guidelines: space and time dependence of virus transmission from speech-driven aerosol transport between two individuals. physical review fluids, 5(12). https://doi.org/10.1103/physrevfluids.5.122501. 54 https://doi.org/10.31920/2050-4306/2019/v8n1a4 https://doi.org/10.1016/s2468-2667(20)30073-6 http://www.fao.org/documents/card/en/c/ca8559en/ https://doi.org/10.1093/jtm/taaa038/5807719 https://doi.org/10.1016/j.epidem.2016.10.001 https://doi.org/10.1186/s12942-020-00225-1 https://www.theigc.org/wp-content/uploads/2020/05/teachout-and-zipfel-2020-policy-brief-.pdf https://www.theigc.org/wp-content/uploads/2020/05/teachout-and-zipfel-2020-policy-brief-.pdf https://doi.org/10.1093/cid/ciaa934 https://www.worldometers.info/world-population/nigeria-population/ https://doi.org/10.1016/j.scitotenv.2020.140929 https://doi.org/10.1103/physrevfluids.5.122501 richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 124 an assessment of the various mitigation strategies to combat desertification in jibia and kaita local government areas of katsina state richard sunday thlakma1*and omale eche john2 1department of geography, faculty of humanities, management and social sciences, federal university of kashere, pmb 0182, gombe state, nigeria 2department of geography, ahmadu bello university zaria, nigeria email: profsrthlakma@gmail.com received 7 april 2019/ revised 19 july 2019/ accepted 22 july 2019/ published 2 august 2019 abstract this study was conducted on an assessment of the various mitigation strategies to combat desertification in jibia and kaita local government areas of katsina state, nigeria. the data use includes satellite imageries for the study such as landsat mss of 1976, landsat tm of 1987, spot xs of 1995 and landsat etm of 2006 as well as structured questionnaires. sixty close ended copies of the questionnaire were administered in the study. purposive sampling method of administering questionnaires was adopted. the percentages land mass covered for each of these variables was determined and estimated in m2. literature was obtained from various agencies which were responsible for desertification control in katsina state. it was found from the reserved forest that in 1976 the percentage of reserved forest was 2.57%. in 1987 however, it increased by 73.9% to 76.47 %. by 1995, it declined by 9.42% to 67.05% and further declined by 0.52% in 2006. effort to combat desertification through the use of reserved forest has been quite significant over the years. also, noticed was a declined in shelter belt from 5.91% in 1987 to 1.097% in 1995 and a shot up to 7.39% in 2006. about 37% of the respondent opined that the deforestation leads to the disappearance of trees while 33% pinioned that it leads to reduction on agricultural productivity. the major strategy adopted to combat desertification is tree planting as supported by 88% of the respondents. it found that desertification as major environmental problem of the study area has reduced drastically from 43.34% in 1976 to 1.29% in 2006. it was also revealed from this study that some organizations such as european economic community/katsina state government eec/ktsg, katsina afforestation project unit ktapu and local government councils are the major agencies that are responsible for mitigating desertification in the study area. keywords: desertification, mitigation, afforestation, shelterbelt and satellite image 1. introduction campbell (1986) and mortimore (1989) define desertification as a process of sustained land degradation (loss of primary production) that results in the inability of the environment to sustain the demands being made upon it by socio – economic systems at geosfera indonesia p-issn 2598-9723, e-issn 2614-8528 vol.4 no. 2 (2019), 124-145, august, 2019 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v4i2.10192 richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 125 existing levels of technology and economic development and under prevailing climatic conditions especially drought. this definition reflects the situation in jibia and kaita local government areas of katsina state which is the study area in which the poor socio-economic activities led the residents to engaged in deforestation which led to the increase in deforestation and drought and also led to increase in micro climate of the study area leading to thereby increasing the rate and magnitudes of deforestation in the jibia and kaita local government areas. deforestation and the degradation of other vegetation, particularly near the margins of deserts have caused once fertile, vegetated land to become barren in a process called desertification. it has been suggested that desertification should be regarded as an extreme form of land degradation occurring where vegetation cover falls below 35 percent on a long term basis (binns 1990). although desert boundaries have shifted over time, deserts have always characterized the earth’s subtropical zones. global patterns of air circulation dictates that the subtropics are regions of subsiding air. when air subsides, it cools down and its capacity to hold moisture decreases, so inhibiting the formation of rain. this accounts for a prevalence of dry climates between latitudes 15o and 30o ‘north and south of the equator. however, these dry climates are extended into other latitudes and their patterns complicated by additional factors, such as distance from the rain-supplying oceans, seasonal high-pressure zones of large continental areas linked with monsoon systems, or the presence of mountain barriers over which air spills on the leeward side. thus creating rain shadows, (united nations conference on desertification nairobi, 1977). the first alarm on the south ward movement of the sahara desert into nigeria was raised by stebbing in (1935) partly in response to other environmental deterioration and desert encroachment. an anglo-french forestry commission in 1937 investigated the evidence of desertification in the northern part of nigeria, and found no evidence to support the report. however, this report created some concern in the colonial emirate. they embarked on tree planting to stop the southward drift of sahara desert, as far back as the 1940s. generally, areas most severely affected by desertification in nigeria are the semi-arid areas lying roughly north of latitude 12on, and these fall within the sudan and sahel zones of nigeria (sagua, et al, 1987). the extent and severity of desertification in nigeria is not fully established neither is the rate of its progression properly documented. however, it is generally agreed that it is by far the most pressing richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 126 environmental problem in the northern part of the country. the visible sign of the phenomenon is the gradual shift in vegetation from bushes and occasional trees to grass and bushes and in the final stages, expansive areas of desert like sand (ariyo, et al 2005). the united nation conventions to combat desertification has called for actions involving international cooperation and a partnership approach. it focuses on improving land productivity, the application of land conservation and sustainable management and water resources, (undpi, 1997). according to the federal government of nigeria (fos, 1999) estimates, between 50% and 75% of bauchi, borno, gombe, jigawa, kano, katsina, kebbi, sokoto, yobe, and zamfara states are being affected by desertification. these states, with a population of about 35 million people account for about 35% of the country’s total land area. table 1. spread of desertification in nigeria (fos, 1999). land area (km2) % of nigeria landmass affected bauchi/gombe 64,605 6.99 borno 70,890 7.67 yobe 45,502 4.93 kano 20,131 2.18 jagawa 23,154 2.51 katsina 24,192 2.62 sokoto/zamfara 65,735 7.12 kebbi 36,800 3.98 in addition, seven adjacent states to the south are reported to have about 10% to 15% of their land areas threatened by processes of desertification. it is estimated that the country is currently losing about 351,000 hectares of its landmass to desert conditions annually, and such conditions are estimated to be advancing south wards at the rate of about 0.6km per year (ariyo et al 2005). in the absence of concrete remedial and mitigative measures, it is estimated that the total cost of environmental degradation in nigeria would amount to about u.s. $5.110 billion per annum, 73% of which will result from land degradation alone (including desertification and erosion (ariyo et al., 2005). njeru, (2009) reported that sand dune stabilization is often done through the use of shelterbelt, woodlot and wind breaks. wind breaks are made from trees and bushes and are used to reduce water erosion and evapotranspiration. they were widely encouraged by development agencies from the middle east. also, the spraying of petroleum and nano clay over semi-arid crop land. this is often done in areas where richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 127 either petroleum or nano clay obtainable for example in iran. in both cases the application of the materials coats seedlings to prevent moisture loss and prevent them from being blown away. community participation in controlling desertification is also encouraged by the united nation convention to combat desertification which calls for a bottom up participatory approach in identifying, implementing, monitoring and evaluating projects that combat desertification and mitigate the effect of drought. project should be initiated and managed with maximum participation of local communities including women, youth, poor farmers, and pastoralists. the unccd offers local communities a very wide range of opportunities based on local communities, local circumstances, tradition, culture, norms, knowledge,and aspiration, (msafirri, 2009). katsina state is one of the states affected by desertification. it has both environmental and socio economic consequences. according to ktarda (2006). all the farmlands within kwangwalam area of katsina state are severely subjected to strong northeasterly harmattan wind, and sheet erosion, which sweeps the top soils, leaving the land bare and unproductive for agricultural purposes. murtala (2003) reported that desertification in katsina state has led to the reduction in the productivity of the land to the bearest minimum, thus eliminating the already meagre food and water resources for livestock. bala (2003) reports that katsina state has many desertification control projects.these include eec/ktsg, katsina arid zone programme, fgn/world bank assisted afforestation project, katsina state committee on drought and desertification control, federal department of forestry, katsina state ministry of agriculture, and katsina state afforestation project unit, ktpu. their major aim is afforestation, and land management to combat the menace of desertification. despite global concern about desertification and other forms of land degradation (soil quality decline, huge acres of land lost due to soil erosion, fertile soils erode away, overgrazing by animals such as cattle, goats, sheep, feed on the grass amongst others) and many years of efforts and investment made for prevention, cure or rehabilitation, the processes of land degradation persist. there are a variety of existing and potential prevention and mitigation strategies for the various desertification and land degradation process, such as physical measures and adaptation strategies integrating changing social, economic, institutional and policy factors. the mitigation of water-related land degradation is crucial in arid and semi-arid areas which are very prone to desertification richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 128 and successful measures invariably depend on an improvement of water storage in the soil (njeru, 2005). an appraisal of desertification strategies is important for the purpose of desert management in katsina state and other parts of nigeria. it is on this background the research aimed at assessment of the various mitigation strategies to combat desertification in jibia and kaita local government areas of katsina state by identify the various mitigation strategies employed in the lgas and assessing the degree of success of the mitigation strategies adopted for mitigating the desertification problem in the lgas. the study area covers jibia and kaita local government areas. these local government areas are located adjacent to each other. they are located approximately between latitudes 12,045n and 13,015n and longitudes 7,000e and 80,30e. the area is located in the extreme north western part of katsina state as show in figure 1. figure 1. map of the study area the climate is seasonally wet and dry. it is dominated by two air masses. these are the rain bearing south-westerly winds and the cold, dry and dusty north-easterly winds, locally known as the “harmattan.” at different times of the year, when one or the other winds prevails, the area experiences either rainfall or dry harmattan depending on the advance or retreat of the other (nyong and kanarogou, 1999). the vegetation here is sudan savannah type.it it is composed of trees scattered over an expanse of grassland. the trees are usually characterized by broad canopies e.g the baobab with large trunk.it is taller and larger than the others species or types found here.it also includes various types of acacias (nilofica, albida and seyal). others are richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 128 and successful measures invariably depend on an improvement of water storage in the soil (njeru, 2005). an appraisal of desertification strategies is important for the purpose of desert management in katsina state and other parts of nigeria. it is on this background the research aimed at assessment of the various mitigation strategies to combat desertification in jibia and kaita local government areas of katsina state by identify the various mitigation strategies employed in the lgas and assessing the degree of success of the mitigation strategies adopted for mitigating the desertification problem in the lgas. the study area covers jibia and kaita local government areas. these local government areas are located adjacent to each other. they are located approximately between latitudes 12,045n and 13,015n and longitudes 7,000e and 80,30e. the area is located in the extreme north western part of katsina state as show in figure 1. figure 1. map of the study area the climate is seasonally wet and dry. it is dominated by two air masses. these are the rain bearing south-westerly winds and the cold, dry and dusty north-easterly winds, locally known as the “harmattan.” at different times of the year, when one or the other winds prevails, the area experiences either rainfall or dry harmattan depending on the advance or retreat of the other (nyong and kanarogou, 1999). the vegetation here is sudan savannah type.it it is composed of trees scattered over an expanse of grassland. the trees are usually characterized by broad canopies e.g the baobab with large trunk.it is taller and larger than the others species or types found here.it also includes various types of acacias (nilofica, albida and seyal). others are richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 128 and successful measures invariably depend on an improvement of water storage in the soil (njeru, 2005). an appraisal of desertification strategies is important for the purpose of desert management in katsina state and other parts of nigeria. it is on this background the research aimed at assessment of the various mitigation strategies to combat desertification in jibia and kaita local government areas of katsina state by identify the various mitigation strategies employed in the lgas and assessing the degree of success of the mitigation strategies adopted for mitigating the desertification problem in the lgas. the study area covers jibia and kaita local government areas. these local government areas are located adjacent to each other. they are located approximately between latitudes 12,045n and 13,015n and longitudes 7,000e and 80,30e. the area is located in the extreme north western part of katsina state as show in figure 1. figure 1. map of the study area the climate is seasonally wet and dry. it is dominated by two air masses. these are the rain bearing south-westerly winds and the cold, dry and dusty north-easterly winds, locally known as the “harmattan.” at different times of the year, when one or the other winds prevails, the area experiences either rainfall or dry harmattan depending on the advance or retreat of the other (nyong and kanarogou, 1999). the vegetation here is sudan savannah type.it it is composed of trees scattered over an expanse of grassland. the trees are usually characterized by broad canopies e.g the baobab with large trunk.it is taller and larger than the others species or types found here.it also includes various types of acacias (nilofica, albida and seyal). others are richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 129 neem (azadrita indica) etc and various types of scattered shrubs such as plastima cassia sigrana, occasionally forming woodlot. most of these trees are xylophytic in nature i.e the ability to resist drought conditions, this is through having longroots, leathery leaves and tiny leaves, as well as development of spines to reduce excessive transpiration. the drought resistant plants do not shed their leaves, but other trees shed their leaves during the dry season. the blades of the grasses wither and dry off during the dry season. their underground stems develop new roots during wet season. grasses in the area hardly grows up to 1m at maturity. pasture in the area is very low and grazing fields are diminishing due to intensive cultivation by the growing population and shrubs usually dry off in the dry season leaving the soil bare, man-made plantations forms a significant portion of vegetation in this area. shelter belts woodlots and wind breaks are all found in the area (shittu, 1999). the study area falls within the chad formation which is made up of sedimentary rocks of cretaceous origin. the relief of the area is of the lower chad plains. the relative relief falls within the average of about 15 to 20 metres, the landform in the study area is sandy plains where great sand sheets occupy most of the area. however, it is not uncommon to find deposits of sand dunes of varying degrees especially when you move towards daddara kusa, kaga villages in jibia local government area and dankama, duma, jiki gisgerewa villages in kaita local government area northwards. the dunes are formed from the deposition of eroded materials carried along by strong winds from across the sahara desert, especially during the dry season, (shittu, 1999). the soil consists of mostly unconsolidated sediments which are predominantly sandy, silt to sandy loam, it is brown or reddish brown in nature. this type of soil is formed from the deposition of eroded materials over sedimentary formation. they are less acidic and well drained, with fairly low clay content. jibia local government area has a population of 167,435 people, while kaita local government area has a population of 182,405(national population commission 2006). the people of jibia and kaita local government areas are predominantly hausa and fulani living together, majority of the inhabitants of the areas are farmers, (subsistence). they also keep animals; almost every family keeps grazing animals like cattle, goats, and sheep. these are usually moved by herders from one place to the other in search of pastures especially during the dry season. some owners of these animals richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 130 move then to areas south of katsina state especially the middle belt region of the country during the dry season in search of past. 2. the methods identify the various mitigation strategies adopted to curb desertification in the study area. this was determined from the satellite imageries, landsat mss of 1976, 50m, landsat tm 30m of 1987, 1995 spot xs20m and 2006 landsat etm. these imageries also covered four decades with one for each decade. the percentages land mass covered for each of these variables was determined and estimated in m2. the imageries used to achieve these objectives, where analyzed by visual interpretation of the imageries. digital image processing of the satellite data and analysis of the stated objectives of the study were carried out using the basic software’s. earth resource data acquisition system, arcgis 10.0 software. these digital data on cd roms were procured from federal department of agriculture and land resources, livestock house mando, in kaduna and the national space research and development agency, abuja. assessing the degree of success for mitigation strategies adopted for mitigating the desertification problem in the lgas. this was determined from structured questionnaire, which was administered to respondents. sixty close ended copies of the questionnaire were administered in the study. purposive sampling method of administering questionnaires was adopted. in addition to the data obtained from the satellite imageries and questionnaire administered, literature was also obtained from various agencies which were responsible for desertification control in katsina state. some of the agencies are european economic community eec/katsina state government (eec/ktgs), katsina state agricultural and rural development authority (ktarda and the afforestation project unit (ktapu). descriptive statistics such as tables, mean, graphs, percentages, frequency, pie chart, bar graphs and line graphs were used for data analysis and presentation. 3. results and discussion 3.1 identify the various strategies adopted to combat desertification available statistics show that the strategies used to combat desertification in the study areas are: through reserved forests, reservoirs for irrigation purposes in other to richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 131 maintain and sustain the soil moisture contents so as to reduced decline in soil fertility, tree planting in the form of shelter belts, woodlots, trees on boundary, the use of vertiver grasses, natural tree regeneration, community and private woodlots (i.e. wind breaker). from the satellite imagery, the estimated land mass (in meter square) was estimated for each of the strategies in four decades: 1976, 1987, 1995 and 2006. percentages and line graphs were then used to analyse the data. the data extracted are presented in table 2. 0 20 40 60 80 100 1970 1980 1990 2000 2010 pe rc en ta g e year rf rf figure 3. trend of reserved forest area in the study area (1976 – 2006) source: author’s analysis richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 132 table 2. afforestation statistics in jibia and kaita lga year reserved forest (rf) in m2 rf (%) reservoir (r) in m2 r (%) tree planting e.g woodlots, natural regeneration, trees on farmland (tp) in m2 tp (%) shelter belt (sb) in m2 sb (%) afforestation total total 1976 1276564.019867 2.57 48331705.96 97.43 49608269.88 100.0 1987 164925844.03604 76.48 49446684.16 22.93 1276564.02 5.91 215649092.2 100.0 1995 193411475.14774 67.05 23312072.56 8.08 68566593.32 23.77 3164720.65 1.10 288454861.6 100.0 2006 199901900.712514 66.54 22958856.09 7.64 75342483.75 25.08 2220605.66 7.39 300423846.2sssssss 100.0 source: author’s analysis richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 133 3.1.1 forest reserve it was found from the reserved forest that in 1976 the percentage of reserved forest was 2.57%. in 1987 however, it increased by 73.9% to 76.47 %. by 1995, it declined by 9.42% to 67.05%. by 2006 it further declined by 0.52%. the decline in forest reserve is small compared to the much increase observed. effort to combat desertification through the use of reserved forest has been quite significant over the years. the line graph, showing the trend of reserved forest from 1976 to 2006 (in decades) is shown above. these changes are seen on figures 3, 4, 5 and 6, this is found in the coloration and sizes of the area coverage of the forest reserve areas, shown on the map, it is from these maps that the statistic for the area coverage of reserved forest where determined. 3.1.2 the use of reservoirs for irrigation purposes the result for reservoirs for irrigation also revealed that it declined in size due to decrease in irrigation farming, high level of deforestation, overgrazing by animals, over cultivation of lands for agricultural purposes and increase in climate change in the study area between 1995 and 2006 by 0.44%. the decline in the area coverage of water bodies (jibia dam in particular) between 1995-2006 was as a result of various purposes, the dam is used for e.g domestic purposes. the shrinking size also is attributable to the building up of sediments in the dam. vegetation is also gradually taking over the immediate surrounding of the dam, thereby reducing its size, the shrinking size of the dam may also be attributed to the amount of rainfall, of the particular year the imagery was taken, the amount of rainfall may be small. the dam however has played a very vital role in the irrigation of the vegetation of its immediate surroundings. this has helped to reduce desertification in the area drastically. figures 5 and 6 shows the dam with a little decrease in the area coverage of the dam, in plates 1 and 2 are different views of the dam which is used for irrigation purposes and also supply water for shelter belts, woodlots, etc. figure 7 shows a decrease in the area coverage of the dam. richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 134 figure 3. afforestation in jibia and kaita local government areas of katsina state (1976) figure 4. afforestation in jabia and kaita lga of kastina state (1987) richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 134 figure 3. afforestation in jibia and kaita local government areas of katsina state (1976) figure 4. afforestation in jabia and kaita lga of kastina state (1987) richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 134 figure 3. afforestation in jibia and kaita local government areas of katsina state (1976) figure 4. afforestation in jabia and kaita lga of kastina state (1987) richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 135 figure 5. afforestation in jabia and kaita lga of kastina state (1995) figure 6. afforestation in jabia and kaita lga of kastina state (2006) richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 135 figure 5. afforestation in jabia and kaita lga of kastina state (1995) figure 6. afforestation in jabia and kaita lga of kastina state (2006) richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 135 figure 5. afforestation in jabia and kaita lga of kastina state (1995) figure 6. afforestation in jabia and kaita lga of kastina state (2006) richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 136 7.4 7.6 7.8 8 8.2 1 2 pe rce nta ge decades 1995 and 2006 r r figure 7. reservoir area. 1995 and 2006 plate 1. different views of jibia dam which is used for irrigation in jibia plate 2. different views of jibia dam which is used for irrigation in jibia 3.1.3 shelter belts the total afforestation area reserved as shelter belt or wind breaker in the study area is small. it declines from 5.91% in 1987 to 1.097% in 1995 however, a shot up to 7.39% in 2006 as it shows in figure 8, the drop in the area coverage of shelter belt in the second decade result from poor maintenance and human interference. figure 8: shelter belt area.1987, 1995 and 2006 this evidence is also shown on the vegetal map, of the study area in figures 4, 5 and 6. figure 3 has no shelter belts on it, this is because shelter belts were not established in the sb(%) 1987 1995 2006 richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 137 first decade of this study. plates 3, 4, 5 and 6, f show evidence of afforestation in the study area. plate 3. grass hedging in kwana shagari in kaita local government area of (evidence of afforestation in the study area) plate 4. grass hedging in gidan saidu buzu in kaita local government area of (evidence of afforestation in the study area). plate 5. afforestation shelter belt in kurya in kaita local government area plate 6. afforestation shelter belt in magama in jibia local government area. 3.1.4 tree planting: woodlots, natural tree regeneration, trees on farmland, trees on farm boundaries etc. richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 138 from the satellite imagery, the data collected to assess tree planting as a strategy to combat desertification shows a decreasing trend in the effort towards tree planting, all through the sampled period. the decline in tree planting may be attributed to poor maintenance of the seedlings to maturity and also human interference and animals browsing. 3.2 determination of the success rate of various strategies used for combating desertification in the study area the third objective of the study is to determine the success rate of the various strategies used in combating desertification in the study area. from the literature, the following strategies were identified: the use of water dams for irrigation purposes, establishment of shelter belts, woodlots, wind breaks, trees on farm lands, natural tree regeneration and the planting of vertiver grasses for sand dune stabilisation. for the study area 60 questionnaires were issued to selected respondents. from the frequency analysis, a 100 percent response rate was recorded. this was because of the interest of the respondents on the issue of desertification. when asked if they knew what desertification was, all the respondents said yes. the frequency table is also shown table 3. table 3. knowledge of the term desertification in the study lgas frequency percent valid percent cumulative percent 60 100.0 100.0 100.0 source: author’s analysis 3.2.1 effects of desertification in the study area they were further asked what effects they would associate with desertification, 33.3% said that it makes productive lands to become unproductive; 10 % said it covers roads with sand; another 10% said it makes agricultural land to become unproductive; yet another 10% said it bring sand to cover their houses, and 36.7% said it leads to the disappearance of trees. the frequency table is shown in table 4. despite the much efforts put in place by various organisations to reduce desertification, there are still some areas where desertification is still a problem as shown in plates 7 and 8. richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 139 plate 7. evidence of desertification in magama, in jibia local government area of katsina state, nigeria. table 4. effects of desertification in the study lgas frequency percentage valid percentage cumulative percentage productive lands become unproductive 20 33.3 33.3 33.3 roads covered by sand 6 10.0 10.0 43.3 poor agricultural productivity 6 10.0 10.0 53.3 bring sand to home 6 10.0 10.0 63.3 trees disappeared 22 36.7 36.7 100.0 total 60 100 100 source: author’s analysis as to whether there are any organizations responsible for the mitigation of desertification in the study areas, all the respondents said yes. the organizations were stated to include: eec, ktapu, local government and state government as displayed in table 5. plate 8. evidence of desertification in dankama, in kaita local government area of katsina state, nigeria. the evidence to prove these are shown on plates 3, 4, display in table 5 and 6 where various organizations that are involves in curbing desertification in ahead their projects on grounds in the study area. richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 140 table 5. organisations involved in mitigating desertification in the study lgas frequency percent valid percent cumulative percent eec 27 45.0 45.0 45.0 state government 2 3.3 33.3 48.3 local government 2 3.3 33.3 51.7 ktapu 29 48.3 48.3 100.0 total 60 100.0 100.0 source: author’s analysis. 3.2.2 strategies for combating desertification on the various strategies adopted to combat desertification, 88.3% said tree planting was the dominant strategy employed e.g woodlots, trees on farmlands and natural tree regeneration. 11.7% said providing shelter belt (another means of tree planting that also serves as wind breaker) was the other strategy employed. the frequency table is shown in table 6. table 6. strategies for combating desertification in the study lgas frequency percent valid percent cumulative percent tree planting e.g woodlots, trees on boundary shelter belt total 53 7 60 88.3 11.7 100.0 88.3 11.7 100.0 88.3 100.0 source: author’s analysis. to determine the success rate of tree planting as a strategy of combating desertification, the respondents were asked to rate the effectiveness of the strategy; 15% said very effective, 50% said effective, 30 percent said fairly effective, while 5% said not effective. the frequency table is presented in table 7. based on this response it can be said that tree planting as a strategy for combating desertification in the study area is effective. there are also efforts been made on the part of the local government, community members to arrest desertification in the study area. some the efforts include: campaigning against tree felling, and by tree planting as well. the frequency table for this response is presented in table 8. richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 141 table 7. effectiveness of the strategies adopted to combat desertification in the study lgas frequency percent valid percent cummulative percent very effective 9 15 15 15 effective 30 50 50 65 fairly effective 18 30 30 95 not effective 3 5 5 100 total 60 100 100 source: author’s analysis table 8. community efforts at combating desertification frequency percent valid percent cummulative percent campaign against 27 45 45 45 tree planting 33 55 55 100 total 60 100 100 source: author’s analysis though results obtained from questionnaires may not be a very empirical judgment in arriving at the results obtained, which says that the strategies adopted for the mitigation of desertification in the study area are effective. the result obtained from the use of questionnaire corroborated with the empirical findings adopted earlier in the first and the second objectives of the study, it shows that desertification as an environmental problem of the study area has reduced drastically from 43.34% in 1976 to 1.29% in 2006. these findings are scientific in nature; it now substantiates the result obtained from the use of questionnaire. personal observation of the researcher, who has lived in study area for close to a decade, also is a proof to substantiate the results obtained from the questionnaires. the observations show that afforestation as a means of curbing desertification, is the major technique adopted in the study area, and it is effective. this evidence is also seen by the various photographs of the projects carried out by various organizations to curb desertification in the study area. from the evidences provided by the satellite imageries used for the research, and the literature review, the results obtained, and personal observations of the researcher, whom has lived in the area on which this research was carried for a decade. it shows that desertification richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 142 as an environmental problem of that region is gradually been phased out totally. this evidence is from various activities put in place to check desertification in this region. there has been an internationally assisted programme from the world bank, katsina arid zone programme which was later called european economic community/katsina state government, (eec/ktsg). agricultural development (ifad). katsina state agricultural and community development project (ksacdp).this are all internationally assisted programmes to combat desertification in katsina state, with the study area inclusive. this evidence also correlated with the (unccd, 1997) which called for international cooperation and a partnership approach in checking desertification. it focuses on improving land productivity, application of land conservation and sustainable management of water resources. the federal government of nigeria also plays a very big role in combating desertification in the study area, through the federal ministry of environment, ecological fund, katsina afforestation project unit/federal government assisted, though ktapu is now funded by the katsina state government. the establishment of the jibia dam which has totally erased desertification on its immediate surroundings are some federal government programmes, currently on ground in katsina state. their major roles are to check desertification. these evidences will support the result gotten in this research that shows that desertification has reduced drastically from about (43.34%) in 1976 to about 1.29% in 2006. all these are efforts put in place by the federal government with the study area inclusive, these programmes have aids in reducing desertification in katsina state drastically. these are practical phenomenon that is observed by the researcher. this relates to the work done by bala (2003), who reported that various agencies have desertification control projects on ground in katsina state. these includes the eec/ktsg, world bank assisted afforestation project. katsina committee on drought and desertification control, katsina afforestation project unit ktapu, katsina state ministry of agriculture etc. their major aim is afforestation and land management to combat the menace of desertification, as an environmental problem. the mitigation strategies used: afforestation is the major mitigation strategy used to combat desertification in the study area. this is done through reserved forests, tree planting in the form of woodlots, natural regeneration, trees on farmland shelter belts establishment, orchards the use of vertiver grasses and trees on farm boundaries. the personal observations of the researcher and the literature reviewed shows that the major mitigation strategies used in the study area by katsina afforestation project unit, (ktapu). katsina agric and rural richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 143 development authority, (ktarda). european economic community / katsina state government,eec/ktsg. adopted afforestation as the major mitigation strategy in checking desertification in the study area. the jibia dam also plays a very important role in checking desertification in the study area. this was supported by findings done by sani (1996). all the evidences given from the results obtained from the research and personal observations of the researcher has justified the findings which said desertification has reduced drastically in the study area. 4. conclusion the assessment of the success of the strategies used in combating desertification in the study area revealed that there has been an appreciable effort put in place to mitigate the effects of desertification. the area coverage of desertification determined from the satellite imageries shows that the spread was much at the first decade of the study, with an area coverage of about 43.34% in 1976 but it shrinked in size overtime to about 1.29% in 2006. afforestation statistics derived from the satellite imageries also shows reserve forest which the area coverage in 1976 was only 2.57% but increased overtime to 66.53% in 2006. the area coverage of shelter belt was 5.91% in 1987 and 7.39 in 2006. with the evidences given, it can now be considered that desertification in the study area has reduced greatly. results obtained from this research, shows that the major strategy for checking desertification in the study area is afforestation and the strategies are effective. this research to some extent has proved such notions not to be totally right, especially the food and agricultural organization (fao) report on desertification in nigeria which shows that desertification is increasing at an alarming rate of about 0.6km into nigeria but in this case it is reducing drastically. references ariyo, j.a, abdullahi, c.j. stigter, o.z onyewotu and i. musa (2005). community participation in planning desertification, control interventions in northern nigeria. lessons from kano state. a paper presented at the conference on prospects and problems of agricultural development in nigeria, held in ahmadu bello university zaria. june 29th – july 2nd, 2005. auwal, u. (2006). an appraisal of desertification in arid zone of bauchi state. unpublished pgdem thesis department of geography, bayero university kano. babura, d.u. (2001). desertifucation in babura local government area. unpublished pgdem thesis. department of geography, bayero university kano. richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 144 bala, a. (2003). an evaluation of drought incidence and hazards in northern nigeria. a paper presented at a seminar on the conservation of the environment. department of geography isa kaita college of education, dutsinma, katsina state 5th november. bins, t. (1990). the threatening deserts: ox blackwell. campbell, d.j (1986). the prospects for desertification in kajiado district kenya, geg. federal office of statistics (1999). annual abstracts of statistics: facts and figures about nigeria. katsina afforestation project unit (2006). a brief on katsina afforestation project unit (fgn assisted) 1st february, 2006. katsina afforestation project unit (2006). annual reports on afforestation projects. katsina afforestation project unit (2009). annual reports on afforestation project. katsina afforestation project unit ktapu. (2009). a brief on katsina afforestation project unit (fgn assited) 5th february, 2009. katsina agricultural and rural development authority ktarda. (2006). land management\ unit, historical background to soil conservation in katsina state. katsina arid zone programme eec/ktsg (1995). soil conservation experience, a paper for ksacdp one day soil conservation strategy workshop 14th december, 1995. katsina arid zone programme eec/ktsg. (1995). historical background of soil conservation efforts in katsina state. msafirri, f. (2009). involving or not involving communities in strategies to combat desertification and drought in kenya. murtala s.k. (2003). an analysis of the problems of desertification in katsina local government area. an nce project, isa kaita college of education dutsin ma, katsina state. nasiru, i.m. (2009). combating desertification and drought in nigeria. daily trust monday, 25th. national population commission (2006). federal republic of nigeria (2009), national population commission official gazette no2, abuja 2nd february 2009 vol 96. njeru, j. (2005). monitoring and modelling crop growth, water use and production crop growth, water use and production under dry land environment, north west of mount kenya. nyong, a.o. and kanoroglou, p.s (1999). the influence of water resources and their locations on rural distribution in north eastern nigeria. a journal of environmental sciences. richard sunday thlakma and omale eche john/geosi vol 4 no 2 (2019) 124-145 145 owen and pickering (1997). global environmental issues. rutledge and new york. sagua, v.o, enabor, e.e, k.o p.r.o, ojanuga a.v, mortimore, m. and kalu, a.e. (1987). ecological disasters in nigeria. drought and desertification federal ministry of science and technology lagos. sani, m. (1996). evaluation of desertification and its effects in jibia local government area of katsina state. unpublished pgdem thesis. department of geography bayero university kano. shittu, k. (1999). an assessment on the socio economic effects of desertification in katsina state. an unpublished bsc project department of geography bayero university kano. stebbing, e.p. (1935). the encroaching sahara: the threat to west africa colonies. a geographical journal. the united nations convention to combat desertification (1997). a new response an age – old problem. united nations conversation on desertification (1977). desertification, it causes and consequences: pergmon press. united nations, department for public information (1997). the united nations convention to combat desertification. whates and jones (1992). land degradation. edward arnold london. 301 research article land cover changes based on cellular automata for land surface temperature in semarang regency fahrudin hanafi1,* , dinda putri rahmadewi1, fajar setiawan2 1department of geography, faculty of social sciences, state university of semarang, sekaran gunungpati, semarang, 50229, indonesia 2limnology research center, indonesian institute of sciences (lipi), bogor, 16911, indonesia received 7 march 2021/revised 17 november 2021/accepted 8 december 2021/published 20 december 2021 abstract land cover changes based on cellular automata for surface temperature in semarang regency has increased significantly due to the continuous rise in its population. therefore, this study aims to identify, analyze and predict multitemporal land cover changes and surface temperature distribution in 2028. data on the land cover map were obtained from landsat 7 and 8 based on supervised classification, while land surface temperature (lst) was calculated from its thermal bands. the collected data were analyzed for accuracy through observation, while cellular automata markov chain was used to predict the associated changes in 2028. the result showed that there are 4 land cover maps with 5-year intervals from 2003 to 2018 at an accuracy of more than 85%. furthermore, the existing land covers were dominated by forest with decreasing trend, while the built-up area continuously increased. the existing land surface temperature range from 20.6°c to 36.6°c, at an average of 28.2°c and a yearly increase of 0.07°c. the temperature changes are positively correlated with the occurrence of land conversion. land cover predictions for 2028 show similar forest dominance, with a 23,4% built-up area at a surface temperature of 28.9°c. keywords: land cover change; cellular automata-markov chain; land surface temperature 1. introduction according to the united nations (2018), population increase is a global problem experienced in every country, with 55% of humans presently living in urban or regional areas, likely increasing by 68% in 2050. these changes tend to affect both local and global climate components, such as the land surface temperature (lst). for example, in nigeria, there was an increase of 19,166.13 ha in urban built areas from 2002 to 2013, with a rise in lst by 6 °c (igun & williams., 2018). geosfera indonesia vol. 6 no. 3, december 2021, 301-318 p-issn 2598-9723, e-issn 2614-8528 https://jurnal.unej.ac.id/index.php/geosi doi : 10.19184/geosi.v6i3.23471 *corresponding author. email address : fahrudin.hanafi@mail.unnes.ac.id (fahrudin hanafi) mailto:fahrudin.hanafi@mail.unnes.ac.id 303 fahrudin hanafi et al. / geosfera indonesia 6 (3), 2021, 301-318 khandelwal et al. (2018) stated that an increase in lst tends to disrupt the climateenergy balance, such as the heat wave phenomenon experienced in 7 major european countries, namely the united kingdom (38.1 °c), germany (41.7 °c), belgium (41.8 °c), france (42.6°c), luxembourg (40.8 °c), scotland (31.6 °c), and the netherlands (40.7 °c) recorded in july 2019. in southeast asia, several major cities also experienced similar conditions. an increase in hotspots was 20% greater than the average lst in hanoi (tran et al., 2017) and at a temperature of 2.9 ° c in jakarta which is higher than in bangkok (estoque et al., 2017). land cover changes also occurred in central java province. moreover, 128 ha of rice fields were converted to settlements or used for other purposes from 2009 to 2010 (bps, 2015). on the contrary, the average air temperature in central java province from 2032 to 2040 was predicted to increase within the range of 0.81 to 0.85 ° c (bmkg, 2019). semarang regency, central java province, indonesia, had a high population growth rate (8.74%) from 2010 to 2016 (bps, 2017a). based on statistical data, in 2016, 1.014 million people with a density of 1.081 people/km2, was recorded. this figure is higher compared to the national average population density of 127 people/km2. however, from 2011 to 2016, agricultural areas were reduced by 0.94% from 60,439.96 to 59,872.49 ha, while land used for other purposes was increased by 1.64%, which is equivalent to 35,148.18 ha (bps, 2017b). this indicates that semarang regency is also susceptible to land and climate changes problems, specifically areas adjacent to the city, which has experienced rapid development and recorded an lst average of 1.32 ° c (nugraha et al., 2016). analysis carried out using past and present spatial data is considered one of the requirements for geographic studies (dadras et al., 2015). cellular automata are the commonly used spatial model of land cover change. it is dynamic and composed of interrelated cells with discrete units (wang et al., 2012). cellular automata are used mainly to generate and predict potential changes (tran et al., 2017). fu & weng (2016) stated that temporal disparities of thermal characteristics due to land cover changes and responses need to be carried out comprehensively. one of the environmental parameters analyzed in this study is lst, estimated from the single thermal channel or split window algorithm method, dependent on the number of bands used (pu et al., 2006). both have a weakness in respect to the atmospheric profile uncertainty, which strongly affects the accuracy of the result (li et al., 2013). however, this is anticipated by inputting the atmospheric profile data into the thermal band spectral radian correction made by the usgs (united states of geological surveys) (coll et al., 2010). 304 fahrudin hanafi et al. / geosfera indonesia 6 (3), 2021, 301-318 study carried out on the land cover change in semarang regency is common for land suitability, flood (susanti et al., 2012), landslide, sedimentation (apriliyana, 2015), carbon stock, and spatial planning review (pangi et al., 2017). these studies were specifically related to land and averaged surface temperatures (kalinda & bandi, 2018). therefore, this study aims to model land cover changes based on raster data using cellular automata related to its surface temperatures in the future. in addition, it also intends to (1) analyze the surface temperature distribution and land cover changes of semarang regency in 2003, 2008, 2013, and 2018, (2) evaluate the relationship between land cover changes and lst, and (3) investigate the distribution of land cover for the following 10 years. 2. methods this field survey was conducted in semarang regency, central java province, from april to june 2019. the area was considered due to the record time of the imagery data input. furthermore, simulation data input only requires 2 land cover imageries, the initial and step year. however, for the sake of detailed information, this study used those acquired in 2003 (initial), 2008 (step 1), 2013 (step 2) and 2018 (step 3), which was compared using population growth and space needed, such as the assumption based on consistent population per built area on initial, and each step. the satellite image data used are (1) landsat 7 path/row 120/65 imagery recorded on may 20, 2003; (2) landsat 7 path/row 120/65 imagery recorded on june 18, 2008; (3) landsat 8 path/row 120/65 imagery recorded on june 24, 2013; (4) landsat 8 path/row 120/65 imagery recorded on august 25, 2018. unfortunately, landsat 7 (2008) had some bad qualities due to the sensor stripping. however, usgs provided corrections using past imagery with a similar location. secondary data used to support the population growth analysis were obtained from the (1) population growth and built area of semarang regency from 2008 to 2018, then (2) slope from astergdem radar image data (usgs). road networks, activity centers, and river patterns were obtained from big and spatial planning of semarang regency data and used as constraint input on land cover simulation. sampling calculation refers to the technical guidelines for collecting and processing spatial data from the geospatial information agency (big). meanwhile, the number of sample points for each land cover type is determined using the proportional stratified sampling method, as shown in figure 1. field surveys are carried out to measure the temperature of the land surface and cover the ground check. this analysis was carried out from april to june 2019, synchronized with 305 fahrudin hanafi et al. / geosfera indonesia 6 (3), 2021, 301-318 the imagery record period. it was assumed that the weather condition, sun duration, and intensity are similar to the imagery and survey data. also, the duration (on distribution sample) is customized from 08.00 to 11.00 am to adjust the recorded time of the imagery. figure 1. location and research sample in semarang regency, central java land cover is classified (maximum likelihood) using envi 5.1 and differentiated according to the method proposed by danoedoro (2006) concerning water, forest, shrub, agricultural, open, and built areas. the classification accuracy threshold used is 85%, thereby determining its mapping by comparing the 2018 image with field observations. data from the previous year's image is compared with the temporal interpretation of google earth. accuracy is determined using a confusion matrix involving the consideration of omission and commission. overall, it indicates the probability that a pixel belongs to a certain class and its representation in the field (lillesand et al., 2004). land cover prediction in 2028 was made with selva's version of idrisi software in accordance with the markov chain method based on cellular automata. meanwhile, markov chain is used to analyze 2 land cover data realized in different years, namely past (2008) and 306 fahrudin hanafi et al. / geosfera indonesia 6 (3), 2021, 301-318 actual information (2019). the transition matrix is focused on the change of any land cover to build area using a 3x3 matrix. built area conditions control this change from the forest, agriculture, or open field, asides the opposite. agriculture is changed from the forest, or open field, besides the opposite including build area, etc. lst is estimated by transforming pixel into spectral radian values using usgs equations to correct surface reflection errors and earth curvature. meanwhile, errors due to atmospheric disturbances, specifically in terms of processing images of land surface temperatures, are determined with the equation proposed by coll et al. (2010). it contains parameters that tend to affect lst, including emissivity, transmission, upwelling and downwelling (kalinda & bandi, 2018). the profile is obtained from the atmospheric correction parameter calculator, modeled according to data input's date, time, and location (tran et al., 2017). the corrected spectral radian values are converted to brightness temperatures using the usgs formula for landsat imagery. this is estimated as lst using the equation proposed by artist & carnahan (1982) and amiri et al. (2009). it is also used to determine an accurate brightness temperature of 8-14 чm wave (artis & carnahan, 1982), while this study utilized bands 6 (10.4 to 12.5 чm) and 10 (10.60 to11.19 чm). the analysis technique is used to determine the effect of each land cover type on lst changes by comparing the t-count with the t-table. the t-test is carried out using a simple linear regression equation while the t-table size is calculated with the formula t (a/2, n-2) = t (0.05/2, 116-2) = t (0.025, 114) = 1.98099, based on the following criteria: (1) assuming the significance value is <0.05 and t-count> t-table, it means that there is an impact, and (2) assuming the significance value is > 0.05 and t-count