‘A Case Study on a Polar Cold Air Outbreak over Fram Strait using a Mesoscale Weather Prediction Model’

The mesoscale weather prediction model ’Lokal-Modell’ (LM) of the Deutscher Wetterdienst is applied to the situation of an Arctic cold air outbreak in the Fram Strait region in April 1998. Observations are available from a flight along 5⁰E carried out during the ARTIST campaign. Initial and time-dependent boundary data for the simulation are taken from a larger scale operational model system. Using the standard configuration of LM, the simulation reproduced the propagation of cold air and the characteristic structure of the atmospheric boundary layer (ABL) in fair agreement with the observations. However, a detailed comparison revealed three basic problems. Firstly, there is evidence that the available data on sea-ice conditions were insufficient approximations to the true state for several reasons. A modification of the sea-ice data towards observations revealed that parts of the discrepancies were due to the original sea-ice data. Secondly, a control run with the model in its standard configuration shows an insufficient warming of the ABL downstream of the ice edge due to underestimation of surface heat fluxes. A simple modification of the approach for the scalar roughness length resulted in the strongest benefit, while comparative studies showed only a slight sensitivity to different types of parametrisation of turbulent mixing or the inclusion of an additional moist convection parametrisation. Thirdly, in all the simulations the deepening of the convective ABL downstream of the ice edge is weaker than observed. This may be partly due to the thermal stratification above the ABL in the analysis data, which is more stable than observed; but it may also be a hint to the fact that processes near the inversion are insufficiently parametrised in mesoscale models with resolutions as used in LM. The simulated cloud layer in the convective ABL is similar to that observed with respect to condensate content, a sharply defined cloud top, a diffuse lower bound, and continuous light precipitation.     

Boundary-Layer Adjustment Over Small-Scale Changes of Surface Heat Flux

Four months of eddy correlation data collected over a grass field and a nearby sage brush community are analyzed to examine the adjustment of the boundary-layer structure as it flows from the heated brush to the snow-covered grass. The grass site includes a 34-m tower with seven levels of eddy correlation data. The midday heat flux over the snow-covered grass and bare ground surfaces is often downward particularly with melting conditions, while the corresponding heat flux over the brush is almost always upward. For most of these cases, a stable internal boundary layer over the snow is well defined in terms of vertical profiles of the buoyancy flux over the snow-covered grass. The stable internal boundary layer is generally embedded within a deeper layer of flux divergence corresponding to increasing upward heat flux with height above the internal boundary layer. With thin snow cover, the surface heat flux over the grass is weak upward due to heating of grass protruding above the snow so that the flow adjusts to a decrease of the upward surface heat flux in the downwind direction. This common case of an adjusting boundary layer contrasts with the formation of an internal boundary layer due to a change of sign of the surface heat in flux the downwind direction. The adjustment of the boundary layer to the decrease of the surface heat flux leads to vertical divergence of the upward heat flux in contrast to the usual heated boundary layer over homogeneous surfaces. The consequences of the cooling due to the vertical divergence of the heat flux are discussed in terms of the heat budget of the adjusting and internal boundary layers.     

Sodar Detected Top-Down Convection in a Nocturnal Cloud-Topped Boundary Layer: A Case Study

Nocturnal convection, originating in a well-mixed marine cloud-topped boundary layer, advected onshore, was observed using a Doppler sodar on the Tyrrhenian coast in Italy. The horizontal and vertical dimensions of the downdrafts were evaluated. The oscillation frequency triggered by the downdrafts at the inversion layer, derived from the harmonic analysis of the sodar measured vertical velocity (w), is compared with the Brunt-Vaisala frequency, obtained from the rawinsonde temperature profile. A similarity function for the ²_w vertical profile was used to fit the sodar experimental data and to retrieve the depth of the mixing layer and the sensible heat flux at the top of the cloud layer. The results are in agreement with the convection layer depth observed in the sodar echoes facsimile record, and with the energy budget evaluated at the top of the cloud layer using the rawinsonde profiles.     

Gas formation mechanism of marine carbonate source rocks in China

It has been proven in exploration practice that thecarbonates in China not only can generate hydrocar-bons, but also form commercial reservoirs. The car-bonates are different from clastic rocks in view of theirdeposition environment as well as their sedimentaryand diagenetic processes. Therefore, the evaluationcriteria and hydrocarbon generation mechanism forcarbonates can not be the same as that for clasticrocks, and it is important to establish a special hydro-carbon generation mechanism and…     

