Turbulence-particle interactions under surface gravity waves

The dispersion and transport of single inertial particles through an oscillatory turbulent aquatic environment are examined numerically by a Lagrangian particle tracking model using a series of idealised test cases. The turbulent mixing is incorporated into the Lagrangian model by the means of a stochastic scheme in which the inhomogeneous turbulent quantities are governed by a one-dimensional k- ε turbulence closure scheme. This vertical mixing model is further modified to include the effects of surface gravity waves including Coriolis-Stokes forcing, wave breaking, and Langmuir circulations. To simplify the complex interactions between the deterministic and the stochastic phases of flow, we assume a time-invariant turbulent flow field and exclude the hydrodynamic biases due to the effects of ambient mean current. The numerical results show that the inertial particles acquire perturbed oscillations traced out as time-varying sinking/rising orbits in the vicinity of the sea surface under linear and cnoidal waves and acquire a non-looping single arc superimposed with the high-frequency fluctuations beneath the nonlinear solitary waves. Furthermore, we briefly summarise some recipes through the course of this paper on the implementation of the stochastic particle tracking models to realistically describe the drift and suspension of inertial particles throughout the water column.     

Global pattern of the ionospheric response to large-scale internal gravity waves

The global pattern of the ionospheric response to large-scale acoustic gravity waves (LS AGW) has been constructed on the basis of an analysis of the large data set available during the 22 March 1979, magnetic storm. Ground-based ionospheric measurements and in-situ satellite measurements from Cosmos-900 were used in this study together with the Joule heating distribution in the high-latitude ionosphere specifically taken at the maxima of two substorms. The characteristics of the reconstructed planetary pattern of the LS AGW have been analysed in detail. It has been established that the LS AGW effects in the ionosphere in terms of both universal and local time were determined by the pattern of high-latitude atmospheric heating, and that the wave front of the LS AGW during both substorms covered practically all local times, i.e. all longitudes. In addition, it was established that one of the sources of the LS AGW was the thermospheric heating in the day-side cusp region. The local time dependence of the amplitude of the AGW effect in both maximum height, h_mF2, and critical frequency, f_OF2, has been reconstructed for the mid-latitude F2 layer. The AGW effects were clearly separated from the electric field effects related to turnings of the interplanetary magnetic field (IMF) B_Z. In the day-time, electric field effects prevailed over the AGW effects, but during the night-time the amplitudes of these two effects were comparable. In contrast to the common view, f_OF2 variations after the AGW passage had a quasi-sinusoidal character both in the day-time and in the night-time. In the night-time ionosphere a high degree of symmetry was observed for the AGW effects in Northern and Southern hemispheres. During the day-time a significant asymmetry was observed in the American longitudinal sector which was related largely to the peculiarities of the heating pattern in the high-latitude ionospheres of the Northern and Southern hemispheres. These observations demonstrate the complexity of the response of the ionosphere at all latitudes to heating of the auroral region.     

Relative phase and amplitude response to gravity waves of minor species in the thermosphere

Abstract

The response of minor species to gravity waves in the thermosphere varies according to the mass of the species. The relative density perturbation of any minor constituent may be related to the relative density perturbation of the atmosphere due to the waves via a complex function of frequency and wavenumber that may be represented as an amplitude and phase response. Peaks and dips in response and large phase shifts are found that are associated with complex poles and zeroes of the response function. These poles and zeroes depend on background quantities, so that the nature of the response is model dependent; collisions are most important. The effects of background and collisions are examined using numerical computations. Relationships to some satellite observations are discussed.     

重力波在中层大气温度波导中的传播模式研究

本文给出了重力波在中层大气温度波导中的导制传播模型,并在此模型的基础上详细讨论了重力波部分导制传播下的对称模式与非对称模式,导出了不同模式下相应的特征函数和色散方程,进一步用离散的方法对两类色散方程进行了求解;同时还利用二维全隐欧拉格式（FICE）对重力波在温度波导中的传播进行了模拟,模拟的结果也成功地展现了对称与非对称两种传播模式.研究表明,下边界的扰动能量在向上传播进入波导区域后被俘获,形成导制传播.不同周期的初始扰动,在波导内均会形成对称与非对称形式两种模式的导制传播,由于两者的行进速度不一致,最终会引起两种不同模式的分离.数值模拟中重力波的水平行进速度与线性模型预测值非常接近.波导中不同模式下重力波的水平波长与初始扰动的水平波长非常一致,然而波导区域内重力波的频率与初始扰动的频率无关,频率不同的初始扰动会激发出相同频率的重力波对称与非对称导制传播模式.这表明在确定的温度波导中,水平波数才是决定重力波传播特性的决定因素.进一步的分析显示,初始扰动的水平波数-频率分布越接近完全导制传播的色散关系时,温度波导中更易于生成以该种模式部分导制传播的重力波.     

