A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.

Factors influencing the climatological mixed layer depth in the South China Sea: numerical simulations

The mixed layer depth (MLD) in the upper ocean is an important physical parameter for describing the upper ocean mixed layer. We analyzed several major factors influencing the climatological mixed layer depth (CMLD), and established a numerical simulation in the South China Sea (SCS) using the Regional Ocean Model System (ROMS) with a high-resolution (1/12A degrees x1/12A degrees) grid nesting method and 50 vertical layers. Several ideal numerical experiments were tested by modifying the existing sea surface boundary conditions. Especially, we analyzed the sensitivity of the results simulated for the CMLD with factors of sea surface wind stress (SSWS), sea surface net heat flux (SSNHF), and the difference between evaporation and precipitation (DEP). The result shows that of the three factors that change the depth of the CMLD, SSWS is in the first place, when ignoring the impact of SSWS, CMLD will change by 26% on average, and its effect is always to deepen the CMLD; the next comes SSNHF (13%) for deepening the CMLD in October to January and shallowing the CMLD in February to September; and the DEP comes in the third (only 2%). Moreover, we analyzed the temporal and spatial characteristics of CMLD and compared the simulation result with the ARGO observational data. The results indicate that ROMS is applicable for studying CMLD in the SCS area.

Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).

Growth of zero-mean-shear turbulent-mixed layer with suspended particles in a 2-layer fluid system

The growth behaviour of zero-mean-shear turbulent-mixed layer containing suspended solid particles has been studied experimentally and analysed theoretically in a two-layer fluid system. The potential model for estimating the turbulent entrainment rate of the mixed layer has also been suggested, including the results of the turbulent entrainment for pure two-layer fluid. The experimental results show that the entrainment behaviour of a mixed layer with the suspended particles is well described by the model. The relationship between the entrainment distance and the time, and the variation of the dimensionless entrainment rate E with the local Richardson number Ri1 for the suspended particles differ from that for the pure two-layer fluid by the factors-eta-1/5 and eta-1, respectively, where eta = 1 + sigma-0-DELTA-rho/DELTA-rho-0.

Interannual horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean: assimilation versus TAO analysis

The interannual anomalies of horizontal heat advection in the surface mixed layer over the equatorial Pacific Ocean in an assimilation experiment are studied and compared with existing observational analyses. The assimilation builds upon a hindcast study that has produced a good simulation of the observed equatorial currents and optimizes the simulation of the Reynolds sea surface temperature (SST) data. The comparison suggests that the assimilation has improved the simulation of the interannual horizontal heat advection of the surface mixed layer significantly. During periods of interrupted current measurements, the assimilation is shown to produce more meaningful anomalies of the heat advection than the interpolation of the observational data does. The assimilation also shows that the eddy heat flux due to the correlation between high-frequency current and SST variations, which is largely overlooked by the existing observational analyses, is important for the interannual SST balance over the equatorial Pacific. The interannual horizontal heat advection anomalies are found to be sensitive to SST errors where oceanic currents are strong, which is a challenge for ENSO prediction. The study further suggests that the observational analyses of the tropical SST balance based on the TAO and the Reynolds SST data contain significant errors due to the large gradient errors in the Reynolds SST data, which are amplified into the advection anomalies by the large equatorial currents.

Physical mechanism and numerical simulations of surface layer temperature inversion in tropical ocean

The one-dimensional Kraus-Turner mixed layer model improved by Liu is developed to consider the effect of salinity and the equations of temperature and salinity under the mixed layer. On this basis, the processes of growth and death of surface layer temperature inversion is numerically simulated under different environmental parameters. At the same time, the physical mechanism is preliminarily discussed combining the observations at the station of TOGA-COARE 0 degrees N, 156 degrees E. The results indicate that temperature inversion sensitively depends on the mixed layer depth, sea surface wind speed and solar shortwave radiation, etc., and appropriately meteorological and hydrological conditions often lead to the similarly periodical occurrence of this inversion phenomenon.

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

Mixed-layer water oscillations in tropical Pacific for ENSO cycle

