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     

Yaxcopoil-1 and the Chicxulub impact

CSDP core Yaxcopoil-1 was drilled to a depth of 1,511 m within the Chicxulub crater. An organic-rich marly limestone near the base of the hole (1,495 to 1,452 m) was deposited in an open marine shelf environment during the latest Cenomanian (uppermost Rotalipora cushmani zone). The overlying sequence of limestones, dolomites and anhydrites (1,495 to 894 m) indicates deposition in various carbonate platform environments (e.g., sabkhas, lagoons). A 100-m-thick suevite breccia (894–794 m) identifies the Chicxulub impact event. Above the suevite breccia is a dolomitic limestone with planktic foraminiferal assemblages indicative of Plummerita hantkeninoides zone CF1, which spans the last 300 ky of the Maastrichtian. An erosional surface 50 cm above the breccia/dolomite contact marks the K/T boundary and a hiatus. Limestones above this contact contain the first Tertiary planktic foraminifera indicative of an upper P. eugubina zone P1a(2) age. Another hiatus 7 cm upsection separates zone P1a(2) and hemipelagic limestones of planktic foraminiferal Zone P1c. Planktic foraminiferal assemblages of Zone Plc to P3b age are present from a depth of 794.04 up to 775 m. The Cretaceous carbonate sequence appears to be autochthonous, with a stratigraphic sequence comparable to late Cretaceous sediments known from outside the Chicxulub crater in northern and southern Yucatan, including the late Cenomanian organic-rich marly limestone. There is no evidence that these sediments represent crater infill due to megablocks sliding into the crater, such as major disruption of sediments, chaotic changes in lithology, overturned or deep dipping megablocks, major mechanical fragmentation, shock or thermal alteration, or ductile deformation. Breccia units that are intercalated in the carbonate platform sequence are intraformational in origin (e.g., dissolution of evaporites) and dykes are rare. Major disturbances of strata by the impact therefore appear to have been confined to within less than 60 km from the proposed impact center. Yaxcopoil-1 may be located outside the collapsed transient crater cavity, either on the upper end of an elevated and tilted horst of the terrace zone, or even outside the annular crater cavity. The Chicxulub site thus records a large impact that predates the K/T boundary impact and mass extinction.     

Crystallographic orientations of high-angle grain boundaries in dynamically recrystallized quartz: First results

The crystallography and geometry of high-angle grain boundaries from dynamically recrystallized quartz have been studied. On the basis of combined electron backscatter diffraction and universal stage measurements, the complete crystallographic orientation of the grain boundaries could be calculated. The u-stage rotation of the grain boundaries to a vertical position reveals that they are never curved but always consist of straight segments. Our results show that these segments preferentially occupy rhombohedral, trapezohedral and bipyramidal orientations, i.e., orientations in a  25–50° girdle to the c-axis. A specific, albeit low, number of segments with special crystallographic orientation, with respect to a neighbouring quartz grain, often shows another special orientation with respect to the other neighbouring grain. Preferred combinations of grain boundary orientations related to both neighbouring grains are (i) low-index rhombohedral and high index trapezohedral, (ii) low-index bipyramidal and low-index trapezohedral or high-index rhombohedral, and (iii) low-index trapezohedral and low or high index trapezohedral. In certain cases, such as at triple junctions, the boundaries occupy specific trapezohedral orientations with a constant angle to the c-axis. This argues for energy isotropy of trapezohedral planes with the same angle to the quartz c-axis. In general, good match coincidence site lattice (CSL) orientations are not preferentially occupied so that most of the studied grain boundaries represent general boundaries. The formation of straight segments in special crystallographic orientations indicates the crystallographic control and implies an energy reduction of certain general boundaries.     

Instabilities of the tidally induced bottom boundary layer in the rotating frame and their mixing effect

To investigate the stability of the bottom boundary layer induced by tidal flow (oscillating flow) in a rotating frame, numerical experiments have been carried out with a two-dimensional non-hydrostatic model. Under homogeneous conditions three types of instability are found depending on the temporal Rossby number Ro_t, the ratio of the inertial and tidal periods. When Ro_t < 0.9 (subinertial range), the Ekman type I instability occurs because the effect of rotation is dominant though the flow becomes more stable than the steady Ekman flow with increasing Ro_t. When Ro_t > 1.1 (superinertial range), the Stokes layer instability is excited as in the absence of rotation. When 0.9 < Ro_t < 1.1 (near-inertial range), the Ekman type I or type II instability appears as in the steady Ekman layer. Being much thickened (100 m), the boundary layer becomes unstable even if tidal flow is weak (5 cm/s). The large vertical scale enhances the contribution of the Coriolis effect to destabilization, so that the type II instability tends to appear when Ro_t > 1.0. However, when Ro_t < 1.0, the type I instability rather than the type II instability appears because the downward phase change of tidal flow acts to suppress the latter. To evaluate the mixing effect of these instabilities, some experiments have been executed under a weak stratification peculiar to polar oceans (the buoyancy frequency N²  10⁻⁶ s⁻²). Strong mixing occurs in the subinertial and near-inertial ranges such that tracer is well mixed in the boundary layer and an apparent diffusivity there is evaluated at 150–300 cm²/s. This suggests that effective mixing due to these instabilities may play an important role in determining the properties of dense shelf water in the polar regions.