Simulation of Atmospheric Boundary Layer Characteristics during Indian Summer Monsoon using Observations from Monsoon Trough Boundary Layer Experiment at Jodhpur, India

The diurnal structure of the boundary layer during Indian summer monsoon period is studied using a one-dimensional meteorological boundary layer model and the observations collected from the Monsoon Trough Boundary Layer Experiment conducted in 1990 at Jodhpur, India. The model was initialized with the observed temperature profiles at 0530 LST on 17 July, 1990 at Jodhpur and was run for 26 hours. The study is carried out with a geostrophic wind speed of 9.5 m s⁻¹ corresponding to the strong wind simulation. The mean thermodynamic and wind structure simulated by the model are in good agreement with those observed from 30 m tower. The computed surface layer characteristics such as the surface fluxes, TKE and standard deviations of velocity components are found to be reasonably in good agreement with those based on turbulence measurements. The shear and buoyancy budget computed from the model are also compared with the turbulence measurements. The integrated cooling budget in the nocturnal boundary layer is examined.     

Spatial transport and spectral transfer of solar wind turbulence composed of Alfvén waves and convective structures II: numerical results

This work follows the paper titled “Spatial transport and spectral transfer of solar wind turbulence composed of Alfvén waves and convective structures I: The theoretical model”, and deals with the detailed physics and numerical solution of a two-component solar wind model, consisting of small-scale Alfvén waves and convected structures. In particular, we present numerical results which qualitatively reflect many of the observed features of the radial and spectral evolution of the turbulent energies, the residual energy, the cross-helicity and Alfvén-ratio in high-speed solar wind streams. These features are the following: the formation of a characteristic “inclined eye”, which evolves between the energy spectra displayed over the frequency axis and tends to close in the radial development of the spectra, a steepening of all spectra towards Kolmogorov-like f^-5/3 spectra, the development of the normalized cross-helicity towards a constant not much less than one and the formation of a “trough” form of the Alfvén ratio with a z-shaped left boundary, By weighting special terms in the equations differently, we can also cast light on the physical role of parametric conversion model terms, wave-structure scattering model terms, nonlinear terms, spherical expansion terms and their effects on the radial evolution of turbulent energies in high-speed solar wind streams.     

Early spring turbulent mixing in an ice-covered Arctic fjord during transition to melting

Observations are presented of currents, hydrography and turbulence in a jet-type tidally forced fjord in Svalbard. The fjord was ice covered at the time of the experiment in early spring 2004. Turbulence measurements were conducted by both moored instruments within the uppermost 5 m below the ice and a microstructure profiler covering 3–60 m at 75 m depth. Tidal choking at the mouth of the fjord induces a tidal jet advecting relatively warmer water past the measurement site and dominating the variability in hydrography. While there was no strong correlation with the observed hydrography or mixing and the phase of the semidiurnal tidal cycle, the mean structure in dissipation of turbulent kinetic energy, work done under the ice and the mixing in the water column correlated with the current when conditionally sampled for tidal jet events. Observed levels of dissipation of turbulent kinetic energy per unit mass, 1.1×10⁻⁷ W kg⁻¹, and eddy diffusivity, 7.3×10⁻⁴ m² s⁻¹, were comparable to direct measurements at other coastal sites and shelves with rough topography and strong forcing. During spring tides, an average upward heat flux of 5 W m⁻² in the under-ice boundary layer was observed. Instantaneous (1 h averaged) large heat flux events were correlated with periods of large inflow, hence elevated heat fluxes were associated with the tidal jet and its heat content. Vertical heat fluxes are derived from shear-probe measurements by employing a novel model for eddy diffusivity [Shih et al., 2005. Parameterization of turbulent fluxes and scales using homogeneous sheared stably stratified turbulence simulations. Journal of Fluid Mechanics 525, 193–214]. When compared to the direct heat flux measurements using the eddy correlation method at 5 m below the ice, the upper 4–6 m averaged heat flux estimates from the microstructure profiler agreed with the direct measurements to within 10%. During the experiment water column was stably, but weakly, stratified. Destabilizing buoyancy fluxes recorded close to the ice were absent at 5 m below the ice, and overall, turbulence production was dominated by shear. A scaling for dissipation employing production by both stress and buoyancy [Lombardo and Gregg, 1989. Similarity scaling of viscous and thermal dissipation in a convecting boundary layer. Journal of Geophysical Research 94, 6273–6284] was found to be appropriate for the under-ice boundary layer.     

