An Intercomparison of Large-Eddy Simulations of the Stable Boundary Layer

Results are presented from the first intercomparison of large-eddy simulation (LES) models for the stable boundary layer (SBL), as part of the Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study initiative. A moderately stable case is used, based on Arctic observations. All models produce successful simulations, in as much as they generate resolved turbulence and reflect many of the results from local scaling theory and observations. Simulations performed at 1-m and 2-m resolution show only small changes in the mean profiles compared to coarser resolutions. Also, sensitivity to subgrid models for individual models highlights their importance in SBL simulation at moderate resolution (6.25 m). Stability functions are derived from the LES using typical mixing lengths used in numerical weather prediction (NWP) and climate models. The functions have smaller values than those used in NWP. There is also support for the use of K-profile similarity in parametrizations. Thus, the results provide improved understanding and motivate future developments of the parametrization of the SBL.     

Influence of the Coriolis force on the formation of a seasonal thermocline

Large eddy simulation (LES) reveals that the Coriolis force plays an important role in seasonal thermocline formation. In the high-latitude ocean, a seasonal thermocline is formed at a certain depth, across which the downward transports of heat and momentum are prohibited. On the other hand, in the equatorial ocean, heat and momentum continue to propagate downward to the deeper ocean without forming a well-defined thermocline. Mechanism to clarify the latitudinal difference is suggested. The depth of a seasonal thermocline h is scaled in terms of both the Ekman length scale λ and the Monin–Obukhov length scale L, as h ??? 0.5(Lλ)^1/2, which is in contrast to the earlier suggestion as h?∝?L.     

Turbulent mixing in a rotating,stratified fluid

Abstract

An experimental study was carried out to investigate the effect of rotation on turbulent mixing in a stratified fluid when the turbulence in the mixed layer is generated by an oscillating grid. Two types of experiments were carried out: one of them is concerned with the deepening of the upper mixed layer in a stable, two-fluid system, and the other deals with the interaction between a stabilizing buoyancy flux and turbulence.

In the first type of experiments, it was found that rotation suppresses entrainment at larger Rossby numbers. As the Rossby number becomes smaller (Ro 0.1), the entrainment rate increases with rotation—the onset of this phenomenon, however, was found to coincide with the appearance of coherent vortices within the mixed layer. The radiation of energy from the mixed layer to the lower non-turbulent layer was found to occur and the magnitude of the energy flux was found to be increased with the rotational frequency. It was also observed that vortices are generated, rather abruptly, in the lower layer as the mixed layer deepens.

In the second set of experiments a quasi-steady mixed layer was found to develop of which the thickness varies with rotation in a fashion that is consistent with the result of the first experiment. Also the rotation was found to delay the formation of a pycnocline.     

Dissolved organic phosphorus utilization and alkaline phosphatase activity of the dinoflagellate <Emphasis Type="Italic">Gymnodinium impudicum</Emphasis> isolated from the South Sea of Korea

This study investigated alkaline phosphatase (APase) activity and dissolved organic and inorganic phosphorus utilization by the harmful dinoflagellate Gymnodinium impudicum (Fraga et Bravo) Hansen et Moestrup isolated from the South Sea of Korea. Under conditions of limited phosphorus, observation of growth kinetics in batch culture yielded a maximum growth rate (μ_max) of 0.41 /day and a half saturation constant (Ks) of 0.71 μM. In time-course experiments, APase was induced as dissolved inorganic phosphorus (DIP) concentrations fell below 0.83 μM, a threshold near the estimated Ks; APase activity increased with further DIP depletion to a maximum of 0.70 pmol/cell/h in the senescent phase. Thus, Ks may be an important index of the threshold DIP concentration for APase induction. G. impudicum utilizes a wide variety of dissolved organic phosphorus compounds in addition to DIP. These results suggest that DIP limitation in the Southern Sea of Korea may have led to the spread of G. impudicum along with the harmful dinoflagellate Cochlodinium polykrikoides in recent years.     

Morphology and Immune-related activities of hemocytes of the mussel Mytilus coruscus (Gould, 1861) from East Sea of Korea

Ocean Science Journal - Mytilus coruscus is a cold water mussel commonly occurring on shallow subtidal rocky substrata along the east coast of Korea. Despite its abundance, there have only been a...     

Changes in mixed layer depth under climate change projections in two CGCMs

Two coupled general circulation models, i.e., the Meteorological Research Institute (MRI) and Geophysical Fluid Dynamics Laboratory (GFDL) models, were chosen to examine changes in mixed layer depth (MLD) in the equatorial tropical Pacific and its relationship with ENSO under climate change projections. The control experiment used pre-industrial greenhouse gas concentrations whereas the 2 × CO₂ experiment used doubled CO₂ levels. In the control experiment, the MLD simulated in the MRI model was shallower than that in the GFDL model. This resulted in the tropical Pacific’s mean sea surface temperature (SST) increasing at different rates under global warming in the two models. The deeper the mean MLD simulated in the control simulation, the lesser the warming rate of the mean SST simulated in the 2 × CO₂ experiment. This demonstrates that the MLD is a key parameter for regulating the response of tropical mean SST to global warming. In particular, in the MRI model, increased stratification associated with global warming amplified wind-driven advection within the mixed layer, leading to greater ENSO variability. On the other hand, in the GFDL model, wind-driven currents were weak, which resulted in mixed-layer dynamics being less sensitive to global warming. The relationship between MLD and ENSO was also examined. Results indicated that the non-linearity between the MLD and ENSO is enhanced from the control run to the 2 × CO₂ run in the MRI model, in contrast, the linear relationship between the MLD index and ENSO is unchanged despite an increase in CO₂ concentrations in the GFDL model.     