Atmospheric boundary layer over steep surface waves

Turbulent air-sea interactions coupled with the surface wave dynamics remain a challenging problem. The needs to include this kind of interaction into the coupled environmental, weather and climate models motivate the development of a simplified approximation of the complex and strongly nonlinear interaction processes. This study proposes a quasi-linear model of wind-wave coupling. It formulates the approach and derives the model equations. The model is verified through a set of laboratory (direct measurements of an airflow by the particle image velocimetry (PIV) technique) and numerical (a direct numerical simulation (DNS) technique) experiments. The experiments support the central model assumption that the flow velocity field averaged over an ensemble of turbulent fluctuations is smooth and does not demonstrate flow separation from the crests of the waves. The proposed quasi-linear model correctly recovers the measured characteristics of the turbulent boundary layer over the waved water surface.     

重力波包在中层大气温度波导中传播的数值模拟

给出了重力波在中层大气温度波导中的反对称形式导制传播的线性理论模型,并采用二维非线性的数值模型对重力波波包在中层大气温度波导中的传播和演变过程进行了模拟研究.模拟的结果表明,下层大气激发的重力波能量进入波导区域后被俘获形成导制传播.重力波在波导内不停地来回反射,垂直方向的自由传播受到限制,能量在波导内沿着水平方向传输,模拟得到的波参数与理论值相近.重力波包在温度波导中传播时伴随着能量泄漏,且能量泄漏的速率随时间变缓,最终总有部分能量被限制在波导区域.重力波在水平方向上传播几百公里后,依然维持着良好的波结构,同时数值模拟也给出了重力波在波导区域内能量密度的时空分布.     

Interpreting the sea surface temperature warming trend in the Yellow Sea and East China Sea

Previous studies have demonstrated that the low-frequency sea surface temperature (SST) variability in the Yellow Sea and East China Sea (YECS) is linked to large-scale climate variability, but explanations on the mechanisms vary. This study examines the low-frequency variability and trends of some atmospheric and oceanic variables to discuss their different effects on the YECS warming. The increasing temperature trend is also observed at a hydrographic section transecting the Kuroshio. The increasing rate of ocean temperature decreases with depth, which might result in an increase in vertical stratification and a decrease in vertical mixing, and thus plays a positive role on the YECS warming. The surface net heat flux (downward positive) displays a decreasing trend, which is possibly a result of the YECS warming, and, in turn, inhibits it. Wind speeds show different trends in different datasets, such that its role in the YECS warming is uncertain. The trends in wind stress divergence and curl have large uncertainties, so their effects on SST warming are still unclear. The Kuroshio heat transport calculated in this study, displays no significantly increasing trend, so is an unlikely explanation for the SST warming. Limited by sparse ocean observations, sophisticated assimilative climate models are still needed to unravel the mechanisms behind the YECS warming.     

Geological structure of the yellow sea area from regional gravity and magnetic interpretation

Regional gravity and aeromagnetic data covering the area of 32°- 38° N, 118°-127° E at the scale of 1：1,000,000 are coordinated and integrated in a synthetic study of the South China Yellow Sea and adjacent areas. Depth to magnetic crystalline basement and its structure are determined by magnetic anomaly inversion. Depth to and thickness of the Paleozoic rock are also revealed by gravity anomaly inversion with constrains of the basement and known seismic information from several profiles. Structure units, main faults, basin boundaries, and sub-suppressions are outlined on the basis of gravity data interpretation.     

Island sheltering of surface gravity waves: model and experiment

A field experiment is used to evaluate a numerical model of the sheltering of gravity waves by islands offshore of the Southern California region. The sheltering model considered here includes only the effects of island blocking and wave refraction over the island bathymetry. Wave frequency and directional spectra measured in the deep ocean (unsheltered region west of the islands) were used as input to the sheltering model and compared with coastal observations. An airborne L-band synthetic aperture radar was used to image the directional properties of the waves in the deep ocean. In addition to the unsmoothed spectra, a unimodal directional spectrum model obtained from fits to the radar spectra was also employed to suppress the high noise level of this system. Coastal measurements were made in about 10 m depth at Torrey Pines Beach with a high resolution array of pressure sensors. The model predictions and data at Torrey Pines Beach agree well in a limited frequency range (0.082 to 0.114 Hz) where the unimodal deep ocean model is appropriate. The prediction that unimodal northern swell in the deep ocean results in a bimodal directional spectrum at Torrey Pines Beach is quantitatively verified. The northern peak of the bimodal spectra is due to waves coming through the window between San Clemente and San Miguel-Santa Rosa Islands. The southerly peak is due to wave refraction over Cortez and Tanner Banks. For lower frequency waves, the effects of strong refraction in the island vicinity are shown qualitatively. Refraction can theoretically supply up to approximately 10% of the deep ocean energy that is otherwise blocked at this site. The modifications of the island shadows due to wave refraction become theoretically negligible for wave frequencies 0.11Hz. Also, local wave generation effects, which are not included in this sheltering model, are shown to be occasionally important for waves with frequencies 0.12Hz.     

