Study on internal waves generated by tidal flow over critical topography

Resonance due to critical slope makes the internal wave generation more effectively than that due to supercritical or subcritical slopes(Zhang et al., 2008). Submarine ridges make a greater contribution to ocean mixing than continental margins in global oceans(Müller, 1977; Bell, 1975; Baines, 1982; Morozov, 1995). In this paper, internal wave generation driven by tidal flow over critical topography is examined in laboratory using Particle Image Velocimetry(PIV) and synthetic schlieren methods in synchrony. Non-tidal baroclinic velocities and vertical isopycnal displacements are observed in three representative regions, i.e., critical, outward-propagating, and reflection regions. Temporal and spatial distributions of internal wave rays are analyzed using the time variations of baroclinic velocities and vertical isopycnal displacement, and the results are consistent with those by the linear internal wave theory. Besides, the width of wave beam changes with the outward propagation of internal waves. Finally, through monitoring the uniformly-spaced 14 vertical profiles in the x-z plane, the internal wave fields of density and velocity fields are constructed. Thus, available potential energy, kinetic energy and energy fluxes are determined quantitatively. The distributions of baroclinic energy and energy fluxes are confined along the internal wave rays. The total depth averaged energy and energy flux of vertical profiles away from a ridge are both larger than those near the ridge.     

A TWO-DIMENSIONAL VORTICITY MODEL OF INTERNAL TIDES GENERATED ON THE CONTINENTAL SHELF/SLOPE

A two-dimensional spectral-difference mode (with vorticity and density equations) of internal tides isdeveloped for studying the genration and propagration of internal tides generated at the continentalshelf/slope. In general, internal tides propagate seaward in deep sea regions and shoreward on the shelf,and are dissipated rapidly. When the Vaisala frequency decreases vertically, waves may be mostly limited to thecontinental slope region. in deep sea region, motions may have either boam-like structure or modal structure,depending on the stratification strerigth and structure, whereas a modal structure may always exist onthe shelf. Waves show strong bottom intensification on the slope when strong stratification exists on thebottom. The barotropic tidal advection may affed the temporal character of internal tides at thecontinental slope, shelf break and shelf regions. but may have little influence on the energy density and energy flux of internal tides. ln the case of strong stratification, waverforms of     

Internal Waves Generated by Tidal Flows over a Triangular Ridge with Critical Slopes

Direct numerical simulations are performed to investigate the generation of internal waves by tide-topography interaction in a lab-scale model. The bottom topography is a triangular ridge with two critical slopes. With increasing tidal forcing, subharmonic instabilities are identified, which cause internal wave beams to become unstable and turbulent. Kinetic energy densities in the upward going beams from the ridge top are stronger than those from the ridge bottom, whereas the reverse is true for the energy flux. This disparity between energy and energy flux is due to the existence of strong pressure disturbances near the ridge bottom. On each side of the critical ridge, there exists an amphidromic point, from which internal wave beams are emitted in opposite directions. The calculated energy conversion rate scales linearly with the square of the forcing amplitude and agrees within 13% of theoretical prediction, even when turbulence occurs. The fraction of radiated baroclinic energy becomes saturated in the range of low excursion parameter considered, which agrees with the behavior in large-scale systems wherein mixing parameterizations must be used. The present work enriches the studies on the generation of internal waves over a critical triangular ridge.     

Inversion analysis of the shear strength parameters for a high loess slope in the limit state

It is difficult to obtain reliable shear strength parameters for the stability analysis and evaluation of high loess slopes.Hence,this paper determines slope elements and physical parameter of 79 slopes with heights of[40,120]m based on the measured loess slopes in the Ganquan and Tonghuang subregions of northern Shaanxi Province,China.In the limit state of the loess slope(stability factor Fs is 1.0),a fitting equation for the slope height and width is established.Then,the model is developed by stability analysis software-SLOPE/W,and the comprehensive shear strength parameters corresponding to different slope heights of the high slope in the study area are obtained by inversion using the Morgenstern-Price method.The results show that when the height of the slope increases,the cohesion c increases in the soil,and the internal friction angle j decreases.This change is consistent with the characteristics reflected in the composition and physical properties of the slopes,and the comprehensive strength parametric curves are very similar between the Ganquan subregion and the Tonghuang subregion.A landslide that occurred in Miaodian-zaitou of Jingyang County,Shaanxi Province,is selected to verify this inversion method,and the results show that the calculated shear strength parameters of the inversion are consistent with the measured value of the actual slope.     

