Grain-size spectra maps of beaches in W. Peloponnese, Greece

Grain-size spectra maps of sands from beaches of the western Peloponnese are related to statistical parameters and can be used to infer characteristics of the sedimentation process involved which were unresolved by previous statistical analysis. On the basis of similarities in spectral patterns, the studied beaches fall into four groups: A, B, C, and D. The periodicity observed on the sinuous pattern of the spectra map of groups A and B is interpreted as indicating edge-wave activity. Regarding group D, the gradual W—E increase of particle size along the coast suggests a W—E increase of wave energy, which is also reflected by the inclined trend of the spectra map. Lack of periodicity on the contour pattern of group C probably reflects intensive sample-to-sample sorting variation.     

A Re-Evaluation of Long-Term Flux Measurement Techniques Part II: Coordinate Systems

To convert measurements of windspeed, eddy flux and scalar concentration into estimates of surface scalar exchange, we implicitly or explicitly assimilate the measurements into mathematical statements of the mass balance in a control volume on a representative patch of the surface. The form of this statement depends on the coordinate system in which it is written and the coordinate system should be chosen so that measurements can be used optimally. This requirement imposes a set of conditions on the coordinates. Here we perform a comparative analysis of some candidate coordinate systems, concentrating on the Cartesian and physical streamline systems. We show that over gentle topography there are definite advantages in working in streamline coordinates. Transforming measurements of vector and tensor quantities measured in the reference frame, s _i, of the anemometer into the reference frame, e _i, of the chosen coordinate system involves using the measured statistics of the wind field to define three Euler rotation angles. We compare the method in most common use, which employs the components of the mean wind vector and the Reynolds stress tensor to define these angles, with the more recent planar-fit method that uses instead an ensemble of mean wind vectors to define the rotations. We find that, in real flows, the standard method has a previously unrecognized closure problem that ensures that the third rotation angle defined using the stress tensor or scalar flux vector will always be in error and often give unphysical results. An alternative procedure is recommended. Finally, the relationships between measurements and model outputs are discussed.     

The Neutral, Barotropic Planetary Boundary Layer, Capped by a Low-Level Inversion