Experimental and analytical investigation of the response of a mud layer to solitary waves

Ocean Dynamics - Although the different aspects of wave-mud interaction have been studied by many researchers, few studies have been conducted on the effect of solitary wave on the particle...     

Simplified longitudinal seismic response of tunnels linings subjected to surface waves

This paper presents an engineering approach for analysing the longitudinal behaviour of tunnels subjected to earthquakes. The tunnel is modelled as a Timoshenko beam connected to the far soil by means of continuous elastic support (Winkler model). Seismic free-field inputs, such as those caused by surface waves travelling parallel to the tunnel axis, were imposed at the base of the springs of the Winkler model, generating bending moments and shear forces on the cross-sections of the tunnel. Closed-form expressions of the tunnel displacements, shear forces, and bending moments were determined at any tunnel section in terms of the seismic excitation, tunnel geometry and material properties, and subgrade reaction modulus of the soil. A dimensional analysis was carried out to ascertain directly the maximum tunnel displacement, bending moment, and shear force.     

Impacts of climate changes on ocean surface gravity waves over the eastern Canadian shelf

A numerical study is conducted to investigate the impact of climate changes on ocean surface gravity waves over the eastern Canadian shelf (ECS). The “business-as-usual” climate scenario known as Representative Concentration Pathway RCP8.5 is considered in this study. Changes in the ocean surface gravity waves over the study region for the period 1979–2100 are examined based on 3 hourly ocean waves simulated by the third-generation ocean wave model known as WAVEWATCHIII. The wave model is driven by surface winds and ice conditions produced by the Canadian Regional Climate Model (CanRCM4). The whole study period is divided into the present (1979–2008), near future (2021–2050) and far future (2071–2100) periods to quantify possible future changes of ocean waves over the ECS. In comparison with the present ocean wave conditions, the time-mean significant wave heights (H _s ) are expected to increase over most of the ECS in the near future and decrease over this region in the far future period. The time-means of the annual 5% largest H _s are projected to increase over the ECS in both near and far future periods due mainly to the changes in surface winds. The future changes in the time-means of the annual 5% largest H _s and 10-m wind speeds are projected to be twice as strong as the changes in annual means. An analysis of inverse wave ages suggests that the occurrence of wind seas is projected to increase over the southern Labrador and central Newfoundland Shelves in the near future period, and occurrence of swells is projected to increase over other areas of the ECS in both the near and far future periods.     

Distribution and pollution assessment of heavy metals in surface sediments in the Yellow Sea

Heavy metal concentrations in surface sediments at 56 stations during two cruises in the Yellow Sea in summer and winter, 2011 were analyzed by inductively coupled plasma-mass spectrometry. The pollution status was assessed via the Geoaccumulation index and Hankanson potential ecological risk index. Higher concentrations of heavy metals (except for Mn) were found in the central Southern Yellow Sea and the western Northern Yellow Sea. The higher contents of Mn were much closer to Shandong Peninsula. Correlation analyses indicated that Pb, Cu, Fe, Ni, Zn and Co probably had the same origin and were controlled by grain size and total organic carbon. Pollution assessment showed that most areas of the Yellow Sea were not or lowly contaminated with the exception of the northwest and south parts of the Southern Yellow Sea showing Cd-contamination. The pollution status of the Yellow Sea in summer was worse than that in winter.     

Mass transport in a thin layer of power-law mud under surface waves

The mass transport velocity in a two-layer system is studied theoretically. The wave motion is driven by a periodic pressure load on the free water surface, and mud in the lower layer is described by a power-law rheological model. Perturbation analysis is performed to the second order to find the mean Eulerian velocity. A numerical iteration method is employed to solve the non-linear governing equation at the leading order. The influence of rheological properties on fluid motion characteristics including the flow field, the surface displacement, the mass transport velocity, and the net discharge rates are investigated based on theoretical results. Theoretical analysis shows that under the action of interfacial shearing, a recirculation structure may appear near the interface in the upper water layer. A higher mass transport velocity at the interface does not necessarily mean a higher discharge rate for a pseudo-plastic fluid mud.     

Modifications of the dispersion relations for surface waves in a layer on a half-space

Summary Dispersion relations for Love and Rayleigh waves in a layer on a half-space are modified by introducing quadratic slownesses instead of velocities. The advantages of this approach are demonstrated on analytical formulae for computing the group velocity.     