Temporal variations in internal tide multimodal structure on the continental shelf,South China Sea

Temporal variations in multimodal structures of diurnal (D ₁) and semidiurnal (D ₂) internal tides were investigated on the continental slope of the Dongsha Plateau, based on 2-month moored acoustic Doppler current profiler observations. Harmonic analysis indicated that the D ₁ components (K ₁ and O ₁) dominated the internal tide field. The vertical structure of the K ₁ constituent presented a first-mode structure while the M ₂ constituent seemed to exhibit a high-mode structure. Amplitude spectra analysis of the current data revealed differences in baroclinic current amplitudes between different water depths. Temporal variations in modal structures ware analyzed, based on the D ₁ and D ₂ baroclinic tides extracted from the baroclinic velocity field with band-pass filters. Analysis showed that the magnitude of the D ₁ internal tide current was much larger than the D ₂ current, and temporal variations in the modal structure of the D ₁ internal tide occurred on an approximately fortnightly cycle. The EOF analyses revealed temporal transformation of multimodal structures for D ₁ and D ₂ internal tides. The enhancement of the D 1 internal tide was mainly due to the superposition of K ₁ and O ₁, according to the temporal variation of coherent kinetic energy.     

浑厚山体随高程增加由表及里地震动响应研究

为探究浑厚山体不同高程由表及里地震动响应规律,以冷竹关山体为例,采用离散元软件建立地形与风化介质组合模型,并从底部边界输入汶川地震波信号,研究该山体两侧边坡的内外动力响应规律。结果表明,随着高程的增加,靠大渡河一侧边坡坡体内与坡表的加速度放大系数均表现出先增大后减小的节律性变化,在近坡顶时增大较快并达到最大值;靠近瓦斯沟一侧边坡坡表受地形起伏的影响,加速度放大系数存在凸坡放大、凹坡减小的特征;相同高程,随边坡由表及里深度的增加,加速度放大系数表现为逐步减小,当距坡表150~200 m时放大曲线趋于平缓;随高程的增加,加速度放大系数由表及里的减小速度变缓,且放大曲线收敛平缓的深度增大;随岩体风化程度的增加,岩体介质波速降低,共振效应使得加速度响应增大,与此同时,斜坡地形与介质组合效应使得坡表峰值加速度放大系数在2.0附近。     

Characteristics of viscous debris flow in a drainage channel with an energy dissipation structure

A new type of drainage channel with an energy dissipation structure has been proposed based on previous engineering experiences and practical requirements for hazard mitigation in earthquakeaffected areas. Experimental studies were performed to determine the characteristics of viscous debris flow in a drainage channel of this type with a slope of 15%. The velocity and depth of the viscous debris flow were measured, processed, and subsequently used to characterize the viscous debris flow in the drainage channel. Observations of this experiment showed that the surface of the viscous debris flow in a smooth drainage channel was smoother than that of a similar debris flow passing through the energy dissipation section in a channel of the new type studied here. However, the flow patterns in the two types of channels were similar at other points. These experimental results show that the depth of the viscous debris flow downstream of the energy dissipation structure increased gradually with the length of the energy dissipation structure. In addition, in the smooth channel, the viscous debris-flow velocity downstream of the energy dissipation structure decreased gradually with the length of the energy dissipation structure. Furthermore, the viscous debris-flow depth and velocity were slightly affected by variations in the width of the energy dissipation structure when the channel slope was 15%. Finally, the energy dissipation ratio increased gradually as the length and width of the energy dissipation structure increased; the maximum energy dissipation ratio observed was 62.9% (where B = 0.6 m and L/w = 6.0).     