The presence of a low-level, capping inversion layer will affect the height and structure of the planetary boundary layer (PBL). Results from models of varying levels of sophistication, including analytical, turbulent kinetic energy (TKE), second-order closure (SOC), large-eddy simulation (LES) and direct numerical simulation (DNS) models, are used to investigate this influence for the neutral, barotropic PBL. Predicted and observed profiles of stress and geostrophic departure components, and integral measures, such as the parameters of Rossby-number similarity theory, are compared for the KONTUR, Marine Stratocumulus, JASIN, Leipzig, Pre-Wangara and Upavon field experiments.Analytical models of the equilibrium value of inversion height z_i, which depend on the surface friction velocity u_*, and both the Coriolis parameter f and the free-flow Brunt-Väisälä frequency N, are found to give reasonable estimates of the PBL height. They also indicate that only the KONTUR and Marine Stratocumulus experiments were strongly influenced by N. More quantitative comparisons would require larger, more comprehensive datasets. The effects of the presence of a capping inversion on the profile structure were found to be insignificant for h_* = |f|z_i/u_* > 0.15.The simple analytical model performed quite well over all values of h_*; it predicted the profiles of the longitudinal stress component (in the direction of the surface stress) better than the lateral component. The more advanced models performed well for small values of h_* (for flow over the sea), but systematically underestimated the cross-isobaric angle for flow over land. These models predicted the profiles of the lateral stress component better than the longitudinal component. The profiles of the analytical model agreed with those of the advanced models when the constant eddy viscosity of the outer layer was increased.Agreement with DNS was achieved by increasing the eddyviscosity of the analytical model by a factor of 5.Zilitinkevich and Esau(2002, Boundary-Layer Meteorology, 104, 371–379)suggest that the neutral, barotropic values of A and B of Rossby-numbersimilarity theory are not universal constants, but depend on the ratio N/|f|. The dependence for A and B is calculated using the analytical model and TKE models. Over the sea (h_* 0.1; N/|f| 100, where we have used the Zilitinkevich-Esau relation to convert between h_* and N/|f|) there is agreement between the model predictions and observations; however over land where the equilibrium boundary-layer height is greater (h_* 0.35; N/|f| 10) the inconsistency between the advanced model predictions (TKE, SOC, LES, and DNS) and observations, as noted previously by Hess and Garratt, still exists. We attribute this disagreement to violations of the strict assumptions of steady, horizontally homogeneous, neutral, barotropic conditions implicit in the observations. At small values of z_i and a strongly stable background stratification (h_* 0.04; N/|f| 1000) both the TKE and analytical models predict that A and B depend significantly on h_*, however observations are unavailable to confirm these predictions. Zilitinkevich and Esau call this case the `long-lived near-neutral PBL', and state that it is found in cold weather at high latitudes.     

One-Equation Turbulence Modelling for Atmospheric and Engineering Applications

A new algebraic turbulent length scale model is developed, based on previous one-equation turbulence modelling experience in atmospheric flow and dispersion calculations. The model is applied to the neutral Ekman layer, as well as to fully-developed pipe and channel flows. For the pipe and channel flows examined the present model results can be considered as nearly equivalent to the results obtained using the standard k– model. For the neutral Ekman layer, the model predicts satisfactorily the near-neutral Cabauw friction velocities and a dependence of the drag coefficient versus Rossby number very close to that derived from published (G. N. Coleman) direct numerical simulations. The model underestimates the Cabauw cross-isobaric angles, but to a less degree than the cross-isobar angle versus Rossby dependence derived from the Coleman simulation. Finally, for the Cabauw data, with a geostrophic wind magnitude of 10 ms^–1, the model predicts an eddy diffusivity distribution in good agreement with semi-empirical distributions used in current operational practice.     

Wind tunnel simulation experiment and investigation on the electrification of sandstorms

~~Wind tunnel simulation experiment and investigation on the electrification of sandstorms~~     

Composition and assemblage characteristics of magnetic minerals in layer S<Subscript>5-1</Subscript> of Xifeng and Duanjiapo sections

Relatively strongly magnetic fine components (< 30μm, XS-4J and DS-4J) which are most environmentally sensitive were separated from layer S_5-1 in the Xifeng and Duanjiapo loess sections and analyzed by MPV-3 for their morphometric characteristics and reflectance, SEM-ESD for their element contents and XRD for their mineral phases, respectively. The results showed that minerals in both samples are dominated by detrial Fe-Ti oxides of aeolian origin. In sample XS-4J the reflectance and iron contents of magnetic minerals are usually high. In addition to magnetite (Fe₃O₄), maghemite (γFe₂O₃) and hematite (Fe₂O₃), some Fe-high oxide (72.25 wt%–86.67 wt%), ilmenite (FeTiO₃), and magnetite-ulvöspinel [Fe(FeCr)O₄, Fe (FeNi)O₄] were also detected. In sample DS-4J obvious negative linear correlations were found between Ti and Fe, and the contents of Mn, Si, Al and Ca are usually high and the minerals are dominated by magnetite (maghemite), goethite (FeOOH) and limonite (containing Si and OH). In addition, the signs of corrosion of magnetic minerals and newly crystallized magnetite (maghemite) were recognized. Differences in the composition and assemblage characteristics of magnetite minerals between XS and DS reflect significant differences in source rocks and preserving conditions.     

黄河三角洲海岸带陆海相互作用概念模式

根据陆海相互作用的观点,黄河流域、黄河干流、河口三角洲及其邻近海区的生态环境相互联系,组成了一个有机的生态系统链,可称为黄河-渤海生态系统。黄河流域的降水量、土壤植被条件使黄河干流具有水少沙多、水沙异源和水资源缺乏的特征。干流入海水、沙通量变化影响了黄河三角洲地区的侵蚀、堆积和发育过程。黄河物质入海后在河口及邻近海域形成了具有高生产力的生态环境和著名的渔场。据此绘出了黄河-渤海生态环境内各环节相互联系的概念模式图,并提出了黄河流域大面积水土保持和南水北调工程逐步展开等新环境下,本区陆海相互作用的研究方向。     

Surface gravitational field and topography changes induced by the Earth’s fluid core motions

Using a Love number formalism, the elastic deformations of the mantle and the mass redistribution gravitational potential within the Earth induced by the fluid pressure acting at the core–mantle boundary (CMB) are computed. This pressure field changes at a decadal time scale and may be estimated from observations of the surface magnetic field and its secular variation. First, using a spherical harmonic expansion, the poloidal and toroidal part of the fluid velocity field at the CMB for the last 40 years is computed, under the hypothesis of tangential geostrophy. Then the associated geostrophic pressure, whose order of magnitude is about 1000 Pa, is computed. The surface topography induced by this pressure field is computed using Love numbers, and is a few millimetres. The mass redistribution gravitational potential induced by these deformations and, in particular, the zonal components of the related surface gravitational potential perturbation (J₂, J₃ and J₄ coefficients), are calculated. Overall perturbations for the J₂ coefficient of about 10^–10, for J₃ of about 10^–11 and for J₄ are found of about 0.3×10^–11. Finally, these theoretical results are compared with recent observations of the decadal variation of J₂ from satellite laser ranging. Results concerning J₂ can be described as follows: first, they are one order of magnitude too small to explain the observed decadal variation of J₂ and, second, they show a significant linear trend over the last 40 years, whose rate of decrease amounts to 7% of the observed value.     

On the sensitivity of computed surges to open-boundary formulation

A coarse-grid (resolution of order 7 km) model of the west coast of Britain is used to examine the sensitivity of computed storm-surge elevations and currents to a range of open-boundary conditions. The storm-surge period 1 to 26 March 1994 is used for this comparison, as it is a time of significant wind activity. Also current measurements in the North Channel of the Irish Sea together with coastal elevation measurements are available for model validation. Elevations and currents previously computed with a coarse-grid shelf-wide model can also be incorporated into the open-boundary condition to examine the influence of far-field effects. Initial model calculations with no far-field input show the importance of including shelf-wide effects from either the external shelf model, or by using observations from coastal gauges interpolated along the open boundary of the west-coast model. Provided the west-coast models open boundary is taken sufficiently far away from the region of interest, in this case the Irish Sea, then either a radiation condition or an elevation-specified condition is appropriate provided far-field effects are taken into account. If these are not included, then neither boundary condition is successful. For the radiation condition it is necessary to include both elevations and currents from a far-field model in order to reproduce the surge. In the case of an elevation-specified boundary condition far-field effects can be incorporated in hindcast calculations by including observed sea-level changes. In a storm-surge prediction calculation the radiation condition with a far-field model is required. Calculations show that computed elevations are spatially more coherent than currents, with flows through the western Irish Sea showing the greatest sensitivity to open-boundary formulation during storm events.Responsible Editor: Phil Dyke     

Tidal currents,energy flux and bottom boundary layer thickness in the Clyde Sea and North Channel of the Irish Sea

A high-resolution three-dimensional model of the Clyde Sea and the adjacent North Channel of the Irish Sea is used to compute the major diurnal and semidiurnal tides in the region, the associated energy fluxes and thickness of the bottom boundary layer. Initially, the accuracy of the model is assessed by performing a detailed comparison of computed tidal elevations and currents in the region, against an extensive database that exists for the M₂, S₂, N₂, K₁ and O₁ tides. Subsequently, the model is used to compute the tidal energy flux vectors in the region. These show that the major energy flux is confined to the North Channel region, with little energy flux into the Clyde Sea. Comparison with the observed energy flux in the North Channel shows that its across-channel distribution and its magnitude are particularly sensitive to the phase difference between elevation and current. Consequently, small changes in the computed values of these parameters due to slight changes of the order of the uncertainty in the open-boundary values to the model, can significantly influence the computed energy flux. The thickness of the bottom boundary layer in the region is computed using a number of formulations. Depending upon the definition adopted, the empirical coefficient C used to determine its thickness varies over the range 0.1 to 0.3, in good agreement with values found in the literature. In the North Channel, the boundary layer thickness occupies the whole water depth, and hence tidal turbulence produced at the sea bed keeps the region well mixed. In the Clyde Sea, the boundary layer thickness is a small fraction of the depth, and hence the region stratifies.Responsible Editor: Phil Dyke