On the entrainment assumption in Schatzmann's integral plume model

The behaviour of stationary, non-passive plumes can be simulated in a reasonably simple and accurate way by integral models. One of the key requirements of these models, but also one of their less well-founded aspects, is the entrainment assumption, which parameterizes turbulent mixing between the plume and the environment. The entrainment assumption developed by Schatzmann and adjusted to a set of experimental results requires four constants and an ad hoc hypothesis to eliminate undesirable terms. With this assumption, Schatzmann’s model exhibits numerical instability for certain cases of plumes with small velocity excesses, due to very fast radius growth. The purpose of this paper is to present an alternative entrainment assumption based on a first-order turbulence closure, which only requires two adjustable constants and seems to solve this problem. The asymptotic behaviour of the new formulation is studied and compared to previous ones. The validation tests presented by Schatzmann are repeated and it is found that the new formulation not only eliminates numerical instability but also predicts more plausible growth rates for jets in co-flowing streams.

The effect of increased convective entrainment on Asian monsoon biases in the MetUM General Circulation Model

We demonstrate that summer precipitation biases in the South Asian monsoon domain are sensitive to increasing the convective parametrisation’s entrainment and detrainment rates in the Met Office Unified Model. We explore this sensitivity to improve our understanding of the biases and inform efforts to improve convective parametrisation. We perform novel targeted experiments in which we increase the entrainment and detrainment rates in regions of especially large precipitation bias. We use these experiments to determine whether the sensitivity at a given location is a consequence of the local change to convection or is a remote response to the change elsewhere. We find that a local change leads to different mean-state responses in comparable regions. When the entrainment and detrainment rates are increased globally, feedbacks between regions usually strengthen the local responses. We choose two regions of tropical ascent that show different mean-state responses, the western equatorial Indian Ocean and western north Pacific, and analyse them as case studies to determine the mechanisms leading to the different responses. Our results indicate that several aspects of a region’s mean-state, including moisture content, sea surface temperature and circulation, play a role in local feedbacks that determine the response to increased entrainment and detrainment.

Characteristics of Langmuir Turbulence in the Ocean Mixed Layer

This study uses large-eddy simulation (LES) to investigate the characteristics of Langmuir turbulence through the turbulent kinetic energy (TKE) budget. Based on an analysis of the TKE budget a velocity scale for Langmuir turbulence is proposed. The velocity scale depends on both the friction velocity and the surface Stokes drift associated with the wave field. The scaling leads to unique profiles of nondimensional dissipation rate and velocity component variances when the Stokes drift of the wave field is sufficiently large compared to the surface friction velocity. The existence of such a scaling shows that Langmuir turbulence can be considered as a turbulence regime in its own right, rather than a modification of shear-driven turbulence. Comparisons are made between the LES results and observations, but the lack of information concerning the wave field means these are mainly restricted to comparing profile shapes. The shapes of the LES profiles are consistent with observed profiles. The dissipation length scale for Langmuir turbulence is found to be similar to the dissipation length scale in the shear-driven boundary layer. Beyond this it is not possible to test the proposed scaling directly using available data. Entrainment at the base of the mixed layer is shown to be significantly enhanced over that due to normal shear turbulence.

The cumulus-capped boundary layer. II: Interface fluxes

This paper considers the relationship between the mean temperature and humidity profiles and the fluxes of heat and moisture at cloud base and the base of the inversion in the cumulus-capped boundary layer. The relationships derived are based on an approximate form of the scalar-flux budget and the scaling properties of the turbulent kinetic energy (TKE) budget. The scalar-flux budget gives a relationship between the change in the virtual potential temperature across either the cloud base transition zone or the inversion and the flux at the base of the layer. The scaling properties of the TKE budget lead to a relationship between the heat and moisture fluxes and the mean subsaturation through the liquid-water flux. The 'jump relation' for the virtual potential temperature at cloud base shows the close connection between the cumulus mass flux in the cumulus-capped boundary layer and the entrainment velocity in the dry-convective boundary layer. Gravity waves are shown to be an important feature of the inversion.

The development of a new dust uplift scheme in the Met Office Unified Model

