1. Home
  2. Archives
  3. Vol 26 (2020) Issue 2
  4. Articles

Analisis Kapasitas Beban Polutan Ruas Hilir Sungai Citarum Menggunakan Qual2kw

Abstract

. The downstream segment of Citarum River flows from Bendung Curug, Karawang Regency and ends at Muara Gembong, Bekasi Regency. The pollution of downstream segment of Citarum River is caused by settlement, agriculture, farm, and industrial activities that dump their waste into river bodies. This causes a decrease in water quality and reduce the carrying capacity of the downstream segment of Citarum River pollutant. The purpose of this research was to analyze the pollutant load capacity (PLC) in downstream segment of Citarum River using QUAL2Kw model. This research was carried out downstream of Citarum River from segments 16 to 20 which were divided into 9 sections with 11 locations monitoring water quality. Sampling was carried out for 3 periods of observation that the sampling is done twice at each location using grab sampling method. Water quality parameters being modeled are BOD, COD, nitrate, and phosphate. The results of the simulation model of PLC for BOD was 14556.67 kg/day-23456.84 kg/day, COD was 24261.12 kg/day-39094.73 kg/day, nitrate was 24261.12 kg/day-39094.73 kg/day, and phosphate was 485.22 kg/day-781.89 kg/day. Overload of BOD was 0%-26%, COD was 61%-68%, phosphate was 78%-7%. Meanwhile, nitrate load in the Citarum River downstream segment doesn

Keywords

INTRODUCTION

The downstream segment of Citarum River starts at Bendung Curug, Karawang Regency and ends at Muara Gembong, Bekasi Regency. Citarum River is included in the category of heavily polluted rivers (KLHK, 2016). Pollution of the river can cause a decrease in water quality. The decline in water quality of Citarum River occurs due to water pollution by industrial and domestic waste discharges that enter the river (KLHK, 2016). Water pollution is the release of all liquids, solids, gases, pathogenic organisms or other substances into the water which will cause the water to be disrupted and become dangerous for the health, safety or welfare of the community and the environment. Water pollution usually comes from human activities such as manufacturing, agriculture, mining, and poor waste management (Masere et al, 2012).

Water pollution occurs when unwanted wastes are spread in the water system so that the water quality changes (Alrumman et al, 2016).

River pollution is defined as a change in the quality of a waters due to human activities, which in turn will disrupt the lives of humans themselves or other living things. These changes can be caused by compounds that enter the river flow that moves downstream along with the water flow or are stored at the bottom, accumulating (especially on sediment) and someday there can also be washing or dilution. These compounds (mainly toxic) accumulate and become a certain concentration which is harmful to the life chain Patricia et al, 2018). Source of river pollution can be classified into point source and non-point source. Point sources refer to sources with discharges entering water bodies in certain locations such as through pipes such as discharges from industry, domestic, and wastewater treatment plants. While non-point sources come from scattered sources such as agricultural activities and runoff water (Afroz et al, 2017).

Pollutant load is the amount of a pollutant contained in water or wastewater. Pollution burden is essentially the amount of pollutant mass in a body of water in a certain period. Pollutant load is expressed in units of pollutant load per unit time (Komarudin, et al, 2015). Pollutant loads entering the downstream segment of Citarum River comes from domestic activities, industry, agriculture, restaurants, hotels, shopping centers, and hospitals. The high load of pollutants entering the river body can reduce the capacity of pollutant load (DTBP) of the river. The purpose of this study is to analyze the pollutant load capacity (PLC) using QUAL2Kw model.

PLC calculation can not be done directly because of the complexity of the factors that influence it, therefore modeling methods are used to simplify this complexity (Komarudin, et al, 2015). The QUAL2Kw model is one method of calculating the pollutant load capacity by making a model in which this model is based on simplifying conditions in the field. (Sara, 2018). QUAL2Kw is a model adapted from the QUAL2K model developed by Dr. Steven C. Chapra from Tufts University (Chapra and Pelletier, 2003). QUAL2Kw is a surface water quality model, one-dimensional, steady state (more specifically a stable flow model because the flow is assumed to be stable but water quality is calculated dynamically with water quality kinetics and heat balance) and implemented in a Microsoft Windows environment (Kannel et al, 2011).

