Generation and evolution of internal solitary waves in the southern Taiwan Strait

ABSTRACT

The generation processes and potential energy sources of internal solitary waves (ISWs) in the southern Taiwan Strait are investigated by driving a high resolution non-hydrostatic numerical model with realistic background conditions. Two main types of ISWs are clarified according to their different energy sources. One is generated by the nonlinear disintegration of remote internal tides emanating from Luzon Strait, and the other type is generated by local tide-topography interaction at the continental slope. The basic properties and evolution processes differ between these two kinds of ISWs. The waves originated from the remote internal tides at Luzon Strait have amplitudes comparable to previous field observations. In contrast, the ISWs generated locally are much weaker than observed waves, even in the presence of a steady offshore background current, which intensifies the generation of onshore ISWs. The ISWs induced by remotely generated M₂ internal tides are stronger than those induced by K₁ internal tides, and the fraction of internal wave energy transmitted onto the shelf is not significantly influenced by the intensity of remotely generated internal tides.     

A numerical study of the barotropic tides and tidal energy distribution in the Indonesian seas with the assimilated finite volume coastal ocean model

The tides and tidal energetics in the Indonesian seas are simulated using a three-dimensional finite volume coastal ocean model. The high-resolution coastline-fitted model is configured to better resolve the hydrodynamic processes around the numerous barrier islands. A large model domain is adopted to minimize the uncertainty adjacent to open boundaries. The model results with elevation assimilation based on a simple nudge scheme faithfully reproduced the general features of the barotropic tides in the Indonesian Seas. The mean root-mean-square errors between the observed and simulated tidal constants are 2.3, 1.1, 2.4, and 1.5 cm for M₂, S₂, K₁, and O₁, respectively. Analysis of the model solutions indicates that the semidiurnal tides in the Indonesian Seas are primarily dominated by the Indian Ocean, whereas the diurnal tides in this region are mainly influenced by the Pacific Ocean, which is consistent with previous studies. Examinations of tidal energy transport reveal that the tidal energy for both of the simulated tidal constituents are transported from the Indian Ocean into the IS mainly through the Lombok Strait and the Timor Sea, whereas only M₂ energy enters the Banda Sea and continues northward. The tidal energy dissipates the most in the passages on both sides of Timor Island, with the maximum M₂ and K₁ tidal energy transport reaching about 750 and 650 kW m^–1, respectively. The total energy losses of the four dominant constituents in the IS are nearly 338 GW, with the M₂ constituent dissipating 240.8 GW. It is also shown that the bottom dissipation rate for the M₂ tide is about 1–2 order of magnitudes larger than that of the other three tidal components in the Indonesian seas.     

Ocean tides and resonance

Tidal currents and surface amplitudes are calculated globally for the dominating diurnal and semidiurnal constituents using an established tidal model under a range of altered bathymetry. The purpose is to evaluate if the well-known amplification of the global tides during the Last Glacial Maximum (LGM) is related to changed propagation properties for the tidal wave or to changed damping due to removal of shelf seas. The response of the tides and tidal dissipation to future sea-level rise is also discussed. The tides in the present and LGM oceans were simulated first, followed by runs where the present day bathymetry was used but the shelf seas removed by the introduction of vertical walls or where sea level is allowed to rise. Previously reported results regarding tidal amplitudes and dissipation rates are reproduced in the control runs. The runs without shelf seas show significantly enhanced tidal amplitudes in the North Atlantic, whereas sea-level rise of 5 m above present levels show a significant shift in the amphidromic points on a local and regional scale but had a limited effect on the open ocean tides. Simulations with very large sea-level rise show a significantly decreased global tidal dissipation, whereas experiments without friction in present-day shallow water display results similar to those with no shelf seas. The results all point towards changing damping properties due to the removal of shelf seas as being the mechanism behind the LGM amplification, and they imply the importance of implementing future sea-level changes properly in tidal simulations.     

The predictability of near-coastal currents using a baroclinic unstructured grid model

