EDDY FORMATION AND INTERANNUAL SEA LEVEL VARIATION

IN THE EAST SEA (SEA OF JAPAN)

USING TOPEX/POSEIDON ALTIMETER DATA
 
 

Byoung-Ju Choi and Fadli Syamsudin

Institute of Marine and Coastal Science, Rutgers University
1 May 2000


[CONTENTS]

     ABSTRACT
    1. INTRODUCTION         2. METHODOLOGY
    3. RESULTS      3.1 Comparison of Altimeter and Tide Gauge data
                              3.2 Eddy Formation and Migration
                              3.3 Interannual variation of sea surface height from 1993 to 1999
    4. DISCUSSION                5. CONCLUSION              6. REFERENCE
    APPENDIX


ABSTRACT

   Eddies formation in 1999 and interannual variation of sea surface height from 1993 to 1999 in the East Sea were  examined using TOPEX/ERS2 and TOPEX/POSEIDON  Altimeter data, respectively. The Eddies were formed in September and continuously developed to become stronger anticyclone Eddies in October, November and December  in 1999 . This finding is a new phenomenon compare to the previous result that Eddies began to form in February and diminished  in May. Interannual Variability of Sea Surface Height in the East Sea showed that Sea Surface Height was higher in 1994, 1998 and 1999, and it was lower in 1993, 1995 and 1996. When the warm water region expands to the North,  the heat content of upper layer increases and sea surface height  rises.  In the contrarily, when cold water region expands to the South,  heat content of upper layer decreases and sea surface height  falls.
 


   1. INTRODUCTION

       Sea level and its horizontal slope and time changes are the surface expression of ocean  processes occurring over large spatial and temporal scales.  Large-scale currents (such as the Kuroshio, extending up to kilometers in depth) and planetary Rossby waves (thousands of km  in horizontal extent) carry the memories of past air-sea exchanges that can affect subsequent weather at great distances from the source of the original air-sea exchange.  Sea level also provides evidence of local heating by solar radiation, the perennial tides, and the addition of  water to the ocean by melting ice caps or glaciers (Ocean ESIP ).

     Altimeters emit a sequence of short pulses at microwave frequencies and then measure the return times to ascertain the instrument-surface distance. The TOPEX/Poseidon altimetric data has been widely used for several applications, such as ocean circulation,  El Nino/La Nina studies, sea level monitoring for environmental safety warning, simulation of larvae dynamics, and so on.  This has been diverging from its first specific mission to study the world's ocean. There are  a number of improvements that has been made on the satellite performances and algorithms  in order to study ocean circulation in the scale of basins and gyres( Fu et al, 1994 ). Nowadays, we can use this altimetric data for smaller scales of physical processes in the ocean.

       East Sea (Sea of Japan) is one of the largest marginal seas. Its mean depth is about 2,700 m. The surface water inflows through a narrow strait with sill depth about 100m and outflows through two narrow straits with sill depths less than 50 m( Tomczak and  Godfrey,1994 ).
 


(a)                                                               (b)
(a) Horizonat

   [Figure 1]  (a)   Horizontal temperature distribution at 100 m depth. (b)   Vertical  structure of temperature along section 133.8'E, 35.4 - 43'N.  From  Altimetry Plots: East/Japan Sea

   We can see a strong horizontal front is formed in the middle of East Sea(Figure 1). It undulates to the North from 38 degree N at the western part and goes down to 37.5 N, 132E and then moves up to 40N again at the eastern part of the East sea (Figure 1 (a) ). The front is separating the Northern Cold Water Region and the Southern Warm Water Region. A Northward western boundary current ( East Korea Warm Current) flows along the coast of Korea and then meanders along the front. The temperature section shows that the the temperature becomes 1 degree C about  300 m depth and almost homogeneous below the depth ( Figure 1 (b) ).

    Our hypothesis are that Eddy formation and migration can be detect by Sea Surface Height Deviation and Mean Sea Surface Height can be computed from climatological hydrographic data. Objectives are to examine Eddy Formation and migration in 1999 and estimate Interannual Variation of Sea Surface Height from 1993 to 1999  in the East Sea..
 

  2. METHODOLOGY

   2.1 TOPEX/POSEIDON Satellite Sensor

  All data taken by satellite are being processed by six science instruments below:

     1. Operational sensors: Dual-frequency radar altimeter (ALT)/NASA for  minimizing the errors caused by the ionospheric free electron ; TOPEX microwave radiometer (TMR)/NASA for estimating total water vapor content in the atmosphere; Laser retroreflector array (LRA)/NASA and Doppler orbitography and radiopositioning integrated by satellite (DORIS); Dual Doppler tracking system receiver/CNES for precision orbit determination.

     2. Experimental sensors: Single-frequency solid state radar altimeter (SSALT)/CNES and Global Positioning System (GPS) demonstration receiver (GPSDR)/NASA.