黄土高原边坡特征与破坏规律的分区研究

为了研究黄土高原自然边坡的特征及破坏规律,按山系与水系或水系的分水岭、地貌单元、地层岩性特征等条件,将黄土高原划分为8个区:临洮—永靖区、天水—通渭区、兰州—会宁区、陇东区、靖边—安塞区、隰县—离石区、甘泉—吉县区和汾渭区.根据极限状态边坡的4个野外判别标准,测量了8个区510个自然极限状态黄土边坡断面,分区采用指数模型回归边坡坡高与坡宽的相关关系,计算各区20、50、100m坡高的边坡稳定系数和失效概率.结果表明:黄土高原的边坡特征与破坏形式具有分区特征,且南北差异性明显.临洮—永靖区边坡坡高与坡宽呈线性关系,表明该区边坡坡度不随坡高变化,边坡稳定性受内摩擦角控制;兰州—会宁区和靖边—安塞区高坡陡,低坡缓,高坡不稳定,易发生错落式滑坡;天水—通渭区、甘泉—吉县区和汾渭区高坡缓,低坡陡,稳定性计算结果显示高坡和低坡都较为稳定,但由于地层结构和地貌的特点,高边坡易发生低速蠕变型滑坡或高速远程滑坡;陇东区边坡整体上较为稳定;隰县—离石区受黏粒含量较高的Q1地层控制,高边坡稳定性较差;50m左右坡高的黄土边坡稳定性对强度指标内聚力、内摩擦角的敏感度都高,易于失稳.     

Numerical investigation of wake-collapse internal waves generated by a submerged moving body

The state-of-the-art OpenFOAM technology is used to develop a numerical model that can be devoted to numerically investigating wake-collapse internal waves generated by a submerged moving body.The model incorporates body geometry,propeller forcing,and stratification magnitude of seawater.The generation mechanism and wave properties are discussed based on model results.It was found that the generation of the wave and its properties depend greatly on the body speed.Only when that speed exceeds some critical value,between 1.5 and 4.5 m/s,can the moving body generate wake-collapse internal waves,and with increases of this speed,the time of generation advances and wave amplitude increases.The generated wake-collapse internal waves are confirmed to have characteristics of the second baroclinic mode.As the body speed increases,wave amplitude and length increase and its waveform tends to take on a regular sinusoidal shape.For three linearly temperature-stratified profiles examined,the weaker the stratification,the stronger the wake-collapse internal wave.     

Influence of mountain-valley morphology on in-situ stress distribution

The in-situ stress field is a key factor controlling the successful construction of a large number of underground structures in mountainous areas, and is intensively affected by the mountainvalley topography. The effects of mountain-valley morphology(the width of the mountain top platform, mountain height, slope angle, and width of the valley bottom) on the distribution of the in-situ stress field were analyzed and interpreted using numerical modeling techniques, where the spatial distribution and maximum values of the horizontal and vertical stresses were analyzed. The results showed that there existed a critical value of the topographic influence depth, where the in-situ stress distribution varied significantly as mountain-valley morphology, after which the influence diminishes. Tectonic action has a more remarkable influence on the in-situ stress distribution than gravitational action under the same mountain-valley morphology. Moreover, the relationships between the magnitudes of these stress components and the morphology variables are described using empirical formulas, which can be directly applied to different topographies to rapidly achieve a rational estimation. The findings of this study can be very useful for quickly understanding the in-situ stress distribution and as stress measurement guidelines.     

Currents and mixing in the northern South China Sea

