Study on the propagation velocity of internal solitary waves in the Andaman Sea using Terra/Aqua-MODIS remote sensing images*

Lina SUN , Jie ZHANG , Junmin MENG ,**

1 First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China

2 Technology Innovation Center for Ocean Telemetry, Ministry of Natural Resources, Qingdao 266061, China

Abstract The Andaman Sea has been a classic study region for internal solitary waves (ISWs) for several decades, and extraordinarily large ISWs are characteristic of the Andaman Sea in the Indian Ocean. This paper presents results on the estimation of the propagation velocity of ISWs in the Andaman Sea that were tracked using 195 image pairs acquired by MODIS National Aeronautics and Space Administration (NASA)Terra/Aqua satellites between January 2014 and December 2018. A total of 562 ISWs were identif ied during the period, and the results of the propagation velocity distribution of ISWs in the Andaman Sea are presented. The estimated propagation velocity of ISWs agrees well with the theoretical results derived from the Korteweg-de Vries (KdV) equation using monthly climatology stratif ication data and local bathymetry.The ISW propagation velocity decreases as they propagate from deep to shallow water; the maximum propagation velocity of 3.27 m/s was estimated on the western side of the Nicobar Islands and minimum speed of 0.54 m/s occurred in the shallow water region of the southeastern Andaman Sea. The results show that the ISW propagation characteristics diff er in the northern, central, and southern regions of the Andaman Sea. In the northern Andaman Sea, the velocity of ISWs propagating westward was greater than that of ISWs propagating eastward at the same water depth. In the central Andaman Sea, the propagation velocity of the ISWs diff ered over a small area at a depth of 2 500 m, and the velocity of ISWs in the deep mixing layer in winter was higher than that in the shallow mixing layer in spring. Monthly variations in ISW propagation velocity were analyzed in the southern Andaman Sea, and the velocity of ISWs diff ered greatly in shallow water and was not signif icantly diff erent in deep water. Water depth and monthly stratif ication play vital roles in controlling the phase speed of ISWs in the Andaman Sea. This study will provide a basis for the propagation and prediction of ISWs in the Andaman Sea.

Keyword: internal solitary waves; propagation velocity; remote sensing; Andaman Sea

1 INTRODUCTION

Internal solitary waves (ISWs) with strong nonlinearity are frequently observed in stratif ied oceans. ISWs are widely distributed in the world’s oceans, especially near rough topography regions in marginal seas and coastal zones. The Andaman Sea in the Indian Ocean has been a classic study region for ISWs for several decades. This region has a special submarine topographic structure and stable stratigraphy, and is considered one of the areas where ISWs are most commonly observed in the world(Jackson, 2007; Magalhaes and Da Silva, 2018).

Internal solitary waves were f irst observed in the Andaman Sea in 1965. Perry and Schimke observed ISWs with amplitudes of 40 m in the southern Andaman Sea, which was the f irst report of ISWs in the Andaman Sea (Perry and Schimke, 1965). Later,papers such as Osborne and Burch (1980) described ISWs that propagated eastwards, separated by distances of approximately 100 km in the southern Andaman Sea, with a maximum amplitude and horizontal velocity of 60 m and 2 m/s, respectively.However, in-situ ISW observations are very sparse in the Andaman Sea. The development of remotesensing technology provides an eff ective means for collecting ISW observational data. Vlasenko and Alpers (2005) and Shimizu and Nakayama (2017)investigated the generation of ISWs and secondary internal waves by using numerical simulations.Studies on the spatial and temporal distribution,generation, propagation, and dissipation mechanisms of ISWs in the Andaman Sea have been reported by using numerical simulations and remote sensing (Sun et al., 2019). Five potential generation locations of mode-1 long living ISWs were revealed in the Andaman Sea using remote sensing observations, and the waves from two potential generation sites between the Nicobar Islands appear to radiate waves in two opposite directions, toward the Andaman Sea and the southern Bay of Bengal (Raju et al., 2019).Unprecedented satellite images showing long-lived short-scale mode-2 internal waves have been documented in the Andaman Sea by nonhydrostatic and fully nonlinear numerical models (Magalhaes et al., 2020). ISWs in the Andaman Sea were remotely generated through energetic surface tide-topography interactions (Apel et al., 1985). Simultaneously, data from geostationary orbit satellites, such as GOCI(Kim et al., 2018), Himawari-8 (Lindsey et al., 2018),and the infrared radiometer VIIRS (Hu et al., 2019),have also been widely used to study ISWs. Satellite remote sensing images can be used not only to study the characteristics of ISWs, but also to invert the amplitude, wavelength, propagation velocity, and other parameters of ISWs (Yaroshchuk et al., 2016;Zhang et al., 2016; Sun et al., 2018).

