Mapping the Coastal Upwelling East of Taiwan Using Geostationary Satellite Data

Coastal upwelling is important for coastal ecosystems and the blue economy because of its large productivity and large potential for catching fish. However, coastal upwelling along the Taiwan east coast has received little attention from the research community. This study used five-year daily Himawari-8 geostationary satellite sea surface temperature data to map the coastal upwelling east of Taiwan during the summer monsoon season. We applied a semi-automatic image process technique based on the topographic position index for the quantitative upwelling mapping. The results show clear evidence of seasonal coastal upwelling along the entire Taiwan east coast, mainly under the influence of upwelling-favorable southwesterly/southerly winds. There are three broad upwelling centers along the Taiwan east coast: north, central, and south. The upwelling around the northern center has the longest upwelling season, lasting from May to September. The upwelling extents are larger between June and August during the height of the summer monsoon.


Introduction
Coastal upwelling is important for coastal ecosystems and the blue economy, because of its elevated productivity and large potential for fish catch resulted from the uplifting of nutrient-rich water towards the sea surface [1,2]. Major coastal upwelling systems are found around the world [1][2][3][4][5][6]. There are numerous upwelling hotspots in these major coastal upwelling systems, a number of which are known to occur in the China seas [3]. In their comprehensive review, ref. [3] identified 12 major upwelling regions in the China seas. Four of these upwelling regions are located in the Taiwan Strait (Southwestern Taiwan Strait, Northwestern Taiwan Strait, Taiwan Bank, and Penghu Islands), and one is located offshore Northeast of Taiwan.
Numerous studies have been focused on the offshore upwelling to the Northeast of Taiwan, which is associated with a cyclonic cold dome [7][8][9][10]. Little literature, however, has indicated the existence of coastal upwelling along the Taiwan east coast.
Taiwan Island is under the subtropical monsoon regime. Hence, the southwesterly/southerly (northeasterly/northerly) winds dominate in summer (winter). Although the southwesterly/southerly summer monsoon (June-August) is usually weaker than the northeasterly/northerly winter monsoon [11], the upwelling-favorable winds are still expected to induce upwelling along the Taiwan east coast, according to the classical Ekman transport theory. A similar mechanism has indeed been identified for the coastal upwelling in the northwestern and southwestern Taiwan Strait along the Fujian coast [3,12,13]. In addition, the Taiwan east coast is also influenced by the poleward flowing Kuroshio Current. However, the Kuroshio Current migrates offshore in summer [9,14]; and as a result, its potential influence on the summer upwelling along the Taiwan east coast is unclear and could be limited or secondary.
This study aimed to investigate the wind-driven coastal upwelling along the east coast of Taiwan using Himawari-8 sea surface temperature (SST) data. The objectives of this study include, (1) to identify significant upwelling-favorable wind events in the summers of recent several years; (2) to map the extent of the coastal upwelling using daily Himawari-8 SST data and examine its relationship with the wind events; and (3) to quantitatively investigate the characteristics of the coastal upwelling.
Remotely sensed SST data have often been used to detect and map upwelling that exhibits local SST anomaly on the sea surface [4,9,10,15,16]. A coastal upwelling event often develops rapidly after a distinct onset and lasts from several days to some weeks, depending on the local condition. The monthly-composited Moderate Resolution Imaging Spectroradiometer (MODIS) SST data used in [4] does not have an adequate temporal resolution to capture the development of an upwelling event. The daily MODIS and Advanced Very High Resolution Radiometer (AVHRR) SST data used in [15,16] are also not suitable for daily upwelling mapping because they often suffer incomplete spatial coverage due to clouds and much lower observation frequency (e.g., twice daily). The Himawari-8 SST data, on the other hand, offer unique advantages in the study of the daily development of upwelling events because of its reliability and high spatial (~2 km) and temporal resolutions (10 min) [17,18]. The merits of the Himawari-8 SST data have been recognized in the two recent upwelling studies [9,10], which used Himawari-8 SST data to investigate the upwelling off Northeastern Taiwan. In [9], the upwelling was only visually examined without attempting to explicitly map the upwelling extent; while in [10], a relatively complex gradient-based edge detection algorithm was used to map the upwelling. In this study, however, we used a scale-independent and relatively straightforward semiautomatic method developed and proved reliable in a recent upwelling study [4]. The near-real-time Himawari-8 SST data and the mapping technique developed in this study have the potential to be operationalized and thus provide a valuable monitoring tool for coastal upwelling.

