In wetlands nutrients mainly originated from anthropogenic factors through domestic waste and housing area development, with the exception inundated with the river flashing from the estuaries. Estuaries are known for its strong environmental gradients and complex physico-chemical process is the main factor for large spatial and temporal variations in the structure and dynamics of any biological communities in that area, including phytoplankton . Phytoplankton is susceptible to fluctuating the productivity therefore leading to variations in distributions and succession. Succession of phytoplankton was affected mostly by combination of biological, such as grazing and competition between species for physical and chemical quality aquatic ecosystem. Physical parameters like, hydrodynamics, rainfall, light and temperature, as well chemical factors, such as availability of dissolved inorganic nutrients, pH and dissolved oxygen in aquatic ecosystem key components in phytoplankton bloom . Phytoplankton community fluctuates due to uptake of high concentration of nutrients, especially nitrogen in the form of nitrate and phosphorus as phosphate are essential phytoplankton growth . Nitrate and ammonia represent as two main sources of dissolved nitrogen as nitrate is derived from anthropogenic origin  while ammonia originates from the discharge of wastewater . Inorganic nitrates and ammonia limit phytoplankton growth  and their ecological significance in rivers and streams has been well documented . There are many factors that affect the relative amount of nitrogen in various forms, which are influences with the discharge of wastewater into water, fertilization and types of wetland, reservoirs or lakes in the path of flowing water . Correlation between nitrogen and phosphorus and abundance further supported the former conclusion that phosphorus is the controlling factor in phytoplankton growth in the estuary where light is not limited  while nitrogen is commonly limited in marine ecosystems . Surface runoff after heavy precipitation initially resulted in increased accumulation of nutrients (nitrate and phosphate) concentration in the estuary. Nitrate is a limited factor for phytoplankton cells density rather than phosphate in Kota Kinabalu coastal water . Abundance and composition of phytoplankton in estuaries are the consequences of nutrients concentration (phosphate, nitrate, ammonia) and might also be the tides as high tides condition allows diatoms from seawater to enter the freshwater and mix with freshwater phytoplankton . On the other hand, the inundated area at KK Wetland is characterized by the environment’s specific properties such as low hydraulic gradient, poor mixing of water columns, with high probability of water retention. Due to the wetlands located in estuary areas, tidal fluctuation plays a vital role in changes in concentration of nutrients. Water flowing in the river channel at KK Wetland is also affected by inputs of high loads of nutrients. The inputs are mainly due to main activities, known as anthropogenic factors, such as domestic waste and development of housing area strongly modify biomass of phytoplankton abundance. Information on the distributed of dissolved inorganic nutrients and the community structure of phytoplankton in wet land is limited, especially characteristics like KK Wetland. This study was attempted to identify the distributed pattern of dissolved inorganic nutrients and community structure of phytoplankton alone the adjacent areas of KK Wetland. It will contribute the scientific data about the concentration of dissolved inorganic nutrients and provide checklist of phytoplankton diversity in Kota Kinabalu Wetland.
2. Materials and Methods
2.1. Description of Study Site
The study was conducted at Kota Kinabalu (KK) Wetland located in west coast of Sabah, Malaysia (Figure 1). Five sampling stations were selected adjacent to KK Wetland as shown in Table 1. Samples were taken for period of period of six months (May 2019 until October 2019).
2.2. Collection of Water Samples for Dissolved Inorganic Nutrients
The water samples from the selected stations were collected by using Van Dorm
Figure 1. Location of sampling stations in Kota Kinabalu Wetland (Sources: QGIS, 2019).
Table 1. Sampling stations with attributions and coordinates.
water sampler within one metre depth in each station. Water sample was then transferred into polyethylene bottles and labelled according to stations. The analysis of dissolved inorganic nutrients was conducted in oceanography laboratory at Borneo Marine Research Institute (BMRI).
2.3. Collection of Phytoplankton for Identification of Species and Cell Density
Phytoplankton net of 20 µm mesh size was towed for 1 minute at the water surface against movement of tide for every station. Phytoplankton was concentrated at the bottom of nets and transferred into polyethylene bottles. They were preserved immediately with Lugol’s iodine. Identification of phytoplankton species was executed thoroughly at Borneo Marine Research Institute (BMRI) based on guidance from .
2.4. Analytical Parameters
2.4.1. Collection of In-Situ Parameters
Water parameters such as salinity (ppt), temperature (˚C), level of DO (mg/L) and pH were measured using YSI multi-parameter probe (Model: Pro Plus, GeoHydrOcean Services) at depths around 1 - 1.5 m.
2.4.2. Laboratory Analysis for Ex-Situ Parameters
The water samples were filtered with 0.45 µm pore-size membrane filter for the determination of nitrate, ammonia and phosphate (mg/L). Nutrients were analyzed with spectrophotometric methods proposed by the American Public Health Association .
2.4.3. Cell Density (Cells/mL)
Phytoplankton cell count was carried out using Sedgewick-Rafter counting chamber under light microscope at 40× magnification. The total number of phytoplankton was counted according to formula by Stirling :
N = Number of plankton cells per litre of original water
A = Total of counted phytoplankton
C = Volume of final concentrate of samples (mL)
V = Volume of field (m3)
F = Number of counted fields
L = Volume of original water (L)
2.4.4. Statistical Analysis
Statistical analyses like one-way ANOVA and Pearson correlation were conducted using the SPSS (Version 19.0). Significant level of p < 0.05 was tested for normality and homogeneity of variances. Diversity indices, such as, species diversity, equitability index and species richness of phytoplankton community were calculated by using PAST software.
3. Results and Discussion
3.1. Spatial In-Situ Parameters
The average salinity during the study period among the stations were recorded in the range of 28 - 32 ppt while the mean of 32 ± 2.10 ppt and the lowest of 28 ± 6.02 ppt observed in the station S1 and S3 respectively. The temperatures were recorded in the range of 30.00˚C - 31.20˚C, fluctuates in less than 1˚C which indicates all five stations have almost similar surface water temperature. The highest dissolved oxygen of 4.46 mg/L was found at the river mouth (S1) while the lowest of 2.68 mg/L was determined at station S4 in the inundated area located inside KK wetland. The value of pH at inundated area of estuary, stations S4 and S5 (7.37 and 7.35) were found lower compared to the pH value of 7.87 in the river mouth, station S1 (Table 2). Significant difference (P = 0.008) was observed in salinity, temperature, dissolved oxygen and values of pH in when compared between months around KK wetland. Temperature, salinity and pH shared similar fluctuating trend throughout the six months while dissolved oxygen fluctuated drastically from June to July (the lowest) and with the increase in the month of August to September (Figure 2).
3.2. Spatial Nutrients Concentration
Concentration (mg/L) of nitrate observed had the lowest values compared to values in concentration determined in ammonia and phosphate. The concentration
Table 2. Values (mean ± SD) of in-situ water parameters (salinity, temperature, dissolved oxygen and pH) collected at the five sampling stations.
of nitrate was observed from 0.005 - 0.012 mg/L (Table 3). The highest concentration of 0.0115 mg/L nitrate was recorded at station S5 while the lowest concentration of 0.0047 mg/L nitrate was found at station S3. However, the concentrations of ammonia were relatively higher at every stations compared to other DIN (Figure 2), especially concentrations of ammonia of 0.2004 mg/L and 2.311 mg/L recorded at the river (stations S2 and S3). The lowest ammonia concentration of 0.0159 mg/L was recorded at station S4 (Table 3).
The concentration of phosphate was found within the range of 0.0089 - 0.0513 mg/L. The lowest of 0.0089 mg/L and the highest of 0.0513 mg/L concentrations of phosphate were determined in station S1 and S4 respectively (Figure 3).
In general, nitrate, ammonia and phosphate showed no significant difference (P = 0.737) in terms among the sampling stations during the study period, except ammonia from station S3. Ammonia and phosphate levels were observed the highest during the month of July and September respectively. Nitrate concentration was observed consistently similar throughout the study period (Figure 3).
The in-situ water parameters such as salinity (ppt), temperature (˚C), dissolved oxygen (mg/L) and pH can greatly affect the dynamics of phytoplankton growth and abundance in the water . There were no significant differences in spatial distribution among the salinity (P = 0.1613), temperature (P = 0.1441),
Figure 2. Trends in Salinity, temperature, dissolved oxygen and pH values from May 2019 to October 2019 (mean ± sd).
Table 3. Nutrients concentration (mean ± SD) in Kota Kinabalu Wetland.
Figure 3. Trends of DIN (Nitrate, Ammonia and Phosphate) concentrations from May 2019 to October 2019 at Kota Kinabalu Wetland area.
dissolved oxygen (P = 0.1128) and values of pH (P = 0.0782) which were recorded in five sampling stations around Kota Kinabalu Wetlands. This might be due to the limitation of study area that were covered in the study, as well as fewer sampling stations which led to less sample variability in analysis. This was the same condition observed for in-situ water parameters between months during the dry season period . In this study, the highest salinity was observed during the pre-monsoon (May) while lowest was at July and August, during the dry period of southwest monsoon, although high salinity was expected in less rainfall period during the dry season  . Dissolved oxygen also plummeted incredibly from 4.34 mg/L on May to 1.34 mg/L on July. This was probably due to the environmental condition as the lowest tide during low tides occurred on July, which promotes decrease in salinity, dissolved oxygen and even values in pH. However, the tidal data was recorded during the study period. Similar conditions was observed recently in tidal lake of Indonesia where the salinity was lower during the lowest tide because saltwater intrusion from the sea did not occur while for dissolved oxygen, diffusion process was restricted hence limiting oxygen supply in water . The same study also suggested the sudden low value of pH level in dry season of low tide was due to the intense decaying process of organic matter as a result of high intensity and frequent exposure of water to sunlight. Salinity and dissolved oxygen concentrations were the lowest in July during the dry season in Marudu Bay . The present study showed that higher salinity was observed at the river mouth near to the bay area where similar finding was found in Brunei Bay during the dry season due to the mixing of seawater . The temperatures were observed similar between sampling stations and months throughout the period of study which was probably due time of sampling therefore dependent to presence of sunlight and heating effect.
Since the samplings were done during the dry season, the concentrations of DIN were not affected by La Nina period. The concentration of nitrate was relatively similar and low compared to ammonia and phosphate in every sampling station and month during the present study. Low concentrations of nitrate were due to longer residence time, thus increases the chances of decoupling between tidal cycle and residual drainage . Decomposition process is known to be higher during the dry season, including the period of residence time . Concentrations of ammonia were the highest at the river channel (S2 and S3) with not much of differences in mean concentration around estuary inundated area (S4 and S5) (Table 3). Across the six months of study, ammonia concentration was the highest in July. Ammonia concentrations were known to be higher than nitrate in estuary and mid-estuary area as a result of regeneration of ammonia from mangrove sediments   with low precipitation of dry season also contributed to the sudden peak of ammonia . The spike increase of ammonia concentration in July was highly due to wastewater effluents from municipal area and use of fertilizer from agricultural activity nearby Kota Kinabalu Wetland area associated with lowest precipitation of the dry period . Similar to trend in ammonia, phosphate mean concentration was observed constant in stations located at the river channel and inundated area of estuary, slightly higher than concentrations at river mouth. Overall phosphate concentration was highest in July. The present average phosphate concentration (0.01 - 0.03 mg/L) is almost similar to Douala estuary  where the concentrations were ranged from 0.03 to 0.04 mg/L in dry season, thus indicating phosphate concentrations were generally low in dry season compared to other monsoon seasons. Phosphate concentrations around Kota Kinabalu wetlands is high around estuary and river channel passing through the housing areas, indicating phosphate levels is induced by sewage and wastewater nearby housing areas . The nutrients level in the present study were lower than the Malaysia NWQS  acceptable limits of 7 mg/L, 0.1 mg/L and 0.2 mg/L for nitrate, ammonia and phosphate respectively which represents the waters around Kota Kinabalu Wetland is suitable for aquatic wildlife.
3.3. Abundance and Diversity of Phytoplankton
A total of 24 genera form 19 families were found in Kota Kinabalu Wetland waters during the period of study. Among them seven were from dinoflagellates (24.95%), 14 from diatoms (55.29%) and three genera, each from Spirotrichea (5.28%), Cyanobacteria (3.33%) and Crysophyta (11.15%). Diatoms were mainly dominated by Chaetoceros sp., Pseduo-nitzschia sp., and Cylindrotheca closterium throughout the study period while Nitzschia sp., Pleurosigma sp., Coscinodiscus sp., Thallassionema sp. and Skeletonema sp. occurred abundantly in certain months from June to September. Dinoflagellates were dominated by only Peridinium quinquecorne and Alexandrium sp. from May until October while massive numbers of Gonyaulax sp. was only observed on August. Cyanobacteria and Codonella sp. of Spirotrichea occurred consistently in small density throughout the June to September.
The highest density of 11.4 × 104 cells/L was recorded at station S1 while the lowest density of 6.29 × 104 cells/L phytoplankton was found at station S3 (Table 4). No significant (P = 0.284) differences were determined with the abundance of phytoplankton among the sampling stations during the study period. Diatoms were found to be dominating every station with the highest percentage of 64.33% at station S1 and of 65.58% at station S2. Abundance of dinoflagellates with highest percentage of 37.56% were observed at station S5 and of 34.48% at station S4. Phytoplankton abundance (cells/L) was the lowest during pre and post-monsoon period (May and October respectively) with a total of 2.685 × 104 cells/L and 3.660× 104 cells/L respectively, while the highest density of 12.115 × 104 cells/L and 10.995 × 104 cells/L of phytoplankton were observed during the dry season, July and August respectively (Table 4).
The Simpson’s Diversity Index (D) for phytoplankton was found to be the highest (0.912) at station S2 and the lowest (0.886) at station S5 (Table 5). However, species evenness is relatively low compared to the diversity as the ranges of the evenness index (J’) is from 0.283 to 0.291, with the lowest inside Kota Kinabalu wetlands inundated area (S5) and highest at the river (S2). Margalef Richness Index (d) was found the highest at station S3 while the lowest was at station S1 with the values of 3.082 and 2.843 respectively (Table 5).
Phytoplankton communities are generally represents as an indicator in determining the changes of nutrients in water and as a vital component for evaluating eutrophication in marine ecosystems . The abundance of phytoplankton was found to be the most abundant (11.4560 × 104 cells/L) at river mouth (S1), as located at the mouth of the Likas Bay. This might be due to the lower turbidity level and mixing of saline and freshwater. Phytoplankton abundance in
Table 4. Phytoplankton abundance (cells/L) in respective months (May 2019 to October 2019).
Table 5. Diversity indices of phytoplankton in Kota Kinabalu Wetland at five sampling sites.
Brahmani estuary  and Perak estuary was observed to be increased with elevating salinity level . The lowest phytoplankton abundance (6.2900 × 104 cells/L) was located at the river channel inside estuary (station S3) of Kota Kinabalu Wetland. The limited abundance of phytoplankton in S3 was due to high turbidity which limits photosynthetic rate for primary productivity. Most estuaries have vast amount of mineral sediments as a result of input from runoffs and suspended by tidal and wave fluctuations . Since the average depth of waters around the study area was within 1 - 1.5 m, process of turbidity driven by mixing of sediments assumed to be high which narrows the photic zone therefore limiting light penetrating into the water column and reducing phytoplankton biomass. However, the fluvial discharges of sediment were not accounted in this study.
Other factor leading to low abundance in S3 is the location was far away from source of nutrients inputs. Phytoplankton abundance was also the highest during July of dry season. There observed huge difference in abundance of phytoplankton between July (12.115 × 104 cells/L) and May (2.685 × 104 cells/L). The phytoplankton abundance was the highest during drier season and the lowest during the inter-monsoon month of May  . The low abundance in May was due to slow current during the inter-monsoon hence limits mixing of water .
Phytoplankton compositions in Kota Kinabalu Wetlands were mainly comprised of diatoms (55.29%), dinoflagellates (24.95%) and cyanobacteria. The species of the Diatom was accounted for 77% of total phytoplankton in marine ecosystem. The diatoms and dinoflagellates were the major dominated group among the phytoplankton in estuaries and marine ecosystem  . Diatoms were observed abundant during the dry season (June, August, and September), pre and post southwest monsoon (May and October). However, densities of diatoms dropped staggeringly in July as the sampling was done during the low tide. The species of Crysophytes were observed as the dominating phytoplankton community at that specific month in low tide estuaries . Dominating dinoflagellates (Peridinium quinquecorne and Alexandrium sp.) which were found in high density in estuary in Kota Kinabalu Wetland during the low tide was due to the tolerance to low salinity and warmer environment . Diatoms such as Chaetoceros sp., Pseduo-nitzschia sp., and Cylindrotheca closterium were most abundant throughout the period of study where these genera represent groups of marine diatoms and can be found in estuary during high tides as the seawater flows into the river and estuary. Nutrient concentrations are not the only factor affecting the phytoplankton blooms because changing in environmental parameter might favor some species to blooms . The environmental parameters such as temperature, salinity, Secchi depth (visibility), pH, and dissolved oxygen are the factors that can affect the algal cells community succession or limit the growth of certain groups of algae . In addition enrichment in nitrogen particularly nitrate lead to an increase in phytoplankton biomass in marine habitat . Tit was documented that salinity is an important factor that influence HABs in Kota Belud’s water other than nutrients concentration . The intrusion of seawater in Kota Kinabalu Wetland might control the composition of phytoplankton community. Pseduo-nitzschia sp., which is also a potential harmful phytoplankton, is capable of producing domoic acid . In Sabah, only four Pseduo-nitzschia sp. had been recorded which are P. pungens, P. calliantha, P. pseudodelicatissima, and P. delicatissima . Genus from diatoms (Chaetoceros, Cylindrotheca, and Pseudo-nitzschia) and dinoflagellates (Peridinium and Alexandrium) persist throughout the entire southwest, pre and post monsoon where their persistence could be stemmed from vast distribution and tolerance to brackish estuarine environmental conditions  . Overall, the total numbers of genera in Kota Kinabalu Wetland (24 genera) were observed relatively smaller compared to studies done by other researchers in tropical estuaries. A total of 37 genera in Philippine mangrove estuary  and 30 genera in Sungai Brunei estuary  indicated the lower diversity of phytoplankton at Kota Kinabalu Wetland.
The phytoplankton abundance in Kota Kinabalu wetlands mainly depends on the nutrient inputs in the estuary with other anthropogenic sources outside of Kota Kinabalu wetlands. The dry season in Kota Kinabalu Wetland resulted in lower concentrations of dissolved inorganic nutrients (nitrate, ammonia and phosphate) while environmental in-situ water parameters were almost similar and did not greatly affect the phytoplankton growth. In Kota Kinabalu Wetland, a total of 24 genera of identified phytoplankton were confirmed with the highest occurrence of dominating phytoplankton diatoms. Among the nutrients nitrate and phosphate are observed with increasing numbers of phytoplankton density. The findings obtained from this study can be used as a guideline and basis for future research, especially during the wet season. Continuous monitoring of environment parameters and succession phytoplankton in different seasons may fulfil the overall goal to achieve sustainable healthy ecosystem of Kota Kinabalu Wetland, which initially cannot be conclude with short period of study.
 Hecky, R.E. and Kilham, P. (1988) Nutrient Limitation of Phytoplankton in Freshwater and Marine Environments: A Review of Recent Evidence on the Effects of Enrichment. Limnology and Oceanography, 33, 796-822.
 Bradley, P.B., Sanderson, M.P., Frischer, M.E., Brofft, J., Booth, M.G., Kerkhof, L.J. and Bronk, D.A. (2010) Inorganic and Organic Nitrogen Uptake by Phytoplankton and Heterotrophic Bacteria in the Stratified Mid-Atalntic Bight. Estuarine, Coastal, and Shelf Science, 88, 429-441.
 Fonge, A.B., Chuyong, B.G., Tening, A.S., Fobid, A.C. and Numbisi, N.F. (2013) Seasonal Occurrence, Distribution and Diversity of Phytoplankton in Douala Estuary, Cameroon. African Journal of Aquatic Science, 38, 123-133.
 Huang, L.M., Jian, W.J., Song, X.Y., Huang, X.P., Liu, S., Qian, P., Yin, K. and Wu, M. (2004) Species Diversity and Distribution for Phytoplankton of the Pearl River Estuary during Rainy and Dry Seasons. Marine Pollution Bulletin, 49, 588-596.
 Anderson, D.M., Gilbert, P.M. and Burkholder, J.M. (2002) Harmful Algal Blooms and Eutrophication: Nutrient Sources, Composition and Consequences. Estuaries, 25, 704-726.
 Azad, S.A., Shaleh, S.R.M. and Soon, T.K. (2016) Temporal and Spatial Distribution of Nutrients and HABs at Coastal Water of Kota Belud, Sabah. Advances in Bioscience and Biotechnology, 7, 233-242.
 Hastuti, A.W., Pancawati, Y. and Surana, I.N. (2018) The Abundance and Spatial Distribution of Plankton Communities in Perancak Estuary, Bali. Earth and Environmental Science, 176, Article ID: 012042.
 Veronica, E., Leksono, A.S., Soemarno and Arfiati, D. (2014) Effect of Water Quality on Phytoplankton Abundance in Hampalan River and Fish Pond of Batajang Village. IOSR Journal of Environmental Science, Toxicology and Food Technology, 8, 15-21.
 Ouyang, Y., Nkedi-kizza, P., Wu, Q.T., Shinde, D. and Huang, C.H. (2006) Assessment of Seasonal Variations in Surface Water Quality. Water Research, 40, 3800-3810.
 Markou, D.A., Sylaois, G.K., Tsihrintzis, V.A., Gikas, G.D. and Haralambidou, K. (2007) Water Quality of Vistonis Lagoon, Northern, Greece: Seasonal Variation and Impact of Bottom Sediments. Desalination, 210, 83-97.
 Leidonald, R., Muhtadi, A., Lesmana, I., Harahap, Z.A. and Rahmadya, A. (2019) Profiles Temperature, Salinity, Dissolved Oxygen and pH in Tidal Lakes. Earth and Environmental Sciences, 260, Article ID: 012075.
 Soon, T.K. and Ransangan, J. (2016) Effects of Environmental Conditions and Nutrients on the Occurrence and Distribution of Potentially Harmful Phytoplankton in Mesotrophic Water. Sains Malaysiana, 45, 865-877.
 Yee, Y.H. and Suratman, S. (2016) Physico-Chemical Parameters Profile during Dry and Wet Seasons in Southern South China Sea: Brunei Bay. Asian Journal of Chemistry, 28, 2146-2152.
 Rivera-Monroy, V.H., de Mustert, K., Twilley, R.R., Castaneda-Moya, E., Romigh, M.M. and Davis, S.E. (2007) Patterns of Nutrient Exchange in a Riverine Mangrove Forest in Shark River Estuary, Florida, USA. Hidrobiologica, 17, 169-178.
 Fernandes, S.O., Eyre, B.D. and Gay, J.M. (2010) Denitrification: An Important Pathway for Nitrous Oxide Production in Tropical Mangrove Sediments (Goa, India). Journal of Environmental Quality, 39, 1507-1516.
 Smith, J., Burford, M.A., Revill, A.T., Haese, R.R. and Fortune, J. (2011) Effect of Nutrient Loading on Biogeochemical Processes in Tropical Tidal Creeks. Biogeochemistry, 108, 359-380.
 Lomoljo, R.M., Ismail, A. and Yap, C.K. (2009) Nitrate, Ammonia and Phosphate Concentrations in the Surface Water of Kuala Gula Bird Sanctuary, West Coast of Peninsular Malaysia. Pertanika Journal of Tropical Agricultural Science, 32, 1-5.
 Reed, M.L, Pinckney, J.L., Keppler, C.J., Brock, L.M., Hogan, S.B. and Greenfield, D.I. (2016) The Influence of Nitrogen and Phosphorus on Phytoplankton Growth and Assemblage Composition in Four Coastal, Southeastern USA Systems. Estuarine, Coastal and Shelf Science, 177, 71-82.
 Garmendia, M., Borja, A., Franco, J. and Revilla, M. (2013) Phytoplankton Composition Indicators for the Assessment of Eutrophication in Marine Waters: Present State and Challenges within the European Diversities. Marine Pollution Bulletin, 66, 7-16.
 Palleyi, S., Kar, R.N. and Panda, C.R. (2008) Seasonal Variability of Phytoplankton Population in the Brahmani Estuary of Orissa, India. Journal of Applied Sciences and Environmental Management, 12, 19-23.
 Nursuhuyati, A.S., Yusoff, F.M. and Shariff, M. (2013) Spatial and Temporal Distribution of Phytoplankton in Perak Estuary, Malaysia during Monsoon Season. Journal of Fisheries and Aquatic Sciences, 8, 480-493.
 Cloern, J.E., Foster, S.Q. and Kleckner, A.E. (2014) Phytoplankton Primary Production in the World’s Estuarine-Coastal Ecosystems. Biogeosciences, 11, 2477-2501.
 Sidabutar, T., Bengen, D.G., Wouthuyzen, S. and Partono, T. (2016) The Abundance of Phytoplankton and Its Relationship to the N/P Ratio in Jakarta Bay, Indonesia. Biodiversitas, 17, 673-678.
 Akhir, M.F.M., Zakaria, N.Z. and Tangang, F. (2014) Inter-Monsoon Variation of Physical Characteristics and Current Circulation along the East Coast of Peninsular Malaysia. International Journal of Oceanography, 2014, Article ID: 527587.
 Majewska, R., Adam, A., Mohammad-Noor, N., Convey, P., Stefano, M.D. and Marshall, D.J. (2017) Spatio-Temporal Variation in Phytoplankton Communities along a Salinity and pH Gradient in a Tropical Estuary (Brunei, Borneo, South East Asia). Tropical Ecology, 58, 251-269.
 Davies, O.A. and Ugwumba, O.A. (2013) Tidal Influence on Nutrients Status and Phytoplankton Population on Okpoka Creek, Upper Bonny Estuary, Nigeria. Journal of Marine Biology, 2013, Article ID: 684739.
 María Trigueros, J. and Orive, E. (2001) Seasonal Variations of Diatoms and Dinoflagellates in a Shallow, Temperate Estuary, with Emphasis on Neritic Assemblages. Hydrobiologia, 444, 119-133.
 Bates, S.S., Hubbard, K.A., Lundholm, N., Montresor, M. and Leaw, C.P. (2018) Pseudo-Nitzschia, Nitzschia, and Domoic Acid: New Research Since 2011. Harmful Algae, 79, 3-43.
 Teng, S.T., Leaw, C.P., Lim, H.C. and Lim, P.T. (2013) The Genus Pseudo-Nitzschia (Bacillariophyceae) in Malaysia, Including New Records and a Key to Species Inferred from Morphology-Based Phylogeny. Botanica Marina, 56, 375-398.
 Snoejis, P. and Weckstrom, K. (2010) Diatoms and Environmental Change in Large Brackish-Water Ecosystems. In: Smol, J.P. and Stoermer, E.F., Eds., The Diatoms: Applications for the Environmental and Earth Sciences, Cambridge University Press, Cambridge, 287-308.
 Aké-Castillo, J.A. and Vázquez, G. (2011) Peridinium quinquecorne var. trispiniferum var. nov. (Dinophyceae) from a Brackish Water Environment. Acta Botanica Mexicana, 94, 125-140.