QUAL2Kw can simulate a number of water parameters such as temperature, pH, BOD, COD, dissolved oxygen, organic nitrogen, ammonia nitrogen, nitrite and nitrate nitrogen, organic phosphorus, inorganic phosphorus, total nitrogen, total phosphorus, total phytoplankton and bottom algae (Kannel et al, 2007).

RESEARCH METHOD

Location and Time of Study

This research was conducted in 6 months from February upto July 2019 in downstream segment of Citarum River starting from segment 16 located at Jl. Demang, Karawang Regency upto segment 20 in Muara Gembong, Bekasi Regency. Sampling is done twice at each location using grab sampling method for 3 observation periods (March, April and June 2019) at 11 water quality monitoring points consisting of 10 main river points and 1 tributary point. The parameter analyzed were BOD, COD, nitrate and phosphate.

Table 1. The Monitoring Locations of The Water Quality

PointLocation of MonitoringCoordinate Point
ELSL
1Jl. Demang, Karawang107°22'06,8"6°24'44,6"
2Jl. Rumambe I, Karawang107°19'28,3"6°21'45,1"
3Jembatan Teluk Jambe,
Karawang
107°18'43,3"6°19'36,8"
4Jembatan Tarumanagara,
Karawang
107°17'16,4"6°17'54,7"
ALabansari, Cikarang Timur,
Bekasi (point source)
107°15'54,1"6°16'56,2"
5Jl. Raya Proklamasi, Karawang107°16'34,7"6°14'57,0"
6Jl. Raya Pebayuran, Bekasi107°16'21,4"6°13'10.9"
7Jl. Raya Tugu Proklamasi,
Rengasdengklok, Karawang
107°17'14,1"6°09'24,9"
8Jl. Teluk Bango, Karawang107°12'31,7"6°05'52,1"
9Jl. Raya Pakisjaya, Karawang107°07'22,3"6°02'50,2"
10Muara Gembong, Bekasi107°02'38,7"5°59'04,5"

Pollutant Load Capacity Analysis

Analysis of pollutant load capacity in downstream segment of Citarum River was carried out using the QUAL2Kw model developed by Dr. Steven C. Chapra. Water quality modeling of downstream segment of Citarum River is divided into 9 reach from the point of 108.5 km at Jl. Demang, Karawang to the point of 0 km in Muara Gembong, Bekasi Regency. The distribution of reach is based on the physical condition of the river, the presence of pollutant sources both point source and non-point source that is suspected to have a major impact on changes in the quality and quantity of the downstream segment of the Citarum River can be seen in Table 2.

In this modeling, the shape of the downstream segment of Citarum River is considered square, the water quality parameters being modeled are BOD, COD, nitrate, and phosphate, input hydrological data comes from field measurements, input meteorological data, input the concentration of point source pollutant from tributaries, activities industry, hotel, shopping center, restaurant, and hospital sourced from the Environmental Agency, input the concentration of non-point source pollutants originate from agriculture and domestic waste of residents in Citarum watershed, water quality monitoring data in March, April and June 2019 are used as data for model verification. Model verification is done by calculating the percentage deviation of the simulation results with the results of field monitoring. Determination of the percentage of permissible deviations is 20% with a confidence level of the model of 80%.

Furthermore, the amount of the pollutant load capacity of the downstream segment of Citarum River is obtained by multiplying the flow rate of each segment by the water quality standard according to its designation. Meanwhile, the amount of pollutant load is obtained through the multiplication of river flow discharges with pollutant concentrations at each water quality monitoring location.

Table 2. Distribution of Reach in Downstream Segment of Citarum River

