Springtime Upwelling and Its Formation Mechanism in Coastal Waters of Manaung Island, Myanmar

: Multisource satellite remote sensing data and the World Ocean Atlas 2018 (WOA18) temperature and salinity dataset have been used to analyze the spatial distribution, variability and possible forcing mechanisms of the upwelling o ﬀ Manaung Island, Myanmar. Signals of upwelling exist o ﬀ the coasts of Manaung Island, in western Myanmar during spring. It appears in February, reaches its peak in March and decays in May. Low-temperature ( < 28.3 ◦ C) and high-salinity ( > 31.8 psu) water at the surface of this upwelling zone is caused by the upwelling of seawater from a depth below 100 m. The impact of the upwelling on temperature is more signiﬁcant in the subsurface layer than that in the surface layer. In contrast, the impact of the upwelling on salinity in the surface layer is more signiﬁcant. Further research reveals that the remote forcing from the equator predominantly induces the evolution of the upwelling, while the local wind forcing also contributes to strengthen the intensity of the upwelling during spring.


Introduction
Upwelling usually refers to the upward movement of water, caused by the divergence of the flow in the surface layer of the ocean [1]. Since deeper water is usually enriched with nutrients, it tends to increase a supply of nutrients to upper oceanic layers and forms the basis for the high productivity of upwelling regions. Consequently, upwelling areas are among the most fertile regions of the global ocean [2]. The production and its variability over this coastal upwelling system are a key concern for the fishing community, since they may affect the day-to-day livelihood of the coastal population and are important for the Indian Ocean rim countries due to their developing country status [3]. Moreover, upwelling is also an important factor modulating regional and global climate. For example, the upwelling along the western coast of Java-Sumatra has changed sea surface temperature (SST) of the warm pool in the eastern Indian Ocean and caused anomalous atmospheric circulation, thereby affecting intraseasonal to decadal variabilities of the tropical climates [4][5][6][7][8][9]. Therefore, the understanding physical processes and their variabilities in the upwelling systems is important because it gives us crucial information regarding variability of a marine ecosystem and its regional climate [10,11].
The Bay of Bengal (BOB) is dominated by the South Asian monsoon. The southwest monsoon prevails in summer (June-August), while the northeast monsoon prevails in winter  Figure 1a indicates the propagation path of the equatorial Kelvin wave and coastal Kelvin wave. The red dots in the blue line denote the selected stations of sea level anomaly (SLA) to track the propagation of Kevin wave along the equator and eastern boundary of the Bay. In Figure 1a, the India coastal upwelling, Sri Lanka coastal upwelling, Sri Lanka Cold Eddy and Sri Lanka are denoted by ICU, SLCU, SLCE and SL, respectively. In Figure 1b

