Waters Of Taiwan Research Articles

Taxonomic review and molecular phylogenetics of the reef gobiid genus Fusigobius from Taiwan with comments on a new species from South China Sea

A taxonomic review is given for the gobiid genus Fusigobius, known for its cryptic and inconspicuous body coloration, from the waters of Taiwan. Six species of Fusigobius have been formally confirmed from Taiwan, including four nominal species (F. duospilus, F. humeralis, F. inframaculatus, and F. neophytus), a newly recorded F. melacron, and a new species, Fusigobius taipinensis n. sp., recently collected from the South China Sea and southern Taiwan. The new species is morphologically close to F. duospilus, but can be clearly distinguished by several important coloration features on the first dorsal fin or body. Phylogenetic analysis and species delimitation inferred from partial sequences of the mtDNA ND5 gene support F. taipinensis as a valid discrete species and F. melacron as the plesiomorphic state for Fusigobius.

First record of Matsubaraea fusiformis (Perciformes: Percophidae) in the waters of Western Taiwan: A taxonomic description

Herein, we present a new record of Matsubaraea fusiformis in western Taiwan waters; the discovery was based on its morphological characters. Five M. fusiformis samples were collected (August 2022) from the offshore waters (depth: 35.3–49.6 m) of Changhua County, western Taiwan, by using commercial bottom trawlers. A total of 5 genera and 10 species of the family Percophidae have been found in the waters of southwestern and northeastern Taiwan, primarily in the sandy deep-sea regions. This paper presents the first record of M. fusiformis in Taiwan. The findings also reveal that M. fusiformis is predominantly found in shallow waters.

Effects of oceanographic conditions on fishery distribution: A case study of chub mackerel (Scomber japonicus) in northeastern Taiwan

Changes in oceanographic conditions can affect species distribution in marine habitats. Global climate change may negatively influence the oceanographic factor–species distribution relationship. Here, we assessed the influence of oceanographic conditions on chub mackerel (Scomber japonicus) distribution in northeastern Taiwan by constructing and using a habitat ensemble model incorporating chub mackerel fishery, climatic oscillation, and oceanography data. Our results indicated that the chub mackerel catch was mainly influenced by the Western Pacific Oscillation. Moreover, sea-surface height and mixed-layer depth exerted the most and least significant effects on chub mackerel distribution, respectively. The chub mackerel catch rate peaked in the study area with a sea-surface temperature of 29 °C, sea-surface chlorophyll of 0.25 mg/m3, sea-surface salinity of 33.7 psμ, and SSH of 0.575 m. Chub mackerel was the most widely distributed in the area between 25°N, 120.5°E and 26.2°N, 121.5°E. Our findings can be used to develop critical adaptation plans for managing chub mackerel fisheries in the northeastern waters of Taiwan. Considering changing climate conditions globally, the incorporation of this knowledge into managerial strategies may aid decision-makers in protecting not only other ocean fisheries but also individuals dependent on them.

Modeling and performance estimation for L-shaped OWC wave energy converters with a theoretical correction for spring-like air compressibility

The investigation of Oscillating Water Columns (OWC) has gained significant attention in recent years thanks to their resilient structural design, allowing them to withstand harsh environmental conditions. This study focuses on the L-shaped OWC (L-OWC) due to its excellent energy capture efficiency. It should be noted that air compressibility plays a significant role in the plenum chamber of the OWC, especially at full scale. The present paper proposes a scaling-rematched approach, facilitating the evaluation of hydrodynamic coefficients and the performance estimation with a theoretical correction for the unmatched scaling problem of spring-like air compressibility. The methodology is applied successfully to a model-scale L-OWC design, employing three-dimensional, incompressible-flow simulations combined with an impeller model. The present study reveals the hydrodynamic characteristics (advantages) of the L-OWC: small fluid damping coefficient and large added mass, and hence the possibility that the spring-like air compressibility can be used to raise the efficiency of power capturing. Furthermore, there exists an interval where the air compressibility has a positive effect on the performance, not only having a significantly longer span than that of the conventional OWC but also much more directly matching the period range of high energy potential found in the wave climate of Northeastern Taiwan waters.

A new species of Pseudodiaptomus (Copepoda, Calanoida, Pseudodiaptomidae) from the coastal waters of northern Taiwan with notes on co-occurring and previously recorded species