Interaction of a negatively buoyant line plume with a density interface

Laboratory experiments were carried out to investigate the interaction between turbulent line buoyant plumes and sharp density interfaces, with the aim of using the results to interpret oceanic observations pertinent to crack openings in the polar ice-cap (leads). These openings take the form of long narrow channels, and are often modeled as line bouyant plumes. The plumes descend as in a homogenoous fluid, impinge on the density interface, and then spread horizontally as gravity currents. Depending on the Richardson number , where Δb is the buoyancy jump across the interface, l_D is the half-width of the plume before the impingement and q₀is the buoyancy flux per unit length of the source, different flow patterns were identified. When Ri < 0.5, the plumes penetrate deep into the bottom layer, deflect horizontally and then spread while showing little vertical rise. When 0.6 < Ri < 5, the penetration is significant, but the fluid bounces back after entraining heavy fluid from the lower layer and then spreads horizontally above the interface as a gravity current. Appreciable mixing between this current and the lower layer was detected when Ri <1. When Ri > 10, the penetration was small and a sharp-nosed gravity current emerged some time after the impact. Measurements were made on the penetration depth, the velocities of the gravity current and the subsurface flow towards the plume, the entrainment rate and other wave parameters. Possible implications of the results for oceanic cases are also discussed.     

The role of molecular diffusion in the deepening of the mixed layer

Turbulent mixing across heat-stratified density interfaces was studied in the laboratory using oscillating-grid generated turbulence. The aim was to study the transition between the entrainment regimes dominated by interfacial wave-breaking and molecular diffusion, and to study the characteristics of the latter. It was observed that, above a critical Richardson number Ri_c, which depends on the Peclet number Pe, the mixing due to wave breaking disappears and that Ri_c Pe⁻ⁿ, where the mean value of the exponent n is approximately . Above Ri_c, the entrainment is molecular-diffusion dominated and takes place through a sequence of events: the buoyancy gradient of the initially sharp density interface is weakened by molecular diffusion until the mixed-layer eddies can engulf a portion of the interfacial layer wherefore the interface sharpens again. Thus, the entrainment events are recurrent with a rate-controlling diffusion stage between them. An entrainment law of the form E Ri⁻²Pe⁻², where E is the entrainment coefficient and Ri is the Richardson number, is suggested for the diffusion-dominated entrainment regime.     

Improvement of the K-profile Model for the Planetary Boundary Layer based on Large Eddy Simulation Data

Modifications of the widely used K-profile model of the planetary boundary layer (PBL), reported by Troen and Mahrt (TM) in 1986, are proposed and their effects examined by comparison with large eddy simulation (LES) data. The modifications involve three parts. First, the heat flux from the entrainment at the inversion layer is incorporated into the heat and momentum profiles, and it is used to predict the growth of the PBL directly. Second, profiles of the velocity scale and the Prandtl number in the PBL are proposed, in contrast to the constant values used in the TM model. Finally, non-local mixing of momentum was included. The results from the new PBL model and the original TM model are compared with LES data. The TM model was found to give too high PBL heights in the PBL with strong shear, and too low heights for the convection-dominated PBL, which causes unrealistic heat flux profiles. The new PBL model improves the predictability of the PBL height and produces profiles that are more realistic. Moreover, the new PBL model produces more realistic profiles of potential temperature and velocity. We also investigated how each of these three modifications affects the results, and found that explicit representation of the entrainment rate is the most critical.     

Retrieval of the single scattering albedo of Asian dust mixed with pollutants using lidar observations

The vertical distribution of single scattering albedos （SSAs） of Asian dust mixed with pollutants was derived using the multi-wavelength Raman lidar observation system at Gwangju （35.10°N,126.53°E）.Vertical profiles of both backscatter and extinction coefficients for dust and non-dust aerosols were extracted from a mixed Asian dust plume using the depolarization ratio from lidar observations.Vertical profiles of backscatter and extinction coefficients of non-dust particles were input into an inversion algorithm to retrieve the SSAs of non-dust aerosols.Atmospheric aerosol layers at different heights had different light-absorbing characteristics.The SSAs of non-dust particles at each height varied with aerosol type,which was either urban/industrial pollutants from China transported over long distances at high altitude,or regional/local pollutants from the Korean peninsula.Taking advantage of independent profiles of SSAs of non-dust particles,vertical profiles of SSAs of Asian dust mixed with pollutants were estimated for the first time,with a new approach suggested in this study using an empirical determination of the SSA of pure dust.The SSAs of the Asian dust-pollutants mixture within the planetary boundary layer （PBL） were in the range 0.88-0.91,while the values above the PBL were in the range 0.76-0.87,with a very low mean value of 0.76 ± 0.05.The total mixed dust plume SSAs in each aerosol layer were integrated over height for comparison with results from the Aerosol Robotics Network （AERONET） measurements.Values of SSA retrieved from lidar observations of 0.92 ± 0.01 were in good agreement with the results from AERONET measurements.