Aeolian mineral dust aerosol is an important consideration in the Earth's radiation budget as well as a source of nutrients to oceanic and land biota. The modelling of aeolian mineral dust has been improving consistently despite the relatively sparse observations to constrain them. This study documents the development of a new dust emissions scheme in the Met Office Unified ModelTM (MetUM) based on the Dust Entrainment and Deposition (DEAD) module. Four separate case studies are used to test and constrain the model output. Initial testing was undertaken on a large dust event over North Africa in March 2006 with the model constrained using AERONET data. The second case study involved testing the capability of the model to represent dust events in the Middle East without being re-tuned from the March 2006 case in the Sahara. While the model is unable to capture some of the daytime variation in AERONET AOD there is good agreement between the model and observed dust events. In the final two case studies new observations from in situ aircraft data during the Dust Outflow and Deposition to the Ocean (DODO) campaigns in February and August 2006 were used. These recent observations provided further data on dust size distributions and vertical profiles to constrain the model. The modelled DODO cases were also compared to AERONET data to make sure the radiative properties of the dust were comparable to observations. Copyright © 2009 Royal Meteorological Society and Crown Copyright

A model study of oceanic mechanisms affecting Equatorial Pacific sea surface temperature during the 1997-98 El Niño

In this study, the processes affecting sea surface temperature variability over the 1992–98 period, encompassing the very strong 1997–98 El Niño event, are analyzed. A tropical Pacific Ocean general circulation model, forced by a combination of weekly ERS1–2 and TAO wind stresses, and climatological heat and freshwater fluxes, is first validated against observations. The model reproduces the main features of the tropical Pacific mean state, despite a weaker than observed thermal stratification, a 0.1 m s−1 too strong (weak) South Equatorial Current (North Equatorial Countercurrent), and a slight underestimate of the Equatorial Undercurrent. Good agreement is found between the model dynamic height and TOPEX/Poseidon sea level variability, with correlation/rms differences of 0.80/4.7 cm on average in the 10°N–10°S band. The model sea surface temperature variability is a bit weak, but reproduces the main features of interannual variability during the 1992–98 period. The model compares well with the TAO current variability at the equator, with correlation/rms differences of 0.81/0.23 m s−1 for surface currents. The model therefore reproduces well the observed interannual variability, with wind stress as the only interannually varying forcing. This good agreement with observations provides confidence in the comprehensive three-dimensional circulation and thermal structure of the model. A close examination of mixed layer heat balance is thus undertaken, contrasting the mean seasonal cycle of the 1993–96 period and the 1997–98 El Niño. In the eastern Pacific, cooling by exchanges with the subsurface (vertical advection, mixing, and entrainment), the atmospheric forcing, and the eddies (mainly the tropical instability waves) are the three main contributors to the heat budget. In the central–western Pacific, the zonal advection by low-frequency currents becomes the main contributor. Westerly wind bursts (in December 1996 and March and June 1997) were found to play a decisive role in the onset of the 1997–98 El Niño. They contributed to the early warming in the eastern Pacific because the downwelling Kelvin waves that they excited diminished subsurface cooling there. But it is mainly through eastward advection of the warm pool that they generated temperature anomalies in the central Pacific. The end of El Niño can be linked to the large-scale easterly anomalies that developed in the western Pacific and spread eastward, from the end of 1997 onward. In the far-western Pacific, because of the shallower than normal thermocline, these easterlies cooled the SST by vertical processes. In the central Pacific, easterlies pushed the warm pool back to the west. In the east, they led to a shallower thermocline, which ultimately allowed subsurface cooling to resume and to quickly cool the surface layer.

Intermittent release of transients in the slow solar wind: 1. Remote sensing observations

The Heliospheric Imager (HI) instruments on board the STEREO spacecraft are used to analyze the solar wind during August and September 2007. We show how HI can be used to image the streamer belt and, in particular, the variability of the slow solar wind which originates inside and in the vicinity of the streamer belt. Intermittent mass flows are observed in HI difference images, streaming out along the extension of helmet streamers. These flows can appear very differently in images: plasma distributed on twisted flux ropes, V‐shaped structures, or “blobs.” The variety of these transient features may highlight the richness of phenomena that could occur near helmet streamers: emergence of flux ropes, reconnection of magnetic field lines at the tip of helmet streamers, or disconnection of open magnetic field lines. The plasma released with these transient events forms part of the solar wind in the higher corona; HI observations show that these transients are frequently entrained by corotating interaction regions (CIRs), leading to the formation of larger, brighter plasma structures in HI images. This entrainment is used to estimate the trajectory of these plasma ejecta. In doing so, we demonstrate that successive transients can be entrained by the same CIR in the high corona if they emanate from the same corotating source. Some parts of the streamers are more effective sources of transients than others. Surprisingly, evidence is given for the outflow of a recurring twisted magnetic structure, suggesting that the emergence of flux ropes can be recurrent.