Abstract We investigated the species diversity of the demersal copepod genus Pseudodiaptomus in the coastal waters of northern Taiwan, and recorded four species new to Taiwanese waters: P. echinatus n. sp., P. ishigakiensis Nishida, 1985, P. philippinensis Walter, 1986, and P. zhengi Wang & Guo, 2019. Pseudodiaptomus echinatus n. sp. is distinguished from the allied species P. serricaudatus (Scott T., 1894) by the female genital double-somite, which is laterally swollen and of subequal length and width, and the male fifth leg, the right endopod of which is oval with a row of spinules increasing in size toward the tip, and the first exopod segment of which has a distomedial digitiform process with rounded tip. In addition, out of the 11 nominal species previously recorded from Taiwanese waters, we found only three species (P. annandalei Sewell, 1919, P. bispinosus Walter, 1984, P. trihamatus Wright, 1937). We were unable to confirm the presence of P. forbesi (Poppe & Richard, 1890), P. incisus Shen & Lee, 1963, P. inopinus Burckhardt, 1913, P. marinus Sato, 1913, P. pacificus Walter, 1986, P. penicillus Li & Huang, 1984, P. poplesia (Shen, 1955), and P. serricaudatus (T. Scott, 1894) by a critical review of available literature based on morphology, locality, habitat, etc. In this regard, it is necessary for a future study to try to verify the presence of these eight species in Taiwanese waters.

Development and Application of an Advanced Automatic Identification System (AIS)-Based Ship Trajectory Extraction Framework for Maritime Traffic Analysis

This study addresses the challenges of maritime traffic management in the western waters of Taiwan, a region characterized by substantial commercial shipping activity and ongoing environmental development. Using 2023 Automatic Identification System (AIS) data, this study develops a robust feature extraction framework involving data cleaning, anomaly trajectory point detection, trajectory compression, and advanced processing techniques. Dynamic Time Warping (DTW) and the Hierarchical Density-Based Spatial Clustering of Applications with Noise (HDBSCAN) algorithms are applied to cluster the trajectory data, revealing 16 distinct maritime traffic patterns, key navigation routes, and intersections. The findings provide fresh perspectives on analyzing maritime traffic, identifying high-risk areas, and informing safety and spatial planning. In practical applications, the results help navigators optimize route planning, improve resource allocation for maritime authorities, and inform the development of infrastructure and navigational aids. Furthermore, these outcomes are essential for detecting abnormal ship behavior, and they highlight the potential of route extraction in maritime surveillance.

Captive spawning, early development, and larviculture of largescale blackfish, Girella punctata (Gray, 1835), and effects of macroalgae and alternative protein sources on its juvenile growth performance

The largescale blackfish, Girella punctata (Gray, 1835), inhabits rocky substrates in the coastal waters of Japan and northeastern Taiwan. Despite its commercial significance, there is limited knowledge regarding its reproductive biology due to the sparse occurrence of broodstocks in its natural environment. This study investigated the reproductive performance and embryonic development of G. punctata, assesses its feeding preferences through controlled feeding trials, and evaluates the impact of various feed types on its growth performance. In total, 124 wild-caught G. punctata were reared for 3–4 years, and their reproductive performance was assessed in 2015 and 2016. Their egg production increased with a sudden water temperature change from 17 ℃ in 2015–22 ℃ and 2016. The total and daily average egg production weights were 3589.22 and 40.79 g, respectively, whereas the total and daily average fertilized egg production weights were 1630.72 and 18.59 g, respectively; moreover, the hatching rate was 42.61 %. During embryonic development, the eggs hatched 45 h post fertilization at 19 ℃. Newly hatched larvae were 2.26 ± 0.04 mm in length with a yolk sac of 1.09 ± 0.02 mm in length. The yolk sac was almost completely absorbed 3 days post hatch (dph), when the larvae grew to a length of 3.42 ± 0.06 mm. We found that fish meal with 40 % protein was the optimal G. punctata fish diet. The weight gain in juveniles fed Cladophora fascicularis ranged from 47 % to 52 %, which was significantly higher than that in those fed Hypnea charoides (p < 0.05). The present study provided a practical guideline on the conditions required for successful fertilization and early development of G. punctata embryos and other larval stages. This information should help ensure greater survival and a higher proportion of healthy G. punctata larvae for mariculture applications.