Analysis of Weak Wind Stable Conditions from the Observations of the Land Surface Processes Experiment at Anand in India

The observations in weak wind stable conditions are scarce. The present study examines the observations from the Land Surface Processes Experiment (LASPEX) conducted at Anand, (Gujarat, India) during the year 1997–1998 to study the characteristics of surface layer under weak wind stable conditions. The observed surface fluxes are compared with those computed using Monin-Obukhov (M-O) similarity theory. The upper air observations and regional climatology are used to justify the persistence of weak wind conditions at Anand. The frequency of occurrence of weak wind stable conditions is observed to be around 67%. In 86% of the cases under weak wind conditions, bulk Richardson number (R_iB) is found to be larger than 0.2. The magnitude of surface fluxes computed from M-O similarity theory is shown to be smaller in comparison to those based on the observations in weak wind stable conditions. Surface fluxes computed using the empirical relations for the eddy diffusivities and drag and heat exchange coefficients are found to be comparable with those based on M-O similarity theory however these fluxes are under-predicted in comparison to the observations. The traditional M-O similarity theory is not able to simulate the observed fluxes well in weak wind stable conditions at Anand.     

南极中山站至Dome A考察断面湍流参数特征近地层湍流参数特征

利用南极中山站至Dome A考察断面上3个自动气象站2005～2007年的观测资料和2008年夏季在中山站附近冰盖获取的湍流观测资料,应用空气动力学方法和涡动相关法计算分析了中山站至Dome A断面上近地层各种湍流参数(感热通量,潜热通量,湍流温度、湿度和速度尺度,地表粗糙度,大气稳定度及动量输送系数)的季节变化、日变...     

近地层能量闭合度对陆面过程模式影响

大量近地层观测试验表明,利用涡动相关法观测的湍流通量小于近地层可利用能量,即近地层能量是不闭合的,这种不闭合度一般为20%甚至更高.而陆面过程模式是基于地气间能量平衡建立,并且模式中的湍流边界层参数化方案通常根据实际观测的湍流通量来确定,因此能量不闭合必将对陆面过程模式造成一定的影响.本文利用2007年春季SACOL站的近地层观测资料,依据能量守恒将能量不闭合中的残余能量通过波文比分配到观测的湍流通量中,即修正涡动相关法观测的湍流通量使得近地层能量达到平衡;之后分别利用观测和修正的湍流通量,建立了能量不闭合和闭合情形下的湍流参数化方案,借助陆面过程模式SHAW,通过数值模拟和对比分析方法考察近地层能量闭合度对陆面过程模式的影响.研究结果表明近地层能量闭合对陆面过程模式有显著的影响:在陆面过程数值模拟中,当应用近地层能量不闭合的湍流通量形成的湍流参数化方案时,陆面过程模式会明显高估地表长波辐射及土壤温度;但当应用修正湍流通量使得近地层能量达到闭合形成的湍流参数化方案后,在不改变任何地表土壤物理生化属性的情况下,陆面过程模式能较好地模拟地表长波辐射和土壤温度.     

Analysis of the effect of the coastal discontinuity on near-surface flow

Conditional sampling is used herein to examine the effect of fetch, stability, and surface roughness changes on wind speeds in the coastal zone. Using data from an offshore wind farm it is shown that at a distance of 1.2–1.7 km from the coast, up to a height of 20 m above the surface, differences in wind speed distributions from onshore and offshore masts are statistically significant for flow moving offshore under all stability conditions. In contrast, differences between the distribution of wind speeds at 38 and 48 m at masts located at the coast and in the coastal zone are not significant for flow moving offshore, indicating that flow at these heights is not fully adjusted to the change in surface roughness (land to sea). These findings are in accordance with calculations of the internal boundary layer (IBL) height which indicate that the IBL would frequently be below the two upper measurement heights at 1.2–1.6 km from the coast. The analyses presented here indicate that the wind speed distribution at a potential offshore wind farm site is not solely dependent on fetch (distance from the coast) but also depends on the stability climate.     

TOGA-COARE IOP期间的海气通量观测结果

在TOGA－COAREIOP期间用涡度相关法对海气热通量进行了船载直接观测．对垂直风速、温度和湿度湍流脉动观测数据的谱分析显示它们在高频区基本满足"－2／3次方律"．对船体简谐震荡影响的讨论从理论上证明该影响在热通量计算中可被忽略．根据以此方法得到的通量求出了中性层结条件下感热和潜热的整体输送系数分别为2．25×10^－3和1．26×10^－3．对海气边界层特性的分析表明该海域的近海层主要呈中性或弱不稳定层结．海气通量的变化与背景环流形势密切相关,潜热通量主要受海面风场强度的影响,而感热通量变化除了风场的影响外,层结变化也是一个重要因素．用整体输送法计算TOGA－COAREIOP期间以及TOGA期间8个航次的通量结果而得到的Bowen比约为0．1,显示潜热通量是暖池大气的主要热源．     