The main modes of interannal variabilities of thermocline and sea surface wind stress in the tropical Pacific and their interactions are investigated, which show the following results. (1) The thermocline anomalies in the tropical Pacific have a zonal dipole pattern with 160 W as its axis and a meridional seesaw pattern with 6-8 degrees N as its transverse axis. The meridional oscillation has a phase lag of about 90 to the zonal oscillation, both oscillations get together to form the El Nino/La Nina cycle, which behaves as a mixed layer water oscillates anticlockwise within the tropical Pacific basin between equator and 12 degrees N. (2) There are two main patterns of wind stress anomalies in the tropical Pacific, of which the first component caused by trade wind anomaly is characterized by the zonal wind stress anomalies and its corresponding divergences field in the equatorial Pacific, and the abnormal cross- equatorial flow wind stress and its corresponding divergence field, which has a sign opposite to that of the equatorial region, in the off-equator of the tropical North Pacific, and the second component represents the wind stress anomalies and corresponding divergences caused by the ITCZ anomaly. (3) The trade winds anomaly plays a decisive role in the strength and phase transition of the ENSO cycle, which results in the sea level tilting, provides an initial potential energy to the mixed layer water oscillation, and causes the opposite thermocline displacement between the west side and east side of the equator and also between the equator and 12 degrees N of the North Pacific basin, therefore determines the amplitude and route for ENSO cycle. The ITCZ anomaly has some effects on the phase transition. (4) The thermal anomaly of the tropical western Pacific causes the wind stress anomaly and extends eastward along the equator accompanied with the mixed layer water oscillation in the equatorial Pacific, which causes the trade winds anomaly and produces the anomalous wind stress and the corresponding divergence in favor to conduce the oscillation, which in turn intensifies the oscillation. The coupled system of ocean-atmosphere interactions and the inertia gravity of the mixed layer water oscillation provide together a phase-switching mechanism and interannual memory for the ENSO cycle. In conclusion, the ENSO cycle essentially is an inertial oscillation of the mixed layer water induced by both the trade winds anomaly and the coupled ocean-atmosphere interaction in the tropical Pacific basin between the equator and 12 degrees N. When the force produced by the coupled ocean-atmosphere interaction is larger than or equal to the resistance caused by the mixed layer water oscillation, the oscillation will be stronger or maintain as it is, while when the force is less than the resistance, the oscillation will be weaker, even break.

Submesoscale activity over the shelf of the northern South China Sea in summer: simulation with an embedded model

We applied a primitive equation ocean model to simulate submesoscale activities and processes over the shelf of the northern South China Sea (NSCS) with a one-way nesting technology for downscaling. The temperature and density fields showed that submesoscale activities were ubiquitous in the NSCS shelf. The vertical velocity was considerably enhanced in submesoscale processes and could reach an average of 58 m per day in the subsurface. At this point, the mixed layer depth also was deepened along the front, and the surface kinetic energy also increased with the intense vertical movement induced by submesoscale activity. Thus, submesoscale stirring/mixing is important for tracers, such as temperature, salinity, nutrients, dissolved organic, and inorganic carbon. This result may have implication for climate and biogeochemical investigations.

多因素对海洋上混合层深度影响的数值模拟

海洋上混合层深度（MLD）是刻画海洋上混合层的一个非常重要的物理参量，MLD的准确确定对人们了解上层海洋动力学、热力学等具有重要的科学意义和实用价值。本文概述了影响MLD的几个主要因素，介绍了MLD的定义方法和主要的两类混合层模型。利用2001年月平均资料，分析了中国南海海表风应力、净热通量、蒸发、降雨及蒸降差的时、空分布变化规律。采用国外先进的MITgcm模式及考虑了非局地作用影响的KPP垂直混合模型，通过改变模式的上边界强迫条件，进行了海表风应力、海表净热通量、蒸降差对MLD影响的敏感性分析，结果表明MLD对海表风应力变化最为敏感。 通过网格嵌套方法，采用MITgcm模式及分辨率为1/4°×1/4°的强迫场资料对南海进行数值模拟，分析了海表温度、盐度、流场及混合层深度的时、空分布特征变化规律，探讨了海表风应力、净热通量、蒸降差对海洋上混合层深度的影响。通过模拟结果与实测资料的比较，同时也表明MITgcm模式应用于中国南海区域范围的MLD数值模拟研究是可行的。 基于考虑了地转和均匀背景流场影响的两层流体界面内波频散关系模型，得到一种新的利用SAR遥感图像反演混合层深度的方法。利用该方法对南海北部SAR内波图像进行了MLD的实例反演研究，并和时空同步的CTD资料进行了对比，结果表明考虑地转和均匀背景流场影响的MLD反演算法更为合理，为更准确地利用SAR图像遥感反演混合层深度奠定了基础。

Instability of Two-Layer Rayleigh-Benard Convection with Interfacial Thermocapillary Effect

The linear instability analysis of the Rayleigh-Allarangoni-Benard convection in a two-layer system of silicon oil 10cS and fluorinert FC70 liquids are performed in a larger range of two-layer depth ratios H, from 0.2 to 5.0 for different total depth H less than or equal to 12 mm. Our results are different from the previous study on the Rayleigh-Benard instability and show strong effects of thermocapillary force at the interface on the time-dependent oscillations arising from the onset of instability convection.