Establishment and analysis of online energy baseline for water purification plants: A case study in Kaohsiung city, Taiwan

Water and energy are closely linked and are crucial for national security and economic development. Most water providers prioritise the stability of water supply and aim to reduce energy consumption under the premise of a stable supply. The average energy required to supply water in Taiwan in one of the lowest worldwide. In the Kaohsiung area, the average energy used by a water purification plant to provide 1 m3 of water is 0.32 kWh/m3, lower than the world average of 0.37 kWh/m3. However, the most energy-consuming plant (Weng Park water purification plant) uses eight times as much energy as the least energy-consuming plant (Pingding water purification plant). Most studies focus on the energy required to provide 1 m3 of water. This study combined attributes of four plants, such as the amount of energy consumed, quantity of water supplied, purified, and collected, and weather data. These data were used to model energy baselines for water providers. Artificial intelligence was imported into Microsoft Azure machine learning to train the model, which was verified using another Kaohsiung plant and one overseas to establish an online energy baseline modelling system that can be applied in various water purification plants.

Poleward migration of tropical corals inhibited by future trends of seawater temperature and calcium carbonate (CaCO3) saturation

Poleward range expansion of marine organisms is commonly attributed to anthropogenic ocean warming. However, the extent to which a single species can migrate poleward remains unclear. In this study, we used molecular data to examine the current distribution of the Pocillopora damicornis species complex in Taiwan waters and applied niche modeling to predict its potential range through the end of the 21st Century. The P. damicornis species complex is widespread across shallow, tropical and subtropical waters of the Indo-Pacific regions. Our results revealed that populations from subtropical nonreefal coral communities are P. damicornis, whose native geographical ranges are approximately between 23°N and 35°N. In contrast, those from tropical reefs are P. acuta. Our analysis of 50 environmental data layers demonstrated that the concentrations of CaCO3 polymorphs had the greatest contributions to the distributions of the two species. Future projections under intermediate shared socioeconomic pathways (SSP) 2–4.5 and very high (SSP5–8.5) scenarios of greenhouse gas emissions showed that while sea surface temperature (SST) isotherms would shift northwards, saturation isolines of two CaCO3 polymorphs, calcite (Ωcal) and aragonite (Ωarag), would shift southwards by 2100. Subsequent predictions of future suitable habitats under those conditions indicated that distinct delimitation of geographical ranges for the two species would persist, and neither would extend beyond its native geographical zones, indicating that tropical Taiwan waters are the northern limit for P. acuta. In contrast, subtropical waters are the southern limit for P. damicornis. We concluded that the decline in CaCO3 saturation would make high latitudes less inhabitable, which could be one of the boundary elements that limit poleward range expansion driven by rising SSTs and preserve the latitudinal diversity gradient (LDG) on Earth. Consequently, poleward migration of tropical reef corals to cope with warming oceans should be reevaluated.

Two new species of the congrid eel genus Ariosoma (Anguilliformes, Congridae, Bathymyrinae) from Indian waters

Two new species have been described from Indian waters, based on the materials collected from Kochi coast, Arabian Sea, Gulf of Mannar and West Bengal coast along the Bay of Bengal. Ariosoma gracilesp. nov. is described, based on five specimens collected from the landings at Kalamukku Fishing Harbour, Arabian Sea. The new species is characterised by longer tail, 55.3–58.7% TL; dorsal-fin origin above pectoral-fin base; no dark or whitish bands on dorsal surface of head, ventral extremities of lower jaw and mid-portion with minute dark pigmentation patch; anus positioned well before the middle of total length; SO canal with 4 pores; 0 or 3 pores on ST canal; pre-dorsal vertebrae 9; pre-anal vertebrae 49–54; total vertebrae 140–142. Ariosoma gracilesp. nov. is closely related to Ariosoma dolichopterum and Ariosoma emmae by sharing similar morphometrics and pre-anal vertebral counts. However, it differs by having more total pores (132–135 vs. 121–129 in A. dolichopterum, 123–126 in A. emmae); fewer pre-anal pores (43–46 vs. 47–51 in A. dolichopterum, 50–53 in A. emmae); more pre-dorsal pores (9 vs. 5–9 in A. dolichopterum, 4–6 in A. emmae). Another new species, Ariosoma kannanisp. nov. is described on the basis of two specimens (157–171 mm TL) from Gulf of Mannar and one specimen (201 mm TL) collected from Shankarpur Fish Landing Centre, West Bengal. This species is similar to Ariosoma megalops, but readily differs by having smaller eyes, smaller interorbital distance and exhibits 10.8% genetic divergence from A. megalops from the Taiwan waters.