Turbulence measurements in a submarine canyon

Profiles of velocity turbulence in Monterey Canyon, made with a recently developed expendable probe, show the existence of a very turbulent bottom boundary layer. The turbulent flow is up to 170 m thick and has peak microscale shears of 1 m s⁻¹ per meter. The rate of dissipation of kinetic energy, based on the observed shear variance, averaged over the depth of the turbulent boundary layer ranged from 70 to 500 × 10⁻⁶W m⁻³. Temperature measurements indicate that the flow was up canyon at a time of low tide. The upper bound for the vertical eddy viscosity is estimated to be17 × 10⁻⁴m²s⁻¹ and for the vertical eddy diffusivity is estimated to be 15 × 10⁻⁴m²s⁻¹. The large vertical scale and the intensity of the observed boundary layer suggest that the flow in Monterey Canyon may be important for the renewal and circulation of water over the continental shelf in the bay area.     

EISCAT observations of unusual flows in the morning sector associated with weak substorm activity

A discussion is given of plasma flows in the dawn and nightside high-latitude ionospheric regions during substorms occurring on a contracted auroral oval, as observed using the EISCAT CP-4-A experiment. Supporting data from the PACE radar, Greenland magnetometer chain, SAMNET magnetometers and geostationary satellites are compared to the EISCAT observations. On 4 October 1989 a weak substorm with initial expansion phase onset signatures at 0030 UT, resulted in the convection reversal boundary observed by EISCAT (at \sim0415 MLT) contracting rapidly poleward, causing a band of elevated ionospheric ion temperatures and a localised plasma density depletion. This polar cap contraction event is shown to be associated with various substorm signatures; Pi2 pulsations at mid-latitudes, magnetic bays in the midnight sector and particle injections at geosynchronous orbit. A similar event was observed on the following day around 0230 UT (\sim0515 MLT) with the unusual and significant difference that two convection reversals were observed, both contracting poleward. We show that this feature is not an ionospheric signature of two active reconnection neutral lines as predicted by the near-Earth neutral model before the plasmoid is “pinched off”, and present two alternative explanations in terms of (1) viscous and lobe circulation cells and (2) polar cap contraction during northward IMF. The voltage associated with the anti-sunward flow between the reversals reaches a maximum of 13 kV during the substorm expansion phase. This suggests it to be associated with the polar cap contraction and caused by the reconnection of open flux in the geomagnetic tail which has mimicked “viscous-like” momentum transfer across the magnetopause.     

Stratosphere-troposphere exchange in an extratropical cyclone, calculated with a Lagrangian method

A Lagrangian technique is developed and applied to calculate stratosphere-troposphere exchange in an extratropical cyclone. This exchange is computed from the potential vorticity or PV along trajectories, calculated from ECMWF circulation data. Special emphasis is put on the statistical significance of the results. The computed field of the cross-tropopause flux is dominated by elongated patterns of statistically significant large downward and small upward fluxes. The downward fluxes mainly occur in the lower part of the considered tropopause folds. The upward fluxes are found near the entrance of the folds, in the tropopause ridges. The ratio between the area averaged downward and upward cross-tropopause fluxes increases with increasing strength of the cyclone. Since the largest fluxes are shown to occur in the regions with the largest wind shear, where PV-mixing is thought to cause large cross-tropopause fluxes, the results are expected to be reliable, at least in a qualitative sense. The position of a tropopause fold along the northwest coast of Africa is confirmed by total ozone observations. The results indicate that the applied Lagrangian technique is an appropriate tool for diagnosing stratosphere-troposphere exchange.     

The effect of flow at Maud Rise on the sea-ice cover – numerical experiments

 The role of seamounts in the formation and evolution of sea ice is investigated in a series of numerical experiments with a coupled sea ice–ocean model. Bottom topography, stratification and forcing are configured for the Maud Rise region in the Weddell Sea. The specific flow regime that develops at the seamount as the combined response to steady and tidal forcing consists of free and trapped waves and a vortex cap, which is caused by mean flow and tidal flow rectification. The enhanced variability through tidal motion in particular modifies the mixed layer above the seamount enough to delay and reduce sea-ice formation throughout the winter. The induced sea-ice anomaly spreads and moves westward and affects an area of several 100 000 km². Process studies reveal the complex interaction between wind, steady and periodic ocean currents: all three are required in the process of generation of the sea ice and mixed layer anomalies (mainly through tidal flow), their detachment from the topography (caused by steady oceanic flow) and the westward translation of the sea-ice anomaly (driven by the time-mean wind).     

