Jurnal Sains Geografi, 3.
, 2024.
DOI: 10.
2210/jsg.
Identifikasi Environment Critical Index (Indeks Kekritisan Lingkunga.
Menggunakan NDVI dan LST di Kota Surabaya Ridho Agusman1,*.
Ilham Badaruddin Mataburu1 1 Geografi.
Universitas Negeri Jakarta.
Jl.
Mangun Muka Raya No.
Jakarta Timur.
DKI Jakarta, 13220 *) Email Korespondensi: ridhoagusman30@gmail.
Abstract Sitasi:
Agusman.
Mataburu.
Identifikasi Environment Critical Indeks (Indeks Kekritisan Lingkunga.
Menggunakan NDVI dan LST di Kota Surabaya.
Jurnal Sains Geografi.
Vol.
1,No.
Sejarah Artikel:
Diterima: 14 Januari 2023 Disetujui: 13 Maret 2023 Publikasi: 30 November 2024 The decline in environmental quality, which can lead to an increase in air temperature, is an adverse impact of the phenomenon of urbanization and industry.
An area will be more comfortable to live in if it is more densely vegetated.
The capital of East Java province is Surabaya, the second largest city in Indonesia.
ECI was analyzed based on LST and NDVI algorithms.
Observations of surface temperature and vegetation density in Surabaya City were conducted in 2013 and 2023 by utilizing Landsat-8 OLI satellite images.
NDVI classification levels are divided into 4, namely non-vegetation, low vegetation, medium vegetation, and high vegetation.
Low vegetation is the most dominating NDVI level in 2013 and 2023 with 80.
6% and 84.
9% respectively.
Surface Temperature Level in Surabaya City has increased in both maximum and minimum temperatures over the last 11 years.
It can be seen that the temperature distribution with a level of >36 degrees Celsius in 2023 experienced a very significant increase, reaching 75.
The Environmental Critical Index in Surabaya City has four classifications, namely non-critical, low, medium, high.
Each level of environmental criticality in Surabaya City was analyzed based on LST and NDVI within 11 years, namely 2013 and 2023.
The level of environmental criticality in Surabaya City has generally decreased due to the increase in vegetation area, although the level of decrease has not been Keyword: NDVI.
LST.
ECI Abstrak Copyright: A 2022 by the authors.
Submitted for possible open access publication under the terms and conditions of the Creative Commons Attribution (CC BY) license .
ttps://creativecommons.
org/license s/by/4.
0/).
Penurunan kualitas lingkungan yang dapat menyebabkan peningkatan suhu udara, merupakan dampak buruk dari adanya fenomena urbanisasi dan Sebuah wilayah akan lebih nyaman untuk dihuni jika semakin lebat vegetasinya.
Ibu kota provinsi Jawa Timur adalah Surabaya, kota terbesar kedua di Indonesia.
ECI dianalisis berdasarkan algoritma LST dan NDVI.
Pengamatan suhu permukaan dan kerapatan vegetasi di Kota Surabaya dilakukan pada tahun 2013 dan 2023 dengan memanfaatkan citra satelit Landsat-8 OLI.
Tingkat klasifikasi NDVI terbagi menjadi 4, yaitu non vegetasi, vegetasi rendah, vegetasi sedang, dan vegetasi tinggi.
Vegetasi rendah menjadi tingkat NDVI yang paling mendominasi pada tahun 2013 dan 2023 dengan masing-masing 80,6% dan 84,9%.
Tingkat Suhu Permukaan di Kota Surabaya mengalami peningkatan pada suhu maksimum maupun minimumnya selama 11 tahun terakhir.
Terlihat distribusi suhu dengan tingkat >36 derajat celcius pada tahun 2023 mengalami peningkatan yang sangat signifikan yaitu mencapai 75,6%.
Indeks kekritisan Lingkungan atau Environmental Critical Index di Kota Surabaya memiliki empat klasifikasi, yaitu non kritis, rendah, sedang.
Masing-masing tingkat kekritisan lingkungan di Kota Surabaya dianalisis berdasarkan LST dan NDVI dalam kurun waktu 11 tahun, yaitu tahun 2013 dan 2023.
Tingkat kekritisan lingkungan di Kota Surabaya secara umum mengalami penurunan karena bertambahnya luasan vegetasi, walaupun tingkat penurunannya belum signifikan.
Kata Kunci: NDVI.
LST.
ECI Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
Pendahuluan Penurunan kualitas lingkungan yang dapat menyebabkan peningkatan suhu udara, merupakan dampak buruk dari adanya fenomena urbanisasi dan industri.
Menurut (Dede et al.
, 2.
, ada beberapa aktivitas manusia yang berkontribusi terhadap pelepasan gas rumah kaca dan peningkatan suhu udara.
Iklim mikro di daerah perkotaan terpengaruh oleh perubahan penggunaan lahan dari vegetasi menjadi area terbangun, yang dibuktikan dengan suhu udara yang lebih tinggi di daerah perkotaan dibandingkan dengan daerah sekitarnya (Erwin.
, 2.
Lingkungan sebagian besar dipelihara oleh vegetasi.
Sebuah wilayah akan lebih nyaman untuk dihuni jika semakin lebat vegetasinya.
Namun, fragmentasi habitat dapat mengubah siklus biologis ekosistem karena perubahan hutan/lahan pembangunan berbagai fasilitas dan kegiatan lain yang menggunakan atau mengubah lanskap (Yanti.
, et al.
, 2.
Temuan ini menunjukkan adanya korelasi substansial antara pembacaan suhu udara permukaan dan tutupan vegetasi.
tinggi suhu udara di wilayah tersebut, semakin rendah tutupan atau kerapatan vegetasi (Zhang et al.
, 2.
Ibu kota provinsi Jawa Timur adalah Surabaya, kota terbesar kedua di Indonesia.
Sebagai salah satu kota besar di Indonesia.
Surabaya memiliki total populasi sebesar 2,87 juta jiwa, menurut hasil sensus BPS Kota Surabaya tahun 2020.
Kota Surabaya memiliki luas wilayah yang cukup besar .
,81 km.
dengan kepadatan penduduk 795 jiwa per km2.
Kepadatan penduduk Surabaya masih membaik Dalam kurun waktu sepuluh tahun, dari tahun 2010 hingga 2020, jumlah penduduk Surabaya bertambah sekitar 108,8 ribu jiwa, atau sebesar 3,94%.
Jika dibandingkan dengan tahun 1980, jumlah penduduk Surabaya telah meningkat hampir 42,47% (BPS Kota Surabaya, 2.
Pemusatan dan penyebaran pelayanan sosial, kegiatan ekonomi, dan kegiatan pelayanan pemerintahan terkonsentrasi di wilayah perkotaan (Dwijayanti.
, 2.
Penggunaan lahan di wilayah perkotaan terus berkembang setiap tahunnya akibat pertambahan penduduk yang sangat pesat.
Kriteria penggunaan lahan ini tentunya perlu mempertimbangkan daya dukung fungsi lahan dengan dukungan dataAedata secara saintifik atau ilmiah.
Surabaya merupakan salah satu wilayah kota yang berfungsi sebagai ibu kota provinsi dan pusat pemerintahan Jawa Timur (Dirk et al.
Adapun hampir dua pertiga dari luas wilayah Kota Surabaya adalah lahan Pembangunan area perumahan dan fasilitas komersial menjadi fokus utama pembangunan kota besar ini.
Menurut (Putra et al.
, 2.
, kondisi tata ruang Kota Surabaya telah mengalami penyusutan secara drastis, yang menunjukkan bahwa area yang dulunya digunakan untuk pertanian, lahan kosong, dan ruang terbuka hijau atau jalur hijau telah berubah menjadi area yang digunakan untuk perumahan, komersial, dan jasa.
Hal ini terjadi sebagai akibat dari perubahan tutupan lahan di Surabaya.
Pertumbuhan bangunan mau tidak mau diperlukan untuk mengimbangi hilangnya ruang terbuka hijau yang disebabkan oleh (Khasanah Priyatmono, 2.
Menurut (Nofrizal, 2.
, jumlah lahan hijau di daerah perkotaan lebih sedikit dibandingkan dengan daerah pinggiran kota, sehingga meningkatkan suhu dan dapat menyebabkan terbentuknya fenomena Urban Heat Island.
Setiap tahun, terjadi peningkatan pembangunan gedung Ae gedung di pusat kota, yang mengakibatkan beberapa daerah yang padat penduduknya memiliki vegetasi yang lebih sedikit.
Keanekaragaman tanaman yang semakin berkurang digantikan oleh permukaan aspal dan beton.
Dibandingkan dengan lokasi yang masih didominasi oleh vegetasi, suhu udara di area yang dibangun akan meningkat secara signifikan.
Emisi panas disebabkan Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
oleh karakteristik material yang digunakan untuk membangun konstruksi perkotaan.
Peningkatan suhu juga disebabkan oleh limbah panas yang dihasilkan oleh perusahaan, pendingin ruangan (AC), dan emisi kendaraan bermotor (Fazeli et al.
, 2016.
Li & Norford, 2016.
Morris et al.
, 2016.
Senanayake et al.
, 2.
Kekritisan lingkungan yang disebabkan oleh efek negatif tersebut mengakibatkan kebutuhan energi yang lebih tinggi, polusi udara dan emisi gas rumah kaca yang lebih banyak, kualitas air yang lebih buruk, dan berdampak pada kenyamanan dan kesehatan (U.
Environment Protection Agency, 2.
Kenaikan suhu udara dapat dikurangi dengan adanya vegetasi.
Mekanisme evapotranspirasi, pembayangan, dan proses fotosintesis dapat menurunkan suhu udara dan mendinginkan atau menurunkan sedikit suhu lingkungan sekitar serta menekan atmosfer (Indrawati et al.
, 2.
Indeks Vegetasi Ternormalisasi (Normalized Difference Vegetation Index.
NDVI) dan (Land Surface Temperatur.
merupakan rasio yang menentukan Indeks Kekritisan Lingkungan (Environmental Criticality Inde.
atau disebut ECI (Senanayake et al.
, 2.
ECI
kemudian dikembangkan menjadi sebuah persamaan yakni suhu udara dan Indeks Area Terbangun (Built-Up Area Inde.
Semakin tinggi suhu udara dan semakin rendah kerapatan vegetasi, maka ECI dapat menjadi pertanda kekritisan lingkungan (Fadlin et al.
, 2.
Lebih lanjut, (Senanayake et al.
, 2013.
menyatakan bahwa nilai ECI yang tinggi dapat digunakan untuk menentukan lokasi UHI.
Oleh karena itu, sangat penting untuk memahami kondisi dan sifat ekosistem di Kota Surabaya.
Dengan menentukan Environmental Critical Index (ECI) atau Indeks kekritisan lingkungan dengan analisis citra pada tahun 2013 dan 2023, analisis dapat dilakukan untuk mengetahui kualitas lingkungan.
Metode Penelitian ECI dianalisis berdasarkan algoritma LST dan NDVI.
Pengamatan dinamika suhu permukaan dan kerapatan vegetasi di Kota Surabaya dilakukan pada dua periode yakni pada tahun 2013 dan 2023 dengan memanfaatkan citra satelit Landsat-8 OLI.
Penelitian ini menggunakan Band 10 untuk mengestimasi nilai brightness temperature dan dikonversi ke suhu permukaan serta Band 4 dan Band 5 untuk menghitung NDVI.
Untuk menghasilkan nilai ECI, langkah yang dilakukan adalah mengidentifikasi LST dan NDVI.
Adapun tahapan identifikasi sebagai berikut:
Penentuan Top of Atmosphere (ToA) Radiance (L).
Tahap awal dalam algoritma ini adalah input Band 10, yaitu konversi digital number (DN) menjadi spectral radians (L) dengan rumusan berikut (Avdan & Jovanovska, yayuI = ycAya O ycEycaycayco yaya.
Adapun:
= Nilai Radianc.
Pancaran (Toa ML = faktor multiplicative rescaling pada spesifik band Qcal = Digital Number (DN) pada piksel disetiap band AL = faktor additive rescaling pada spesifik Band Penentuan Brightness Temperature (BT).
Tahapan selanjutnya adalah konversi nilai spectral radians yang diperoleh dari persamaan .
menjadi atSatellite Brightness Temperature dengan menggunakan Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
algoritma sebagai berikut (Weng et al.
yaAycN = .
acuan untuk proporsi vegetasi dan tanah Nilai Pv dapat diperoleh (Digavinti et al.
, 2017.
Sobrino et al.
ycEyc = ( Adapun:
BT = Brightness Temperature dalam satuan kelvin.
K1 = Konstanta dari kalibrasi suhu absolut (K) K2 = Konstanta dari kalibrasi spektral radian pada spesifik Band yang terdapat pada metadata citra.
L = Nilai Pancaran.
Adapun:
NDVImin = Nilai minimum pada pixel
NDVI
NDVImax = Nilai maksimum pada pixel
NDVI
Nilai Pv yang telah diperoleh akan digunakan untuk menghitung nilai LSE, dengan persamaan berikut (Sobrino et , 2.
Nilai BT kemudian dikonversi ke dalam satuan derajat Celsius dengan persamaan:
yaAycN (EE) = yaAycN .
Oe 273,1.
ya = 0,004 ycU ycEycO 0,986.
Penentuan Indeks Kerapatan Vegetasi (NDVI).
NDVI merupakan pengukuran kehijauan vegetasi yang berguna untuk memahami kerapatan pada vegetasi dan menilai perubahan kesehatan pada NDVI dihitung dengan persamaan berikut:
0,004 = nilai rata-rata emisivitas vegetasi rapat ycAyaycOya = ( .
Adapun:
0,986 = nilai emisivitas standar lahan Estimasi Land Surface Temperature (LST).
Perhitungan LST (Digavinti et al.
, 2.
Adapun:
yaycIycN = DVI = (Band 5 Ae Band .
/ (Band 5 Band .
, .
RED = nilai Digital Number pada RED Band Adapun:
NIR = nilai Digital Number pada NIR Band = panjang gelombang tengah Band 10 Penentuan Land Surface Emissivity (LSE).
Untuk menghitung emisivitas permukaan, diperlukan menghitung proportional vegetation (P.
, yaitu proporsi vegetasi yang memberikan estimasi setiap area di bawah tutupan Nilai pixel dari NDVI menjadi ( ))] BT = Brightness Temperature .
C).
E = emisivitas = 14.
380 mK Penentuan Indeks Kekritisan Lingkungan (ECI).
ECI merupakan perbandingan nilai LST dengan nilai NDVI yang menunjukkan tingkat Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
kekritisan lahan pada suatu wilayah Jika semakin tinggi nilai ECI, maka semakin tinggi kekritisan lahan di wilayah tersebut.
stretching nilai piksel .
dilakukan untuk menyamakan kedua nilai piksel antara LST dan NDVI.
Formula ECI pada penelitian (Fitriani et al.
, 2.
dijabarkan sebagai persamaan berikut:
yayaya = Adapun tingkat indeks kekritisan lingkungan ini dibagi menjadi 4 kelas, yaitu tidak kritis, rendah, sedang, dan Hasil dan Pembahasan 1 Hasil NDVI Kota Surabaya Tahun 2013 dan .
Gambar 1.
NDVI Kota Surabaya Tahun 2013 Gambar 2.
NDVI Kota Surabaya Tahun 2023 Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
NDVI
Normalized Difference Vegetation Index merupakan salah satu menganalisis ECI.
Berdasarkan hasil analisis NDVI Kota Surabaya tahun 2013 dan 2023 tidak ada perubahan yang signifikan pada kedua tahun tersebut.
Adapun terlihat pada peta.
Kota Surabaya masih didominasi oleh tingkat kerapatan vegetasi yang rendah yang berwarna oranye dalam kurun waktu 2013-2023.
LST Kota Surabaya 2013-2023 Gambar 3.
Land Surface Temperatur (LST) Kota Surabaya Tahun 2013 Gambar 4.
Land Surface Temperatur (LST) Kota Surabaya Tahun 2023 Selain NDVI sebagai parameter dalam identifikasi ECI.
LST atau Land Surface Temperature juga menjadi salah satu parameter dalam menganalisis ECI.
Berdasarkan hasil analisis LST Kota Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
Surabaya tahun 2013 dan 2023 terdapat perubahan yang sangat signifikan terutama tingkat suhu >36 derajat celcius dari distribusi suhu daratan Kota Surabaya ini.
Terlihat pada peta bahwa tingkat suhu > 36 elcius yang berwarna merah semakin meluas dalam kurun waktu 2013 Ae 2023.
ECI Kota Surabaya 2013-2023 Gambar 5.
Environmental Critically Index (ECI) Kota Surabaya Tahun 2013 Gambar 6.
Environmental Critically Index (ECI) Kota Surabaya Tahun 2023 Berikut Environmental Critical Index (ECI) Kota Surabaya tahun 2013 dan 2023.
Pada peta terlihat bahwa distribusi ECI di Kota Surabaya digambarkan dengan warna merah Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
sebagai tingkat ECI tinggi, sedangkan yang berwarna biru merupakan tingkat ECI yang tidak kritis.
2 Pembahasan NDVI Kota Surabaya 2013-2023 Distribusi nilai kerapatan vegetasi Kota Surabaya diperoleh dengan analisis NDVI.
Klasifikasi NDVI pada tahun 2013 dan 2023 agar mengetahui tingkat perubahan kerapatan vegetasi.
Menurut Nurgiantoro & Aris .
dalam Fitriani et.
area yang memiliki nilai NDVI >0,50 merupakan kategori bervegetasi penuh.
Nilai NDVI tertinggi pada tahun 2013 yaitu sebesar 0,598 dan terendah ada pada nilai -0,17.
Sedangkan nilai NDVI tertinggi pada tahun 2023 yaitu sebesar 0,58 dan terendah -0,18.
Berdasarkan pengklasifikasian NDVI ini ditentukan dengan 4 kelas yaitu non vegetasi (-1 Ae .
, lalu vegetasi rendah .
Ae 0,.
, vegetasi sedang .
,25 Ae 0,.
dan vegetasi tinggi (>0,.
Table 1.
Klasifikasi NDVI Kota Surabaya
Tahun 2013
NDVI 2023
Luas Ha Persen
Non Vegetasi
1129,49
3,44
Rendah
26476,83
80,63
Sedang 5204,32 15,85 Tinggi 25,48 0,08 yaitu non vegetasi, vegetasi rendah, vegetasi sedang, dan vegetasi tinggi.
Vegetasi rendah menjadi tingkat NDVI yang paling mendominasi pada tahun 2013 dan 2023 dengan masing-masing 80,6% dan 84,9%.
Ada peningkatan tingkat vegetasi rendah kurang lebih 4% selama 11 tahun.
Lalu tingkat kerapatan dengan non vegetasi mengalami penurunan sekitar 1% dalam kurun waktu 11 tahun.
Untuk tingkat kerapatan sedang mengalami penurunan sekitar 3% dan tingkat kerapatan tinggi juga mengalami penurunan yang tidak terlalu signifikan sebanyak 0,01% dalam kurun waktu 11 tahun terakhir.
Land Surface Temperatur (LST) Kota Surabaya 2013-2023 Analisis Land Surface Temperature (LST) menunjukkan adanya distribusi serta tingkatan suhu permukaan yang dilakukan pada 11 tahun terakhir yaitu tahun 2013 dan 2023 di Kota Surabaya.
Terlihat pada peta bahwa distribusi suhu maksimum yang awalnya terkonsentrasi pada bagian tengah Kota Surabaya meluas ke segala arah.
Nilai suhu maksimum pada tahun 2013 yaitu 38,8 derajat celcius dan nilai suhu minimum sebesar 25,2 derajat Sedangkan pada tahun 2023, nilai suhu maksimum sebesar 46,3 derajat celcius dan nilai suhu minimum sebesar 26,8 derajat Table 3.
Klasifikasi Land Surface Temperature (LST) Kota Surabaya Tahun 2013 Table 2.
Klasifikasi NDVI Kota Surabaya
Tahun 2023
LST 2023
Luas Ha Persen
NDVI 2023
Luas Ha Persen
1327,76
4,04
Non Vegetasi
883,30
2,69
1214,03
3,70
Rendah
27885,76
84,93
2051,27
6,25
Sedang 4043,55 12,31 3396,48 10,34 Tinggi 22,45 0,07 >36 24849,02 75,67 Table 4.
Klasifikasi Land Surface Temperatur Tingkat klasifikasi NDVI terbagi menjadi 4.
Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
(LST) Kota Surabaya Tahun 2023
LST 2023
Luas Ha Persen
6225,41
18,96
5055,00
15,40
9652,49
29,39
10465,40
>36
1439,82
klasifikasi Quantile.
Table 3.
Klasifikasi Environmental Critically Index (ECI) Kota Surabaya Tahun 2013
ECI 2023
Luas Ha Persen
31,87
Tidak Kritis 7887,81 24,02 4,38 Rendah
8485,07
25,84
Sedang 25,23 Tinggi 8180,03 24,91 Tingkat Suhu Permukaan di Kota Surabaya maksimum maupun minimumnya selama 11 tahun terakhir.
Terlihat bahwa distribusi suhu dengan tingkat >36 derajat celcius pada tahun 2023 mengalami peningkatan yang sangat signifikan yaitu mencapai 75,6% dengan luasan sekitar 24 ribu hektar.
Ini dapat dibandingkan jika dilihat distribusi suhu dengan tingkat >36 pada tahun 2013 hanya sebesar 4% dengan luasan hanya sekitar 1400-an hektar.
Artinya terjadi peningkatan suhu maksimum sekitar 71% selama 11 tahun di Kota Surabaya.
Peningkatan suhu yang signifikan ini dapat diasumsikan karena dalam 11 tahun terakhir terjadi peningkatan aktivitas manusia sebagai bagian dari urbanisasi seperti meningkatnya jumlah polusi akibat kendaraan bermotor dan kegiatan industri, meningkatnya area pemukiman.
Environmental Critically Indeks (ECI) Kota Surabaya 2013Ae2023 Analisis ECI di Kota Surabaya terlihat pada peta menunjukkan tidak terdapat perubahan yang signifikan antara 11 tahun terakhir, yaitu tahun 2013 dan 2023.
Tingkat kekritisan lingkungan tinggi digambarkan berwarna merah pada peta, sedangkan yang berwarna biru merupakan tingkat kekritisan lingkungan yang tidak kritis.
Indeks kekritisan lingkungan ini terbagi menjadi 4 klasifikasi, yaitu tidak kritis, rendah, sedang, dan tinggi.
Masing-masing pembagian kelas Table 4.
Klasifikasi Environmental Critically Index (ECI) Kota Surabaya Tahun 2023
ECI 2023
Luas Ha Persen
Tidak Kritis 7618,15 23,20 Rendah
8459,92
25,76
Sedang 8393,05 25,56 Tinggi 88366,08 25,48 Tidak terdapat perubahan ECI yang signifikan pada 11 tahun terakhir.
Adapun peningkatan terjadi justru pada tingkat tidak kritis sebanyak hampir 1%.
Selain itu tingkat kekritisan rendah juga mengalami kenaikan walau hanya sebesar 0,1%.
Sedangkan tingkat sedang dan rendah mengalami penurunan masing-masing 0,3% dan 0,5%.
Tingkat kekritisan lingkungan di Kota Surabaya ini mengalami penurunan dapat disebabkan karena adanya pertumbuhan vegetasi selama 11 tahun terakhir.
Terlihat bahwa luasan non vegetasi yang sebelumnya mempunyai luas sekitar 1129 hektar pada tahun 2013, berkurang menjadi 800-an hektar pada tahun 2023.
Oleh karena itu, peningkatan vegetasi di suatu wilayah dapat mengurangi tingkat kekritisan lingkungan.
Kesimpulan Indeks Lingkungan Environmental Critical Index di Kota Surabaya memiliki empat klasifikasi, yaitu non kritis, rendah, sedang, tinggi.
Masing22 Jurnal Sains Geografi, 1.
, 2023.
DOI: 10.
2210/jsg.
masing tingkat kekritisan lingkungan di Kota Surabaya dianalisis berdasarkan Land
Surface Temperature (LST) dan Normalized Difference Vegetation Index (NDVI) dalam kurun waktu 11 tahun, yaitu tahun 2013 dan Tingkat kekritisan lingkungan di Kota
Surabaya secara umum mengalami penurunan karena bertambahnya luasan vegetasi, walaupun tingkat penurunannya belum signifikan Daftar Pustaka