Stratospheric water vapour and high climate sensitivity in a version of the HadSM3 climate model

It has been shown previously that one member of the Met Office Hadley Centre single-parameter perturbed physics ensemble – the so-called "low entrainment parameter" member – has a much higher climate sensitivity than other individual parameter perturbations. Here we show that the concentration of stratospheric water vapour in this member is over three times higher than observations, and, more importantly for climate sensitivity, increases significantly when climate warms. The large surface temperature response of this ensemble member is more consistent with stratospheric humidity change, rather than upper tropospheric clouds as has been previously suggested. The direct relationship between the bias in the control state (elevated stratospheric humidity) and the cause of the high climate sensitivity (a further increase in stratospheric humidity) lends further doubt as to the realism of this particular integration. This, together with other evidence, lowers the likelihood that the climate system's physical sensitivity is significantly higher than the likely upper range quoted in the Intergovernmental Panel on Climate Change's Fourth Assessment Report.

Changes in light intensity can influence age at sexual maturity in domestic pullets

1. Shaver White and ISA Brown pullets were reared to 140 d in groups of 8 in cages on a 10-h photoperiod of incandescent light and maintained at an illuminance of 3 or 25 lux, or transferred from 3 to 25 lux or from 25 to 3 lux at 63 or 112 d of age. 2. There was no significant difference in sexual maturity, measured as eggs per 100 bird.d at 139 and 140 d, for ISA Brown maintained on 3 or 25 lux, but Shaver White pullets exposed to constant 3 lux matured significantly later than those maintained on 25 lux. 3. In Shaver Whites, sexual maturity was significantly delayed by an increase from 3 to 25 lux at 63 and 112 d, and advanced by a decrease from 25 to 3 lux at 112 d. Sexual maturity of ISA Browns was not significantly affected by a change in illuminance at 63 or 112 d, though responses were in the same direction as for Shaver Whites. 4. In both breeds, total feed consumed to 112 d was higher for birds on 3 lux than 25 lux, but lower between 112 d and 140 d when birds on 25 lux underwent rapid sexual development. In both breeds, body weight at 63 d was higher for birds exposed to 3 lux than 25 lux, but body weight gain thereafter was similar for the two light intensities. 5. In both breeds, plasma luteinising hormone (LH) concentration at 63 and 112 d was lower in birds maintained on 3 lux than 25 lux. At 63 and 112 d, transfers from 25 to 3 lux depressed, whereas transfers from 3 to 25 lux at 63 d, but not at 112 d, increased plasma LH. 6. Advances or delays in sexual maturity induced by changes in illuminance were not correlated with differences in feed intake, body weight gain, or with changes in plasma LH. 7. One possible explanation for the inverse relationship between the direction of change in illuminance at 63 and 112 d in pullets exposed to a 10-h photoperiod and the age at which they became sexually mature is that changes in light intensity and/or spectral composition affect the entrainment of the circadian rhythm of photoinducibility, to effect a phase shift in the photoinducible phase and/or the responsiveness of phototransduction pathways.

Volatile flavour profile of goat meat extracted by three widely used techniques

Three procedures for the isolation of volatiles from grilled goat meat were compared: dynamic headspace entrainment on Tenax TA, simultaneous steam distillation-extraction, and solid-phase microextraction. Headspace entrainment on Tenax TA extracted the highest number of Maillard-derived volatile compounds. Two hundred and three volatile components were identified: 159 are reported for the first time in goat meat. Most of the volatiles detected (155) were lipid oxidation products, such as hydrocarbons, aldehydes, alcohols, ketones, carboxylic acids and esters. Forty-eight Maillard-derived compounds were identified. comprising pyrazines, pyrroles, thiophenes, furanthiol derivatives, alkyl and alicyclic sulphides, pyridines, and thiazoles. Some reported character impact compounds of cooked meat, e.g., 12-methyltridecanal, (EE)-2,4-decadienal, methional, and dimethyl trisulphide were identified in the volatile profile of goat meat, together with a series of C-2 to C-5 alkylformylcyclopentenes, which have been reported in cooked chicken, pork, beef and lamb, as being important for the characteristic flavour impression of different animal species. (C) 2009 Published by Elsevier Ltd.