A limited domain, coastal ocean forecast system consisting of an unstructured grid model, a meteorological model, a regional ocean model, and a global tidal database is designed to be globally relocatable. For such a system to be viable, the predictability of coastal currents must be well understood with error sources clearly identified. To this end, the coastal forecast system is applied at the mouth of Chesapeake Bay in response to a Navy exercise. Two-day forecasts are produced for a 10-day period from 4 to 14 June 2010 and compared to real-time observations. Interplay between the temporal frequency of the regional model boundary forcing and the application of external tides to the coastal model impacts the tidal characteristics of the coastal current, even contributing a small phase error. Frequencies of at least 3 h are needed to resolve the tidal signal within the regional model; otherwise, externally applied tides from a database are needed to capture the tidal variability. Spatial resolution of the regional model (3 vs 1 km) does not impact skill of the current prediction. Tidal response of the system indicates excellent representation of the dominant M ₂ tide for water level and currents. Diurnal tides, especially K ₁, are amplified unrealistically with the application of coarse 27-km winds. Higher-resolution winds reduce current forecast error with the exception of wind originating from the SSW, SSE, and E. These winds run shore parallel and are subject to strong interaction with the shoreline that is poorly represented even by the 3-km wind fields. The vertical distribution of currents is also well predicted by the coastal model. Spatial and temporal resolution of the wind forcing including areas close to the shoreline is the most critical component for accurate current forecasts. Additionally, it is demonstrated that wind resolution plays a large role in establishing realistic thermal and density structures in upwelling prone regions.     

Excess of bottom-released methane in an Arctic shelf sea polynya in winter

Latent heat polynyas are regions generating strong ice formation, convection and extensive water mass formation. Here we report on the effects of these processes on resuspension of sediments and subsequent methane release from the seafloor and on the resulting excess methane concentration in surface water on a polar shelf during winter. The study is based on measurements of concentration and δ¹³C values of methane, water temperature, salinity, light transmission and sea ice data collected in March 2003 in Storfjorden, southern Svalbard. In winter, strong and persistent northeasterly winds create polynyas in eastern Storfjorden and cause ice formation. The resulting brine-enriched water cascades from the Storfjordbanken into the central depression thereby enhancing the turbulence near the seafloor. A distinct benthic nepheloid layer was observed reflecting the resuspension of sediments by the cascading dense bottom water. High concentrations of ¹³C-depleted methane suggest submarine discharge of methane with the resuspended sediments. As the source of the submarine methane, we propose recent bacterial methanogenesis near the sediment surface because of extremely high accumulation rates of organic carbon in Storfjorden. Convective mixing transports newly released methane from the bottom to the sea surface. This eventually results in an excess concentration in surface water with respect to the atmospheric equilibrium, and a sea-air flux of methane during periods of open water. When a new ice cover is formed, methane becomes trapped in the water column and subsequently oxidized. Thus, the residual methane is strongly enriched in ¹³C in relation to the δ¹³_CCH4${\delta }^{13}{\mathrm{C}}_{{\mathrm{CH}}_4}$ signature of atmospheric methane. Our results show that latent heat polynyas may induce a direct pathway for biogases like methane from sediments to the atmosphere through coupling of biogeochemical and oceanographic processes. Extrapolating these processes to all Arctic ocean polynyas, we estimate a transfer of CH₄ between 0.005 and 0.02 Tg yr⁻¹. This is not a large contribution but the fluxes from the polynyas are 20–200 times larger than the ocean average and the methane evasion process in polynyas is certainly one that can be altered under climate change.     

Inversion for tides in the Eastern North Pacific Ocean