Modeling of swordtip squid (Uroteuthis edulis) monthly habitat preference using remote sensing environmental data and climate indices

Understanding the spatial arrangement of species in maritime settings necessitates the study of oceanography. Hence, doing a study on the correlation between oceanography and species dispersion is imperative, considering the impacts of global climate change. The study used a generalized additive modeling approach to analyze the influence of oceanographic conditions on the distribution of swordtip squid in northeastern Taiwan, integrating fishing data, climatic oscillation and oceanography. Among seven oceanographic characteristics, bottom sea temperature (SSTB), sea surface height (SSH), sea surface chlorophyll (SSC), and sea surface temperature (SST) showed significant influence in generalized additive model (GAM) analysis (combined deviance explained: 40.30%). The monthly catch rate of swordtip squid is influenced by six climatic oscillations, with the Pacific Decadal Oscillation having the most significant impact, accounting for 31% of the distribution, followed by the North Pacific Gyre Oscillation at 10.8% and the Western Pacific Oscillation at 6.05%. From 2015 to 2019, the main areas for squid fishing were situated in the northeastern waters of Taiwan, precisely within the geographical coordinates of 25°N to 28°N and 121.5°E to 125°E. This study provides crucial insights for managing swordtip squid fisheries in Taiwan's northwest waters, highlighting the importance of incorporating oceanographic conditions relating to climate change information into decision-making to protect global ocean fisheries and their dependent communities.

Oil Spill Sensitivity Analysis of the Coastal Waters of Taiwan Using an Integrated Modelling Approach

Pollution caused by marine oil spills can lead to persistent ecological disasters and severe social and economic damages. Numerical simulations are useful and essential tools for accurate decision making during emergencies and planning response actions. In this study, we applied the Princeton Ocean Model (POM) to determine current data, including seawater velocity, salinity, and temperature, and we obtained the fate and trajectory of spilled oil using OpenOil. Several probable oil slicks around Taiwan were simulated over time (12 months) and space (four spill locations in the marine area of each coastal city or county) using the model. The percentage risk under the effect of an oil spill is estimated. The risk zone of the coastal waters of Taiwan was identified based on the frequency of simulated oil slicks hitting the coast and sensitive resources. This information not only helps authorities guide the preparation of effective plans to minimise the impacts of oil spill incidents but could also be used to improve regulations related to shipping and vessel navigation in regional seas.

Assessment of Tide Model Prediction and Discrepancy in Shallow Waters of Taiwan to Improve Data Fusion Methods

Assessment of Tide Model Prediction and Discrepancy in Shallow Waters of Taiwan to Improve Data Fusion Methods

A trawl collected dataset of Johnius (Actinopterygii, Sciaenidae) species in central-western Taiwanese waters.

Sciaenidae is one of the most important coastal fisheries in Taiwan, both in production and economic value. It is also significant as the main targetted diet of Chinese white dolphins, Sousachinensis, especially for the genus Johnius, such as J.taiwanensis, J.belangerii and J.distinctus, which is primarily found in central-western Taiwan coastal waters. Despite an abundance of Johnius species occurrences reported in the Global Biodiversity Information Facility (GBIF) and the Taiwan Biodiversity Information Facility (TaiBIF) data portals (Mozambique, Australia, Taiwan, Korea, India, Indonesia, South Africa, Pakistan, Vietnam and China), there are no specific datasets that properly document the regional distribution of this genus, especially in Taiwanese waters. Thus, this paper describes a dataset of genus Johnius occurrences in waters on the central-western coast of Taiwan. The data collection for the present study was conducted from 2009 until 2020 and comprised 62 sampling events and 133 occurrence records. All fish specimens were collected by trawling in Miaoli, Changhwa and Yunlin Counties, Taiwan and brought back to the lab for identification, individual number count and body weight measurement. These processing data have been integrated and established in the Taiwan Fish Database and published in GBIF. This dataset contains six Johnius species and 2,566 specimens, making it comprehensive Johnius fish fauna and spatial distributional data on the coastal habitat in central-western Taiwanese waters. This dataset contains 133 occurrence records of Johnius species (Sciaenidae) with 2,566 specimens, making it the most extensive public dataset of Johnius distribution records in Taiwan. The publication of this dataset through the TaiBIF and GBIF dataset platforms demonstrated that the number of Johnius spatial and temporal records in Taiwan waters is influenced by the topographical structure of the Changyun Rise (CYR) in combination with the cold current of the China Coastal currents and bound with the warm currents of the Kuroshio and the South China Sea on the central-western coast of Taiwan. The data serve as the foundation for understanding the biogeography and Johnius species ecology in Taiwan's coastal waters, which present a 2°C water temperature difference split at the CYR.