Vacuum-induced bubble formation in liquid-tempered chocolate

Bubble inclusion is one of the fastest growing operations practiced in the food industry. A variety of aerated foods is currently available in supermarkets, and newer products are emerging all the time. This paper aims to combine knowledge on chocolate aeration with studies performed on bubble formation and dispersion characteristics. More specifically, we have investigated bubble formation induced by applying vacuum. Experimental methods to determine gas hold-up (volume fraction of air), bubble section distributions along specific planes, and chocolate rheological properties are presented. This study concludes that decreasing pressures elevate gas hold-up values due to an increase in the number of bubble nuclei being formed and release of a greater volume of dissolved gases. Furthermore, bubbles are observed to be larger at lower pressures for a set amount of gas because the internal pressure needs to be in equilibrium with the surrounding pressures. Temperature-induced changes to the properties of the chocolate have less of an effect on bubble formation. On the other hand, when different fats and emulsifiers are added to a standard chocolate recipe, milk fat was found to increase, significantly, the gas hold-up values and the mean bubble-section diameters. It is hypothesized that this behavior is related to the way milk fats, which contain different fatty acids to cocoa butter, crystallize and influence the setting properties of the final product. It is highlighted that apparent viscosity values at low shear rate, as well as setting behavior, play an important role in terms of bubble formation and entrainment.

Evaluating forecasts of the evolution of the cloudy boundary layer using diurnal composites of radar and lidar observations

Observations of boundary-layer cloud have been made using radar and lidar at Chilbolton, Hampshire, UK. These have been compared with output from 7 different global and regional models. Fifty-five cloudy days have been composited to reveal the mean diurnal variation of cloud top and base heights, cloud thickness and liquid water path of the clouds. To enable like-for-like comparison between model and observations, the observations have been averaged on to the grid of each model. The composites show a distinct diurnal cycle in observed cloud; the cloud height exhibits a sinusoidal variation throughout the day with a maximum at around 1600 and a minimum at around 0700 UTC. This diurnal cycle is captured by six of the seven models analysed, although the models generally under-predict both cloud top and cloud base heights throughout the day. The two worst performing models in terms of cloud boundaries also have biases of around a factor of two in liquid water path; these were the only two models that did not include an explicit formulation for cloud-top entrainment.

Using Doppler radar with a simple explicit microphysics model to diagnose problems with ice sublimation depth-scales in forecast models

Several previous studies have attempted to assess the sublimation depth-scales of ice particles from clouds into clear air. Upon examining the sublimation depth-scales in the Met Office Unified Model (MetUM), it was found that the MetUM has evaporation depth-scales 2–3 times larger than radar observations. Similar results can be seen in the European Centre for Medium-Range Weather Forecasts (ECMWF), Regional Atmospheric Climate Model (RACMO) and Météo-France models. In this study, we use radar simulation (converting model variables into radar observations) and one-dimensional explicit microphysics numerical modelling to test and diagnose the cause of the deep sublimation depth-scales in the forecast model. The MetUM data and parametrization scheme are used to predict terminal velocity, which can be compared with the observed Doppler velocity. This can then be used to test the hypothesis as to why the sublimation depth-scale is too large within the MetUM. Turbulence could lead to dry air entrainment and higher evaporation rates; particle density may be wrong, particle capacitance may be too high and lead to incorrect evaporation rates or the humidity within the sublimating layer may be incorrectly represented. We show that the most likely cause of deep sublimation zones is an incorrect representation of model humidity in the layer. This is tested further by using a one-dimensional explicit microphysics model, which tests the sensitivity of ice sublimation to key atmospheric variables and is capable of including sonde and radar measurements to simulate real cases. Results suggest that the MetUM grid resolution at ice cloud altitudes is not sufficient enough to maintain the sharp drop in humidity that is observed in the sublimation zone.