Data
The World Ocean Atlas 2018 (WOA18) has been newly released by the Ocean Climate Laboratory of the National Centers for Environmental Information (NCEI) and the National Oceanic and Atmospheric Administration (NOAA) of the U.S. The WOA18 collected temperature and salinity samples obtained by all ship-deployed Conductivity-Temperature-Depth (CTD) packages, profiling floats, moored and drifting buoys, gliders and undulating oceanographic recorder profiles. Using an  Figure 1a indicates the propagation path of the equatorial Kelvin wave and coastal Kelvin wave. The red dots in the blue line denote the selected stations of sea level anomaly (SLA) to track the propagation of Kevin wave along the equator and eastern boundary of the Bay. In Figure 1a, the India coastal upwelling, Sri Lanka coastal upwelling, Sri Lanka Cold Eddy and Sri Lanka are denoted by ICU, SLCU, SLCE and SL, respectively. In Figure 1b The strong upwelling around the Madras coast was first reported in 1964 from in-situ data [25], and the upwelled water was found to come from a 30 m layer induced by the strong winds along the coast of Madras. Based on hydrographic data collected during the summer monsoon of 1989, Remote Sens. 2020, 12, 3777 3 of 15 the upwelling zone off the eastern Indian Peninsula (Figure 1a) is found to extend in coastal waters from Madras to Visakhapatnam in summer [26,27], within a range of approximately 40 km offshore. Alongshore winds are an important dynamic mechanism for its formation [15]. Below this upwelling band, downwelling is often observed suggesting the presence of an undercurrent [15].
As mentioned above, previous studies mainly focused on the upwelling in the western boundary of the Bay [19][20][21][22][23][24][25][26][27][28]. In contrast, there are very few reports on the upwelling in its eastern boundary. Yesaki and Jantarapagdee [29] noted wind-induced upwelling over the continental shelf off the west coast of Thailand. The observations of La Fond [30] indicated that during April-May 1963, SST in the northeastern Bay exceeded 28 • C while it was relatively low (<27 • C) in the northern coasts off Myanmar, reflecting the possible impact of cold water upwelling from the subsurface. Follow-up studies also demonstrated the low-temperature zone located at the northwestern coasts of Myanmar 15-20 • N during the winter monsoon, which was considered a sign of upwelling [12,31]. This conclusion was also supported by the satellite-derived Chlorophyll-a (Chl-a) data [32]. Akester [33] inferred that the high primary productivity in the continental shelf off Myanmar 16-20 • N was caused by the upwelling there in April 2015, and the further evidence of upwelling, including low-temperature, high-salinity and low-dissolved oxygen, was shown in the coastal region.
In summary, although the existing studies have discovered signs of upwelling in the northwestern coasts of Myanmar in winter and spring, most studies were only based on scattered observations from a single cruise. Until now, little has been known about the spatial distribution and seasonal evolution of the upwelling, and there is a lack of research on its dynamic mechanisms. Therefore, in this paper, satellite observations and the World Ocean Atlas 2018 (WOA18) temperature and salinity data were used to analyze the characteristics and causes of the upwelling in the northwestern coast of Myanmar during the winter monsoon. The remainder of this paper is organized as follows: Section 2 introduces the data and method used in this study; Section 3 analyzes the spatial distribution and evolution process of the upwelling in the study area; Section 4 presents the dynamic mechanisms of the upwelling; discussions and, finally, conclusions are stated in Sections 5 and 6, respectively.

Data
The World Ocean Atlas 2018 (WOA18) has been newly released by the Ocean Climate Laboratory of the National Centers for Environmental Information (NCEI) and the National Oceanic and Atmospheric Administration (NOAA) of the U.S. The WOA18 collected temperature and salinity samples obtained by all ship-deployed Conductivity-Temperature-Depth (CTD) packages, profiling floats, moored and drifting buoys, gliders and undulating oceanographic recorder profiles. Using an objective analysis technique, the raw data are processed into gridded data with a horizontal resolution of 0.25 • and 102 vertical levels above 5500 m (including 21 layers, at 5 m spacing, for the top 100 m). Data are presented for climatological composite periods (annual, seasonal, monthly, seasonal and monthly difference fields from the annual mean field, and the number of observations) at 102 standard depths. Standard error of the mean fields is binned into several ranges depending on the depth level. In our study region (90.5-97.5 • E, 15.5-21.5 • N), the error of temperature in upper 100 m layer ranges from 0.00 to 0.89 • C, and salinity error ranges between 0.00 and 0.33 psu [34,35]. In this study, climatological monthly temperature and salinity were used in the analysis of the upwelling distribution.
The Moderate Resolution Imaging Spectroradiometer (MODIS) daily Chl-a data were provided by the U.S. National Aeronautics and Space Administration (NASA). The chlorophyll concentration was derived using the OC3 algorithm [36]. There is a root mean squared error (RMSE) of 1.228 mg m −3 for MODIS Chl-a against in situ observations [37]. MODIS Chl-a with a horizontal resolution of 4 km [38] spanning between 1 January 2003 and 31 March 2020, were used to study seasonal evolution of the upwelling within the study area.
Sea surface wind data were from the cross-calibrated multi-platform (CCMP) product provided by the NASA's Physical Oceanography Distributed Active Archive Center (PODAAC). CCMP is a multisource fusion surface wind data derived from satellite observations, including measurements of Special Sensor Microwave Imager (SSM/I), Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) and Tropical Rainfall Measuring Mission Microwave Imager (TMI). To create the CCMP winds, an enhanced variational analysis method (VAM) performs quality control and combines all available Remote Sensing Systems (RSS) cross-calibrated wind data with available conventional ship and buoy data and European Centre for Medium-Range Weather Forecasts (ECMWF) analyses. Overall, the CCMP analysis has the best overall fit to the in-situ observations with an error speed difference ranging from 1.6 m s −1 versus ships to 0.6 m s −1 versus the higher-quality TAO buoys. This is also seen in the RMSE direction fit, which ranges from 11.5 • to 7.0 • . This dataset has a horizontal resolution of 0.25 • × 0.25 • and a temporal resolution of 6 h from the period 1993 to 2016 [39].
Daily sea level anomaly (SLA) data were provided by Copernicus Marine Environment Monitoring Service (CMEMS), with a horizontal resolution of 0.25 • × 0.25 • . This gridded data is derived from the multi-satellite measurements of TOPEX/Poseidon (T/P), European Remote Sensing Satellite-1 (ERS-1) and European Remote Sensing Satellite-2 (ERS-2). There is an error of~0.02 m for the SLA data [40]. The SLA data we used are from between 1 January 1993 and 31 March 2020. These two data sets of CCMP wind and CMEMS SLA were averaged to obtain weekly and monthly values and then were used to analyze the dynamic mechanisms of the upwelling.