A number of improvements of those sensors performance make it possible to study Ocean Circulation in the basin scale such as Gyres and Eddies.

    2.2 Relation between sea surface height, heat content of upper layer, and upper layer thickness

    The East Sea has well defined 2 layers in vertical temperature structure (Figure 1-b).
The sea  surface height is a good measure of heat content of upper ocean and upper layer thickness. Since 8 degree C isotherm is located at the top of thermocline, the depth of 8 degree C can be a upper layer thickness. Warm surface layer has about 150 m layer thickness in Southern part of the East Sea. Lower layer is colder than 1 degree C.

Since heat content is given as Heat content = c x density x Temperature,
heat content of upper 300 m layer is the integration of heat content over depth.

Heat content of upper layer = c x average density x average Temperature x 300 m.

 where c is a specific heat,  4,200 J / ( Kg degree C ).

Average temperature over 300 m can represents the heat content of upper 300 m layer.
 
 


(a)                                                (b)

 [Figure 2] a. Comparison  between annual mean adjusted sea level and heat content of upper 300m layer at Ulleung Island (South Western East Sea). b. Comparison between the depth of 8 degree C isotherm and mean adjusted sea level. The depth of 8 degree C isotherm represents a thermocline depth or upper layer thickness. From Choi (1994).
 

   Heat content of upper layer and sea level variation are highly correlated (Figure 2-a). The correlation coefficient is 0.84. When the average temperature increases, the sea level rises. The thickness of upper layer and sea level are linearly correlated (Figure 2-b). When the thermocline goes down, upper layer thickness increases, heat content of upper 300 layer increase, sea surface height increases. When the sea surface height is high, the heat content of upper ocean is high and the upper layer thickness is thick.

   2.3 Geostrophic Calculation

             We used UNESCO Training Module for application of TOPEX/POSEIDON satellite Altimetry to produce maps of eddies and other geostrophic features developed by Donald R. Kobayashi, et al (1999). The method can compute sea level height from measured TOPEX altimeter data. We create a simple fortran program to calculate its East-West and North-South velocity components from where we made our current vector plot using Generic Mapping Tool (GMT).  We are also comparing our result with observational data and sea level height anomaly and its geostrophic current taken from interactive map in the internet produced by Colorado Center for  Astrodynamics Research (CCAR), University of Colorado.

                The dz/dx and dz/dy directional gradients from the resulting sea surface topography are used in the following equations to estimate the u (East-West) and v (North-South) components of the geostrophic current in units of centimeters per second:

                 u = -(g/f)dz/dy       and        v = (g/f)dz/dx

                     where

                      g = 980 cm per second (gravitational acceleration)

                      f = 2 W sin ( phi )

                      W = 7.29x10-5 radians per second (earth angular rotational velocity)

                      phi = latitude
 

3. RESULTS

   3.1 Comparison of Altimeter and Tide Gauge data
 

Ulleung Island is located at the South Eastern part of the East Sea(131'E, 37.5'N).
Because it is a small island and the water depth around it is deeper than 1,000 m,
it is a good place to measure sea level and compare its sea level data from tide gauge with the air-borne altimeter data. We compared the annual mean deviation from 1993 to 1998. They have a good linear relation and their ranges are from -4 cm to +8 cm (Figure 3). Ulleung island has higher sea level in 1995 and 1998.
 


  [Figure 3]. Comparison of Altimeter data and Tide Gauge data
 
 
 

   3.2 Eddy Formation and Migration

        Sea surface height deviation data are downloaded from  Colorado University . After we download data, we examine the sea surface height deviation at the middle of each month.

   [Figure 4] Formation and Migration of eddies at each month in 1999.

January 15, 1999:    We have two  strong warm anticyclonic eddies at 137E, 39N and 134E, 37-38N and a very weak  anticyclonic eddy at 130E, 37-37.5N in the southern part of East Sea.

 February 15, 1999:  The two strong Eddies found in January is decaying but The  weaker Eddy is increasing its intense.