As one of the most important parameters of ISWs,the propagation velocity is aff ected by various factors,such as the topography, pycnocline, and tidal currents,which can directly ref lect the propagation characteristics of ISWs. Semidiurnal tides are the main driving force of ISWs on the continental shelf.Based on the distance between two or more packets of ISWs in a single image, the propagation velocity of ISWs can be calculated by assuming that the ISW group has the same period (12.42 h) as the semidiurnal tide (Porter and Thompson, 1999; Li et al., 2000;Zhao et al., 2004). However, this method ignores the propagation direction of ISWs, and ISWs will travel for several hundred kilometers during the long halfday tidal period. With the increasing enrichment of satellite remote sensing data, ISWs can be captured from two adjacent remote sensing images of the same location, and the propagation velocity of the ISWs can be calculated from the space displacement and time interval (Liu et al., 2014; Hong et al., 2015). This method has higher digital accuracy than the singleimage method (Jackson et al., 2013). Although the existing satellite data are abundant, it is not easy to acquire and capture two sets of ISW images in the same scene acquired in a short time. This method is restricted by the remote sensing image resolution,satellite observation geometry, solar azimuth, altitude,and data acquisition and many other conditions play an important role (Jackson and Alpers, 2010). The time lag between the two sets of images used to calculate the propagation velocity should not be too short. There are many mixed pixels in remote sensing images with low spatial resolution due to the limitation of the resolution of the satellite sensor.This paper presents results on the estimation of the propagation velocity of ISWs in the Andaman Sea tracked using 195 image pairs of MODIS imagery over the period from January 2014 to December 2018.A total of 562 ISWs were identif ied during the period and the results of the propagation velocity distribution of ISWs in the Andaman Sea are presented. This paper introduces satellite remote sensing imagery covering 5 years between January 2014 and December 2018, which is used to study the propagation velocity of ISWs in the Andaman Sea (Section 2). The method for calculating the propagation velocity of ISWs by remote sensing imagery and the method for theoretically calculating the propagation velocity of ISWs theoretically, are then presented (Section 3).Section 4 presents the analysis of the propagation velocity characteristics of ISWs in the Andaman Sea.Section 5 is the discussion. The last section presents the conclusion.

2 SATELLITE OBSERVATION DATA

MODIS is one of the primary sensors onboard the Terra and Aqua satellites, which were launched in 1999 and 2002, respectively. MODIS is a passive moderate resolution imaging spectroradiometer with 490 detectors in 36 spectral bands ranging from 0.4 microns (visible light) to 14.4 microns (thermal infrared) with a maximum spatial resolution of 250 m and a scanning width of 2 330 km. The Terra and Aqua satellites pass through the equator at approximately 10꞉30 and 13꞉30 local time each day,respectively. The two satellites pass through the same position on the same day with an approximately 3-h time interval. In view of the unique orbital and band characteristics of Terra/Aqua MODIS, the advantages of obtaining MODIS images of the same location more than once a day are presented, and the propagation velocity of ISWs in the Andaman Sea is studied in this paper.

Fig.1 MODIS image of ISWs in the Andaman Sea

To estimate the propagation velocity of ISWs, we need a pair of images that capture the same packet at two diff erent times. A total of 195 such image pairs with suitable combinations are selected in the Andaman Sea spanning 5 years between 2014 and 2018. The geo-locations of the leading wavefront of the ISW packet are traced by taking the central pixels of bright patches in the imagery. An ensemble of all such leading wavefronts is combined to produce a complete picture of ISWs over the Andaman Sea, as shown in Fig.1a. The red curve represents the ISW wavefront from Terra-MODIS imagery, and the blue curve represents the ISW wavefront from Aqua-MODIS imagery. The ISWs are mainly distributed in the southern, central, and northern parts of the Andaman Sea. The propagation characteristics of ISWs are observed diff erently in the northern, central,and southern regions of the Andaman Sea.