Oscillatory Instability of Rayleigh-Marangoni-Benard Convection in Two-Layer System

The convective instabilities in two or more superposed layers heated from below were studied extensively by many scientists due to several interfacial phenomena in nature and crystal growth application. Most works of them were performed mainly on the instability behaviors induced only by buoyancy force, especially on the oscillatory behavior at onset of convection (see Gershuni et. Al.(1982), Renardy et. Al. (1985,2000), Rasenat et. Al. (1989), and Colinet et. Al.(1994)) . But the unstable situations of multi-layer liquid convection will become more complicated and interesting while considering at the same time the buoyancy effect combined with thermocapillary effect. This is the case in the gravity reduced field or thin liquid layer where the thermocapillary effect is as important as buoyancy effect. The objective of this study was to investigate theoretically the interaction between Rayleigh-Bénard instability and pure Marangoni instability in a two-layer system, and more attention focus on the oscillatory instability both at the onset of convection and with increasing supercriticality. Oscillatory behavious of Rayleigh-Marangoni-Bénard convective instability (R-M-B instability) and flow patterns are presented in the two-layer system of Silicon Oil (10cSt) over Fluorinert (FC70) for a larger various range of two-layer depth ratios (Hr=Hupper/Hdown) from 0.2 to 5.0. Both linear instability analysis and 2D numerical simulation (A=L/H=10) show that the instability of the system depends strongly on the depth ratio of two-layer liquids. The oscillatory instability regime at the onset of R-M-B convection are found theoretically in different regions of layer thickness ratio for different two-layer depth H=12,6,4,3mm. The neutral stability curve of the system displaces to right while we consider the Marangoni effect at the interface in comparison with the Rayleigh-Bénard instability of the system without the Marangoni effect (Ma=0). The numerical results show different regimes of the developing of convection in the two-layer system for different thickness ratios and some differences at the onset of pure Marangoni convection and the onset of Rayleigh-Bénard convections in two-layer liquids. Both traveling wave and standing wave were detected in the oscillatory instability regime due to the competition between Rayleigh-Bénard instability and Marangoni effect. The mechanism of the standing wave formation in the system is presented numerically in this paper. The oscillating standing wave results in the competition of the intermediate Marangoni cell and the Rayleigh convective rolls. In the two-layer system of 47v2 silicone oil over water, a transition form the steady instability to the oscillatory instability of the Rayleigh-Marangoni-Bénard Convection was found numerically above the onset of convection for ε=0.9 and Hr=0.5. We propose that this oscillatory mechanism is possible to explain the experimental observation of Degen et. Al.(1998). Experimental work in comparison with our theoretical findings on the two-layer Rayleigh-Marangoni-Bénard convection with thinner depth for H<6mm will be carried out in the near future, and more attention will be paid to new oscillatory instability regimes possible in the influence of thermocapillary effects on the competition of two-layer liquids

Improved performance in organic light emitting diodes with a mixed electron donor-acceptor film involved in hole injection

Organic light emitting diodes using a mixed layer of electron acceptor 3, 4, 9, 10 perylenetetracarboxylic dianhydride and electron donor copper phthalocyanine (PTCDA:CuPc) on indium tin oxide (ITO) anodes were fabricated. The device properties were found to be strongly dependent on the thickness of the PTCDA:CuPc film: both the power efficiency and the driving voltage of the device were optimized with a thickness of PTCDA:CuPc ranging from 10 to 20 nm. As compared to the conventional ITO/CuPc hole injection structure, the ITO/PTCDA:CuPc hole injection structure could remarkably enhance both the luminance and the power efficiencies of devices. A mechanism of static-induced, very efficient hole-electron pairs generation in mixed PTCDA:CuPc films was proposed to explain the experimental phenomena. The structural and optical properties of PTCDA:CuPc film were examined as well. (c) 2007 American Institute of Physics.

Numerical calculation of the response of the south China sea to typhoon imbudo

The response of the South China Sea (SCS) to Typhoon Imbudo was examined using POM model. The results indicated that SST decreased by 2-6 degrees C with a rightward-biased response as Typhoon Imbudo passed across the SCS. Due to a strong mixing process, the mixed layer (ML) depth deepened as much as 10-60 m and ML heat budget lost 824.78 W/m(2), which was OF dominated by the vertical mixing. By the response of upper ML heat transport, the temperature below the ML increased and oscillated near the inertial period. Furthermore, strong inertial currents were generated by the storm with the max currents up to 1.4 m/s in the upper ML.

The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42,134-149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313-351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313-351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971-5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift. (C) 2008 Elsevier Ltd. All rights reserved.