Method to Identify Area of the Upwelling and Its Intensity
Water that rises to the surface as a result of upwelling is characterized by low-temperature and high-salinity. Both SST and SSS have been extensively used as an indicator of upwelling in many previous studies (e.g., [41,42]). Correspondingly, upwelling areas were identified by a consideration of both SST and SSS from WOA18 in our study. As we shall see below that the region around Myanmar coasts with SST below 28.3 • C has a relatively high salinity. Thus, the regions around the coast of Sri Lanka with SST less than 28.3 • C are used to describe the regions affected by the coastal upwelling.
In addition, the upwelled water from subsurface is typically rich in nutrients. In the upwelling region, surface water has relatively high Chl-a concentration (e.g., [19,41]), besides its typical relatively cold and saline nature. Thus, Chl-a concentration is used as an indicator for the intensity of upwelling in previous works (e.g., [43]). In this study, climatological weekly Chl-a concentration averaged in the upwelling region during January-April is constructed from the daily MODIS Chl-a data and then it is used to explore the evolution processes of the intensity of the upwelling. A larger magnitude of Chl-a concentration denotes a stronger upwelling.

Spatial-Temporal Distribution of the Coastal Upwelling off Manaung Island Using WOA18 Data
The sample amounts of temperature and salinity in the surface layer for each 1 • × 1 • grid during spring (March to May) are shown in Figure 2a,b, respectively. Generally speaking, observations of temperature and salinity have a good coverage over the study region. In the coastal areas, the temperature and salinity samples around Manaung Island are 7-18 and 6-16, respectively, and both temperature and salinity samples in open ocean are far larger than that in the nearshore areas. Thus, WOA18 provides adequate data for this study to investigate the spatial structure of the upwelling in the northwestern coasts of Myanmar. 31.8 psu, and the high-salinity center is located in the northwestern coast of Manaung Island. The high-salinity (>31.8 psu) and low-temperature (<28.3 °C) zones do not completely overlap, with the high-salinity zone occurring slightly to the north of the low-temperature zone. This low-temperature (<28.3 °C) and high-salinity (>31.8 psu) water at the surface may reflect the occurrence of upwelling near the Manaung Island in spring. This is broadly consistent with the location of upwelling signals revealed by Akester [33] using in-situ data obtained in April 2015.   (30-70 m). This means that the impact of upwelling on temperature is more pronounced in the subsurface than that in the surface layer. This may be caused by two factors. Firstly, the subsurface layer lies within the seasonal thermocline and the vertical temperature gradient there is larger than that in the near-surface layer, and thus the impact of upwelling on Distributions of surface temperature and salinity in spring are shown in Figure 2c,d, respectively. SST in the study region is relatively high (>28 • C), owing to strong shortwave radiation during spring [44]. The highest SST (29.5 • C) appears in the coastal areas around the estuary of the Irrawaddy River. A low-temperature zone is found distributed along the northwestern coasts of Myanmar, with a cold center located off the coast of Manaung Island. Taken the 28.3 • C isotherm as the boundary of the upwelling, the upwelling zone appears as "tongue-shaped", extending across 17.7-19.3 • N offshore of the Manaung Island, to the east of 92.3 • E. Affected by a large amount of runoff and precipitation in the northern Bay, salinity of the study area is relatively low (<33.0 psu). The lowest salinity (~30 psu) occurs offshore of the northwestern Myanmar ( Figure 2d). Similar to the low-temperature zone, a high-salinity zone is also formed in the waters around the Manaung Island, with salinity exceeding 31.8 psu, and the high-salinity center is located in the northwestern coast of Manaung Island. The high-salinity (>31.8 psu) and low-temperature (<28.3 • C) zones do not completely overlap, with the high-salinity zone occurring slightly to the north of the low-temperature zone. This low-temperature (<28.3 • C) and high-salinity (>31.8 psu) water at the surface may reflect the occurrence of upwelling near the Manaung Island in spring. This is broadly consistent with the location of upwelling signals revealed by Akester [33] using in-situ data obtained in April 2015. The amplitude of the isotherm uplift is significantly smaller in the near-surface layer (above 30 m layer) than that in the subsurface layer (30-70 m). This means that the impact of upwelling on temperature is more pronounced in the subsurface than that in the surface layer. This may be caused by two factors. Firstly, the subsurface layer lies within the seasonal thermocline and the vertical temperature gradient there is larger than that in the near-surface layer, and thus the impact of upwelling on temperature is more significant in the subsurface layer. Secondly, the intensity of upwelling may be partially suppressed by strong salinity stratification in the near-surface layer [16], which further weakens upwelling near the surface and thus results in weaker uplift of isotherm there. The sectional distribution of salinity also reveals upwelling of high-salinity water in the subsurface layer off the coast of Manaung Island (Figure 3b). However, in contrast to temperature, the upwelling of high-salinity water is more significant in the upper layer above 30 m than that in the subsurface layer. This is mainly because both the freshwater input from river discharge and the precipitation into the surface make the vertical salinity gradient larger in the upper layer than that in the subsurface layer [45]. Therefore, even under the condition of a weaker uplift of isotherms in the near-surface layer (Figure 3a), the uplift of isohalines is much more significant in the upper layer.
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 16 temperature is more significant in the subsurface layer. Secondly, the intensity of upwelling may be partially suppressed by strong salinity stratification in the near-surface layer [16], which further weakens upwelling near the surface and thus results in weaker uplift of isotherm there. The sectional distribution of salinity also reveals upwelling of high-salinity water in the subsurface layer off the coast of Manaung Island (Figure 3b). However, in contrast to temperature, the upwelling of highsalinity water is more significant in the upper layer above 30 m than that in the subsurface layer. This is mainly because both the freshwater input from river discharge and the precipitation into the surface make the vertical salinity gradient larger in the upper layer than that in the subsurface layer [45]. Therefore, even under the condition of a weaker uplift of isotherms in the near-surface layer (Figure 3a), the uplift of isohalines is much more significant in the upper layer. Seasonal evolution for the coastal upwelling off Manaung Island is further explored using WOA18 climatologically monthly temperature and salinity data ( Figure 4). In February, a highsalinity zone starts to form in the surface off the coasts of Manaung Island, while the low-temperature signal is not significant. This is mainly because the southward flowing low-temperature water from the northern Bay along the coasts during that period [46] covers the upwelled cold water from subsurface. In fact, both the low-temperature and high-salinity signals of the upwelling are clearly presented in the subsurface layers between 40 and 80 m around the coastal areas of Manaung Island in February ( Figure 5). In March, the extent of low-temperature and high-salinity water is the largest, indicating that the coastal upwelling reaches its peak in this month. Subsequently, the extent of lowtemperature and high-salinity water significantly shrinks in April and disappears in May. Thus, April and May are periods for the decay and dissipation of the coastal upwelling of Manaung Island. Seasonal evolution for the coastal upwelling off Manaung Island is further explored using WOA18 climatologically monthly temperature and salinity data (Figure 4). In February, a high-salinity zone starts to form in the surface off the coasts of Manaung Island, while the low-temperature signal is not significant. This is mainly because the southward flowing low-temperature water from the northern Bay along the coasts during that period [46] covers the upwelled cold water from subsurface. In fact, both the low-temperature and high-salinity signals of the upwelling are clearly presented in the subsurface layers between 40 and 80 m around the coastal areas of Manaung Island in February ( Figure 5). In March, the extent of low-temperature and high-salinity water is the largest, indicating that the coastal upwelling reaches its peak in this month. Subsequently, the extent of low-temperature and high-salinity water significantly shrinks in April and disappears in May. Thus, April and May are periods for the decay and dissipation of the coastal upwelling of Manaung Island.