Reach NameLocationElevation (usl)Downstream Coordinate
NoUpstreamDownstreamReach
Length
of
Reach
Up
strea
m
Down
stream
SLEL
Initial ReachFinal Reach(Km)(Km)(m)(m)°'"°'"
0Headwater108,52062444,
6
107226,8
1HeadwaterJl. Rumambe13,994,6201962145,1071928,
Reach NameLocationElevation (usl)Downstream Coordinate
NoUpstreamDownstreamReach
Length
of
Reach
Up
strea
m
Down
stream
SLEL
Initial ReachFinal Reach(Km)(Km)(m)(m)°'"°'"
I, Karawang13
2Jl.
Rumambe
I,
Karawang
Jembatan
Teluk Jambe
5,684,0191861936,
8
1071843,
3
3Jembatan
Teluk
Jambe
Jembatan
Tarumanagar
a
10,674,8181661754,
7
1071716,
4
4Jembatan
Tarumanag
ara
Jl. Raya
Proklamasi,
Karawang
9,269,21613614571071634,
7
5Jl. Raya
Proklamasi
, Karawang
Jl. Raya
Pebayuran,
Bekasi
5,657,6131261310,
9
1071621,
4
6Jl. Raya
Pebayuran,
Bekasi
Jl. Raya
Tugu
Proklamasi,
Karawang
11,639,312116924,
9
1071714,
1
7Jl. Raya
Tugu
Proklamasi
, Karawang
Jl. Teluk
Bango,
Karawang
18,322,91176552,
1
1071231,
7
8Jl. Teluk
Bango,
Karawang
Jl. Raya
Pakisjaya,
Karawang
16,45,6746250,
2
107722,
3
9Jl. Raya
Pakisjaya,
Karawang
Muara
Gembong,
Bekasi
17,30415594,5107238,
7

RESULTS AND DISCUSSION

Pollutant Load of Point Source

Point source pollutants in the downstream segment of Citarum River comes from tributaries, industrial activities, offices, shopping centers, hotels, hospitals and hospitals. Calculation of point source pollutant load can be seen in Table 3.

Table 3. Pollutant Load of Point Source

Pollutant Load(kg/day)
ReachLocationBODCODNitratePhos
phate
1Jl. Demang,
Karawang
Jl. Rumambe
I, Karawang
3800.029899.540.111.80
2Jl. Rumambe I,
Karawang
Jembatan
Teluk Jambe,
Karawang
1027.974437.450.305.16
3Jembatan Teluk
Jambe,
Karawang
Jembatan
Tarumanagara,
Karawang
565.32800.940.111.91
4Jembatan
Tarumanagara,
Karawang
Jl. Raya
Proklamasi,
Karawang
15317.7481650.790.030.49
5Jl. Raya
Proklamasi,
Karawang
Jl. Raya
Pebayuran,
Bekasi
0.460.540.0010.01
6Jl. Raya
Pebayuran,
Bekasi
Jl. Raya Tugu
Proklamasi,
Rengas
dengklok,
Karawang
1.502.600.0010.02
7Jl. Raya Tugu
Proklamasi,
Rengasdengklok,
Karawang
Jl. Teluk
Bango,
Karawang
----
8Jl. Teluk Bango,
Karawang
Jl. Raya
Pakisjaya,
Karawang
----
9Jl. Raya
Pakisjaya,
Karawang
Muara
Gembong,
Bekasi
----

Based on Table 3., the highest pollutant load of point source was in reach 4, because in addition to the pollutant load originate from industrial activities, offices and hospitals, there is also a tributary (Cibeet River) which greatly contributes to the increase pollutant load of point source in that reach.

Pollutant Load of Non-Point Source

Pollutant load of non-point source in downstream segment of Citarum River comes from domestic activities and agricultural activities in watershed. Calculation of pollutant load of non-point source can be seen in Table 4.

Table 4. Pollutant Load of Non-Point Source