The response of ionospheric convection in the polar cap to substorm activity

We report multi-instrument observations during an isolated substorm on 17 October 1989. The EISCAT radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71°-78°. SAMNET and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. IMP-8 magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux to the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71° by the time of the expansion phase onset. A westward electrojet, centred at 65.4°, occurred at the onset of the expansion phase. This electrojet subsequently moved poleward to a maximum of 68.1° at 2000 UT and also widened. During the expansion phase, there is evidence of bursts of plasma flow which are spatially localised at longitudes within the substorm current wedge and which occurred well poleward of the westward electrojet. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the “distant” neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase. It is not obvious whether the electrojet mapped to a near-Earth neutral line, but at its most poleward, the expanded electrojet does not reach the estimated latitude of the polar cap boundary.     

Comparison of a simple planetary boundary layer model with measurements of a turbulent boundary layer under pack ice

A simple analytical theory which describes the motion in a turbulent planetary boundary layer near a rough sea bed by using a two-layer eddy viscosity model is presented. The vertical structure of the current in the boundary layer is presented, and comparisons are made with data fromMcPhee and Smith (1976, Journal of Physical Oceanography,6, 696–711) obtained from measurements of the turbulent boundary layer under drifting Arctic ice.     

Stable Atmospheric Surface Layer in Modelling with Constant Ground Temperature

Influence of the surface roughness, geostrophic wind speed and initial stable stratification on surface layer scales, namely the friction velocity, temperature scale and the Monin-Obukhov length, and on non-dimensional profiles has been studied in quasi-stationary PIAPBLM (Prague Institute of Atmospheric Physics Boundary Layer Model) runs. Modifications of the stability functions, of the non-local mixing length concept, and of surface layer parametrization have improved the model results. With the ground temperature kept constant in time, the boundary layer is subjected to a turbulent cooling induced by surface roughness and to a counteracting turbulent heat transfer directed downwards. If the wind speed is lower, a rather mixed layer develops capped by a more stable layer. The estimated Monin-Obukhov length then slightly increases and the non-dimensional gradients overpredict the empirical values even more.     

Altitude integration effects in the skewness of type-2 coherent echoes

Possible effects of signal reception from different electrojet heights in a skewness of auroral coherent spectra are studied assuming that the “inherent” spectral line due to plasma turbulence is of type-2 and symmetrical. For reasonable ionospheric parameters, the altitude integrated spectra are expected to be skewed negatively for positive mean Doppler shift, in agreement with radar observations at small aspect angles. However, the spectra could be skewed positively if the turbulent layer or the electron density profile is shifted to high altitudes of \sim120 km. This change of spectral shape will not be observed experimentally if, at the same time, either the electron collision frequency is enhanced or the “inherent” spectral width is increased. Observational results are discussed in view of the predictions given.     

边界层参数化方案对边界层热力和动力结构特征影响的比较

本文利用高分辨率中尺度WRF模式,通过改变边界层参数化方案进行多组试验,评估该模式对美国北部森林地区边界层结构的模拟能力,同时比较了五种不同边界层参数化方案模拟得出的边界层热力和动力结构.结果表明：除个别方案外,配合不同边界层方案的WRF模式都能成功模拟出白天对流边界层强湍流混合特征和夜间稳定边界层内强逆温、逆湿和低空急流等热力和动力结构.非局地YSU、ACM2方案在白天表现出强的湍流混合和卷夹,相比于局地MYJ、UW方案,模拟的对流边界层温度更高、湿度更低、混合层高度更高、感热通量更大,更接近实际观测,这表明在不稳定层结下考虑非局地大涡输送更为合理,但局地方案在风速和风向的预报上存在一定优势.TEMF方案得到的白天局地湍流混合强度为所有方案中最弱,混合层难以发展,无法体现对流边界层内气象要素垂直分布均匀的特点.对于夜间稳定边界层的模拟,不同参数化方案之间的差异较小,但是YSU方案在一定程度上高估了机械湍流,导致局地湍流混合偏强,从而影响了其对稳定边界层的模拟能力.     

Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.