A regional model of tides in the Eastern North Pacific Ocean is developed through the use of inversion with two-dimensional finite element codes. Since global tide models are least accurate in coastal environments, modeling tides on a regional scale allows tidal propagation and interaction along the coast to be more accurately represented. In this respect, a regional model can act as a liaison between open ocean dynamics and physical processes more pertinent to coastal systems. The region of interest in this study extends from the Aleutian Islands to Southern California and includes deep ocean, continental shelf, and shallow water features. Boundary conditions are determined from nonlinear inversion of harmonic data from both shallow water and deep ocean tide gauges. Spatial patterns of amplitudes and phases from the model are examined for major constituents. Results are also compared to global tide models at selected stations.     

Internal tides and waves near the continental shelf edge

Abstract

The subject is reviewed from the viewpoints of theory, internal tide and wave structure and their implications.

A wider theoretical context suggests scope for further investigation of natural or nearly-trapped forms above the inertial frequency.

Although internal tides in many locations are observed to have first-mode vertical structure, higher modes are seen offshore from shallow shelf-break forcing and for particular Froude numbers, and may be expected locally near generation. Bottom intensification is often observed where the sea floor matches the characteristic slope. Solitons form from internal tides of large amplitude or at large changes of depth.

Internal tides and solitons are observed also at many sills and in straits, and to intensify in canyons.

Non-linear effects of the waves, especially solitons, include the conveyance of water, nutrients, ‘‘mixing potential'’ etc. away from their source to other locations, and the generation of mean currents. The waves transfer energy and possibly heat between the ocean and shelf, may be a source of medium frequency waves on the shelf (periods of minutes) and can contribute to interior mixing and overturning, bottom stirring and sediment movement.     

Quarter diurnal shelf resonances and tidal bed stress in the English Channel

The mean sea level and mean bed stress due to tidal co-oscillations in the presence of quadratic friction is examined analytically and numerically. In some idealised situations under conditions of near M₄ shelf resonances it is shown that phase relationships tend to exist between the M₂ tidal currents and the M₄ tidal currents which combine to give maximum currents in the flood or ebb direction. In the presence of quadratic friction these effects result in a mean bed stress and affect mean sea level. It is shown that these idealised responses are in part relevant to the sea level differences and sand transport paths due to tides around the British Isles.     

Tidal friction for semidiurnal tides

The quadratic law of bottom friction demands an increased frictional coefficient for the S₂ and N₂ tides with respect to a dominant M₂ tidal signal. A numerical model of the semidiurnal tides in the northeast Atlantic gives an increase in friction of ∼35% for the N₂, S₂ and K₂ tides with respect to the M₂ tide and this value is close to a theoretically derived estimate for the region, providing confirmation of the general widescale applicability of the quadratic friction law. Small differences in friction also occur for the N₂, S₂ and K₂ tides and this is attributed to the interaction of the L₂ and μ₂ tides with the M₂ tide in the presence of quadratic friction. Energy dissipation relationships, anomalous K₂/S₂ amplitude ratios and the role of quadratic friction on 18.6 year variations of semidiurnal forcing are examined.     

Sediment mass and distribution in the South China Sea since the Oligocene

The sediment macro-distribution patterns and their evolutionary characteristics in the South China Sea (SCS) are discussed based on a quantification of the sediment mass from the be- ginning of seafloor spreading in the Oligocene to the Present. Above the pre-Oligocene base, the total sediment mass for the whole SCS is estimated to be 1.44×1016 t, with the highest average accumulation rate of ~22 g·cm-2·ka-1 in the Oligocene. Having no large abyssal fans but fast accumulation in sedimentary basins on the continental shelf and slope, the SCS shows quite different sedimentary characters not only from the open ocean but also from small backarc basins along the marginal West Pacific, apparently controlled by the coupling between local tectonics and global climate changes.     

On sediment accumulation rates and stratigraphic completeness: Lessons from Holocene ocean margins