March 15, 1999: We found a new weak cyclonic Eddy at 130E, 39N.

April 15, 1999: The weak cyclonic Eddy found in March is decay.

May 15, 1999: The Anticyclonic Eddies are beginning to form again at 131E, 39N; 137-138E, 39N; and 134E, 37-38N.

June 15, 1999: Depression of the sea level in the eastern part of East sea is increasing and make the strong southward geostrophic current moving farther into western part of the region. The Anticyclonic Eddy of 131E, 39N and 134E, 37-38N are decay, but  Eddy of 137-138E, 39N becomes stronger.
 

July 15, 1999: It has relatively similar pattern with the previous month but here it is  giving an increase in the movement of warm water to the west. All Eddies are decay.

August 15, 1999:  It is beginning to form a large area of Cyclonic Eddy at 132-134E, 39N.

September 15, 1999: The cyclonic Eddy in August is becoming stronger and shifting to 135E, 39N. We found also strong anticyclonic Eddy at 131E, 38N.

October 15, 1999:  Strong anticyclonic Eddies are developed at 131-132E, 39-40N; 133-134E, 37-38N; 136-137E, 38-39N.
 

November 15, 1999: The strong geostrophic pattern of Eddies are similar than previous month but a little bit shifting to the south .

December 15, 1999:  We found more energetic eddies in the southern part of East sea.
 

    3.3 Interannual variation of sea surface height from 1993 to 1999
 

       Sea surface height deviation data are downloaded from  Ocean ESIP . After we got the data, we subset the data for the East Sea using FORTRAN program and plot the data using GMT plotting software. Sea  surface height deviation are computed at 1 x 1 degree intervals with respect to the mean of 4-year period, 1993-1996.
 
 

  [Figure 5] Sea level annual mean deviation at each year from 1993 to 1999.

1993                                              1994
   1993: negative deviation in western and southern part of the East Sea
   1994: strong positive deviation in western and southern part of the East Sea

1995                                              1996

   1995: negative deviation in eastern part of the East Sea
   1996: strong negative deviation in southern part of the East Sea

 1997                                                        1998

   1997: positive deviation in western and southern part of the East Sea
   1998: strong positive deviation in southern part of the East Sea

       1999
           
   1999: strongest positive deviation in southern part of the East Sea. It exceeds our color index.
 

    Because the strong horizontal front is formed along  40'N and it confines warm water region  below 40'N, the sea surface height, heat content, and upper layer thickness highly change below 40'N in East Sea.  To summarize our Figure 5 we computed the average of sea surface height deviations in Southern East Sea (35' - 40'N) in each year.
 
 

[Table 1]  The average of sea surface height deviations
in Southern East Sea (unit: cm)
  Year 1993  1994   1995   1996   1997   1998   1999 
 AVG.  -0.5    2.0 -0.4   -1.4   0.8     1.5    3.5

      We can find that the sea surface height ( or heat content of upper layer ) was higher in 1994, 1998 and 1999, and it was lower in 1993, 1995 and 1996.
 

   4. Discussion

      4.1 Eddy Formation

From historical Eddy formation in the East Sea, we found that The Eddy began to form at February where strong colder northerly winter wind system blew in the region and diminished in May as the onset of weaker summer wind system. In contrast, we found different condition of Eddies formation in 1999. The anticyclonic Eddies already formed in January and became to decay in February, March, April and diminished in May.

We found a new result that has never been observed before of the early anticyclonic Eddy formation  in September and continued until October, November and December in 1999, and even in January 2000. We don't have enough wind data and supporting analysis to explain this event yet and we need to analyze our data after we add Mean Sea Surface Height to the sea surface height deviation.
 

     4. 2 Interannual Sea Surface Height Variation

     What is the  oceanographic process that changes sea surface height ( or heat content of upper layer ) in Southern East Sea? We examined the deviation of sea surface height. we cannot understand the oceanographic process from the deviation. We need mean sea surface height.
 


Real Sea Surface Height = Mean Sea Surface Height  + Sea Surface Height Deviation


 


 We can make the mean sea surface height from long term measurement of altimetry data, numerical model, and observation data. US Navy generates Mean Sea Surface Height from long term altimeter data, climatological data and Numerical model.

