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.     

Summer boundary-layer height at the plateau site of Dome’C,antarctica

Measurements of the mean and turbulent structure of the planetary boundary layer using a sodar and a sonic anemometer, and radiative measurements using a radiometer, were carried out in the summer of 1999–2000 at the Antarctic plateau station of Dome C during a two-month period. At Dome C strong ground-based inversions dominate for most of the year. However, in spite of the low surface temperatures (between −50 and −20 °C), and the surface always covered by snow and ice, a regular daytime boundary-layer evolution, similar to that observed at mid-latitudes, was observed during summertime. The mixed-layer height generally reaches 200–300 m at 1300–1400 LST in high summer (late December, early January); late in the summer (end of January to February), as the solar elevation decreases, it reduces to 100–200 m. A comparison between the mixed-layer height estimated from sodar measurements and that calculated using a mixed-layer growth model shows a rather satisfactory agreement if we assign a value of 0.01–0.02 m s⁻¹ to the subsidence velocity at the top of the mixed layer, and a value of 0.003–0.004 K m⁻¹ to the potential temperature gradient above the mixed layer.     

Measurements of Flows Over Isolated Valleys

Wind-tunnel measurements of the flow over an isolated valley both normal and at an angle (45°) to a simulated neutrally stable atmospheric boundary layer are presented. Attention is concentrated on the nature of the flow within the valley itself. The work formed part of a wider study that included detailed field measurements around an African desert valley and some limited comparisons with that work are included. A scale of about 1:1000 was used for the laboratory work, in which an appropriate combination of hot wire and particle image velocimetry was employed. For a valley normal to the upwind flow, it is shown that the upstream influence of the valley extends to a distance of at least one half of the axial valley width upstream of the leading edge, whereas differences in mean flow and turbulence could be identified well beyond two valley widths from the downwind edge. Non-normal wind angles lead to significant along-valley flows within the valley and, even at two valley heights above the valley ridge level, there remains a significant spanwise flow component. Downwind turbulence levels are somewhat lower in this case, but are still considerably higher than in the undisturbed boundary layer. At both flow angles, there are significant recirculation regions within the valleys, starting from mean separation just beyond the leading edge, but the strong spanwise flow in the 45° case reduces the axial extent of the separated zone. The flow is shown to be in some ways analogous to flow over an isolated hill. Our results usefully enhance the field data and could be used to improve modelling of saltation processes in the field.     

Variance Method to Determine Turbulent Fluxes of Momentum And Sensible Heat in The Stable Atmospheric Surface Layer

Evidence is presented that in the stable atmospheric surface layer turbulent fluxes of heat and momentum can be determined from the standard deviations of longitudinal wind velocity and temperature, σ _u and σ _T respectively, measured at a single level. An attractive aspect of this method is that it yields fluxes from measurements that can be obtained with two-dimensional sonic anemometers. These instruments are increasingly being used at official weather stations, where they replace the standard cup anemometer–wind vane system. With methods such as the one described in this note, a widespread, good quality, flux network can be established, which would greatly benefit the modelling community. It is shown that a ‘variance’ dimensionless height (ζ _σ) defined from σ _u and σ _T is highly related to the ‘conventional’ dimensionless stability parameter ζ=z/L, where z is height and L is the Obukhov length. Empirical functions for ζ _σ are proposed that allow direct calculation of heat and momentum fluxes from σ _u and σ _T. The method performs fairly well also during a night of intermittent turbulence.     

Second Tesseral Harmonic Torque Due to the DynamicS of the Oceanic Surface Layer as Detected by TOPEX/POSEIDON Altimetry 1993-2000

Estimates of the second tesseral torque due to the variations in the radial space position of the mean ocean surface as monitored by TOPEX/POSEIDON altimeter system are derived. The magnitude of the studied torque may be compared to the tidal torque and to the tesseral torque caused by deformations due to the Earths rotation. However, such torque estimates strongly depend on the thickness of the ocean surface layer adopted in the spherical model of which the dynamics is believed to be responsible for the derived torque. The dependence on the thickness is discussed.     

Impact of copper smelting on lakes in the southern Ural Mountains,Russia, inferred from chironomids

Karabash (52°2 N, 60°10 E) is a copper smelting town in the southern Ural Mountains of Russia. The town is affected by sulphur dioxide emissions and deposition of metal-rich particulates from the smelter, acid drainage from old mine workings, and leachates from disused waste dumps and tailings dams. The close proximity of houses to these sources of pollution is of concern to human health and has devastated terrestrial vegetation in the environs. The environmental impact of the smelter on lakes in the area has been assessed using chironomids. Short sediment cores were taken from 16 lakes within a 50 km radius of the smelter and the composition of the chironomid fauna from the bottom of each core, representing conditions prior to the commissioning of the smelter in 1910, was compared with the present chironomid fauna in the surface sediments. Redundancy analysis (RDA) showed that changes in the chironomid fauna of most lakes were driven by trophic change, independent of the industrial activity. Lakes and ponds adjacent to the smelter and waste dumps, which directly receive contaminated waters were devoid of macro- and mesofauna and flora, but there is no evidence that other lakes have been severely impacted by smelter emissions. Local geology ensures that the lakes are well-buffered to the effects of acid deposition which will limit the bioavailability of metals in the water column and sediment.     

Lithological and geochemical record of anthropogenic changes in recent sediments of a small and shallow lake (Lake Pusty Staw,northern Poland)

The article presents the results of lithological and geochemical investigations of recent sediments of Lake Pusty Staw. The analysed sediments document about 300 years of the history of this lake. Historical materials indicate that significant changes in the catchment of the lake took place from the beginning of the 18th century (deforestation and afforestation) followed by 19th century tourist development, and 20th century industrialisation. The sediments were dated using the ²¹⁰Pb and ¹³⁷Cs method and core lithology. These made it possible to establish a reliable chronology to the year 1730. The calculated sedimentation rates ranged between 0.17 and 0.83 cm year^–1. There was a period of intensified erosion caused by land clearance between 1734 and 1810, which resulted in an acceleration in sedimentation rate (0.36 cm year^–1) and a change of lithological type of sediment from detritus gyttja to clayey gyttja. On the basis of Cu, Zn, Cd and Pb content in the sediments, it was determined that the beginning of pollution of the lake dates back to the 19th century and was caused probably by the existence of a health resort. A systematic increase in pollution occurred in the 20th century as a result of industrial plants. Normalised with respect to Al, the content of heavy metals in polluted sediments was from several to twenty times higher than in sediments of the preindustrial period, and a comparison of historical materials with changes of sediments in most cases made it possible to identify the direct causes of the increase in pollution.     

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.