Stratigraphic completeness is a fundamental consideration when deciphering the mass accumulation history of sediments and the geologic record of earth and ocean processes. In this study, stratigraphic completeness was examined in the context of late Holocene sedimentary successions using published sediment accumulation rates for five ocean margin systems (Amazon shelf, Hudson estuary, northern California shelf, Mid-Atlantic slope, Santa Monica Bay). Plots of mass accumulation rate versus time span of averaging were used to determine how rates scale with measurement period, and to estimate levels of stratigraphic completeness for comparison within and among margin systems. Statistically significant inverse correlations between accumulation rate and time span of averaging are indicated for all but one of these systems—most of the sedimentary records examined are stratigraphically incomplete. At the 10³-yr level of resolution, completeness is 20–48% for strongly tidal estuarine (Hudson estuary) and deltaic shelf (Amazon shelf) sites, 51–91% for accretionary shelves (northern California shelf) and slopes (Mid-Atlantic slope), and 85–100% for a sediment-starved slope (Santa Monica bay). Mass accumulation rates converge to a relatively narrow range (0.01–0.1 g cm⁻² yr⁻¹) at the 10⁴-yr level of resolution, consistent with the notion that there are universal controls on sediment accumulation rate, i.e., rate of sea-level rise and sediment supply. Among sites on the Amazon and northern California shelves, within-system completeness varies by ∼10–20% on account of site-specific sedimentary processes that preferentially trap or disperse suspended sediment. Overall, stratigraphic completeness increases with water depth shelf-to-slope, yet depth is not a robust predictor of completeness in general owing to differences in strata-forming processes among shallow-marine environments. Significantly, completeness varies inversely with instantaneous deposition rate as the most sediment-rich systems tend to exhibit the most incomplete sedimentary records. The findings of this study emphasize the importance of considering time span and fidelity when interpreting the accumulation history of modern and Holocene sedimentary strata.     

A study of tides and currents in Cook Strait, New Zealand

Greater Cook Strait (GCS) lies between the North and the South Islands of New Zealand. Its location at the convergence of the Pacific and Indo-Australian tectonic plates leads to interesting bathymetry with an adjacent shallow shelf and deep ocean trench as well as numerous crossing faults and complex shoreline geometry. Our purpose in this study is to examine tides and currents in GCS and, in particular, identify the major forcing mechanisms for the residual currents. Toward this end, we use an unstructured-grid numerical model to reproduce the tides and currents, verify these results with observations and then use the model to separate the various forcing mechanisms. The physical forcing includes nonlinear generation from tides and tidal currents, differences in sea level between the Pacific Ocean and Tasman Sea boundaries, density variations, wind stress and river discharge into GCS. Each of these mechanisms is important in different areas.     

Sunlight and tidal interaction as a mechanism of carbon transport in an ice-covered region: Results of a coupled ice–ocean ecosystem model

In order to clarify the mechanism of carbon transport in an ice-covered ecosystem in Lake Saroma (44^°N${44}^{°}\mathrm{N}$, 143^°E${143}^{°}\mathrm{E}$, Hokkaido, Japan), a three-dimensional numerical calculation using a coupled ice–ocean ecosystem model was conducted. This model comprises an ocean ecosystem model, an ice ecosystem model, and equations for the coupling between ice and ocean. Comparisons of calculated results with observational data confirm that the calculation well reproduced the in situ phenomena with respect to tides, tidal currents, concentrations of POC and chlorophyll a in ice and in water, and sinking fluxes beneath the ice. The analysis of the organic carbon budget based on the calculation reveals that tide-induced transport, the enhancement of biological production in a pelagic system, and the physical release of organic matter from ice associated with ice-melting are important factors affecting the carbon transport during the ice-melting season. The carbon transport has a one-day time cycle. This is because principal driving forces are sunlight, and diurnal tides. The described mechanism of “sunlight and tidal pumping” is one of the most important features of carbon transport in a coupled ice–water ecosystem.     

Distribution,characteristics and potential impacts of chromophoric dissolved organic matter (CDOM) in Hudson Strait and Hudson Bay,Canada

