| ||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Mesoscale studies in the Indian Ocean using altimeter data
ABSTRACT
1.0 INTRODUCTIONThe circulation patterns in the Indian Ocean, forced by the reversing southwest (May to August) and northeast (October to January) monsoonal winds, represent a unique oceanographic problem that has been the subject of a number of investigations [e.g. Cox, 1976, Anderson et. al., 1976, Rao et. al., 1989, Bruce et. al., 1994 Schott et. al., 1994 Evensen et. al., 1994 Shetye et. al., 1994 ]. The South West (SW) monsoon is normally observed over the Arabian Sea from May to August and the North East (NE) monsoon from October to January Indian Meteorological~Society, 1994. The surface circulation in the Indian Ocean is schematicly represented in Figure 1 for the two monsoon periods. During the NE monsoon there is a strong southward flowing current off the northeast African coast called the Somali Current. During the SW monsoon the westward flowing Southern Equatorial Current (SEC) feeds the Somali Current (SC) which then flows northward with speeds up to 200~cm/s Pickard et. al., 1990, and turning offshore at 2 -3 deg. S. The SC develops two anticyclonic gyres, the Great Whirl (GW) (centered at about 8 deg. N) and the Southern Gyre (SG) (located in the area 0 - 5 deg. N, 53 deg. E) [ Perigaud et. al., 1988 Shott et. al., 1994 ]. Significant mesoscale variability has also been observed in the northern Indian Ocean and a number of eddies have been examined by previous investigators. Satellite altimetry is particularly suited for studying mesoscale variability and has been used in the Southern Indian Ocean by investigators such as Wakker et. al., 1990, Snaith et. al., 1996, Grundling, 1995. Evensen et. al., 1995 assimilated Geosat altimeter data in an ocean circulation monsoon to study the shedding of Agulhas rings. Further altimeter studies in the northern Indian Ocean have been done by Perigaud et. al., 1988 who studied the variability of the Somali Current using Seasat altimeter data. Perigaud et. al., 1992, Perigaud et. al., 1993 studied the annual/interannual variability of currents in the area using Geosat altimeter data. Bruce et. al., 1994 studied the formation of a large anticyclonic eddy in the eastern Arabian Sea during the NE monsoon using Geosat altimeter data, and Greiner et. al., 1993 assimilated Geosat sea level variations into a nonlinear reduced gravity shallow-water model. However, so far the extensive data set from ERS-1 and Topex/Poseidon altimeters have not been used jointly in the Indian Ocean circulation. Therefore the main objective of our investigation was to gain an improved understanding of the oceanic variability in response to the monsoon system in the northern Indian Ocean using altimeter data. The altimeter data used and the method of processing are described in Section 2. In Section 3 the details of the SSH variability field from the altimeters is examined, and the temporal and spatial characteristics obtained from Complex Principal Component Analysis (CPCA) are presented. In Section 4 results pertaining to the propagation and life time of eddies are presented, and in Section 5 a discusion and conclusion is given.
(a) SW monsoon (August) (b) NE monsoon (January) Figure 1. Surface circulation during the SW and NE monsoon. Modified from Molinari et. al., 1990. Somali Current (SC), Great Whirl (GW), Indian Monsoon Current (IMC), Equatorial Jet (EJ), Southern Equatorial Current (SEC), Northern Equatorial Current (NEC), Southern Equatorial Counter Current (SECC), East African Counter Current (EACC), Coastal Current (CC) (Units in cm). 2.0 ALTIMETER DATA PROCESSINGThe altimeter data used in the study are from the Geosat (17 days repeat period, 01.11.86 - 01.09.89), the ERS-1 (35 days repeat period 01.04.92 - 31.12.93) and the Topex/Poseidon (10 days repeat period, 01.01.93 - 10.08.95) missions. The first 10 cycles of Topex/Poseidon data are not used in this study, because of attitude problems with the satellite AVISO, 1992 . The data were processed using standard repeat track analysis to obtain SSH anomalies [Cheney et. al., 1991, Samuel, 1993]. For each data point, Root Mean Square (RMS) of the SSH were computed. The SSH anomalies are also interpolated to a regular space-time grid using a simple exponential weighting scheme [Samuel, 1993]. Geosat data are gridded on a 1.0 x 1.0 deg. spatial grid at 8 days intervals, ERS-1 on a 0.5 deg. x 0.5 deg. spatial grid at 15 days time intervals and Topex/Poseidon on a 2.0 x 2.0 deg. spatial grid at 5 days intervals. These grid spacings and intervals were choses in order to properly present the spatial and temporal sampling characteristics of the different satellite missions. The spatial/temporal modes of the data sets were extracted using Complex Principal Component Analysis (CPCA) on the gridded data sets. CPCA is a variant of Principal Component Analysis which is used to decompose the variability in the data sets into orthogonal modes, and is designed specifically to handle propagating features [Horel, 1984, Preisendorfer, 1988]. CPCA is done only for Topex/Poseidon and not for ERS-1 since the duration is too short (one and a half years). (a) Geosat (b) ERS-1 (c)Topex/Poseidon Figure 2. RMS Sea Surface height anomaly distribution in the Northern Indian Ocean from Geosat, ERS-1 and Topex/Poseidon. (Units in cm) 3.0 SPATIAL AND TEMPORAL VARIABILITY3.1 Root Mean Square variabilityRoot Mean Square (RMS) of the SSH provides an indication of
the distribution of mesoscale variability in the study area. Over
most of the study area the RMS variability is low, less than 7 cm
(Figure 2). High RMS variability is observed in highly energetic
areas such as the Somali Current area (15 - 17 cm), and along the
western part of the Bay of Bengal (13 - 15 cm). High variability
is also observed along 10 deg. S at 45 - 50 deg. E and at 75 - 82
deg. E (11 - 13 cm) and in the eastern part of the Bay of Bengal
(9 - 11 cm). The high RMS values are about half of those observed
in highly energetic areas such as the Gulf Stream or the Agulhas,
where similar analysis using altimeter data have shown RMS values
up to 30 cm [Le Traon et. al., 1990]
and 32 cm [Wakker et. al., 1990].
In the Nordic Sea, [Samuel et. al.,
1994] observed values with a maximum of 15 cm. There is a
close relationship between the areas of greatest SSH variability
and the major currents, the Somali Current, East Indian Coastal
Current [McCreary et. al., 1996;
Shankar et. al., 1996] and
Southern Equatorial Current, and higher variability is most
prevalent in the western parts of the oceans. As an overall
assessment the distribution of RMS height in the Indian Ocean
seems to be consistent for the three data sets. In general the
ERS-1 shows more extended areas with high variability (9 - 11
cm), which is particularly seen over the Chagos Laccadive Ridge
(north-south ridge, extending between 73 deg. - 78 deg. E). Other
areas with high RMS distribution are around 78 deg. E, 10 deg. S
where the Southern Equatorial Current is found and around 47 deg.
E, 10 deg. S which is in the area of the East African Counter
Current . Geosat indicates values of 11 - 13 cm in these areas,
while the other two altimeters show values of 7 - 9 cm. In the
eastern part of the Bay of Bengal values of 13 - 15 cm are shown
for all three altimeters. 3.2 CPCAComplex principle Component Analysis (CPCA) decomposes the data set into a number of orthogonal modes, each of which are represented by a complex time series and complex spatial eigenvector. The time series indicate the temporal evolution of the spatial averaged variability, while the spatial patterns (eigenvectors) indicate the distribution of temporaly averaged variability. The modes are ordered according to the associated eigenvalues which give an indication of the fraction of the total variability that is explained by the respective modes. Usually major portion of the variability will be contained in the first modes, and therefore only these modes need to be explained. In our case, the first two modes of the Geosat and Topex/Poseidon data set accounted for 28.3 % and 41.0 % respectively, and these will be discussed. Figure 3 shows the phase and the amplitude time series of the first two modes for Geosat and Topex/Poseidon. Dominant modes for both the data sets have a clear annual period as seen in the phase plots Figure 3a and 3c. The phase plot for the Geosat first mode shows a clear annual signal while for Topex/Poseidon there is some leakage from shorter periods. Both the data sets indicate that the annual signal has maxima around January and July, with small variations from year to year, possibly related to variability in the onset of the monsoon.
Figure 3. Temporal characteristics for first and second Geosat and Topex/Poseidon SSH CPCA mode. The amplitude is in cm, and the phase is from -pi to +pi.} The amplitude time series (Figure 3b and 3d) indicate when the annual signal is strongest. For Geosat the strongest is January 1989 and also reluctivly strong in January 1987 and 1988, May 1988 and 1989. The amplitude time series for Topex/Poseidon mode one (Figure 3d) suggests that this mode actually contains two signals, the annual signal and a signal with a period of approximately two months. As mentioned earlier, this is also evident in the corresponding phase plot (Figure 3c). The fact that the two signals are combined in one mode suggests that they are correlated. Several investigators have observed oscillations with a 40--60 days period in the western Indian Ocean using in situ data and numerical models [Mysak et. al., 1984; Schott et. al., 1988; Swallow et. al., 1988; Woodberry et. al., 1989] and explain these oscillations as caused by horizontal shear (barotropic instability) which generate eddies at a period of 40--50 days. This signal is possibly not resolved by the 17 days repeat period of Geosat, explaining the absence from Figure 3a and 3b. The phase time series of the second mode for Geosat (Figure 3e) show a clear semi-annual period with maxima in May, August and November, February. This mode therefore appears to be related to the onset and decay phase of the SW and NE monsoons. For the Topex/Poseidon phase no clear semi annual period is observed (Figure 3g). There is a strong indication of a half yearly period from May'93 to November'93, but for the rest of the sampling period, another phase with a period of approximately 2 months similar to that above, overrides the semi-annual signal. The amplitude (Figure 3h) show that the period when the semi-annual phase is observed is weeker than the 2 months period observed. In order to examine the significance of the peaks and lows in the temporal characteristics, and to describe the varying physical attributes of propagating features in traveling wave fields, a truncated data set may be reconstructed using the first two modes. The results will be discussed in section 3.3. Figure 4. SSH anomaly field reconstructed from the first and second Geosat CPCA mode during the NE monsoon. 3.3 Large scale features.A reconstruction using the annual and semi-annual mode for Geosat and Topex/Poseidon was done to describe the varying physical attributes of propagating features. Only the anomaly fields during the NE monsoon in 1988 for Geosat (Figure 4) are shown, since the propagation is most easily seen here. The annual mode shows a feature, (A) which is observed in the Somali Current area (Figure 4) during the 1987/88 NE monsoon from the Geosat data. No clear evidence of direction of propagation is observed. In the Bay of Bengal a rotational movement (B) is observed, but it is not possible to determine the direction, but we assume that this is an anticyclonic gyre due to analysis of results from a multi-layer, adiabatic model, driven by climatological monthly mean winds winds [Poterma et. al., 1991] and ship drift data [Defant, 1961; Cutler et. al., 1984] which have observed a large anticyclonic flow in the Bay of Bengal. For the semi-annual mode it appears that the main phenomena is Rossby waves. Our investigation identifies Rossby waves along 7 deg. N (C), along equator (D) and between 5 deg. S and the equator (E). The westward propagating waves have a speed up to 6 km/d along 7 deg. N (C), up to 13 km/d along equator (D) and 7 km/d between 5 deg. S and the equator (E), indicating a decrease in speed away from the equator. Our results show generally lower speeds than Perigaud et. al., 1993 interannual sea level studies, using Geosat data based on crossover analysis and shallow water simulation. Perigaud et. al., 1993 observed Rossby waves which takes from 4 months to cross the basin at equator, and have a speed up to 12 km/d along 15 deg. S. Although Topex/Poseidon accounts for more variability than the two Geosat modes, the propagating waves are easier observed using Geosat altimeter data, mainly due to a better spatial sampling. The same features where also observed from the Topex/Poseidon data (not shown here) with the same characteristics. 3.4 MESOSCALE EDDIESCloser inspection of the gridded SSH anomaly data set showed evidence of several mesoscale eddies. In this section some of the observed eddies will be examined. Figure 5 shows the location and rotation of the most persistent eddies observed from the three data sets. We observed eddies to be more dominant during the monsoon periods, and the two areas, the Somali Current area in the Arabian Sea (30 deg. - 61 deg. E, 0 - 20 deg. N) and the western part of the Bay of Bengal (75 - 95 deg. E, 0 - 20 deg. N) examined, indicate eddies with a maximum amplitude up to 45 cm in July/August and December/January, with a spatial scale ranging between 250--1000 km. The typical life cycle of these eddies is illustrated by examining consecutive time steps of the gridded Geosat, ERS-1 and Topex/Poseidon data in the area of the Somali Current. The growth of an anticyclonic eddy E1, located approximately at 8.5 deg. N, 53 deg. E, is observed from June to September 1989 from the data sets. This eddy is most intense during the first half of August, with amplitudes of 30 - 35 cm (some of the individual passes over this eddy show amplitudes up to 45 cm, before interpolation), with a spatial extent up to 750 km. This compares well with the in situ data from Fischer et. al., 1996. This eddy propagates northeastward with a speed up to 17 km/d. The anticyclonic eddy E2, located approximately at 13.5 deg. N, 56.5 deg. E, is also observed from June to September 1989, and is more intense at the end of August, with amplitudes of 20 - 25 cm, and a spatial extent ranging from 250 - 550 km. This eddy propagated southeastward with a speed up to 7 km/d, and in September it seems to coalence with E1. The third eddy observed is the cyclonic eddy E3, located at approximately 9 deg. N, 51 deg. E, with maxima of 25 - 30 cm at the end of June and during the first half of August. This cold wedge occurs as a result of the two anticyclonic eddies E1 and E2 and defines the so called northern wall [Tsai et. al., 1992; Luther et. al., 1985]. No clear evidence of propagation is seen until August where a slight northward propagation, with a speed up to 17 km/d is observed. The life time of these eddies is observed to be more than three months. Eddy E1 is assumed to be the eddy called the Great Whirl [Perigaud et. al., 1988; Fischer et. al., 1996] and E2 the northern Socotra eddy [Ali, 1990]. The anticyclonic eddy E4 and the cyclonic eddy E6 are observed during the NE monsoon in the same position as E1 and E3. Geosat show that E4 is most intense (30 - 35 cm) during the first half of December in 1988, while the ERS-1 and Topex/Poseidon data show that this eddy is most intense during the beginning of October 1993 (25 - 30 cm) and 1994 (16 - 20 cm) respectively. Geosat show a southward propagation up to 6 km/d, while there is no clear evidence of propagation from the other two data sets. Geosat show that the cyclonic eddy E6 is most intense at the end of December 1993 (25 - 30 cm). The ERS-1 data show that this eddy is most intense during October (25 - 30 cm) and in November 1994 (16 - 20 cm) for Topex/Poseidon. All three show a southward propagation up to 6 km/d. In the Bay of Bengal the anticyclonic eddies E10, E11 and a cyclonic eddy E13 were observed during the SW monsoon. These are most intense during the months May and June, with maxima of 20 - 25 cm. E11 has a southward propagation with a speed up to 7 km/d, E13 has a westward propagation with a speed up to 20 km/d, and the northward propagating (11 km/d) anticyclonic eddy E12, was observed only from the Geosat data. During the NE monsoon an anticyclonic northward propagating eddy E14 with a speed up to 27 km/d, and a cyclonic westward propagating eddy E17, with a speed up to 10 km/d is observed. These eddies are most intense during the second half of December 1988, and 1994.
Figure 5. Overview of the location and rotation for eddies observed in the Indian Ocean using altimetry data. 4.0 DISCUSSION AND CONCLUSIONFor the first time SSH data from three altimeter missions are used to study mesoscale variability in the Arabian Sea and the Bay of Bengal. The correlation between the three altimeters is very good. The Root Mean Square of the SSH anomalies indicates where the high energetic areas are. This formed the basis for the further investigation. Using CPCA analysis we were able to study the temporal and spatial characteristics of the monsoon periods. The annual period based on Geosat data indicate maxima in January and June, indicating the time when the monsoon is most intense. The semi-annual period indicate maxima in May, August and November, February, and appears to be related to the onset and decay phase of the SW and NE monsoons. Topex/Poseidon indicate oscillations caused by horizontal shear which generate eddies at a period of 40--50 days [Mysak et. al., 1984; Schott et. al., 1988; Swallow et. al., 1988; Woodberry et. al., 1989]. Westward propagating features, possibly Rossby waves are also identified from the CPCA, along equator with speeds of 13 km/d. Along 4 deg. - 10 deg. N speeds up to 6~km/d are observed, and between 5 deg. S and equator they travel with a speed of 7~km/d. This indicate a decrease in speed away from equator, which compare well with Perigaud et. al., 1992 results. Mesoscale eddies are located, and the propagation, spatial range, amplitude and life time of these are observed during the monsoon periods. These results compare well with previous investigators who have used altimeter measurements [Ali, 1990; Perigaud et. al., 1988; Perigaud et. al., 1992; Greiner et. al., 1991]. The dimensions of the eddies and the horizontal movement are within the expected range of the observations and model simulations. The correlation between the three altimeter data sets used in this investigation is also very good. Small variations occur, due to the satellites different temporal and spatial sampling and due to different errors and corrections. The Great Whirl, the Socotra eddy and a cyclonic eddy are identified in the Arabian Sea during the SW monsoon. During the NE monsoon an anticyclonic eddy, E4 located in the same place as the Great Whirl was observed, and another Socotra eddy. A cyclonic eddy was also observed north of the E4. A total of three anticyclonic eddies and one cyclonic eddy are observed in the Bay of Bengal, and their propagation speed and direction is investigated. The propagation of eddies during the SW monsoon seems to be eastward and throughout the months of the NE monsoon the eddies have a westward propagation, following the atmospheric circulation. The Southern Gyre, which together with the Great Whirl forms the well-known two gyre system in the Somali Current, was the only major eddy not observed in this investigation. Inspite of few investigations in especially the Bay of Bengal investigating the sea surface variability, our results confirms and strengthen previous results using other techniques [Potemra et. al., 1991;Shetye et. al., 1991; Shetye et. al., 1996]. 5.0 References
Keywords: ESA European Space Agency - Agence spatiale europeenne, observation de la terre, earth observation, satellite remote sensing, teledetection, geophysique, altimetrie, radar, chimique atmospherique, geophysics, altimetry, radar, atmospheric chemistry |
|||||||||||||||||||||||||||||||||||||||||||||||||||||
| Copyright 2000 - European Space Agency. All rights reserved. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||