Study Area, Study Period, and Its Climate
The climate of Taiwan Island is predominantly affected by the East Asian monsoon. Hence, the weak southwesterly/southerly (strong northeasterly/northerly) winds usually prevail in summer (winter) [11]. The waters surrounding Taiwan Island are influenced by three major ocean currents [19]. Taiwan's east coast is mainly affected by the northwardflowing Kuroshio Current; while the Taiwan Strait is also affected by the intrusion of the Kuroshio Current, in addition to the southward flowing China Coastal Current in winter and the northward-flowing extension of the South China Sea Warm Current (or the Taiwan Strait Warm Current) throughout the year. In general, the SST is lower in the Taiwan Strait than that in the east of Taiwan, with an annual mean of >26 • C [20]. The salinity ranges 31-34 psu and has a relatively large spatial-temporal variability in the Taiwan Strait; while, it remains relatively stable east of Taiwan [20].
In this study, the study area was located to the east of Taiwan, with two sub-areas: one 20 km buffer and the other 20-70 km buffer from the east coast ( Figure 1). The coastal upwelling is defined as the upwelling area mainly within the 20 km buffer, but may extend to the 20-70 km buffer. The study period is from May to September each year, between 2015 and 2019. Among these months, June to August is the main summer monsoon season; while, May and September are the transition period.