The characteristics of chromophoric dissolved organic matter (CDOM) were studied in Hudson Bay and Hudson Strait in the Canadian Arctic. Hudson Bay receives a disproportionately large influx of river runoff. With high dissolved organic matter (DOM) concentrations in Arctic rivers the influence of CDOM on coastal and ocean systems can be significant, yet the distribution, characteristics and potential consequences of CDOM in these waters remain unknown. We collected 470 discrete water samples in offshore, coastal, estuarine and river waters in the region during September and October 2005. Mixing of CDOM appeared conservative with salinity, although regional differences exist due to variable DOM composition in the rivers discharging to the Bay and the presence of sea-ice melt, which has low CDOM concentrations and low salinity. There were higher concentrations of CDOM in Hudson Bay, especially in coastal waters with salinities <28$<28$, due to river runoff. Using CDOM composition of water masses as a tracer for the freshwater components revealed that river runoff is largely constrained to nearshore waters in Hudson Bay, while sea-ice melt is distributed more evenly in the Bay. Strong inshore–offshore gradients in the bio-optical properties of the surface waters in the Hudson Bay cause large variation in penetration of ultraviolet radiation and the photic depth within the bay, potentially controlling the vertical distribution of biomass and occurrence of deep chlorophyll maxima which are prevalent only in the more transparent offshore waters of the bay. The CDOM distribution and associated photoprocesses may influence the thermodynamics and stratification of the coastal waters, through trapping of radiant heating within the top few meters of the water column. Photoproduction of biologically labile substrates from CDOM could potentially stimulate the growth of biomass in Hudson Bay coastal waters. Further studies are needed to investigate the importance of terrestrial DOM in the Hudson Bay region, and the impact of hydroelectric development and climate change on these processes.     

Multi-scale modeling of Puget Sound using an unstructured-grid coastal ocean model: from tide flats to estuaries and coastal waters

Water circulation in Puget Sound, a large complex estuary system in the Pacific Northwest coastal ocean of the United States, is governed by multiple spatially and temporally varying forcings from tides, atmosphere (wind, heating/cooling, precipitation/evaporation, pressure), and river inflows. In addition, the hydrodynamic response is affected strongly by geomorphic features, such as fjord-like bathymetry and complex shoreline features, resulting in many distinguishing characteristics in its main and sub-basins. To better understand the details of circulation features in Puget Sound and to assist with proposed nearshore restoration actions for improving water quality and the ecological health of Puget Sound, a high-resolution (around 50 m in estuaries and tide flats) hydrodynamic model for the entire Puget Sound was needed. Here, a three-dimensional circulation model of Puget Sound using an unstructured-grid finite volume coastal ocean model is presented. The model was constructed with sufficient resolution in the nearshore region to address the complex coastline, multi-tidal channels, and tide flats. Model open boundaries were extended to the entrance of the Strait of Juan de Fuca and the northern end of the Strait of Georgia to account for the influences of ocean water intrusion from the Strait of Juan de Fuca and the Fraser River plume from the Strait of Georgia, respectively. Comparisons of model results, observed data, and associated error statistics for tidal elevation, velocity, temperature, and salinity indicate that the model is capable of simulating the general circulation patterns on the scale of a large estuarine system as well as detailed hydrodynamics in the nearshore tide flats. Tidal characteristics, temperature/salinity stratification, mean circulation, and river plumes in estuaries with tide flats are discussed.     

The oceanic tides in the South Atlantic Ocean