We investigated the vertical distribution of current velocity data of the entire water column at a site on the continental shelf of the northern South China Sea (SCS) from August 4 to September 6, 2007, and found that the characteristics of barotropic and baroclinic tides are mainly diurnal. During the observation period, we also estimated the mixing before and after the passage of Typhoon Pabuk. We found that the internal-wave-scale dissipation rate, the turbulent dissipation rate, and the mixing rate in every water layer increased by about an order of magnitude after the typhoon passage. We analyzed a case of abrupt strong current and calculated the mixing rate before, during, and after the typhoon event. The results show that the internal-wave-scale dissipation rate and the mixing rate in every water layer increased by about two orders of magnitude during the event, while the turbulent dissipation rate increased by about an order of magnitude. Passage of the abrupt strong current could also have increased the mixing rate of affected seawater by more than an order of magnitude. However, the passage of the typhoon differed in that there was an increase in mixing only in the lower layer where the abrupt strong current was particularly strong. The variation of the mixing rate may help us to understand the effects of typhoons and abrupt strong currents on the mixing of seawater.     

AN EXACT SOLUTION OF SUB-INERTIAL FREQUENCY BAROCLINIC WAVES

By using a coordinate transformation, an exact solution of internal tides with sub-inertial frequency isobtained when the bottom slope is linear and the Vaisala frequency is constant. Accordingly thedispersion relations of free waves are presented. This solution is suitable for general coastal low-frequencybaroclinic waves with zero alongshore wavenumber.     

Spurious dianeutral mixing in a global ocean model using spherical centroidal voronoi tessellations

In order to quantitatively evaluate the spurious dianeutral mixing in a global ocean model MPAS-Ocean(Model for Prediction Across Scales) using a spherical centroidal voronoi tessellations developed jointly by the National Center for Atmospheric Research and the Los Alamos National Laboratory in the United States, we choose z* vertical coordinate system in MPAS-Ocean, in which all physical mixing processes, such as convection adjustment and explicit diffusion parameter schemes, are omitted, using a linear equation of state. By calculating the Reference Potential Energy(RPE), front revolution position, time rate of RPE change, probability density function distribution and dimensionless parameter χ, from the perspectives of resolution, viscosity, Horizontal Grid Reynolds Number(HGRN), Re?, and momentum transmission scheme, using two ideal cases, overflow and baroclinic eddy channel, we qualitatively analyze the simulation results by comparison with the three non-isopycnal models in Ilicak et al.(2012), i.e., MITGCM, MOM, and ROMS. The results show that the spurious dianeutral mixing in the MPAS-Ocean increases over time. The spurious dianeutral transport is proportional to the HGRN directly and is reduced by increasing the lateral viscosity or using a finer resolution to control HGRN. When the HGRN is less than 10, spurious transport is reduced significantly. When using the proper viscosity closure, MPAS-Ocean performs better than MITGCM and MOM, closely to ROMS, in the 2D case without rotation, and much better than the above-mentioned three ocean models under the condition of 3D space with rotation due to the cell area difference between the hexagon cell and the quadrilateral cell with the same resolution. Both the Zalesak(1979) flux corrected transport scheme and Leith closure in MPAS-Ocean play an excellent role in reducing spurious dianeutral mixing. The performance of Leith scheme is preferable to the condition of three-dimensional baroclinic eddy.     

Propagation of internal waves up continental slope and shelf

In a two-dimensional and linear framework, a transformation was developed to derive eigensolutions of internal waves over a subcritical hyperbolic slope and to approximate the continental slope and shelf. The transformation converts a hyperbolic slope in physical space into a flat bottom in transform space while the governing equations of internal waves remain hyperbolic. The eigensolutions are further used to study the evolution of linear internal waves as it propagates to subcritical continental slope and shelf. The stream function, velocity, and vertical shear of velocity induced by internal wave at the hyperbolic slope are analytically expressed by superposition of the obtained eigensolutions. The velocity and velocity shear increase as the internal wave propagates to a hyperbolic slope. They become very large especially when the slope of internal wave rays approaches the topographic slope, which is consistent with the previous studies.     

Rotary spectrum analysis of tidal current in the southern Yellow Sea