On the propagation of certain types of surface waves in a non-homogeneous elastic layer

Summary The propagation of a certain type of surface waves in a non-homogeneous elastic layer of finite thickness lying in welded contact with a semi-infinite homogeneous elastic medium has been investigated in this paper. The surface wave is characterised by the fact that the dilatation and the vertical displacement component are both zero. It has been proved that such disturbance can propagate only when the phase velocity of the wave is greater than the shear wave velocity and the group velocity is inversely proportional to the phase velocity.     

南海三维热岩石圈结构及地表构造响应

为系统了解南海海底热流特征和岩石圈热结构与地表构造的关系,本文首先以最新海底热流为约束,采用三维数值模拟方法计算得到整个南海的"热"岩石圈底部边界（底界）的埋深,然后开展详细分析壳、幔热流贡献以及岩石圈热结构与地表不同构造单元之间的对应关系.结果显示南海岩石圈底部边界埋深为30~100 km,最浅埋深集中分布在南海海盆水深大于3000 m的区域（如南海西北次海盆、西南次海盆和东部次海盆西侧）,以及南海北部陆缘中部、莺歌海盆地、中建南盆地和南薇西盆地,体现了南海深部结构对地表不同构造单元的控制作用;靠近海盆东侧的马尼拉海沟具有较低的海底热流,对应岩石圈底界埋藏也较深,可能与南海和古南海洋壳俯冲有关;吕宋岛弧（119°E-122°E和13°N-17°N）表现出较高的海底热流,同时具有较浅的热岩石圈底界埋深,可能与南海板块和菲律宾海板块双向俯冲有关.     

Variations in the characteristics of acoustic gravity waves according to simulation data

The characteristics of different-scale acoustic gravity waves (wavelengths of 100–1200 km, periods of 10–50 min) under different geophysical conditions have been studied using a numerical model for calculating the vertical structure of these waves in a nonisothermal atmosphere in the presence of an altitudedependent background wind and in a situation when molecular dissipation is taken into account. It has been established that all considered acoustic gravity waves (AGWs) effectively reach altitudes of the thermosphere. The character of the amplitude vertical profile depends on the AGW scales. The seasonal and latitudinal differences in the AGW vertical structure depend on the background wind and temperature. A strong thermospheric wind causes the rapid damping of medium-scale AGWs propagating along the wind. Waves with long periods to a lesser degree depend on dissipation in the thermosphere and can penetrate to high altitudes. A change in the geomagnetic activity level affects the background wind vertical distribution at high latitudes, as a result of which the AGW vertical structure varies.     

PCBs and its coupling with eco-environments in Southern Yellow Sea surface sediments

In this paper, the spatial distribution and source of the PCBs in surface sediments of the Southern Yellow Sea (SYS) and influencing factors, such as the sediment characteristics (components, relative proportions and total organic carbon contents), and hydrodynamic conditions were analyzed. PCB concentrations in the surface sediments ranged from 518-5848 pg/g, with average values of 1715 pg/g decreasing sharply compared to last year. In the study area, the PCB pollution level in the middle area was the highest, followed by that of the east coast and the west coast, respectively. Although the PCB level in the coastal areas was lower than that in the middle areas, it was proven in our study that the Yellow Sea obtained PCBs by virtue of river inputs. There was a positive and pertinent correlation between the clay proportion and PCB concentrations, and the increase of the PCB concentrations was directly proportional to the increase of TOC contents, with r=0.61, but it was contrary to the sediment grain size. Consequently, the factors controlling PCB distribution had direct or indirect relationships with sediment grain size; moreover, the hydrodynamic conditions determined the sediment components and grain size. In conclusion, hydrodynamic conditions of the Yellow Sea were the most important influencing factors effecting the distribution of PCBs in the surface sediments of the SYS.     

The instability of internal gravity waves to localised disturbances

The instability of an internal gravity wave due to nonlinear wave-wave interaction is studied theoretically and numerically. Three different aspects of this phenomenon are examined. 1. The influence of dissipation on both the resonant and the nonresonant interactions is analysed using a normal mode expansion of the basic equations. In particular, the modifications induced in the interaction domain are calculated and as a result some modes are shown to be destabilised by dissipation. 2. The evolution of an initial unstable disturbance of finite vertical extent is described as the growth of two secondary wave packets travelling at the same group velocity. A quasi-linear correction to the basic primary wave is calculated, corresponding to a localised amplitude decrease due to the disturbance growth. 3. Numerical experiments are carried out to study the effect of a basic shear on wave instability. It appears that the growing secondary waves can have a frequency larger than that of the primary wave, provided that the shear is sufficient. The instability of waves with large amplitude and long period, such as tides or planetary waves, could therefore be invoked as a possible mechanism for the generation of gravity waves with shorter period in the middle atmosphere.