Dynamic Mechanisms for Seasonal Evolution of the Coastal Upwelling of Manaung Island Using Satellite Measurements
As is well known, local surface wind is an important factor influencing the generation and evolution of coastal upwelling [28,47]. Therefore, the role of local wind on the formation and evolution of the coastal upwelling of Manaung Island is explored here.
Affected by the South Asian monsoon, the predominant wind in the northern Bay is northwesterly/northerly during winter and spring [12]. The direction of wind stress is roughly parallel to the NW-SE orientation of the coastline of northwestern Myanmar (Figure not shown). According to the classic Ekman transport theory, i.e., in the Northern Hemisphere, the wind blows in a direction parallel to the coast (on the left side of the wind direction), and it causes net movement of surface water at about 90 degrees to the right of the wind direction. Because the surface water flows

Dynamic Mechanisms for Seasonal Evolution of the Coastal Upwelling of Manaung Island Using Satellite Measurements
As is well known, local surface wind is an important factor influencing the generation and evolution of coastal upwelling [28,47]. Therefore, the role of local wind on the formation and evolution of the coastal upwelling of Manaung Island is explored here.
Affected by the South Asian monsoon, the predominant wind in the northern Bay is northwesterly/northerly during winter and spring [12]. The direction of wind stress is roughly parallel to the NW-SE orientation of the coastline of northwestern Myanmar (Figure not shown). According to the classic Ekman transport theory, i.e., in the Northern Hemisphere, the wind blows in a direction parallel to the coast (on the left side of the wind direction), and it causes net movement of surface water at about 90 degrees to the right of the wind direction. Because the surface water flows