Wind Data, Upwelling Index and Significant Upwelling-favorable Wind Events
We used wind data to calculate an upwelling index (UI). The Ekman-based upwelling index has been widely used for upwelling studies [2,[4][5][6]. The wind data were obtained from CFSv2 (Climate Forecast System version 2) data of NCEP (National Center for Environmental Prediction) [21]. Based on ground observation, satellite data assimilation and reanalysis, the Climate Forecast System (CFS) is a global high-resolution coupling system of atmospheric-sea-land. The wind data have a spatial resolution of 0.5° × 0.5° and a temporal resolution of 6 h.
The daily composited wind data were first calculated from the 6-hourly data. The daily UI (in m 2 s −1 ) was then calculated using the following equations:

Wind Data, Upwelling Index and Significant Upwelling-favorable Wind Events
We used wind data to calculate an upwelling index (UI). The Ekman-based upwelling index has been widely used for upwelling studies [2,[4][5][6]. The wind data were obtained from CFSv2 (Climate Forecast System version 2) data of NCEP (National Center for Environmental Prediction) [21]. Based on ground observation, satellite data assimilation and reanalysis, the Climate Forecast System (CFS) is a global high-resolution coupling system of atmospheric-sea-land. The wind data have a spatial resolution of 0.5 • × 0.5 • and a temporal resolution of 6 h.
The daily composited wind data were first calculated from the 6-hourly data. The daily UI (in m 2 s −1 ) was then calculated using the following equations: Remote Sens. 2021, 13,170 4 of 20 where τ is the wind stress, ρ a is the air density (1.22 kg m −3 ), ρ w is the seawater density (1026 kg m −3 ), C d is the drag coefficient calculated using the empirical equation of [22], f is the Coriolis parameter (s −1 ), α is the wind direction, β is the general shoreline orientation, ω is the rotation rate of the Earth (7.2921 × 10 −5 rads −1 ), ϕ is the latitude, and V is the wind speed (in m s −1 ). We calculated UI at three coastal locations (Figure 1), representing the northern, central and southern sections of the Taiwan east coast. The β values were set at 18 • , 15 • , and 20 • from the northern aspect, for the northern, central, and southern locations, respectively. The significant upwelling-favorable wind events from May to September each year were then identified using the following criteria: 1.
≥5 days of positive UI, satisfying one of the following conditions: a. ≥5 days of consecutive positive UI; b.
≥3 days of consecutive positive UI, followed by 1 day of negative UI, then followed by ≥2 days of consecutive positive UI; c.
≥4 days of consecutive positive UI, followed by 2 days of negative UI, then followed by ≥3 days of consecutive positive UI; d.
≥5 days of consecutive positive UI, followed by 3 days of negative UI, then followed by ≥5 days of consecutive positive UI.

2.
The cumulative UI for the first five days ≥2.0 m 2 s −1 These criteria consider both the duration and intensity in the definition of a significant upwelling-favorable wind event. Note that a significant wind event is allowed to expand a few days into April or October to maintain the integrity of the event.

Himawari-8 SST Data and Upwelling Mapping
The Himawari-8 meteorological satellite was launched by the Japan Meteorological Agency (JMA) in October 2014 [17]. The Himawari-8 is a new generation geostationary satellite carrying an Advanced Himawari Imager (AHI), capable of providing geophysical data at a spatial resolution of~2 km and a temporal resolution of 10-min full-disk frequency [17]. The Himawari-8 (H-8) SST data used in this study were processed by the Japan Aerospace Exploration Agency (JAXA) based on a quasi-physical SST algorithm and a Bayesian cloud screen method [18]. In [18], being evaluated against the buoy data which were measured at 20-30 cm, the ocean surface skin temperature of the H-8 SST data have a bias of −0.16 • C and a Root Mean Square Difference of~0.59 • C. After considering the "cool skin effect" of the skin temperate as compared with the buoy temperature, the accuracy of the H-8 SST data is deemed reliable.
Specifically, we used cloud-free daily H-8 SST data processed to Level 3, with the highest quality-levels of 4 and 5, representing skin SST according to the GHRSST standard and format. Note that in 2015 the time-series starts from 13 July when the H-8 data became available.
Ref. [4] demonstrated that the topographic position index (TPI) [29] was an appropriate method to identify and map the negative local SST anomaly in the upwelling study. In this study, a similar semi-automatic technique based on the TPI method was developed for mapping the coastal upwelling off the Taiwan east coast from the daily H-8 SST data.
The TPI is a local-based image processing algorithm [29] that has been successfully used to map the ocean currents [30,31] and the coastal upwelling [4]. In this study, the TPI was calculated from the daily H-8 SST data using a nominated circular window with a radius of 50 cells (~100 km). This window size was deemed large enough to capture the coastal upwelling east of Taiwan.
To map upwelling areas, firstly we selected areas of negative TPI values using a pre-defined threshold, satisfying the following condition: where the S_SD(TPI) is the spatial standard deviation of the TPI image; θ is a threshold which was set at 0.75 in this study. Next, we further selected only those areas that have an areal extent greater than 100 km 2 as suitable upwelling candidates. It was assumed that an upwelling area should have a reasonably large local SST anomaly (SST_A). The local SST anomaly of an upwelling candidate was calculated using the equation below: where SST c is the mean SST value of an upwelling candidate and SST b is the mean SST value of the 20 km buffer around (but excluding) the upwelling candidate. An upwelling candidate was removed from the candidate set if its SST anomaly is less than a threshold, for example: where δ was set at 0.5 • C in this study.
The threshold values of θ and δ were selected through a trial-and-error process and the experience gained from a previous upwelling study [4]. In the final step, only the upwelling candidates that are entirely or mainly within the 20 km buffer from the coast were regarded as coastal upwelling.
To examine the relationship between the upwelling-favorable winds and the coastal upwelling mapped from the H-8 SST data, a 20 km buffer was generated for each of the three wind locations ( Figure 1). Within the 20 km buffer, one of three upwelling status was determined from the SST mapping results: (1) detected, when upwelling was detected within the buffer; (2) not-detected, when non-upwelling was detected within the buffer; (3) uncertain, when more than 50% of the buffer has no data in H-8 SST. Because the surface upwelling signature could lag behind a favorable wind event due to the travel time taken to uplift the deep water, the relationship was examined with three additional days extended to the end of each significant upwelling-favorable wind event.

Analyzing Upwelling Characteristics
After mapping the coastal upwelling east of Taiwan, we were able to quantitatively analyze several upwelling characteristics. As in [4], we calculated the area of influence (AoI), SST anomaly (SST_A), and chlorophyll-a ratio (Chla_R) of the coastal upwelling. These three upwelling properties represent different proxies of upwelling strength.
The AoI was calculated as the areal extent of the mapped upwelling areas. The SST_A was calculated using Equation (6) as the difference between the mean SST value of the upwelling areas and the mean SST value of the study area (i.e., the 70 km buffer from the coast). Similarly, Chla_R was calculated as the ratio between the mean chlorophyll-a concentration of the upwelling area and the mean chlorophyll-a concentration of the study area. The data used to calculate Chla_R were the daily composites of H-8 chlorophyll-a data processed by the JAXA to Level 3 using the algorithm developed in [32]. The spatial resolution of the chlorophyll-a data is~5 km. Note that these three characteristics were not calculated for those days when upwelling could only be partially mapped due to the cloud coverage.

Results
As examples, the H-8 SST image on 19 July 2016, the TPI image calculated from the SST image, and the corresponding chlorophyll-a image are displayed in Figures 1 and 2. Two upwelling areas were identified along the central and southern sections of the Taiwan east coast, with lower SST values, negative TPI values, and higher chlorophyll-a concentrations than the offshore area.
Two upwelling areas were identified along the central and southern sections of the Taiwan east coast, with lower SST values, negative TPI values, and higher chlorophyll-a concentrations than the offshore area.

Coastal Upwelling with Significant Upwelling-Favorable Wind Events
Each year between May and September during 2015-2019, several significant upwelling-favorable wind events were identified with variable durations along the east coast of Taiwan (Figures 3-5; Table 1). The wind events with the longest durations occur between June and August, lasting for more than 40 days (Figures 3-5; Table 1). At the northern location, the significant wind events usually start in late April and last into early October; while, at the southern and central locations, the wind events usually start in May or June and last only to mid-September (Figures 3-5). During the entire study period, the northern location experienced more such events (21 events) and lasted more days (340 days) than the central (14 events and 245 days) and southern (15 events and 265 days) locations (Table 2). On average, the southern and northern locations were influenced by much stronger upwelling-favorable winds than the central location, with the mean UI values of 1.37, 1.21, and 0.56 m 2 s −1 , respectively ( Table 2).

Coastal Upwelling with Significant Upwelling-Favorable Wind Events
Each year between May and September during 2015-2019, several significant upwellingfavorable wind events were identified with variable durations along the east coast of Taiwan (Figures 3-5; Table 1). The wind events with the longest durations occur between June and August, lasting for more than 40 days (Figures 3-5; Table 1). At the northern location, the significant wind events usually start in late April and last into early October; while, at the southern and central locations, the wind events usually start in May or June and last only to mid-September (Figures 3-5). During the entire study period, the northern location experienced more such events (21 events) and lasted more days (340 days) than the central (14 events and 245 days) and southern (15 events and 265 days) locations (Table 2). On average, the southern and northern locations were influenced by much stronger upwellingfavorable winds than the central location, with the mean UI values of 1.37, 1.21, and 0.56 m 2 s −1 , respectively ( Table 2).
The results show clear co-occurrence between the significant upwelling-favorable wind events and the detection of surface upwelling signature from the H-8 SST data, with an overall detection rate of 86.6% ± 17.1% (Figures 3-5; Table 1). Specifically, the upwelling detection rates are 92.9% ± 9.8%, 77.0% ± 22.9%, and 86.7% ± 15.7% for the north, the central, and the south locations, respectively (Table 1). Indeed, the surface upwelling signature often persisted several days after the end of the wind events (Figures 3-5). Uncertainties often occurred on the dates of strong winds (Figures 3-5), which were likely associated with overcast or stormy weather conditions preventing the H-8 sensor from acquiring cloud-free data. Remote Sens. 2021, 13, x FOR PEER REVIEW 7 of 25

Upwelling Maps
To examine the overall spatial pattern of the coastal upwelling east of Taiwan, we generated an upwelling frequency map by combining the daily upwelling maps of those upwelling dates that have complete H-8 SST data (n = 291) during the entire study period (Figure 6a). The overall spatial pattern shows that during the significant wind events, upwelling can occur along the entire east coast of Taiwan (Figure 6a). In general, further away from the coast, less frequent coastal upwelling occurs. Although upwelling could extend more than 20 km offshore, most of the upwelling occurs within 20 km off the coast (Figure 6a). In particular, the northern and southernmost coastal areas experience upwelling more often, with a frequency greater than 70% (Figure 6a). In terms of spatial extent, the areas with high (70-90%), moderate (30-70%) and low (10-30%) upwelling frequency are 15% (~1250 km 2 ), 44% (~3820 km 2 ) and 41% (~3570 km 2 ), respectively (Figure 7a).  There are substantial intra-seasonal variations in upwelling extent (Figure 8). The upwelling area increases substantially from May (~7648 km 2 ) to June (~8976 km 2 ) and July (~9516 km 2 ), especially at the northern section of the Taiwan east coast, then decreases slightly in August (~9176 km 2 ) and further reduces in September (~6968 km 2 ). Three upwelling centers (north, central, and south), with greater than 70% of upwelling frequency, are clearly visible from May to July (Figure 8a-c). The central upwelling center disappears in August; and by September, only the northern upwelling center still remains (Figure 8c,d). In May, the coastal upwelling almost always (>90% frequency) occurs at the northern center, as well as at a tiny area near the southern tip of the southern upwelling center (Figure 8a). In June, the extents of the three upwelling centers expand, with an increased and a reduced area experiencing a very-high frequency of upwelling at the southern center and the northern center, respectively (Figure 8b). In terms of spatial extent, the area experiencing very-high frequency upwelling occupies ~500 km 2 and ~830 km 2 in May and June, respectively (Figure 7b). The largest area with a high frequency of upwelling occurs in August (~1870 km 2 ), followed by June (~1620 km 2 ). There are substantial intra-seasonal variations in upwelling extent (Figure 8). The upwelling area increases substantially from May (~7648 km 2 ) to June (~8976 km 2 ) and July (~9516 km 2 ), especially at the northern section of the Taiwan east coast, then decreases slightly in August (~9176 km 2 ) and further reduces in September (~6968 km 2 ). Three upwelling centers (north, central, and south), with greater than 70% of upwelling frequency, are clearly visible from May to July (Figure 8a-c). The central upwelling center disappears in August; and by September, only the northern upwelling center still remains (Figure 8d,e). In May, the coastal upwelling almost always (>90% frequency) occurs at the northern center, as well as at a tiny area near the southern tip of the southern upwelling center (Figure 8a). In June, the extents of the three upwelling centers expand, with an increased and a reduced area experiencing a very-high frequency of upwelling at the southern center and the northern center, respectively (Figure 8b). In terms of spatial extent, the area experiencing very-high frequency upwelling occupies~500 km 2 and~830 km 2 in May and June, respectively (Figure 7b). The largest area with a high frequency of upwelling occurs in August (~1870 km 2 ), followed by June (~1620 km 2 ). Remote Sens. 2021, 13, x FOR PEER REVIEW 14 of 22 There are also substantial interannual variations in upwelling extent (Figures 7a and  9). The total upwelling area was the largest in 2016 (~10330 km 2 ), followed by in 2015 (~9120 km 2 ), and the smallest in 2019 (~7260 km 2 ). The three upwelling centers were clearly visible in 2016, 2017, and 2018 (Figure 9b-d). The very-high frequency of upwelling occurred at part of the southern section in 2015 and 2016, and at part of the northern section in 2018 (Figure 9a,b,d). Overall, the coastal upwelling occurred less frequently in 2019 than that in other years. In total, the year 2016 had the largest upwelling area with high and very-high upwelling frequency (~1900 km 2 ), followed by the year 2015 (~1600 km 2 ) (Figure 7a). There are also substantial interannual variations in upwelling extent (Figures 7a and 9). The total upwelling area was the largest in 2016 (~10330 km 2 ), followed by in 2015 (~9120 km 2 ), and the smallest in 2019 (~7260 km 2 ). The three upwelling centers were clearly visible in 2016, 2017, and 2018 (Figure 9b-d). The very-high frequency of upwelling occurred at part of the southern section in 2015 and 2016, and at part of the northern section in 2018 (Figure 9a,b,d). Overall, the coastal upwelling occurred less frequently in 2019 than that in other years. In total, the year 2016 had the largest upwelling area with high and very-high upwelling frequency (~1900 km 2 ), followed by the year 2015 (~1600 km 2 ) (Figure 7a).

Upwelling Characteristics
The monthly mean upwelling area is highest in June (~5200 km 2 ) and lowest in September (~3100 km 2 ) (Figure 10a). Statistically, the monthly mean upwelling area is significantly higher in June than that in September or in May; while, it is significantly higher in July and August than that in September. The monthly mean SST anomaly between the upwelling area and the study area (the 70 km buffer) is largest (0.99 °C) in May but smallest (0.76 °C) in August (Figure 10b). Statistically, the monthly mean upwelling's SST_A is significantly higher in May and September than that in August. In terms of chlorophyll-a ratio between the upwelling area and the study area, the monthly mean Chla_R is the largest (3.1) in September, which is significantly higher than that (2.1-2.2) in all other months statistically (Figure 10c).

Upwelling Characteristics
The monthly mean upwelling area is highest in June (~5200 km 2 ) and lowest in September (~3100 km 2 ) (Figure 10a). Statistically, the monthly mean upwelling area is significantly higher in June than that in September or in May; while, it is significantly higher in July and August than that in September. The monthly mean SST anomaly between the upwelling area and the study area (the 70 km buffer) is largest (0.99 • C) in May but smallest (0.76 • C) in August (Figure 10b). Statistically, the monthly mean upwelling's SST_A is significantly higher in May and September than that in August. In terms of chlorophyll-a ratio between the upwelling area and the study area, the monthly mean Chla_R is the largest (3.1) in September, which is significantly higher than that (2.1-2.2) in all other months statistically (Figure 10c). In terms of interannual variation, the annual mean upwelling area was highest in 2016 (~5125 km 2 ) and continued to reduce in the following three years (Figure 11a). Statistically, the annual mean upwelling area in 2016 was significantly higher than that in all of the following three years; while, the annual mean upwelling area (~4270 km 2 ) in 2017 was also significantly higher than that (~3250 km 2 ) in 2019. The annual mean SST anomaly was the largest (1.0 °C) in 2018 and the lowest (0.81 °C) in 2017 (Figure 11b). Their difference was statistically significant. The annual mean Chla_R was the highest (2.5) in 2018, and the lowest (2.2) in 2016 (Figure 11c). However, the differences between these four years were not statistically significant.

Discussion
This study shows clear evidence of seasonal coastal upwelling along the Taiwan east coast (Figures 6a and 8). The coastal upwelling is mainly detected during the summer monsoon season when the southwesterly/southerly winds are upwelling-favorable (Figures 3-5). The coastal upwelling is hence believed to be wind-driven, with a similar mechanism as the coastal upwelling along the Fujian coast on the western side of the Taiwan Strait [3,12,13]. There are two pieces of additional evidence to support this finding. Firstly, a greater number of days with significant upwelling-favorable winds is associated with a greater number of upwelling days, both at the level of individual wind events (r = 0.96; n = 50; p <0.001; Figure 12a) and at the level of the upwelling centers (r = 0.98; n=3; p <0.01; In terms of interannual variation, the annual mean upwelling area was highest in 2016 (~5125 km 2 ) and continued to reduce in the following three years (Figure 11a). Statistically, the annual mean upwelling area in 2016 was significantly higher than that in all of the following three years; while, the annual mean upwelling area (~4270 km 2 ) in 2017 was also significantly higher than that (~3250 km 2 ) in 2019. The annual mean SST anomaly was the largest (1.0 • C) in 2018 and the lowest (0.81 • C) in 2017 (Figure 11b). Their difference was statistically significant. The annual mean Chla_R was the highest (2.5) in 2018, and the lowest (2.2) in 2016 (Figure 11c). However, the differences between these four years were not statistically significant. In terms of interannual variation, the annual mean upwelling area was highest in 2016 (~5125 km 2 ) and continued to reduce in the following three years (Figure 11a). Statistically, the annual mean upwelling area in 2016 was significantly higher than that in all of the following three years; while, the annual mean upwelling area (~4270 km 2 ) in 2017 was also significantly higher than that (~3250 km 2 ) in 2019. The annual mean SST anomaly was the largest (1.0 °C) in 2018 and the lowest (0.81 °C) in 2017 (Figure 11b). Their difference was statistically significant. The annual mean Chla_R was the highest (2.5) in 2018, and the lowest (2.2) in 2016 (Figure 11c). However, the differences between these four years were not statistically significant.

Discussion
This study shows clear evidence of seasonal coastal upwelling along the Taiwan east coast (Figures 6a and 8). The coastal upwelling is mainly detected during the summer monsoon season when the southwesterly/southerly winds are upwelling-favorable (Figures 3-5). The coastal upwelling is hence believed to be wind-driven, with a similar mechanism as the coastal upwelling along the Fujian coast on the western side of the Taiwan Strait [3,12,13]. There are two pieces of additional evidence to support this finding. Firstly, a greater number of days with significant upwelling-favorable winds is associated with a greater number of upwelling days, both at the level of individual wind events (r = 0.96; n = 50; p <0.001; Figure 12a) and at the level of the upwelling centers (r = 0.98; n=3; p <0.01;

Discussion
This study shows clear evidence of seasonal coastal upwelling along the Taiwan east coast (Figures 6a and 8). The coastal upwelling is mainly detected during the summer monsoon season when the southwesterly/southerly winds are upwelling-favorable (Figures 3-5). The coastal upwelling is hence believed to be wind-driven, with a similar mechanism as the coastal upwelling along the Fujian coast on the western side of the Taiwan Strait [3,12,13]. There are two pieces of additional evidence to support this finding. Firstly, a greater number of days with significant upwelling-favorable winds is associated with a greater number of upwelling days, both at the level of individual wind events (r = 0.96; n = 50; p < 0.001; Figure 12a) and at the level of the upwelling centers (r = 0.98; n = 3; p < 0.01; Table 2). Secondly, a stronger wind-driven upwelling index is also associated with a higher upwelling detection rate (r = 0.80; n = 3; p < 0.1; Table 2), although the positive correlation at the level of individual wind events is relatively weak (r = 0.40; n = 50; p < 0.01; Figure 12b).
Remote Sens. 2021, 13, x FOR PEER REVIEW 17 of 22 Table 2). Secondly, a stronger wind-driven upwelling index is also associated with a higher upwelling detection rate (r=0.80; n=3; p<0.1; Table 2), although the positive correlation at the level of individual wind events is relatively weak (r = 0.40; n = 50; p <0.01; Figure 12b). Although upwelling can be detected along the entire east coast of Taiwan, it generally occurs around three upwelling centers (Figures 6a, 8 and 9). However, there are clear spatial and temporal variabilities among the three upwelling centers. The northern and southern upwelling centers are persistent upwelling centers during the summer monsoon season; while, the central upwelling center, which can be absent in some months and years, is much weaker in intensity (Figures 8 and 9; Tables 1 and 2). Upwelling occurs more often at the northern and southern sections of the coast, most likely due to a larger number of days with significant upwelling-favorable winds and an overall stronger upwelling index (Tables 1 and 2). In particular, the northern upwelling center usually has the longest upwelling season, lasting through the entire extended summer season from May to September (Figure 8), which is most likely because of the prolonged summer monsoon season at this location ( Figure 3). The upwelling's spatial patterns mapped in this study (Figures 6a,8 and 9) generally agree with those of the temporally-averaged SST maps (Figures 6b, 13 and 14). In particular, the three upwelling centers, which have lower SST than the offshore area, are also clearly visible in the temporally-averaged SST maps. This confirms the reliability of the semi-automatic upwelling mapping technique used in this study. Although upwelling can be detected along the entire east coast of Taiwan, it generally occurs around three upwelling centers (Figures 6a, 8 and 9). However, there are clear spatial and temporal variabilities among the three upwelling centers. The northern and southern upwelling centers are persistent upwelling centers during the summer monsoon season; while, the central upwelling center, which can be absent in some months and years, is much weaker in intensity (Figures 8 and 9; Tables 1 and 2). Upwelling occurs more often at the northern and southern sections of the coast, most likely due to a larger number of days with significant upwelling-favorable winds and an overall stronger upwelling index (Tables 1 and 2). In particular, the northern upwelling center usually has the longest upwelling season, lasting through the entire extended summer season from May to September (Figure 8), which is most likely because of the prolonged summer monsoon season at this location ( Figure 3). The upwelling's spatial patterns mapped in this study (Figures 6a, 8 and 9) generally agree with those of the temporally-averaged SST maps (Figures 6b, 13 and 14). In particular, the three upwelling centers, which have lower SST than the offshore area, are also clearly visible in the temporally-averaged SST maps. This confirms the reliability of the semi-automatic upwelling mapping technique used in this study.
The upwelling areas are larger in the main summer monsoon season between June and August than those in the transition period of May and September (Figures 8 and 10a). This is likely due to the higher upwelling index in these months ( Figure 15). Larger upwelling areas in June-August coincide with lower upwelling SST anomaly (Figure 10a,b). This might be because upwelling in June-August could be sourced from relatively shallower (and warmer) water than that in May and September, due to a stronger summer thermocline barrier [33]. A much higher chlorophyll-a ratio in September than that in other months (Figure 10c) is due to higher chlorophyll-a concentrations (up to 0.7 mg m −3 ) in the upwelling areas. This could be due to, at the beginning of autumn, the strong summer thermocline barrier becomes weaker, and surface cooling leads to deeper mixed layer depth [33], and as a result, nutrients from deeper waters can be uplifted to the surface during an upwelling event in early autumn. It should be noted that, however, this study could not separate chlorophyll-a of the terrestrial source from that of the upwelling source [9].  The upwelling areas are larger in the main summer monsoon season between June and August than those in the transition period of May and September (Figures 8 and 10a). This is likely due to the higher upwelling index in these months ( Figure 15). Larger upwelling areas in June-August coincide with lower upwelling SST anomaly ( Figure  10a,b). This might be because upwelling in June-August could be sourced from relatively The time-series of this study is not sufficient to examine the interannual variability of the upwelling characteristics. According to the Ocean Niño Index, in the summer of 2016, the strong El Niño condition was quickly transitioned into the La Niña condition. This might have something to do with the strongest upwelling activities detected in the summer season of 2016 measured by the upwelling area (Figure 11a). In contrast, weak El Niño conditions during late spring and early summer of 2019 might have contributed to the weakest upwelling activities in the summer season of 2019 (Figure 11a). Similar ENSO influence on the upwelling of the Taiwan Strait has been documented in previous studies [3].
The role of the Kuroshio Current in the summer coastal upwelling east of Taiwan remains unclear. The incursion of the Kuroshio Current to the East China Sea plays an important role in the upwelling off Northeastern Taiwan [7][8][9][10]. A similar mechanism could help form and enhance the northern upwelling center of east Taiwan. In addition, the intrusion of the Kuroshio Current into the South China Sea through the Luzon Strait may form an anticyclonic loop current off Southern Taiwan [34][35][36][37]. The influence of the Kuroshio Current to the southern upwelling center of east Taiwan, however, is likely to be limited because this loop current more likely occurs in winter [36].

Conclusions
Coastal upwelling is important for coastal ecosystems and the blue economy. This study demonstrates the advantages of using Himawari-8 SST data to monitor the daily development of upwelling events because of its reliability and high spatial (~2 km) and temporal (10 min) resolutions. The semi-automatic, topographic position index based image processing method is simple, intuitive, and scale-independent. The method is able to map the spatial extent of the upwelling SST signature and, as a result, to quantitatively investigate the upwelling characteristics such as the area of influence, SST anomaly, and chlorophyll-a concentrations. Most importantly, this study offers the pioneering insight into the summer upwelling along the entire Taiwan east coast. The coastal upwelling east of Taiwan is believed to be driven by the southwesterly/southerly summer monsoon winds, with a similar mechanism as the coastal upwelling along the Fujian coast on the western side of the Taiwan Strait.
In summary, the key findings of this study are: • Wind-driven upwelling occurs along the entire Taiwan east coast during the summer monsoon season; • There are three board upwelling centers along the Taiwan east coast: north, central, and south; The time-series of this study is not sufficient to examine the interannual variability of the upwelling characteristics. According to the Ocean Niño Index, in the summer of 2016, the strong El Niño condition was quickly transitioned into the La Niña condition. This might have something to do with the strongest upwelling activities detected in the summer season of 2016 measured by the upwelling area (Figure 11a). In contrast, weak El Niño conditions during late spring and early summer of 2019 might have contributed to the weakest upwelling activities in the summer season of 2019 (Figure 11a). Similar ENSO influence on the upwelling of the Taiwan Strait has been documented in previous studies [3].
The role of the Kuroshio Current in the summer coastal upwelling east of Taiwan remains unclear. The incursion of the Kuroshio Current to the East China Sea plays an important role in the upwelling off Northeastern Taiwan [7][8][9][10]. A similar mechanism could help form and enhance the northern upwelling center of east Taiwan. In addition, the intrusion of the Kuroshio Current into the South China Sea through the Luzon Strait may form an anticyclonic loop current off Southern Taiwan [34][35][36][37]. The influence of the Kuroshio Current to the southern upwelling center of east Taiwan, however, is likely to be limited because this loop current more likely occurs in winter [36].

Conclusions
Coastal upwelling is important for coastal ecosystems and the blue economy. This study demonstrates the advantages of using Himawari-8 SST data to monitor the daily development of upwelling events because of its reliability and high spatial (~2 km) and temporal (10 min) resolutions. The semi-automatic, topographic position index based image processing method is simple, intuitive, and scale-independent. The method is able to map the spatial extent of the upwelling SST signature and, as a result, to quantitatively investigate the upwelling characteristics such as the area of influence, SST anomaly, and chlorophyll-a concentrations. Most importantly, this study offers the pioneering insight into the summer upwelling along the entire Taiwan east coast. The coastal upwelling east of Taiwan is believed to be driven by the southwesterly/southerly summer monsoon winds, with a similar mechanism as the coastal upwelling along the Fujian coast on the western side of the Taiwan Strait.
In summary, the key findings of this study are: • Wind-driven upwelling occurs along the entire Taiwan east coast during the summer monsoon season; • There are three board upwelling centers along the Taiwan east coast: north, central, and south; • The upwelling around the northern center has the longest upwelling season, lasting from May to September; • The upwelling extents are larger between June and August during the height of the summer monsoon.