Fig.2 ISW images acquired by Terra-MODIS (13 Mar. 2017 at 04∶00 UTC) and Aqua-MODIS (13 Mar. 2017 at 07∶00 UTC)in the southern Andaman Sea

3 METHOD

3.1 Multitemporal image comparison method

The propagation velocity of ISWs was determined by using two diff erent satellites for the multitemporal image (MTI) method, and the leading crests in two diff erent groups of ISWs were used for the tidal period image (TPI) method (Hong et al., 2015). The TPI method is mainly applicable to tidal ISWs, but has a slightly decreased accuracy. This result is because ISWs will change greatly under the inf luence of the external environment, such as the terrain and f low f ield, during long-distance propagation during the long tidal cycle. This paper studies the propagation velocity of ISWs in the Andaman Sea based on the MTI method using MODIS onboard National Aeronautics and Space Administration (NASA) Terra/Aqua satellites.

The time interval between Terra and Aqua MODIS images was approximately three hours. During this interval, the ISWs could travel distances ranging from several kilometers to tens of kilometers. Therefore,the propagation velocity of ISWs isv=l/t, wherelis the propagation distance of the ISW andtis the propagation time. Figure 2 shows the Terra and Aqua MODIS remote sensing images acquired on 13 March 2017, showing ISW signatures in the southern Andaman Sea. Only one of the wave packets was intercepted to clearly show the ISW bands in the remote sensing images. The coverage areas of the two images are displayed as the solid black box in Fig.2a.Figure 2b & c shows the MODIS images acquired on 13 March 2017 at 04꞉00 and 07꞉00 UTC, respectively.To accurately obtain the propagation distance of ISWs acquired on Terra and Aqua MODIS images, the coverage areas used to calculate the propagation distance of ISWs are displayed as solid red boxes.Figure 3 shows the intensity prof iles across the transect lines in the area displayed as the solid red box in Fig.2b & c. The blue curve represents the ISWs in the Terra-MODIS image, and the red curve represents the ISWs in the Aqua-MODIS image. The distances between the Terra and Aqua MODIS images were measured, as shown in Fig.3, and the propagation velocity of ISWs was computed using a data acquisition interval of 3 h.

3.2 Theoretical phase velocity of ISWs

The famous nonlinear equation (Korteweg-de Vries (KdV) equation) describing the solitary wave characteristics was derived by Korteweg and de Vries(1895). Thus, the characteristics of a single soliton in the two-layer ocean model can be described by the KdV equation in dimensional form. The propagation velocity of ISWs can be determined theoretically by solving the KdV equation using water depth, ocean stratif ication, and background current information.The formula is as follows:

Fig.3 Cross sections of the Terra-MODIS image (blue line) and the Aqua-MODIS image (red line) brightness (ref lectance)

where,tis time,ηis the amplitude of ISWs, andC0is the linear internal wave velocity, which is also a linear term.αandβare the f irst order nonlinear term coeffi cient and the dissipative term coeffi cient,respectively.

The shallow water term and dissipative term can be ignored if the viscous eff ect of seawater and the frictional eff ect of sea bottom are not considered. The solution for the above KdV equation is as follows:

whereη0is the amplitude of ISWs,Cpis the phase velocity of ISWs, andlis the half-wavelength of the ISWs.

For a two-layer mixed ocean system, the phase velocityCpand half-wavelengthlare:

4 RESULT

4.1 Propagation velocity distribution of ISWs