The finite element ocean tide model of Le Provost and Vincent (1986) has been applied to the simulation of the M₂ and K₁ components over the South Atlantic Ocean. The discretisation of the domain, of the order of 200 km over the deep ocean, is refined down to 15 km along the coasts, such refinement enables wave propagation and damping over the continental shelves to be correctly solved. The marine boundary conditions, from Dakar to Natal, through the Drake passage and from South Africa to Antarctica, are deduced from in situ data and from Schwiderski’s solution and then optimised following a procedure previously developed by the authors. The solutions presented are in very good agreement with in situ data: the root mean square deviations from a standard subset of 13 pelagic stations are 1.4 cm for M₂ and 0.45 cm for K₁, which is significantly better overall than solutions published to date in the literature. Zooms of the M₂ solution are presented for the Falkland Archipelago, the Weddell Sea and the Patagonian Shelf. The first zoom allows detailing of the tidal structure around the Falklands and its interpretation in terms of a stationary trapped Kelvin wave system. The second zoom, over the Weddell Sea, reveals for the first time what must be the tidal signal under the permanent ice shelf and gives a solution over that sea which is generally in agreement with observations. The third zoom is over the complex Patagonian Shelf. This zoom illustrates the ability of the model to simulate the tides, even over this area, with a surprising level of realism, following purely hydrodynamic modelling procedures, within a global ocean tide model. Maps of maximum associated tidal currents are also given, as a first illustration of a by-product of these simulations.     

Form drag is a major component of bed shear stress associated with tidal flow in the vicinity of an isolated sand bank,Torres Strait,northern Australia

Tidal current and elevation data were collected from five oceanographic moorings during October 2004 in Torres Strait, northern Australia, to assess the effects of large bedforms (i.e., sand banks) on the drag coefficient (C_D) used for estimating bed shear stress in complex shallow shelf environments. Ten minute averages of tidal current speed and elevation data were collected for 18 days at an on-bank site (<7 m water depth) and an off-bank site (<10 m). These data were compared to data collected simultaneously from two shelf locations (<11 m) occupied to measure regional tidal behaviour. Overall C_D estimates at the on- and off-bank sites attained 7.0±0.1×10⁻³ and 6.6±0.1×10⁻³, respectively. On-bank C_D estimates also differed between the predominant east–west tidal streams, with easterly directed flows experiencing C_D=7.8±0.18×10⁻³ and westerly directed flows C_D=6.4±0.12×10⁻³. Statistically significant differences between the off-bank and on-bank sites are attributed to the large form drag exerted by the sand banks on the regional tidal currents, and statistically significant differences between the westward and eastward flows is ascribed to bedform asymmetry. Form drag from the large bedforms in Torres Strait comprises up to 65% of the total drag coefficient. When constructing sediment transport models, different C_D estimates must therefore be applied to shelf regions containing steep bedforms compared to regions that do not. Our results extend the limited inventory of seabed drag coefficients for shallow shelf environments, and can be used to improve existing regional seabed mobilisation models, which have direct application to environmental management in Torres Strait.     

Tide model comparison over the Southwestern Atlantic Shelf

Sea surface height (SSH) as measured by satellites has become a powerful tool for oceanographic and climate related studies. Whereas in the open ocean good accuracy has been achieved, more energetic dynamics and a number of calibration problems have limited applications over continental shelves and near the coast. Tidal ranges in the Southwestern Atlantic (SWA) continental shelf are among the highest in the world ocean, reaching up to 12 m at specific locations. This fact highlights the relevance of the accuracy of the tidal correction that must be applied to the satellite data to be useful in the region. In this work, amplitudes and phases of tidal constituents are extracted from five global tide models and three regional models and compared to the corresponding harmonics estimated from coastal tide gauges (TGs) and satellite altimetry data. The Root Sum Square (RSS) of the misfit of the common set of the five tidal constituents solved by the models (M₂, N₂, S₂, K₁ and O₁) is higher than 18 cm close to the coast for two of the regional models and higher than 24.5 cm for the rest of the models considered. Both values are too high to provide an accurate estimation of geostrophic non-tidal currents from satellite altimetry in the coastal region. On the other hand, the global model with the highest spatial resolution has a RSS lower than 4.5 cm over the continental shelf even when the non-linear M₄ overtide is considered. Comparison with in-situ current measurements suggests that this model can be used to de-tide altimetry data to compute large-scale patterns of SSH and associated geostrophic velocities. It is suggested that a local tide model with very high resolution that assimilates in-situ and satellite data should meet the precision needed to estimate geostrophic velocities at a higher resolution both close to the coast and over the Patagonian shelf.