Assessing the Fishing Impact on the Marine Ecosystem of Guishan Island in the Northeastern Waters of Taiwan Using Ecopath and Ecosim

The northeastern waters of Guishan Island constitute one of the crucial fishing grounds for coastal trawl fishery in Taiwan and have been exploited for many decades. To construct the marine ecosystem and to examine the interactions among trophic levels of fisheries resources in the waters of Guishan Island, historical catch, catch composition, biological information, fishing effort, environmental data such as sea surface temperature, salinity, and nutrients were analyzed using Ecopath with Ecosim. The results indicated that the longline and drift net fisheries have a very minor incidental catch of cetaceans, with a fishing mortality (F) of 0.01 year−1 and an exploitation rate (E) of 0.03. The F and E were 0.308 year−1 and 0.617 for small skates and rays, and were 0.261 year−1 and 0.580, respectively, for small sharks. The F and E of the dolphinfish, Coryphaena hippurus, an important pelagic species, were 0.411 year−1 and 0.245, respectively. Fisheries had negative impact on major commercial species except the dolphinfish and the oil fish, Lepidocybium spp., which benefited from the reduction of their predators or competitors. The keystone species of the Guishan Island marine ecosystem is phytoplankton, which has the lowest trophic level and great biomass, and is an important energy source of the ecosystem. The influences of zooplankton and anchovy rank as second and third, respectively, with regard to the keystone species in the ecosystem due to their great biomass. Regarding the biomass of less abundant species, carangids had the highest influence followed by hairtail due to their feeding habits. The results of simulations using Ecosim indicated that the hairtail, small sharks, skates and rays, mackerels, and marine eels will benefit if fishing efforts are reduced by 30%. On the other hand, the biomass of phytoplankton, zooplankton, demersal benthivores, and shrimps will decrease due to the increase in the biomass of their predators.

Analysis of Unidirectional Wave Spectral Characteristics in the Northeastern Waters of Taiwan

Marine energy development has been actively promoted in recent years. Analyzing wave data in the surrounding sea areas is crucial for improving wave power generator efficiency. This study collected unidirectional wave spectrum data from the northeast coastal area of Taiwan (e.g., Pengjia Islet, Fugui Cape, and NTOU). It aimed to analyze the main wave types (i.e., wind waves, swell, or mixed sea) and their proportions during the northeast and non-northeast monsoon periods. The results indicate that swell waves dominate at the NTOU buoy, while wind waves dominate at the other two stations. In addition, the analysis of wave spectra ratios revealed a low percentage of unimodal spectra (only about 0.5%), and bimodal and unformed spectra presented as predominant. Furthermore, this study also introduced a brand new JW-J wave spectrum, which combines the JW and JONSWAP spectra to describe low-frequency (swell) and high-frequency (wind wave) systems, showing superior performance compared to the Torsethaugen spectrum. Finally, this research analyzed the JW-J spectrum parameters from the Pengjia Islet and Fugui Cape stations and applied these results to the NTOU bouy for effectively describing wind/swell wave variations.

Utilizing Numerical Models and GIS to Enhance Information Management for Oil Spill Emergency Response and Resource Allocation in the Taiwan Waters