A total of 195 MODIS image pairs were acquired in the Andaman Sea in January 2014 and December 2018 from the NASA satellites Terra and Aqua. The morphology of 562 ISWs in the 195 image pairs was added to GIS software as shown in Fig.4a. The crest lines of the diff erent colors in Fig.4a represent the ISWs detected by the diff erent satellites. The red curves show the ISWs acquired by Terra, and the blue curves show the ISWs acquired by Aqua. The ISWs detected by Terra and Aqua on the same day were mainly distributed in the southern, central, and northern parts of the Andaman Sea, and most ISW pairs were detected in the southern Andaman Sea.The entire research area of the Andaman Sea was divided into 0.5°×0.5° cells, and the average propagation velocity of ISWs in each cell was calculated. The spatial distribution of the ISW propagation velocity is shown in Fig.4b. The propagation velocity of ISWs in the deep sea of the central Andaman Sea was relatively fast, with a rate of approximately 3 m/s. In the southern and central parts of the Andaman Sea, the propagation velocity of ISWs decreases gradually in the process of transmission from west to east, and the velocities of ISWs in the shallow waters of the eastern coast are less than 1 m/s. The ISWs in the northern Andaman Sea are divided into two types. The f irst type of ISW was generated in the waters around the northern Andaman Sea islands and propagated eastward into the Andaman Sea. The ISWs of this type were signif icantly aff ected by topographic changes. Along the propagation direction of the ISWs, the water depth was shallow in the north and relatively deep in the south, which resulted in the propagation velocity of ISWs in the south being obviously higher than that in the north. Another type of ISW originates at the continental slope of the northern Andaman Sea and propagates southwest toward the Andaman Islands.The propagation velocity of ISWs was signif icantly reduced when it reached the sea area around the Andaman Sea islands.

Fig.4 Distribution of the ISW crest line of the leading wave (a), propagation velocity (b), and theoretical propagation velocity (c)

Figure 4c shows the diagrams of propagation velocity calculated by Eq.5. The propagation velocity of ISWs obtained by remote sensing corresponded with the results of the theoretical calculations. The propagation velocity of ISWs was almost consistent with the variation in water depth. The ISWs were distributed in water depths ranging from tens of meters to thousands of meters in this paper, and the propagation velocity of the ISWs increased gradually with the depth of water, as shown in Fig.5, in which the propagation velocity was calculated by the 562 ISWs in the 195 image pairs in January 2014 and December 2018 from Terra and Aqua. The propagation velocity of ISWs was below 1 m/s when the water depth was less than 100 m, and the propagation velocity of the ISWs increased signif icantly with increasing water depth within the range of 200 m. When the water depth was below 500 m, the propagation speed of ISWs was less than 2 m/s. The propagation velocity of ISWs can reach more than 3 m/s when the depth is greater than 3 000 m. Since the upper layer depth does not vary greatly, the total water depth can show a strong correlation with ISW velocity. In general, the representative water depth is an essential factor.

Fig.5 Relationship between the velocity of ISWs and the depth of water in the Andaman Sea

Fig.6 Remote sensing images of ISWs matched by Terra and Aqua MODIS in the study area

4.2 Propagation velocity of ISWs at the same depth

The propagation velocity of ISWs was aff ected by various factors; diff erent ISWs may not have the same velocity at the same location, and the propagation velocity of ISWs at diff erent times at the same location may not be diff erent either. The propagation velocity of ISWs constantly changed throughout the life cycle.The propagation velocity of ISWs at the local scale was studied, where the central longitude and latitude were 95.3°E and 9.65°N, respectively. The water depth in the study area was approximately 2 500 m.Ten pairs of images were obtained that simultaneously detected ISWs spanning 5 years between 2014 and 2018 (Fig.6b-k). The coverage areas of the study area are displayed as solid black boxes in Fig.6a.

The propagation velocity of ISWs in the study area can be obtained by using the above 10 pairs of remote sensing images, as shown in Table 1. Due to the depth of water in this region, the propagation velocity of ISWs was relatively fast, ranging from 2.31 to 3.27 m/s. The propagation velocity of ISWs varies greatly in areas with the same water depth. The maximum propagation velocity of 3.27 m/s was estimated on 24 January 2015, and the minimum propagation velocity of the ISWs of 2.31 m/s was obtained on 16 March 2018. Figure 7 shows the average monthly propagation velocity of the ISWs in the study area. The propagation velocity of ISWs varies at the same water depths and shows certain seasonal variations, in which the maximum propagation velocity of ISWs occurs in January.

Fig.7 Average monthly velocity of ISWs in a small region of the central Andaman Sea