Dynamic Mechanisms for Seasonal Evolution of the Coastal Upwelling of Manaung Island Using Satellite Measurements
As is well known, local surface wind is an important factor influencing the generation and evolution of coastal upwelling [28,47]. Therefore, the role of local wind on the formation and evolution of the coastal upwelling of Manaung Island is explored here.
Affected by the South Asian monsoon, the predominant wind in the northern Bay is northwesterly/northerly during winter and spring [12]. The direction of wind stress is roughly parallel to the NW-SE orientation of the coastline of northwestern Myanmar (Figure not shown). According to the classic Ekman transport theory, i.e., in the Northern Hemisphere, the wind blows in a direction parallel to the coast (on the left side of the wind direction), and it causes net movement of surface water at about 90 degrees to the right of the wind direction. Because the surface water flows away from the coast, the water is replaced with water from below [48]. Thus, the northwesterly wind stress will cause an offshore Ekman transport there, thereby triggering the upwelled subsurface water to the surface layer as a compensatory for the offshore transport, which induces the formation of coastal upwelling. To facilitate further analysis, the coordinate axes are rotated clockwise by 45 • , and the wind stress is decomposed into directions parallel to and perpendicular to the coasts, with the southeast direction along the coastline being positive (as shown in Figure 1b). During the northeast monsoon (December to April of the following year), the coast of Manaung Island experiences the northwesterly wind, favorable to the formation of upwelling. The wind stress is generally weak (0.01-0.02 N m −2 ) and reaches its maximum (approximately 0.02 N m −2 ) in January (Figure 6a). At the same time, the wind stress curl reaches the maximum (approximately 2 × 10 −7 N m −3 ) (Figure 6b). Then, the wind direction reverses to a southeasterly one during the monsoon transition period in May, which is not favorable to the formation of upwelling. Thus, local wind stress may contribute to the formation of the upwelling in spring, although its intensity is relatively weak, and the timing of its maximum is two months out of phase with the peak of upwelling (in March). That is to say that although the local wind stress contributes to the development of upwelling, there may be other factors (i.e., remote forcing from the equator) that have important contributions to the development of the coastal upwelling of Manaung Island, as we shall see below.
Remote Sens. 2020, 12, x FOR PEER REVIEW 8 of 16 away from the coast, the water is replaced with water from below [48]. Thus, the northwesterly wind stress will cause an offshore Ekman transport there, thereby triggering the upwelled subsurface water to the surface layer as a compensatory for the offshore transport, which induces the formation of coastal upwelling. To facilitate further analysis, the coordinate axes are rotated clockwise by 45°, and the wind stress is decomposed into directions parallel to and perpendicular to the coasts, with the southeast direction along the coastline being positive (as shown in Figure 1b). During the northeast monsoon (December to April of the following year), the coast of Manaung Island experiences the northwesterly wind, favorable to the formation of upwelling. The wind stress is generally weak (0.01-0.02 N m −2 ) and reaches its maximum (approximately 0.02 N m −2 ) in January (Figure 6a). At the same time, the wind stress curl reaches the maximum (approximately 2 × 10 −7 N m −3 ) (Figure 6b). Then, the wind direction reverses to a southeasterly one during the monsoon transition period in May, which is not favorable to the formation of upwelling. Thus, local wind stress may contribute to the formation of the upwelling in spring, although its intensity is relatively weak, and the timing of its maximum is two months out of phase with the peak of upwelling (in March). That is to say that although the local wind stress contributes to the development of upwelling, there may be other factors (i.e., remote forcing from the equator) that have important contributions to the development of the coastal upwelling of Manaung Island, as we shall see below.  Figure 1a and 1b, respectively. All climatological monthly mean data is averaged from the daily/6 h data for the period 1993 to 2016.
Besides local wind, remote forcing from the equator is also the main forcing that induces seasonal variabilities of the BOB. In particular, it plays a dominant role at the eastern and northern boundary of the Bay [49][50][51]. It affects the interior Bay via the planetary waves (including coastal Kelvin waves and the westward propagating Rossby waves triggered by the equatorial processes) [52]. Referring to Cheng et al. [51], the propagation path of the planetary waves from the equator to  Figure 1a,b, respectively. All climatological monthly mean data is averaged from the daily/6 h data for the period 1993 to 2016.
Besides local wind, remote forcing from the equator is also the main forcing that induces seasonal variabilities of the BOB. In particular, it plays a dominant role at the eastern and northern boundary of Remote Sens. 2020, 12, 3777 9 of 15 the Bay [49][50][51]. It affects the interior Bay via the planetary waves (including coastal Kelvin waves and the westward propagating Rossby waves triggered by the equatorial processes) [52]. Referring to Cheng et al. [51], the propagation path of the planetary waves from the equator to the interior Bay along the 200 m isobath (the blue line in Figure 1a) is selected, and the SLA evolutions at each station of December and May of the following year are illustrated (Figure 6c). The upwelling Kelvin wave originates from the equator, propagates into the eastern boundary of the BOB in early February and reaches the coastal region of Manaung Island in early March, resulting in a remarkable decrease in SLA there (Figure 6c). This further amplifies the upwelling induced by the local wind stress (Figure 6a,b), so that the upwelling peaks in March. Subsequently, the downwelling Kelvin wave, triggered by equatorial westerlies in early April, reaches the Manaung Island in early May, thereby making an important contribution to the decay of the upwelling.
Rao et al. [53] demonstrated that the aforementioned equatorial upwelling (downwelling) Kelvin wave is driven by equatorial easterlies (westerlies) wind. The lead-lag correlation analysis between SLA off the Manaung Island (Station 26 shown in Figure 1a) and that at the central equator (Station 1 shown in Figure 1a) shows that the correlation coefficient is positive and reaches maximum of 0.64 (significant at 99% confidence level) when variation of SLA at Manaung Island (SLA_n) lags behind SLA in the central equator (SLA_s) by 32 days (Figure 7 Figure 1a) is selected, and the SLA evolutions at each station of December and May of the following year are illustrated (Figure 6c). The upwelling Kelvin wave originates from the equator, propagates into the eastern boundary of the BOB in early February and reaches the coastal region of Manaung Island in early March, resulting in a remarkable decrease in SLA there (Figure 6c). This further amplifies the upwelling induced by the local wind stress (Figure 6a,b), so that the upwelling peaks in March. Subsequently, the downwelling Kelvin wave, triggered by equatorial westerlies in early April, reaches the Manaung Island in early May, thereby making an important contribution to the decay of the upwelling. Rao et al. [53] demonstrated that the aforementioned equatorial upwelling (downwelling) Kelvin wave is driven by equatorial easterlies (westerlies) wind. The lead-lag correlation analysis between SLA off the Manaung Island (Station 26 shown in Figure 1a) and that at the central equator (Station 1 shown in Figure 1a  The above analysis has demonstrated that both local wind and equatorial forcing have important impacts on the formation and evolution of the coastal upwelling of Manaung Island. The concentration of Chl-a is another important indicator of the upwelling intensity [43]. The MODIS images of Chl-a for the upwelling period between February and March in 2020 (Figure 8) demonstrate that the short-term evolution of Chl-a distribution in the upwelling zone off Manaung Island is consistent with that of SLA. This consistency underpins that both Chl-a and SLA can well characterize in the intensity of the coastal upwelling of Manaung Island. The above analysis has demonstrated that both local wind and equatorial forcing have important impacts on the formation and evolution of the coastal upwelling of Manaung Island. The concentration of Chl-a is another important indicator of the upwelling intensity [43]. The MODIS images of Chl-a for the upwelling period between February and March in 2020 (Figure 8) demonstrate that the short-term evolution of Chl-a distribution in the upwelling zone off Manaung Island is consistent with that of SLA. This consistency underpins that both Chl-a and SLA can well characterize in the intensity of the coastal upwelling of Manaung Island. Remote Sens. 2020, 12, x FOR PEER REVIEW 10 of 16 The daily MODIS Chl-a data in the core area of upwelling (the red box in Figure 2c), from January to April, between 2003 and 2020, were composited and then bin-averaged to obtain the climatologically weekly averaged Chl-a time series from January to April. The SLA_n, SLA_s, alongshore wind stress, and wind stress curl were also analyzed using the same method. Then, all these time series were normalized by 1 standard deviation of each ( Figure 9). The concentration of Chl-a in the upwelling zone (blue line) continues to increase from January, reaching its peak in early March, and then it decreases during late March and throughout April (Figure 9). This evolution pattern is consistent with variations in the upwelling intensity revealed by the WOA18 temperature and salinity data. As there is a significant correlation between Chl-a concentration and SLA_n with a correlation coefficient of −0.76, SLA_n serves as a good indicator for the intensity of the upwelling at Manaung Island. Thus, the relative contribution of each influencing factor to variations in the upwelling was evaluated by calculating its correlation with SLA_n. The results of correlation analysis (Table 1) demonstrate that there is a high positive correlation between weekly SLA_n and the weekly SLA_s with SLA_s leading SLA_n by 35 days, and their correlation coefficient is 0.60, significant at 95% confidence level. In contrast, SLA_n is negatively correlated with alongshore wind stress and local wind stress curl, with correlation coefficients of −0.35 and −0.33, respectively. These two correlations are not statistically significant (p > 0.05), which indicates that the local wind stress in the upwelling region of Manaung Island has a relatively small impact on the variation in SLA_n. Therefore, the correlation between SLA_n and the equatorial forcing is much stronger than the local wind (0.60 vs. −0.35/−0.33), underpinning that although both local wind and equatorial forcing have impacts on evolution of the upwelling at Manaung Island, the remote forcing from the equator plays a more important role than the local wind does. The daily MODIS Chl-a data in the core area of upwelling (the red box in Figure 2c), from January to April, between 2003 and 2020, were composited and then bin-averaged to obtain the climatologically weekly averaged Chl-a time series from January to April. The SLA_n, SLA_s, alongshore wind stress, and wind stress curl were also analyzed using the same method. Then, all these time series were normalized by 1 standard deviation of each ( Figure 9). The concentration of Chl-a in the upwelling zone (blue line) continues to increase from January, reaching its peak in early March, and then it decreases during late March and throughout April (Figure 9). This evolution pattern is consistent with variations in the upwelling intensity revealed by the WOA18 temperature and salinity data. As there is a significant correlation between Chl-a concentration and SLA_n with a correlation coefficient of −0.76, SLA_n serves as a good indicator for the intensity of the upwelling at Manaung Island. Thus, the relative contribution of each influencing factor to variations in the upwelling was evaluated by calculating its correlation with SLA_n. The results of correlation analysis (Table 1) demonstrate that there is a high positive correlation between weekly SLA_n and the weekly SLA_s with SLA_s leading SLA_n by 35 days, and their correlation coefficient is 0.60, significant at 95% confidence level. In contrast, SLA_n is negatively correlated with alongshore wind stress and local wind stress curl, with correlation coefficients of −0.35 and −0.33, respectively. These two correlations are not statistically significant (p > 0.05), which indicates that the local wind stress in the upwelling region of Manaung Island has a relatively small impact on the variation in SLA_n. Therefore, the correlation between SLA_n and the equatorial forcing is much stronger than the local wind (0.60 vs. −0.35/−0.33), underpinning that although both local wind and equatorial forcing have impacts on evolution of the upwelling at Manaung Island, the remote forcing from the equator plays a more important role than the local wind does.

Discussions
In the upwelling region, upwelled supplies of nitrate and the tremendous blooms of phytoplankton they support, render these regimes the "marine ranchings" [55]. This phytoplankton biomass then feeds into productive food chains which support a significant share of the biological resources that humans harvest from the ocean, and indeed, attract commercial fishers and fisheries [56]. Therefore, the coastal upwelling has a profound impact on coastal populations, and it is particularly important for developing countries (e.g., Myanmar). Despite of the importance of the upwelling, until now, little has been known about the coastal upwelling off Myanmar. This is mainly due to the sparseness of observational data and to the multi-scale aspects of upwelling that are difficult to measure and simulate in climate/ocean models.
Our analysis demonstrates that the center of the upwelling is located at the coastal region off Manaung Island during spring. This is generally consistent with the location of the upwelling reported in the previous studies [30,32,33]. However, seasonal evolution characteristics were not discussed in these previous works due to the scatted measurements from a single cruise used in these studies. In contrast, our study indicates that the upwelling generates in February, then peaks in March and finally decays in May. Based on the satellite observations, the influence of local winds and remote equatorial forcing on the evolution of the upwelling is further discussed by the correlation analysis for the first time. Nevertheless, the physical mechanisms of the upwelling are preliminary. Models will be helpful to isolate the dynamics and assess the relative contribution of various processes in the future study. Our results would contribute to the understanding of coastal processes (including the coastal upwelling) off Myanmar and regional climate variability and their prediction.
Note that the Indian Ocean Dipole (IOD) and El Niño and Southern Oscillation (ENSO) are major climate modes that affect the interannual variability in the northern Indian Ocean [57][58][59]. They can cause interannual variabilities in the regional upwelling in the Indian Ocean by inducing variations in the sea surface wind field [59,60]. For example, significantly interannual variability of the Java-Sumatra upwelling is primarily driven by local winds and remote equatorial Indian Ocean winds associated with IOD/ENSO events [9,18,61]. As our analysis has shown, the coastal upwelling off Myanmar is mainly caused by the remote equatorial forcing and to a lesser degree by the local winds. Both these two forcing indicate significant interannual variabilities caused by IOD and ENSO events [58]. However, how and to what extent the IOD and ENSO events affect the interannual variabilities of the upwelling off Myanmar are still unknown and await further study. A series of cruises are suggested to be carried out in the upwelling region off Myanmar coasts during normal years and IOD/ENSO years to address these important scientific issues.

Conclusions
In this paper, WOA18 temperature and salinity datasets, combined with the MODIS Chl-a data, SLA and the CCMP wind data, were used to analyze the spatial structure and seasonal evolution of the coastal upwelling around the Manaung Island. Possible dynamic mechanisms leading to the evolution process of the upwelling were also explored for the first time. The upwelling near the Manaung Island has been identified, which begins in February, reaches its peak in March and disappears in May. Low-temperature (<28.3 • C) and high-salinity (>31.8 psu) water at the surface of this upwelling zone mainly upwells from a depth of 100 m. The impact of upwelling on temperature is found to be more significant in the subsurface layer than that in the surface layer. In contrast, the impact of upwelling on salinity in the surface layer is more significant. This different effect of the upwelling on the temperature and salinity in upper layer is possibly caused by the strong salinity stratification there.
During winter and spring, the local wind near the Manaung Island is dominated by northwesterly/northerly wind. It tends to induce offshore transport in the upper water volume, and consequently is favorable to the formation of coastal upwelling. In addition, the upwelling Kelvin wave, driven by the equatorial easterly wind stress during January-March, propagates along the eastern boundary of the Bay of Bengal, which promotes uplifting of the thermocline and halocline in the coastal regions around the Manaung Island, further enhancing the upwelling. In April, the downwelling Kelvin wave, propagating into the eastern boundary of the Bay from the equator, deepens the thermocline in the waters around the Manaung Island, and simultaneously the local wind is weakening. Both these factors may play an important role in the dissipation of upwelling in this area. Correlation analysis further demonstrates that the remote equatorial forcing exerts a more significant effect on the upwelling variability than the local wind does. Although our study is qualitative, the physical processes are robust features of this region and their quantification will be a major challenge for future observational as well as modeling studies.
Despite the joint effects of remote equatorial forcing and local wind can well explain the generation and evolution of the coastal upwelling around the Manaung Island, we note that the center of the upwelling is located at the Manaung Island where the offshore bottom topography near the Manaung Island has the largest gradient towards the coast. This sharp change in topography may have a significant impact on the distribution of the upwelling center [62], which awaits further analysis in the future.