Magma underplating and Hannuoba present crust-mantle transitional zone composition: Xenolith petrological and geochemical evidence

Studies on the deep-seated xenoliths from global volcanoes reveal that the present petrological crust-mantle boundary between the lower crust and the upper mantle actually is a transitional layer from mainly mafic granulites to ultramafic spinel lher-zolites[1,2], i.e. a transitional zone distinctive from the seismological Moho[3]. Oceanic lithosphere crust- mantle transitional zone can be established from the study on the exposed ophiolites. However, as for the continental lithosphere, since …     

Sudden enhancement of sedimentation flux of <Superscript>210</Superscript>Pb<Subscript>ex</Subscript> as an indicator of lake productivity as exemplified by Lake Chenghai

The fundamental assumption of ²¹⁰Pb sediment dating is the stable flux of ²¹⁰Pb_ex, which was derived from atmosphere and then transferred into sediments via lake water. When the sedimentation rate is relatively constant, the ²¹⁰Pb_ex activity in sediments will be exponentially reduced with sedimentation age. ²¹⁰Pb_ex in lake water is incorporated into sediments mainly via organic particulates. If the sedimentation flux of organic matter in lake water is suddenly increased, ²¹⁰Pb_ex will be significantly deposited and then transferred into sediments. On the one hand such sudden purification effect is obviously unfit for the fundamental assumption of ²¹⁰Pb dating; on the other hand, the sudden enhancement of ²¹⁰Pb_ex flux would be indicative of the conspicuous variation of primary productivity of lake water. This problem will be discussed in accordance with the variation trend of ²¹⁰Pb_ex in the vertical profile of recent sediments of Lake Chenghai, Yunnan Province. The sediment core was collected from the deep-water area of Lake Chenghai in June 1997. The vertical profile of ¹³⁷Cs activity is characterized by a tree-peak pattern. This profile gave reliable ages, and also showed the stability of sediment accumulation in the recent ten years. The vertical profile of ²¹⁰Pb_ex activity displays a specific distribution of peaks, and is similar to the vertical profile of C_org. This phenomenon seems to be related to the mechanism of constraining the transfer of ²¹⁰Pb_ex into lake sediments. The average atomic ratios of H_org/C_org and C_org/N_org in Lake Chenghai sediments are 5.51 and 7.04, respectively, indicating that the organic matter was predominantly derived from the remains of endogenic algae. In terms of the three-stage evolutionary characteristics of organic matter in sediments, i.e., “deposition-de-composition-accumulation”, the sedimentation fluxes (F(C_org)) of organic carbon (C_org) since 1970 were calculated by modeling. The sedimentation fluxes of ²¹⁰Pb_ex (F(²¹⁰Pb_ex)) in different years display good synchronous relations with the sedimentation flux of organic carbon (F(C_org)), especially in the years of 1972–1974 and 1986–1989. The variation of F(C_org) led to the variation of F(²¹⁰Pb_ex); the variation of F(²¹⁰Pb_ex) reflects, to some extent, the historical variation of lake productivity.     

Pliocene—pleistocene boundary in the Kleszczów Graben at Belchatów,central Poland

The Pliocene-Pleistocene boundary, defined using pollen analysis, has been recorded in sections at the Belchatów outcrop, Kleszczów Graben, central Poland. The boundary is located at the top of the green clay subunit of the Tertiary Clayey-sandy unit. A pollen flora of Reuverian C type, dominated by pine, 10–20% of ‘Reuverian elements’ and scarce thermophilous Tertiary trees, was found directly below this subunit. Cold stage floras of Praetiglian type were, in turn, recorded in the ?ekińsko Formation, lying directly on the green clays. The Pliocene—Pleistocene transition at Belchatów demonstrates the characteristic features of this boundary in western Europe, that is, the disappearance of Tertiary thermophilous tree pollen taxa, followed by a cooling to open-forest or forest-steppe conditions. Both pollen analysis and geological data suggest continuous sedimentation at least from the Upper Pliocene to the Prae-Tiglian, with no change of sediment sources and no marked hiatuses. The fluvial deposits of the Praetiglian ?ekińsko Formation, although Pleistocene in age, mark a final period in the Tertiary evolution of the Kleszczów Graben. Later, there is a hiatus to the deposition of the first glacigenic sediments in the region possibly of Elsterian age.     

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.