The propagation velocity of ISWs was inf luenced not only by topographic depth but also by the stratif ied structure of the ocean. The change in season causes a change in sea temperature, followed by a change in the depth of the ocean mixing layer, which aff ects the propagation velocity of ISWs. Figure 8 shows the seasonal variation curves of temperature, salinity,density, and buoyancy frequency in the study area.Diff erent colors represent the temperature, salinity,density, and buoyancy frequency curves of diff erent months. The propagation velocity of ISWs in the deep mixing layer in January was higher than that in the shallow mixing layer in April. This result shows that monthly stratif ication plays a vital role in controlling the ISWs in the Andaman Sea.

4.3 Propagation velocity of ISWs at diff erent depths

4.3.1 The propagation velocity of ISWs in the northern Andaman Sea

There were two types of ISWs in the northern Andaman Sea. The f irst type of ISW was generated in the northern Andaman Sea islands and propagated eastward into the Andaman Sea, eventually dissipating on the eastern continental shelf of the Andaman Sea. Another type of ISW originates at the continental slope of the northern Andaman Sea and propagates southwest toward the Andaman Islands,eventually dissipating on the continental shelf near the northern Andaman Islands. Figure 9a shows a MODIS image acquired on 19 March 2014 containing

Table 1 Propagation velocity of ISWs as measured at the same locations from 10 pairs of Terra and Aqua MODIS images

ISWs with two propagation directions. Three groups of ISW packets that were seaward propagate toward the northern Andaman Sea. Three groups of ISW packets propagating southwest toward the Andaman Islands were visible in the sea. Twenty pairs of remote sensing images of ISW in the northern Andaman Sea were used to calculate the propagation velocity of ISW in both directions, as shown in Fig.9b. The red dot indicates the propagation velocity of ISWs traveling southwest, and the blue box indicates the propagation velocity of ISWs traveling southeast in the northern Andaman Sea. The propagating velocity of the ISW in the west direction is signif icantly higher than that in the east direction where the water depth is less than 500 m. The propagation velocity of the ISW is similar at water depths of 1 000-2 500 m, but overall, the velocity of the ISW propagating westward was higher than that of the ISW propagating eastward at the same depth in the northern Andaman Sea.

Fig.8 Vertical temperature, salinity, density, and buoyancy frequency prof iles for the diff erent seasons

Fig.9 The MODIS image of the northern Andaman Sea acquired on 19 March 2014 (a); the velocity of ISWs in the northern Andaman Sea (b)

The propagation velocity of ISWs increased rapidly at depths below 1 000 m, while the propagation velocity did not increase signif icantly at depths between 1 000 m and 2 500 m. The propagation velocity was mainly concentrated in the range of 2.0-2.5 m/s. The velocity of the ISWs propagating westward was obviously larger than that of the ISWs propagating eastward at the same depth in the northern Andaman Sea. We believe this diff erence was probably caused by the change in topography.Figure 10 displays the change in water depth along the direction of the principal axis of ISW propagation.The blue curves show the ISWs propagating southeast toward the continental shelf, and the red curves show the ISWs propagating southwest toward Andaman Island. ISWs propagating southwestward were generated at the slope break of the northern Andaman Sea and then propagated to the deep sea in the southeast direction; the depth did not change signif icantly along the crestline of ISWs. ISWs propagating eastward were generated in the North Andaman Islands and propagated southeast into the Andaman Sea, and the topography along the crestline of ISWs changed signif icantly. The depth in the northern sea was shallow, and the water depth in the southern sea was relatively deep. This diff erence leads to a large diff erence in the propagation velocity of the same ISWs, which was fast in the south and slow in the north, thus aff ecting the overall propagation velocity of the ISWs.

Fig.10 Topography prof iles along the propagation direction of ISWs in the northern Andaman Sea

Fig.11 Variation in propagation velocity and water depth of ISWs in the southern Andaman Sea

4.3.2 The velocity of ISWs in the south of the Andaman Sea