(a)   Sea Surface Height Deviation    (b)   Sea Surface Height Residual

      [Figure 6]  (a) Sea Surface Height Deviation (b) Sea Surface Height Residual on April 21, 2000. From  Altimetry Plots: East/Japan Sea
 

  A mean steric heights of surface relative to 500 m were calculated using long term hydrographic observation data from KODC and JODC (Figure 7-a ) and we used the mean steric height as a mean sea surface height (Figure 7-b).
 
 

  [Figure 7] (a) Climatological mean temperature at 100 m depth (b) Estimated Mean Surface Height from climatological density data (c) Actual sea surface height distribution from 1993  to 1999.

(a)                                                                      (b)

We can see a strong horizontal front(5 degree isotherm, green color) is formed in the middle of East Sea(Figure 7- a, b). It undulates to the North from 38 degree N at the western part and goes down to 37.5 N, 132E and then moves up to 40N again at the eastern part of the East sea. The front is separating the Northern Cold Water Region and the Southern Warm Water Region. A Northward western boundary current ( East Korea Warm Current) flows along the coast of Korea and then meanders along the front.
            (c) 1993
         
   1993: The front(green color) is fluctuating.

1994                                                             1995

   1994: The front(green color) move to the North and Warm Water Region expands.
   1995: The front(green color) move to the North and Warm Water Region expands only in Central part of the East Sea.

 1996                                                             1997

   1996: The front(green color) move to the South and Cold Water Region expands.
   1997: The front(green color) is formed in the same position as climatology but sea surface height in southern part of the East Sea increases.

1998                                                          1999

  1998: The front(green color) move to the North and Warm Water Region expands. And sea surface height in southern part of the East Sea increases.
  1999: The front(green color) moves to the North and Warm Water Region expands. And sea surface height in southern part of the East Sea increases.

     We can notice that the movement of the front in the middle of the East Sea changes the oceanographic conditions such as upper layer thickness, heat content, and sea surface height. When Warm Water Region expanded to the North in 1994, 1998 and 1999, upper layer thickness increased, heat content of upper layer increased, and sea surface height rose. When Cold Water Region expanded to the South in 1996, upper layer thickness decreased, heat content of upper layer decreased, and sea surface height fell.

  We need to examine what cause the interannual variation of sea surface height in Southern East Sea in the future.  There are several factors that can affect sea surface heights of the East Sea : First one is the change of inflow and outflow through the straits. If the amount of warm water inflow through Korea and Sushima straits increase, the sea surface height in the East Sea will increase. Second one is wind. If a strong winter monsoon blow from North West, the warm water region becomes smaller and sea surface height will decrease.  Air-sea interaction is the third one. If the air cools the sea surface in cold winter, the water temperature become lower and sea surface height will decrease.
 

    5. CONCLUSION

       We compared the annual mean deviations of sea level data from tide gauge with the air-borne TOPEX/POSEIDON altimeter data from 1993 to 1998 at Ulleung Island.
They have a good linear relation and their ranges are from -4 cm to +8 cm. Ulleung Island has higher sea level in 1995 and 1998.

   Eddies formation in 1999 was examined using TOPEX/ERS2 Altimeter data. The Eddies were formed in September and continuously developed to stronger anticyclonic Eddies in October, November and December. This is a new result compare to the previous years when Eddies began to form in February and diminished  in May.

   Interannual variation of sea surface height from 1993 to 1999 in the East Sea  was studied by TOPEX/POSEIDON Altimeter data. Interannual Variability of Sea Surface Height in the East Sea showed that Sea Surface Height was higher in 1994, 1998 and 1999, and it was lower in 1993, 1995 and 1996. When the warm water region expands to the North,  heat content of upper layer increases , and sea surface height  rises.  When cold water region expands to the South,  heat content of upper layer decreases , and sea surface height  falls.
 


  Acknowledgment

        We would like to express our sincere thank to Hank Statscewich for helping us to use GMT plotting. Prof. Scott Glenn, Prof. Jim Miller, and Prof. Jennifer Francis for continuing support during the project. All friends that make this project finish.


    6.  REFERENCE

1. Byoung Ju Choi (1994): Relations between Sea Level Variations and Oceanographic Conditions around Ulleung Island, Master Thesis in Seoul National University, 54pp.

2. Delcroix, T., C., Gautier, (1987): Estimates of Heat Content Variations from Sea Level Measurements in the central and western Tropical Pacific from 1979 to 1985, Journal of Physical Oceanography, 17, 725-734.