Pollutant Load
(kg/day)
ReachParameterDomesticAgricultureTotal
BOD6412.923.306416.22
COD8817.764.958822.71
1Total N-3.673.67
Nitrat1.60-1.60
Fosfat27.251.8329.09
BOD5671.722.765674.48
COD7798.624.137802.75
2Total N-3.063.06
Nitrat1.42-1.42
Fosfat24.101.5325.64
BOD5288.393.425291.81
COD7271.535.137276.66
3Total N-3.803.8
Nitrat1.32-1.32
Fosfat22.481.9024.38
BOD8318.213.798322.00
COD11437.545.6811443.22
4Total N-4.214.21
Nitrat2.08-2.08
Fosfat35.352.1137.46
BOD2699.562.112701.67
COD3711.903.173715.06
5Total N-2.342.34
Nitrat0.67-0.67
Fosfat11.471.1712.65
BOD5160.502.155162.65
COD7095.683.237098.91
6Total N-2.392.39
Nitrat1.29-1.29
Fosfat21.931.1923.13
Pollutant Load
(kg/day)
ReachParameterDomesticAgricultureTotal
BOD4266.524.004270.52
COD5866.466.005872.46
7Total N-4.444.44
Nitrat1.07-1.07
Fosfat18.132.2220.36
BOD6004.284.676008.94
COD8255.887.008262.88
8Total N-5.195.19
Nitrat1.50-1.50
Fosfat25.522.5928.11
BOD5893.374.435897.80
COD8103.386.658110.03
9Total N-4.924.92
Nitrat1.47-1.47
Fosfat25.052.4627.51

Based on Table 4, the highest pollutant load of non-point source comes from domestic in watershed. The largest pollutant load of non-point source was in reach 4, because in that reach the population density is quite high compared to other reach and there was also a large enough agricultural land so that runoff that is generated can increase the burden of pollutants entering the river.

The results of Existing Condition simulation using QUAL2Kw

Model simulations are carried out to represent the actual river conditions. The model that has been run will produce a simulation graph of the parameters modeled from the downstream segment of Citarum River.

Based on Figure 1, it can be seen that the results of simulation models for BOD almost all approach the results of analysis. Model simulation results show that BOD concentrations almost entirely exceed quality standards according to KEPGUB West Java No. 39 of 2000 for class B, C, D. The simulation results of the model are influenced by the increase in debit and levels of pollutants entering downstream segment of Citarum River. The high BOD concentration in the downstream segment of Citarum River is caused by domestic activities that dump waste directly into river

bodies. Industrial activities around the river also contribute to increasing the BOD level in the river.

Large industries located in the upstream and downstream of Citarum watershed in the form of the textile industry, the food processing industry, and the footwear/shoes manufacturing industry act as contributors to waste in Citarum River sub-watershed. In addition, micro and small business activities also increase the burden of waste entering Citarum River bodies (Kusuma et al, 2018). A high BOD value indicates a decrease in DO value because the oxygen available in the water is being consumed by bacteria, causing fish and other aquatic organisms to be unable to survive in the river (Pathak and Limaye, 2011).

2

Figure 1. The results of BOD Model Simulation in Downstream Segment of Citarum River

Model simulation results and analysis data for COD can be seen in Figure 2. COD concentration both model simulation results and the results of analysis data have

far exceeded the quality standard based on KEPGUB West Java No. 39 of 2000 for classes B, C, D which is 10 mg/L. Similar to BOD, high levels of COD are also caused by waste originate from domestic and industrial activities in watershed. The COD reflects the amount of chemically oxidized organic compound or the amount of oxygen needed to oxidize organic substances to CO2 and H2O. Certain organic materials present in wastewater are immune to biological degradation and some are toxic even at low concentrations. COD levels in wastewater will decrease along with the reduction in the concentration of organic matter contained in wastewater (Metcalf and Eddy, 2003).

1

Figure 2. The results of COD Model Simulation in Downstream Segment of Citarum River

Based on Figure 3, it can be seen that the results of the nitrate model simulation have approached the analysis data. Both the model simulation results and analysis data results of nitrate concentrations in downstream segment of Citarum River were still below the quality standard according to KEPGUB West Java No. 39 of 2000 for classes B, C, D (10 mg/L). Sources of nitrate pollutants in downstream segment of Citarum