The ISWs in the southern Andaman Sea were mainly generated in the ridge near the Gretel Strait and propagated eastward into the Andaman Sea. The black rectangle shown in Fig.11a is the study area.Along the propagation direction of ISWs, the water depth gradually decreases, and the propagation velocity of ISWs gradually increases with increasing water depth. When the water depth was more than 1 500 m, the propagation velocity of ISWs did not increase with increasing water depth but tended to remain stable and even decreased. This f inding was mainly due to the ISWs produced in the gray ridge near the strait, which spread eastward into the Andaman Sea. ISWs f irst spread into the sea area that was more than 2 000-m deep; the ISWs were still in the process of growth in this phase, and the growth process of ISWs strengthened. However, the propagation velocity was not the largest in the deepsea portion of the study area. As the ISWs continue to travel eastward, their strength gradually increases,and their propagation velocity reaches a maximum of 2.8 m/s at depths of 1 200-1 500 m. ISWs continue to travel to the eastern coast, their velocity gradually decreases with decreasing water depth, and the minimum value is approximately 0.54 m/s.

Fig.12 The propagation velocity of the ISWs in the southern Andaman Sea varies with the water depth from 1 to 12 months

The propagation velocity of the ISWs was calculated based on the Terra/Aqua MODIS remote sensing images covering the southern Andaman Sea from January 2014 to December 2018. The average propagation velocity of ISWs was calculated by month and water depth (Fig.12). The diff erent colored columns represent the propagation velocity of the ISWs in diff erent months. The propagation velocity of ISWs in diff erent months varied greatly with the change in water depth, and the propagation velocity of ISWs in shallow water areas varied greatly. For water depths below 100 m, the velocity of ISWs can reach 1.9 m/s in August, while the average velocity in May is only 0.54 m/s. At water depths of 1 300-1 700 m, the maximum velocity of ISWs was 2.8 m/s in May. In general, the average velocity of ISWs increases with increasing water depth, as shown in the yellow column in Fig.12.

5 DISCUSSION

In this paper, the propagation velocity of ISWs in the Andaman Sea was demonstrated in detail, rather than just individual cases in a typical region. Although the existing satellite data are abundant, it is not easy to acquire and capture two sets of ISW images in the same scene. This is restricted by the remote sensing image resolution, satellite observation geometry,solar azimuth, altitude, and data acquisition and many other conditions play an important role. There are many mixed pixels in remote sensing images with low spatial resolution due to the limitation of the resolution of the satellite sensor. The time lag between the two sets of images used to calculate the propagation velocity should not be too long or too short. If the time is too long, internal solitary waves will be aff ected by the external environment such as topography during their propagation and the propagation direction will also change, so that the propagation velocity calculated is not accurate. If the time lag is too short, the ISWs travel a distance of only a few pixels, this will cause large errors to the calculation of propagation velocity. Terra and Aqua are two satellites that cooperate with each other in time and space, and two images were eff ectively obtained at the same point on the same day, with a time interval of approximately 3 h. The remote sensing images used in this paper cover the Andaman Sea every month and the typical occurrence area of ISWs, which allowed the propagation velocity of ISWs to be detected.

However, only 195 pairs of MODIS image were matched in the Andaman Sea in January 2014 and December 2018. Thus it can be seen the data acquisition is valuable. Scanning width of the data used in this paper is 2 330 km, which can obtain a complete wave crest line with a length of hundreds of kilometers. The spatial resolution of the data used is 250 m, the position error of the extracted internal solitary wave in the remote sensing image is about 1 pixels. This distance is miniscule compared to the long distance of the internal solitary wave in three hours. Therefore, the position error of the extracted internal solitary wave has little inf luence on the calculation results of propagation velocity.

The propagation velocity of ISWs has a certain regional character and aff ected by various factors,such as the topography, pycnocline, and tidal currents.The results in this paper were consistent with existing reports that water depth was one of the important factors aff ecting the propagation of ISWs. Overall,the propagation velocity of ISWs increases with increasing depth. It was found that the propagation velocity of ISWs diff ered throughout the same range.A 0.7°×0.7° area in the central Andaman Sea was taken as the study area, where the water depth was approximately 2 500 m. The research shows that the propagation velocity of ISWs within this range exhibited great diff erences. In the season with a deep mixing layer, the propagation velocity of ISWs within this range was relatively large. This result indicates that the depth of the pycnocline was another important factor aff ecting the propagation velocity of the ISWs.Under the inf luence of the external environment, the propagation velocity of the ISWs varies seasonally,where the propagation velocity of ISWs in shallow water areas exhibits signif icant seasonal variation mainly because the shallow water area is greatly aff ected by the light and the external environment,and the seasonal thermocline changes signif icantly,while the thermocline in the deep water area is relatively stable.