This study aims to establish a comprehensive workflow for developing emergency response plans for both actual and scenario oil spill incidents in the Taiwan waters while addressing the resource allocation for oil spill containment as well. This workflow comprises two vital components. The first component involves the integration of numerical tools and observational data, which includes the incorporation of wind data from sources such as the National Centers for Environmental Prediction (NCEP) or meteorological stations. Additionally, it incorporates ocean current data simulated by the semi-implicit cross-scale hydroscience integrated system model (SCHISM) into the general NOAA operational modeling environment (GNOME) model, which is a new approach for this purpose. In order to assess the efficacy of this component, two distinct case studies were conducted. The first case study focused on an incident in a northern coastal area of Taiwan under open sea conditions, whereas the second case study examined an incident within a major commercial harbor in central Taiwan. The second component of this workflow involves creating oil risk maps by integrating the results from the first component with specific geographical factors into Google Earth. These oil risk maps serve multiple purposes. They offer real-time information to emergency response commanders regarding oil spill hazard prediction, and they also enable the effective development of emergency response strategies and disposal plans for potential oil spill incidents. This is achieved by generating risk maps for various scenarios using the approach outlined in the first component. Additionally, these maps assist in the assessment and planning of resource allocation for oil containment.

Study of the underwater soundscape in coastal waters of Hsinta Power Plant

This study investigates the acoustic environment surrounding the Hsinta power plant, situated in the southwest coastal waters of Taiwan, in proximity to Xingda Port, Yongan Port, Yongan Liquefied Natural Gas Plant, and an artificial reef area via passive acoustics monitoring (PAM) and Automatic Identification System (AIS). Shipping noise (including 45% fishing boats, 20.8% large ships, and 6% tug) generated by inbound and outbound ships significantly affects areas within a range of 2 to 10 km from the study site. Notably, mechanical signals at 60 and 120 Hz, not originating from PAM, have been observed, likely attributed to the Hsinta power plant. Moreover, a distinct series of periodic biological signals within the frequency range of 500 to 800 Hz is consistently recorded near the artificial reef area between 8 p.m. and 10 p.m. Assessing the impact of local ship noise by comparing ambient noise levels below the 100 Hz band before and after the coal shipping seasons (April to September), revealing a difference of 5 to 10 dB. These findings provide valuable insights into the acoustic characteristics of the study area, emphasizing the importance of understanding soundscape dynamics amidst industrial activities and maritime traffic. [Work supported in part by the National Science and Technology Council, Taiwan under Grant NSTC 112-2222-E-130-002-.]

The observation records from whale and dolphin watching inshore of Hualien, eastern Taiwan.

The eastern waters of Taiwan have been lacking baseline and research data for several years. This study was initiated by Kuroshio Ocean Education Foundation (KOEF) in collaboration with the Turumoan whale-watching company since 1998, collecting long-term ecological data for cetaceans in the inshore of Hualien, eastern Taiwan. This dataset includes 10,675 records of cetacean sightings from June 1998 to December 2021. Collection of cetacean sighting records was paused for one year in 2001 due to budgetary reasons. All of the sighting records were collected by whale-watching boat guides that were trained by KOEF. Following a standardised protocol, guides used a handheld GPS device and cetacean sighting record sheets to document information about the cetacean species identification, location, time, number of individuals, the presence of mother-calf pairs and mixed-species groups and other states of each sighting during a whale-watching tour. The collection of citizen-science data during this period has significantly advanced Taiwan's cetacean baseline data in the study area. Additionally, we make data available to the public in the form of citizen-science, making a substantial contribution to the advancement of ocean scientific research. We have published the dataset on Global Biodiversity Information Facility, allowing users around the world to download the dataset. This is currently the largest dataset of cetacean sighting records in Taiwan (last updated on 2023-09-05). We have also recorded several species on the International Union for Conservation of Nature (IUCN) Red List categorised as "Data Deficient" during our survey inshore of Hualien, eastern Taiwan, including Ginkgo-toothed beaked whales (Mesoplodonginkgodens Nishiwaki & Kamiya, 1958), Omura's whales (Balaenopteraomurai Wada, Oishi & Yamada, 2003) and killer whales (Orcinusorca (Linnaeus, 1758)). There are also sperm whales (Physetermacrocephalus Linnaeus, 1758), categorised as "Vulnerable" in the IUCN Red List and false killer whales (Pseudorcacrassidens (Owen, 1846)), categorised as "Near Threatened". This study is also the first and only long-term study that has documented cetaceans in the study area.

Occurrence, Health Risk Assessment, and Management of PPCPs in Drinking Water of Taiwan

Occurrence, Health Risk Assessment, and Management of PPCPs in Drinking Water of Taiwan