40 Jurnal Teknik Lingkungan Vol.26 No. 2 Ayu Widya Utami, Pramiati Purwaningrum, dan Diana Irvindiaty Hendrawan.

River comes from agricultural activities, the use of fertilizers, and settlement activities. Nitrates derived from ammonium that enter the river body mainly through domestic waste concentrations in the river will decrease when getting further from the point of disposal due to the activity of microorganisms in the water for example bacteria nitrosumonas (Mustofa, 2015). Nitrate-nitrogen levels in natural waters are almost never more than 0.1 mg/L, but if nitrate levels are greater than 0.2 mg/L will lead to eutrophication which further stimulates the growth of algae and aquatic plants rapidly (Effendi, 2003).

1

Figure 3. The results of Nitrate Model Simulation in Downstream Segment of Citarum River

Based on Figure 4., the results of the phosphate model simulation almost all approach the analysis data. The model simulation results and analysis results of phosohate in downstream segment of Citarum River have far exceeded the quality standard for class 1 based on PP No. 82 of 2001 which is equal to 0.2 mg/L.

0

Figure 4. The results of Phosphate Model Simulation in Downstream Segment of Citarum River

High levels of phosphate in downstream segment of Citarum River occur due to agricultural runoff and domestic activity waste, especially human waste and laundry waste containing detergents that are discharged directly into river bodies. In agricultural areas orthophosphate is derived from fertilizer that enters rivers or lakes through drainage and rainwater flow. Polyphosphates can enter rivers through wastewater from residents and industries that use phosphate-containing detergents, such as the metal industry and so on. Organic phosphate is found in wastewater from residents (feces) and food scraps. Organic phosphate can also occur from orthophosphate which is dissolved through biological processes because both bacteria and plants absorb phosphate for its growth (Amirullah and Ismail, 2009).

Based on the model simulation results of the existing downstream segment of Citarum River, the percentage deviations from the BOD, COD, nitrate, and phosphate parameters were 5%, 8%, 13%, and 13%. Deviation values of less than 20% are 42 Jurnal Teknik Lingkungan Vol.26 No. 2 Ayu Widya Utami, Pramiati Purwaningrum, dan Diana Irvindiaty Hendrawan.

considered acceptable in predicting such concentrations in downstream segment of Citarum River and can be used to calculate river pollutant loads. The percentage of errors in modeling can occur due to water sampling activities at each monitoring point (Zhang et al, 2012). The model will be more accurate if adding various input variables including basic algae, sediment oxygen demand, organic nitrogen, total and organic phosphorus, etc. in the monitoring program (Hossain et al, 2014).

Pollutant Load Capacity Analysis

Determination of the pollutant load capacity is based on the value of the difference between the existing pollutant load in the river with the permissible pollutant load that is the concentration of each water quality parameter according to the water class in the quality standard. The difference can be negative or positive. If the difference between pollutant load and capacity was negative, it is necessary to reduce pollutant load. If the difference between pollutant load and capacity is positive, it can be said that the river is still able to accommodate the amount of loads before it exceeds its capacity.

Table 5. Pollutant Load Capacity