For the ISWs in the northern Andaman Sea, there are two types of ISWs. The f irst type of ISW was generated in the northern Andaman Sea islands and propagated eastward into the Andaman Sea, eventually dissipating on the eastern continental shelf of the Andaman Sea. Another type of ISW originates at the continental slope of the northern Andaman Sea and propagates southwest toward the Andaman Islands,eventually dissipating on the continental shelf near the Northern Andaman Island. The velocity of the ISWs propagating in the west direction was obviously higher than that of the ISWs propagating in the east direction at the same water depth, which was mainly aff ected by the change in topography. This diff erence was mainly because the ISWs propagating to the west were located in the deep sea along the direction of the wave crest line of the ISWs, where there were limited changes in water depth. For the ISWs propagating eastward, along the direction of the wave crest line of the ISWs, the shallow sea continental shelf area was located to the north, and the deep-water area was located to the south. This condition results in an ISW traveling at diff erent speeds from north to south, and then the whole propagation velocity of the ISWs is aff ected.

The propagation velocity of ISWs plays a very important role in the study of their propagation characteristics and prediction, and is also one of the key parameters in the numerical simulation of ISWs.Therefore, the propagation velocity of ISWs has important academic research value. The results of the study on the propagation velocity of ISWs will provide a valuable scientif ic basis for the study of the propagation and prediction of ISWs in the Andaman Sea. Unfortunately, due to the limitation of satellite remote sensing data and the propagation characteristics of the ISWs themselves, the propagation velocity of the ISWs on the continental shelf of the shallow sea could not be detected. The detection of the velocity in these areas requires satellite remote sensing data with higher resolution or results obtained by numerical inversion.

6 CONCLUSION

The Andaman Sea has been a classic study region for ISWs because of its special geographical location and climate. As one of the important parameters of ISWs, propagation velocity has important research signif icance. The propagation velocity of the ISWs was calculated based on 195 MODIS image pairs acquired in the Andaman Sea in January 2014 and December 2018 from the NASA satellites Terra and Aqua. A total of 562 ISWs were identif ied during the period, and the results of the propagation velocity distribution of ISWs in the Andaman Sea were presented. The velocity distribution of ISWs in the Andaman Sea was plotted using a 0.5°×0.5° grid, and the distribution of the propagation velocity of ISWs in the Andaman Sea was obtained. The propagation velocity distribution of ISWs in the Andaman Sea was consistent with the distribution of water depth. The propagation velocity of ISWs in the deep area in the central and western regions was relatively high. The maximum propagation velocity of 3.27 m/s was estimated on the western side of the Nicobar Islands,and the minimum speed was 0.54 m/s in the shallow water region of the southeastern Andaman Sea. ISW propagation characteristics were observed to diff er in the northern, central, and southern regions of the Andaman Sea. The propagation velocity in this paper focuses on the three main occurrence regions of ISWs in the Andaman Sea.In the northern Andaman Sea, the velocity of the ISWs propagating westward was greater than that of the ISWs propagating eastward at the same water depth in the northern Andaman Sea. In the central part of the Andaman Sea, the propagation velocity of ISWs was studied in a 0.7°×0.7° grid area with a water depth of 2 500 m. The propagation velocity of ISWs in the deep mixing layer in winter was higher than that in the shallow mixing layer in summer. In the southern Andaman Sea, the monthly variations in ISW propagation velocity were analyzed. The velocity of ISWs diff ered greatly in shallow water and was not signif icantly diff erent in deep water. In general, the average velocity of ISWs increases with increasing water depth. The results indicate that water depth and monthly stratif ication play a vital role in controlling the propagation velocity of ISWs in the Andaman Sea.

7 DATA AVAILABILITY STATEMENT

Sequence data that support the f indings of this study have been deposited in the NASA Worldview application (https://worldview.earthdata.nasa.gov).

8 ACKNOWLEDGMENT

We thank the NASA MODIS Rapid Response System website announced that swath images.