3.  D.R. Kobayashi, J.J. Polovina (1999): UNESCO Training Module for application of TOPEX/POSEIDON satellite Altimetry to produce maps of eddies and other geostrophic features for fisheries research, National Marine Fisheries Service Honolulu Laboratory, NOAA. (http://www.nmfs.hawaii.edu/eco/unesco1.html).

4. Hahn, S.B., (1991): Tendency of Sea Water Temperature Rise and Mean Sea Level Rise in the Korean Water (unpublished).

5.  Kang, Y.Q. , B.D., Lee, (1985): On the Annual Variation of Mean Sea Level along the coast of Korea, Journal Oceanography Society Korea, 20 (1), 22-30.

6.  Kang, Y. , H., KAng, (1991): Long-Term Fluctuations of Water Temperature in the Upper 200m off the Southeast Coast of Korea, Bull. Korean Fisheries Society, 24 (6), 450-458.

7. L.L. Fu, E.J. Christensen, C.A. Yamarone Jr., et al (1994):   TOPEX/POSEIDON mission overview,
Journal of Geophysical  Research, 99, no. C12, 24369- 24381.

8. Oh, I.S., A.B., Rabinovich, M.S. Park, and R.N., Mansurov, (1993): Seasonal Sea Level
Oscillation in the East Sea, Journal Oceanography Society Korea, 28, 1-16.

9. Tomcazak, M. and J. Godfrey (1994) : Regional Oceanography : an introduction Pergamon, 422pp.


    APPENDIX

 Available TOPEX/POSEIDON DATA

   I.  NOAA Topex/Poseidon Analysis   The joint NASA/CNES satellite altimeter, Topex/Poseidon, has been operating since September 1992.  At NOAA these data are are used to form three different types of sea level time series:

        1. Topex/Poseidon Along-track Deviations : Global Topex/Poseidon sea level deviations (relative to the 3-year mean 1993-95) are computed at 1-degree latitude intervals along the satellite track. These data are available via anonymous ftp as either yearly files or one file for the entire T/P mission.

       2. Topex/Poseidon Gridded Deviations: Sea level deviations, computed along the Topex/Poseidon track as described above, are averaged in 4x1-degree (longitude x latitude) cells to construct regular grids of sea level deviation with respect to the 3-year period, 1993-1995. The data are averaged over monthly and seasonal periods. The gridded data are also available via anonymous ftp for monthly and seasonal periods.

      3. Topex/Poseidon Gridded Anomalies: The gridded deviations described above are converted to anomalies by removing the annual and semiannual harmonics. These maps show interannual changes relative to 1993-95 for monthly, seasonal, and annual time periods. The gridded data are also available via anonymous ftp for monthly, seasonal, and annual periods.

  II.  Ocean ESIP  is a member of NASA's Earth Science Information Partner (ESIP) Federation.

   1.  Uniform grid (1° x 1° x 5-day averages)   SeaSea Surface Height Residuals
   2.  Along-track grid (6.2 km spacing)    SeaSea Surface Height Residuals
 

    Sea Surface Height (SSH): Sea surface height (SSH) is the height of the sea surface above the reference ellipsoid. It is calculated by subtracting the corrected altimeter range (see above) from the altitude.

        SSH = Altitude - Corrected Range

   Residual Sea Surface: The residual sea surface will be defined here as the sea surface height minus the mean sea surface or the geoid and minus known effects, i.e., tides and inverse barometer.

    Residual Height = SSH - Geophysical Surface
                          - Tide Effects
                          - Inverse Barometer

    [ Geophysical Surface = Mean Sea Surface or Geoid ]
 

 III.  University of Texas, Center for Space Research (UTCSR)   TOPEX/POSEIDON Dynamic Ocean Topography and Sea Level Anomalies  (1992-Present)

  1. Gridded and Smoothed Data Files (1x1 deg resolution) sea level anomalies
  2. Sea Level Anomaly Images
 

 IV. Colorado University  Colorado University Near real time data Host

   1.   data host  Near Real Time SSH DATA
   2.  Gostrophic velocity viewer Goestrophic velocity
 

 Application of Altimeter data

 Aplication of TOPEX/POSEIDON Satellite Altimetry to  Simulate dynamics of larvae of the spiny lobster  in the NorthWestern Hawaiian Islands, 1995-96.     In Press, Fisheries Bulletin