ReachParameterPollutant
Load
Capacity
(kg/day)
Pollutant
Load
(kg/day)
Difference
(kg/day)
Percentage of
Overload
(%)
BOD14556.6716012.34-1455.679%
0COD24261.1269556.63-45295.5165%
Nitrat24261.12800.6223460.500%
Fosfat485.223906.04-3420.8288%
BOD15709.7918664.66-2954.8716%
COD26182.9873820.25-47637.2765%
1Nitrat26182.98913.1825269.800%
Fosfat523.663234.86-2711.2084%
ReachParameterPollutant
Load
Capacity
(kg/day)
Pollutant
Load
(kg/day)
Difference
(kg/day)
Percentage of
Overload
(%)
BOD16596.5216974.06-377.542%
COD27660.8774935.95-47275.0763%
2Nitrat27660.871080.9626579.910%
Fosfat553.222852.50-2299.2881%
BOD17187.9117043.88144.040%
COD28646.5276917.13-48270.6163%
3Nitrat28646.521243.5427402.990%
Fosfat572.932658.53-2085.6078%
BOD20905.7116885.164020.560%
4COD34842.8688669.93-53827.0861%
Nitrat34842.861517.5833325.280%
Fosfat696.864623.32-3926.4685%
BOD21217.1117995.953221.160%
COD35361.8690502.14-55140.2961%
5Nitrat35361.861645.6133716.250%
Fosfat707.244690.38-3983.1485%
BOD21569.8623702.44-2132.589%
COD35949.7797865.45-61915.6863%
6Nitrat35949.771987.4533962.320%
Fosfat719.005129.15-4410.1586%
BOD22140.1028233.06-6092.9622%
7COD36900.17105702.47-68802.3065%
ReachParameterPollutant
Load
Capacity
(kg/day)
Pollutant
Load
(kg/day)
Difference
(kg/day)
Percentage of
Overload
(%)
Nitrat36900.172298.8634601.310%
Fosfat738.005445.75-4707.7586%
BOD22814.0226082.87-3268.8513%
COD38023.37114634.33-76610.9667%
8Nitrat38023.372438.5435584.830%
Fosfat760.475686.16-4925.6987%
BOD23456.8431866.22-8409.3926%
COD39094.73122959.31-83864.5868%
9Nitrat39094.732473.3736621.360%
Fosfat781.895784.25-5002.3686.5%

Based on Table 4., pollutant load of BOD, COD, and phosphate in all reach exceeds pollutant load capacity (PLC), while for nitrate pollutant load does not exceed its PLC. The percentage of BOD overload is 0%-26%, COD overload was 61%-68%, and phosphate overload was 78%-87%.

CONCLUSION

Simulation results of existing conditions for pollutant load of BOD, COD, nitrate, and phosphate in downstream segment of Citarum River almost all have exceeded its pollutant load capacity (PLC). The results of the PLC for BOD model simulation was 14556.67 kg/day-23456.84 kg/day while pollutant load of BOD was 16012.34 kg/day-23456.84 kg/day so that the BOD overload was 0%-6%. The results of PLC for COD model simulation was 24261.12 kg/day-39094.73 kg/day while pollutant load of COD was 69556.63 kg/day-122959.31kg/day so that the COD overload was 61%-68%. The result of PLC for nitrate model simulation was 24261.12 kg/day-39094.73 kg/day while pollutant load of nitrate was 800.62kg/day-2473.37kg/day so that the nitrate overload was 0%. The results of PLC for phosphate model simulation was 485.22 kg/day-781.89 kg/day while pollutant load of phosphate was 3906.04 kg/day-5784.25 kg/day so that the phosphate overload was 78% -87%.

Research Intelligence

Data from OpenAlex ↗

Metrics

2
Citations
0.18
FWCIfield-weighted
56th
Percentilevs same year + field
Article
Work type
Open Access

Citation Trend

Citation Timeline

YearCitations
20232

Institution Network