A complete set of one-month Acoustic Doppler Profiler (ADP) current data at a station in the southern Yellow Sea (SYS) is analyzed using the rotary spectrum method. The results revealed different rotary properties between barotropic and baroclinic tidal currents. The barotropic and baroclinic tidal currents rotate elliptically counter-clockwise and clockwise, respectively. Meanwhile, baroclinic bottom tidal currents are almost along-isobath. The baroclinic cross-isobath velocities attenuate quickly at the bottom, implying important effects of bottom topography on the cross-isobath motions.     

Observation and numerical experiments for drag coefficient under typhoon wind forcing

This paper presents a study on drag coefficients under typhoon wind forcing based on observations and numerical experiments. The friction velocity and wind speed are measured at a marine observation platform in the South China Sea. Three typhoons: SOULIK(2013), TRAMI(2013) and FITOW(2013) are observed at a buoy station in the northeast sea area of Pingtan Island. A new parameterization is formulated for the wind drag coefficient as a function of wind speed. It is found that the drag coefficient(Cd) increases linearly with the slope of 0.083′10~(-3) for wind speed less than 24 m s~(-1). To investigate the drag coefficient under higher wind conditions, three numerical experiments are implemented for these three typhoons using SWAN wave model. The wind input data are objective reanalysis datasets, which are assimilated with many sources and provided every six hours with the resolution of 0.125?×0.125?. The numerical simulation results show a good agreement with wave observation data under typhoon wind forcing. The results indicate that the drag coefficient levels off with the linear slope of 0.012′10~(-3) for higher wind speeds(less than 34 m s~(-1)) and the new parameterization improvese the simulation accuracy compared with the Wu(1982) default used in SWAN.     

A study on the baroclinic structure of the abyssal circulation

