Simulations of NLIW in the SCS using SUNTANS

1
Simulations of NLIW in the SCS using SUNTANS

• Z. Zhang, O. B. Fringer, R. L. Street
• Environmental Fluid Mechanics Laboratory
• Stanford University
• 28 August 2007

2
Generation, Propagation, Conversion, Dissipation
Nonhydrostatic SUNTANS simulation with 11 million
grid cells Horizontal resolution 500 m near
Dongsha and Luzon and 5 km near
boundaries Vertical resolution 10 m 150 m (100
vertical levels) Forcing first 8 tidal
constituents using OTIS at boundaries with 50 km
thick sponge layers Stratification Horizontally
uniform data from model of C. Wu.
3
(No Transcript)
4
(No Transcript)
5
(No Transcript)
6
(No Transcript)
7
(No Transcript)
8
(No Transcript)
9
(No Transcript)
10
(No Transcript)
11
(No Transcript)
12
(No Transcript)
13
(No Transcript)
14
(No Transcript)
15
(No Transcript)
16
(No Transcript)
17
(No Transcript)
18
(No Transcript)
19
(No Transcript)
20
(No Transcript)
21
(No Transcript)
22
(No Transcript)
23
(No Transcript)
24
(No Transcript)
25
(No Transcript)
26
(No Transcript)
27
(No Transcript)
28
(No Transcript)
29
(No Transcript)
30
(No Transcript)
31
Comparisons to observations
32
SPRING TIDE JUNE 2005
Observations (S. Ramp, NPS)
SUNTANS
33
SSH, 25 June 2005
SAR image from internalwaveatlas.com
34
Generation mechanism
Isotherm displacements at mooring B1 are due to
nonlinear waves generated at L1, 31 hours
earlier, which is the time it takes for the
first-mode internal tide to propagate from
station L1 to B1. Lagged peak displacements
correspond to generation during peak ebb (ebb is
defined as positive currents in the figure).
Data courtesy Steve Ramp, NPS
35
2D simulations with SUNTANS
Levitus stratification
28 C
Depth at sill DS 200 m
3 C
Ocean depth D0 3500 m
Barotropic forcing at diurnal frequency
Radiation of first-mode baroclinic wave. Sponge
layers are also employed to damp internal waves
at the boundaries.
36
Case with Frsill0.27
Isotherms 16, 20, 24, 28 degrees C
37
Case with Frsill1.60
Isotherms 16, 20, 24, 28 degrees C
38
Beginning of ebb tide
t10
End of ebb tide
t2T/2
LT/2c1
Beginning of flood Lee wave release
t3T/2
Arrival at B1
t4t3LB1/c1
DtLB1/c1
39
Lee wave vs. flood wave
15o Isotherms after 2 tidal periods
First flood
First ebb
Fr0.25
Increasing Fr
Stronger ebb pushes peak of depression farther
from the crest of the sill.
Fr3.00
After 2 tidal periods, the flood wave has
propagated roughly 2.5 mode-1 wavelengths,
while the ebb wave has propagated roughly 1
mode-1 wavelength. The crests of the flood
wave with increasing amplitude have propagated
further because of amplitude dispersion, while
the crests of the ebb wave have been delayed due
to the excursion of the peak of the depression
during the ebb.
40
Nonhydrostatic effects (Frsill1.60)
Nonhydrostatic code
Isotherms 16, 20, 24, 28 degrees C
Hydrostatic code
41
Nonhydrostatic-Hydrostatic Comparison
Dispersion in the hydrostatic model is purely
numerical. The numerical dispersion is much
smaller than the physical, nonhydrostatic
dispersion, which leads to excessive steepening
of the front. The oversteepened front is
diffused due to grid-scale numerical diffusion,
and this causes a reduction in the wave
amplitude. Reduction in the wave
amplitude reduces the amplitude dispersion and
thereby reduces the speed of propagation of the
wavetrain.
15o C isotherm after 3 tidal periods.
42
Using the KdV equation as a model
Initial wave Half-sine wave (first-mode internal
tidal wavelength) with amplitude 70 m.
43
Computed vs. Modeled (KdV) Results
SUNTANS
KdV
44
Effects of rotation
16 oC isotherm after 2 K1 tides
(1) Reduction of along-sill tidal excursion
distance in the ebb direction causing wave with
rotation to appear much faster. (2) Larger
amplitude wave without rotation eventually
overtakes wave with rotation due to
amplitude dispersion.
(2)
(1)
45
Conclusions

• Generation mechanism lee-wave release
• Nonlinear wavetrains in the SCS are likely
  generated as a result of the formation of large
  lee waves of depression in the Luzon Strait which
  are "released" at the end of the ebb tide and
  subsequently propagate into the SCS.
• Nonhydrostatic effects
• The hydrostatic model has insufficient dispersion
  (roughly ¼ that of the nonhydrostatic model), and
  this dispersion is purely numerical.
• The lack of sufficient dispersion leads to
  oversteepening and excessive numerical diffusion,
  which reduces the amplitude of the wavetrain, and
  this in turn causes a reduction in the
  propagation speed of the wavetrain.
• Effects of rotation
• Rotation leads to generation of smaller
  wavetrains which travel more slowly as a result
  of weaker amplitude dispersion.

46
(No Transcript)
47
(No Transcript)
48
(No Transcript)
49
(No Transcript)
50
28 June 2005, 1220
27 June 2005, 1130
28 June 2005
30 June 2005, 1401
29 June 2005, 1311
Generation is strongly influenced by
three- dimensional topography.