References

  1. A. Mustofa (2015). Kandungan Nitrat dan Pospat sebagai Faktor Tingkat Kesuburan Perairan Pantai. Jurnal DISPROTEK, 6(1) DOI: 10.34001/jdpt.v6i1.193
  2. Amirullah dan Ismail, M. (2009). Biokimia Kesehatan Edisi I. Pustaka As-Salam: Makassar.
  3. C. Patricia (2018). Analisis Daya Tampung Beban Pencemar Nitrat dan Fosfat di Sungai Ciliwung dengan Model QUAL2Kw. Tugas Akhir. Program Studi Teknik Lingkungan. Fakultas Arsitektur Lanskap dan Teknologi Lingkungan. Universitas Trisakti.
  4. H. Effendi (2003). Telaah Kualitas Air Bagi Pengelolaan Sumber Daya dan Lingkungan Perairan. Jakarta: Kanisius (Anggota IKAPI)
  5. H. Pathak and S. N. Limaye (2011). Interdependency between Physicochemical Water Pollution indicators: A Case Study of River Babus, Sagar, M. P., India, Analele UniverstaNii din Oradea-Seria Geografie (June) Article No. 211103-515, 23-29
  6. Kementerian Lingkungan Hidup dan Kehutanan Republik Indonesia. Alokasi Beban Pencemaran Sungai Citarum. (2016).
  7. M. Komarudin, S. Hariyadi and B.Kurniawan (2015). Analisis Daya Tampung Beban Pencemar Sungai Pesanggrahan (Segmen Kota Depok) dengan Menggunakan Model Numerik dan Spasial. Jurnal Pengelolaan Sumberdaya Alam dan Lingkungan, 5(2), 121-132. DOI: 10.19081/jpsl.5.2.121 DOI: 10.19081/jpsl.5.2.121
  8. M. A. Hossain, I. M. Sujaul, and M. A. Nasly (2014). Application of QUAL2Kw for Water Quality Modeling in The Tunggak River, Kuantan, Pahang, Malaysia. Research Journal of Recent Sciences, 3(6), 6-14
  9. M. S. B. Kusuma, L. M. Hutasoit, B. Kombaitan, M. Syafila, T. Setiadi, I. N. P. Aryantha, and Y. A. Piliang (2018). Strategi Pengelolaan Terpadu Penyelesaian Permasalahan DAS Citarum-Forum Guru Besar ITB 2018. Pidato Forum Guru Besar ITB. Forum Guru Besar Institut Teknologi Bandung, Bandung.
  10. Metcalf and Eddy (2003). Wastewater Engineering: Treatment and Reuse, 4th ed. Newyork: McGraw Hill Companies, Inc.
  11. P. S. Sara (2018). Analisis Daya Tampung Beban Pencemar BOD dan COD di Sungai Ciliwung dengan Model QUAL2Kw. Tugas Akhir. Program Studi Teknik Lingkungan. Fakultas Arsitektur Lanskap dan Teknologi Lingkungan. Universitas Trisakti.
  12. P. R. Kannel,. S. R. Kanel, S.Lee, Y-S.Lee, T.Y. Gan (2011). A Review of Public Domain Water Quality Models for Simulating Dissolved Oxygen in Rivers and Streams. Environ Model Assess, 16(2), 183-204. DOI: 10.1007/s10666-010-9235-1 DOI: 10.1007/s10666-010-9235-1
  13. P. R.Kannel, S. Lee, Y-S Lee, S. R Kanel and G. J. Pelletier (2007). Application of Automated QUAL2Kw for Water Quality Modeling and Management in the Bagmati River, Nepal. Ecological Modelling, 202, 503-517.
  14. R. Afroz, and A. Rahman (2017). Health Impact of River Water Pollution in Malaysia. International Journal of Advanced and Applied Sciences, 4(5), 78-85
  15. R. Zhang, X. Qian, X.Yuan, R. Ye, B. Xia and Y. Wang (2012). Simulation of Water Environmental Capacity and Pollution Load Reduction Using QUAL2K for Water Environmental Management. Int. J. Environ. Res. Public Health, 9, 4504-4521. DOI: 10.3390/ijerph9124504 DOI: 10.3390/ijerph9124504
  16. S. A. Alrumman, A. F. El-Kott and S. M. A. S. Keshk (2016). Water Pollution: Source & Treatment. American Journal of Environmental Engineering, 6(3), 88-98. DOI: 10.5923/j.ajee.20160603.02. DOI: 10.5923/j.ajee.20160603.02
  17. S.C. Chapra, and G.J. Pelletier (2003). QUAL2K: A Modeling Framework for Simulating River and Stream Water Quality (Beta Version): Documentation and Users Manual. Civil and Environmental Engineering Dept., Tufts University, Medford, MA.Churchill, M.A., Elmore, H.L., and Buckingham, R.A. 1962. The prediction of stream reaeration rates. J. Sanit. Engrg. Div., ASCE, 88(4),1-46
  18. T.P.Masere, A. Munodawafa, and T.Chitata (2012). Assessment of human impact on water quality along Manyame River. International Journal of Development and Sustainability, 1(3), 754-765