In this paper, the linear continuously stratified ntxiel of the abyssal circulation proposed by Pedlosky (1992) was emended to include the second order term - (γθm) in the vertical turbulent mixing pararneterization of - (ω‘θ‘）x = kvθxx- θxxx, in which kv is a vertical diffusion coefficient,and γ is the second order coefficient of turbulent mixing (or simply called γ-term and γ＜0 is only allcmed). The influence of the γ-tema on the baroclinic structure of the abyssal circulation driven by up-welling out of the abyss was investigated, It was found that the γ-term has a noticeable influence on thebaroclinic structure of the upwelling driven abyssal circulation. For uniform upwelling, it favors the baroelinic layering of the abyssal circulation in the eastern part of the basin, but prevents the layering in thewest. In addition, this parameter was found to decrease the vertically averaging meridional velocity of theabyssal circulation frcrn the west to the east on the southern boundary. For upwelling localized near theeastern boundary, the γ - term favors baroclinic layering of the abyssal circulation in the whole basin.Especially, on the southern boundary the γ - term could strengthen the vertically averaging meridional ve-loeity in the west, but gready weaken it in the east. The ntntel presented here might be considered as anextension of the Pedlosky baroclinic ntntel of the abyssal circulation.     

Current and thermohaline characteristics of the Arabian Sea during January 1998

Based on a ship survey during January 1998, the characteristics of the flow, the thermohaline properties and the volume transport of the Arabian Sea are discussed. A strong westward flow exists between 10.5?N and 11?N, part of which turns to the south as the Somali current near the coast at about 10?N and the rest turns north. At the passage between the African continent and the So- cotra Island, the northern branch separates into two flows: the left one enters the passage and the right one flows eastward along the southern slope of the island. Off the island the flow separates once more, most of it meandering northeast and a small fraction flow- ing southeast. Volume transport calculation suggests that the tidal transport is one or two orders of magnitude smaller than the total transport in this region and it becomes more important near the coast. The average velocity of the flow in the upper layer (0-150 m) is about 20 cm s-1, with a maximum of 53 cm s-1 appearing east of the Socotra Island, and the subsurface layer (200-800 m) has an aver- age velocity of 8.6 cm s-1; the velocity becomes smaller at greater depths. The depth of the seasonal thermocline is about 100 m, above which there is a layer with well mixed temperature and dissolved oxygen. High-salinity and oxygen-rich water appears near the surface of the northern Arabian Sea; a salinity maximum and oxygen minimum at 100 m depth along 8?N testifies the subduction of surface water from the northern Arabian Sea. Waters from the Red Sea and the Persian Gulf also influence the salinity of the area.     

Effect of environmental gradients on the quantity and quality of fallen logs in Tsuga longibracteata forest in Tianbaoyan National Nature Reserve, Fujian province, China

We investigated the quantity and quality 0f fallen l0gs in different Tsuga l0ngibracteata f0rest c0mmunities in the Tianba0yan Nati0nal Nature Reserve. We used redundancy analysis t0 determine the spatial distributi0n 0f fallen l0gs in the different f0rest c0mmunities and t0 analyze the relati0nships am0ng stand structure, t0p0graphic fact0rs and human disturbance. The v0lume, c0vered area, mean l0g length and number 0f fallen l0gs differed significantly am0ng f0rest types （P 〈 0.05）, but mean diameter at breast height sh0wed n0 significant difference （P 〉 0.05）. The l0g v0lume and c0vered area in different f0rest types sh0wed the f0ll0wing trend： T. l0ngibracteata pure f0rest 〈 T. l0ngibracteata ＋ Olig0staehyum scabrifl0rur 〈 T. l0ngibraeteata ＋ hardw00d 〈 Rh0d0dendr0n simiarum ＋ T. l0ngibraeteata 〈 T. l0ngibraeteata ＋ Phyll0stachys heter0cycla pubescens. The spatial distributi0n patterns 0f l0gs quantity and quality indicated that l0g v0lume and c0vered area were str0ngly affected by envir0nmental fact0rs in the f0ll0wing 0rder： human disturbance 〉 elevati0n 〉 sl0pe p0siti0n 〉 b0le height 〉 tree height 〉 sl0pe aspect 〉 density 〉 basal area 〉 sl0pe gradient. The relative c0ntributi0n 0f envir0nmental variables 0n the t0tal variance was t0p0graphy （76%） 〉 disturbance （42%） 〉 stand structure （35%）. T0p0graphy and disturbance c0mbined explained 8.2% 0f the variance. Fallen l0~s auantitv and aualitvwere negatively related t0 elevati0n and sl0pe p0siti0n, and p0sitively ass0ciated t0 human disturbance. The l0g v0lume decreased fr0m n0rthern t0 s0uthern sl0pes. Envir0nmental fact0rs had the highest impact 0n class I （slightly decayed）, and l0west impact 0n class V （highly decayed）.     

Numerical simulation of wave field in the South China Sea using WAVEWATCH III

Wave fi elds of the South China Sea(SCS) from 1976 to 2005 were simulated using WAVEWATCH III by inputting high-resolution reanalysis wind fi eld datasets assimilated from several meteorological data sources. Comparisons of wave heights between WAVEWATCH III and TOPEX/Poseidon altimeter and buoy data show a good agreement. Our results show seasonal variation of wave direction as follows: 1. During the summer monsoon(April–September), waves from south occur from April through September in the southern SCS region, which prevail taking about 40% of the time; 2. During the winter monsoon(December–March), waves from northeast prevail throughout the SCS for 56% of the period; 3. The dominant wave direction in SCS is NE. The seasonal variation of wave height H s in SCS shows that in spring, H s ≥1 m in the central SCS region and is less than 1 m in other areas. In summer, H s is higher than in spring. During September–November, infl uenced by tropical cyclones, H s is mostly higher than 1 m. East of Hainan Island, H s 2 m. In winter, H s reaches its maximum value infl uenced by the north-east monsoon, and heights over 2 m are found over a large part of SCS. Finally, we calculated the extreme wave parameters in SCS and found that the extreme wind speed and wave height for the 100-year return period for SCS peaked at 45 m/s and 19 m, respectively, SE of Hainan Island and decreased from north to south.