Perturbation growth at the convective scale for CSIP IOP18

The Met Office Unified Model is run for a case observed during Intensive Observation Period 18 (IOP18) of the Convective Storms Initiation Project (CSIP). The aims are to identify the physical processes that lead to perturbation growth at the convective scale in response to model-state perturbations and to determine their sensitivity to the character of the perturbations. The case is strongly upper-level forced but with detailed mesoscale/convective-scale evolution that is dependent on smaller-scale processes. Potential temperature is perturbed within the boundary layer. The effects on perturbation growth of both the amplitude and typical scalelength of the perturbations are investigated and perturbations are applied either sequentially (every 30 min throughout the simulation) or at specific times. The direct effects (within one timestep) of the perturbations are to generate propagating Lamb and acoustic waves and produce generally small changes in cloud parameters and convective instability. In exceptional cases a perturbation at a specific gridpoint leads to switching of the diagnosed boundary-layer type or discontinuous changes in convective instability, through the generation or removal of a lid. The indirect effects (during the entire simulation) are changes in the intensity and location of precipitation and in the cloud size distribution. Qualitatively different behaviour is found for strong (1K amplitude) and weak (0.01K amplitude) perturbations, with faster growth after sunrise found only for the weaker perturbations. However, the overall perturbation growth (as measured by the root-mean-square error of accumulated precipitation) reaches similar values at saturation, regardless of the perturbation characterisation.

How well do high resolution models reproduce tropical convection?

Cascade is a multi-institution project studying the temporal and spatial organization of tropical convective systems. While cloud resolving numerical models can reproduce the observed diurnal cycle of such systems they are sensitive to the chosen resolution. As part of this effort, we are comparing results from the Met. Office Unified Model to data from the Global Earth Radiation Budget satellite instrument over the African Monsoon Interdisciplinary Analyses region of North Africa. We use a variety of mathematical techniques to study the outgoing radiation and the evolution of properties such as the cloud size distribution. The effectiveness of various model resolutions is tested with a view to determining the optimum balance between resolution and the need to reproduce the observations.

Estimating the spatial scale of herbicide and soil interactions by nested sampling, hierarchical analysis of variance and residual maximum likelihood

An unbalanced nested sampling design was used to investigate the spatial scale of soil and herbicide interactions at the field scale. A hierarchical analysis of variance based on residual maximum likelihood (REML) was used to analyse the data and provide a first estimate of the variogram. Soil samples were taken at 108 locations at a range of separating distances in a 9 ha field to explore small and medium scale spatial variation. Soil organic matter content, pH, particle size distribution, microbial biomass and the degradation and sorption of the herbicide, isoproturon, were determined for each soil sample. A large proportion of the spatial variation in isoproturon degradation and sorption occurred at sampling intervals less than 60 m, however, the sampling design did not resolve the variation present at scales greater than this. A sampling interval of 20-25 m should ensure that the main spatial structures are identified for isoproturon degradation rate and sorption without too great a loss of information in this field.

Independent evaluation of the ability of spaceborne radar and lidar to retrieve the microphysical and radiative properties of ice clouds

The combination of radar and lidar in space offers the unique potential to retrieve vertical profiles of ice water content and particle size globally, and two algorithms developed recently claim to have overcome the principal difficulty with this approach-that of correcting the lidar signal for extinction. In this paper "blind tests" of these algorithms are carried out, using realistic 94-GHz radar and 355-nm lidar backscatter profiles simulated from aircraft-measured size spectra, and including the effects of molecular scattering, multiple scattering, and instrument noise. Radiation calculations are performed on the true and retrieved microphysical profiles to estimate the accuracy with which radiative flux profiles could be inferred remotely. It is found that the visible extinction profile can be retrieved independent of assumptions on the nature of the size distribution, the habit of the particles, the mean extinction-to-backscatter ratio, or errors in instrument calibration. Local errors in retrieved extinction can occur in proportion to local fluctuations in the extinction-to-backscatter ratio, but down to 400 m above the height of the lowest lidar return, optical depth is typically retrieved to better than 0.2. Retrieval uncertainties are greater at the far end of the profile, and errors in total optical depth can exceed 1, which changes the shortwave radiative effect of the cloud by around 20%. Longwave fluxes are much less sensitive to errors in total optical depth, and may generally be calculated to better than 2 W m(-2) throughout the profile. It is important for retrieval algorithms to account for the effects of lidar multiple scattering, because if this is neglected, then optical depth is underestimated by approximately 35%, resulting in cloud radiative effects being underestimated by around 30% in the shortwave and 15% in the longwave. Unlike the extinction coefficient, the inferred ice water content and particle size can vary by 30%, depending on the assumed mass-size relationship (a problem common to all remote retrieval algorithms). However, radiative fluxes are almost completely determined by the extinction profile, and if this is correct, then errors in these other parameters have only a small effect in the shortwave (around 6%, compared to that of clear sky) and a negligible effect in the longwave.

Using aircraft measurements to determine the refractive index of Saharan dust during the DODO Experiments

Much uncertainty in the value of the imaginary part of the refractive index of mineral dust contributes to uncertainty in the radiative effect of mineral dust in the atmosphere. A synthesis of optical, chemical and physical in-situ aircraft measurements from the DODO experiments during February and August 2006 are used to calculate the refractive index mineral dust encountered over West Africa. Radiative transfer modeling and measurements of broadband shortwave irradiance at a range of altitudes are used to test and validate these calculations for a specific dust event on 23 August 2006 over Mauritania. Two techniques are used to determine the refractive index: firstly a method combining measurements of scattering, absorption, size distributions and Mie code simulations, and secondly a method using composition measured on filter samples to apportion the content of internally mixed quartz, calcite and iron oxide-clay aggregates, where the iron oxide is represented by either hematite or goethite and clay by either illite or kaolinite. The imaginary part of the refractive index at 550 nm (ni550) is found to range between 0.0001 i to 0.0046 i, and where filter samples are available, agreement between methods is found depending on mineral combination assumed. The refractive indices are also found to agree well with AERONET data where comparisons are possible. ni550 is found to vary with dust source, which is investigated with the NAME model for each case. The relationship between both size distribution and ni550 on the accumulation mode single scattering albedo at 550 nm (ω0550) are examined and size distribution is found to have no correlation to ω0550, while ni550 shows a strong linear relationship with ω0550. Radiative transfer modeling was performed with different models (Mie-derived refractive indices, but also filter sampling composition assuming both internal and external mixing). Our calculations indicate that Mie-derived values of ni550 and the externally mixed dust where the iron oxide-clay aggregate corresponds to the goethite-kaolinite combination result in the best agreement with irradiance measurements. The radiative effect of the dust is found to be very sensitive to the mineral combination (and hence refractive index) assumed, and to whether the dust is assumed to be internally or externally mixed.

Hydrolysis of N2O5 on sub-micron mineral salt aerosols

The kinetics of reactive uptake of gaseous N2O5 on sub-micron aerosol particles composed of aqueous ammonium sulfate, ammonium hydrogensulfate and sodium nitrate has been investigated. Uptake was measured in a laminar flow reactor, coupled with a differential mobility analyser (DMA) to obtain the aerosol size distribution, with N2O5 detection using NO chemiluminescence. FTIR spectroscopy was used to obtain information about the composition and water content of the aerosol particles under the conditions used in the kinetic measurements. The aerosols were generated by the nebulisation of aqueous salt solutions. The uptake coefficient on the sulfate salts was in the range [gamma]=0.0015 to 0.033 depending on temperature, humidity and phase of the aerosol. On sodium nitrate aerosols the values were much lower, [gamma]<0.001, confirming the inhibition of N2O5 hydrolysis by nitrate ions. At high humidity (>50% r.h.) the uptake coefficient on liquid sulfate aerosols is independent of water content, but at lower humidity, especially below the efflorescence point, the reactivity of the aerosol declines, correlating with the lower water content. The lower uptake rate on solid aerosols may be due to limitations imposed by the liquid volume in the particles. Uptake on sulfate aerosols showed a negative temperature dependence at T>290 K but no significant temperature dependence at lower temperatures. The results are generally consistent with previous models of N2O5 hydrolysis where the reactive intermediate is NO2+ produced by autoionisation of nitrogen pentoxide in the condensed phase.

Black carbon measurements in the boundary layer over western and northern Europe

Europe is a densely populated region that is a significant global source of black carbon (BC) aerosol, but there is a lack of information regarding the physical properties and spatial/vertical distribution of rBC in the region. We present the first aircraft observations of sub-micron refractory BC (rBC) aerosol concentrations and physical properties measured by a single particle soot photometer (SP2) in the lower troposphere over Europe. The observations spanned a region roughly bounded by 50° to 60° N and from 15° W to 30° E. The measurements, made between April and September 2008, showed that average rBC mass concentrations ranged from about 300 ng m−3 near urban areas to approximately 50 ng m−3 in remote continental regions, lower than previous surface-based measurements. rBC represented between 0.5 and 3% of the sub-micron aerosol mass. Black carbon mass size distributions were log-normally distributed and peaked at approximately 180 nm, but shifted to smaller diameters (~160 nm) near source regions. rBC was correlated with carbon monoxide (CO) but had different ratios to CO depending on location and air mass. Light absorption coefficients were measured by particle soot absorption photometers on two separate aircraft and showed similar geographic patterns to rBC mass measured by the SP2. We summarize the rBC and light absorption measurements as a function of longitude and air mass age and also provide profiles of rBC mass concentrations and size distribution statistics. Our results will help evaluate model-predicted regional rBC concentrations and properties and determine regional and global climate impacts from rBC due to atmospheric heating and surface dimming.

Evaluation of the model representation of the evolution of convective systems using satellite observations of outgoing longwave radiation

We introduce a technique for assessing the diurnal development of convective storm systems based on outgoing longwave radiation fields. Using the size distribution of the storms measured from a series of images, we generate an array in the lengthscale-time domain based on the standard score statistic. It demonstrates succinctly the size evolution of storms as well as the dissipation kinematics. It also provides evidence related to the temperature evolution of the cloud tops. We apply this approach to a test case comparing observations made by the Geostationary Earth Radiation Budget instrument to output from the Met Office Unified Model run at two resolutions. The 12km resolution model produces peak convective activity on all lengthscales significantly earlier in the day than shown by the observations and no evidence for storms growing in size. The 4km resolution model shows realistic timing and growth evolution although the dissipation mechanism still differs from the observed data.

Joint effects of population density and toxicant exposure on population dynamics of Capitella sp I

Very few studies have analyzed the dependence of population growth rate on population density, and even fewer have considered interaction effects of density and other stresses, such as exposure to toxic chemicals. Yet without such studies we cannot know whether chemicals harmful at low density have effects on carrying capacity or, conversely, whether chemicals reducing carrying capacity are also harmful at low density, impeding a population's capacity to recover from disturbance. This study examines the combined effects of population density and a toxicant (fluoranthene) on population growth rate (pgr) and carrying capacity using the deposit-feeding polychaete Capitella sp. I as a test organism. Populations were initiated with a stable age distribution, and population density and age/size distribution were followed during a period of 28 wk. Fluoranthene (FLU), population density, and their interaction influenced population growth rate. Population growth rate declined linearly with the logarithm of population biomass, but the slope of the relationship was steeper for the control populations than for populations exposed to 50 mug FLU/(g sediment dry mass). Populations exposed to 150 mug FLU/(g sediment dry mass) went extinct after 8 wk of exposure. Despite concerns that toxicant effects would be exacerbated at high density, we found the reverse to be the case, and effects of fluoranthene on population growth rate were much reduced in the region of carrying capacity. Fluoranthene did. reduce carrying capacity by 46%, and this could haven important implications for interacting species and/or sediment biogeochemical processes.

Low-temperature sol-gel preparation of ordered nanoparticles of tungsten carbide/oxide

Tungsten carbide/oxide particles have been prepared by the gel precipitation of tungstic acid in the presence of an organic gelling agent [10% ammonium poly(acrylic acid) in water, supplied by Ciba Specialty Chemicals]. The feed solution; a homogeneous mixture of sodium tungstate and ammonium poly(acrylic acid) in water, was dropped from a 1-mm jet into hydrochloric acid saturated hexanol/concentrated hydrochloric acid to give particles of a mixture of tungstic acid and poly(acrylic acid), which, after drying in air at 100 degrees C and heating to 900 degrees C in argon for 2 h, followed by heating in carbon dioxide for a further 2 h and cooling, gives a mixture of WO, WC, and a trace of NaxWO3, with the carbon for the formation of WC being provided by the thermal carbonization of poly(acrylic acid). The pyrolyzed product is friable and easily broken down in a pestle and mortar to a fine powder or by ultrasonics, in water, to form a stable colloid. The temperature of carbide formation by this process is significantly lower (900 degrees C) than that reported for the commercial preparation of tungsten carbide, typically > 1400 degrees C. In addition, the need for prolonged grinding of the constituents is obviated because the reacting moieties are already in intimate contact on a molecular basis. X-ray diffraction, particle sizing, transmission electron microscopy, surface area, and pore size distribution studies have been carried out, and possible uses are suggested. A flow diagram for the process is described.

Relationships between gluten rheological properties and hearth loaf characteristics

The rheological properties of fresh gluten in small amplitude oscillation in shear (SAOS) and creep recovery after short application of stress was related to the hearth breadbaking performance of wheat flours using the multivariate statistics partial least squares (PLS) regression. The picture was completed by dough mixing and extensional properties, flour protein size distribution determined by SE-HPLC, and high molecular weight glutenin subunit (HMW-GS) composition. The sample set comprised 20 wheat cultivars grown at two different levels of nitrogen fertilizer in one location. Flours yielding stiffer and more elastic glutens, with higher elastic and viscous moduli (G' and G") and lower tan 8 values in SAOS, gave doughs that were better able to retain their shape during proving and baking, resulting in breads of high form ratios. Creep recovery measurements after short application of stress showed that glutens from flours of good breadmaking quality had high relative elastic recovery. The nitrogen fertilizer level affected the protein size distribution by an increase in monomeric proteins (gliadins), which gave glutens of higher tan delta and flatter bread loaves (lower form ratio).

Comparison of small and large deformation rheological properties of wheat dough and gluten

The rheological properties of dough and gluten are important for end-use quality of flour but there is a lack of knowledge of the relationships between fundamental and empirical tests and how they relate to flour composition and gluten quality. Dough and gluten from six breadmaking wheat qualities were subjected to a range of rheological tests. Fundamental (small-deformation) rheological characterizations (dynamic oscillatory shear and creep recovery) were performed on gluten to avoid the nonlinear influence of the starch component, whereas large deformation tests were conducted on both dough and gluten. A number of variables from the various curves were considered and subjected to a principal component analysis (PCA) to get an overview of relationships between the various variables. The first component represented variability in protein quality, associated with elasticity and tenacity in large deformation (large positive loadings for resistance to extension and initial slope of dough and gluten extension curves recorded by the SMS/Kieffer dough and gluten extensibility rig, and the tenacity and strain hardening index of dough measured by the Dobraszczyk/Roberts dough inflation system), the elastic character of the hydrated gluten proteins (large positive loading for elastic modulus [G'], large negative loadings for tan delta and steady state compliance [J(e)(0)]), the presence of high molecular weight glutenin subunits (HMW-GS) 5+10 vs. 2+12, and a size distribution of glutenin polymers shifted toward the high-end range. The second principal component was associated with flour protein content. Certain rheological data were influenced by protein content in addition to protein quality (area under dough extension curves and dough inflation curves [W]). The approach made it possible to bridge the gap between fundamental rheological properties, empirical measurements of physical properties, protein composition, and size distribution. The interpretation of this study gave indications of the molecular basis for differences in breadmaking performance.

Stress relaxation behavior of wheat dough, gluten, and gluten protein fractions

Relaxation behavior was measured for dough, gluten and gluten protein fractions obtained from the U.K. biscuitmaking flour, Riband, and the U.K. breadmaking flour, Hereward. The relaxation spectrum, in which relaxation times (tau) are related to polymer molecular size, for dough showed a broad molecular size distribution, with two relaxation processes: a major peak at short times and a second peak at times longer than 10 sec, which is thought to correspond to network structure, and which may be attributed to entanglements and physical cross-links of polymers. Relaxation spectra of glutens were similar to those for the corresponding doughs from both flours. Hereward gluten clearly showed a much more pronounced second peak in relaxation spectrum and higher relaxation modulus than Riband gluten at the same water content. In the gluten protein fractions, gliadin and acetic acid soluble glutenin only showed the first relaxation process, but gel protein clearly showed both the first and second relaxation processes. The results show that the relaxation properties of dough depend on its gluten protein and that gel protein is responsible for the network structure for dough and gluten.

Transient development of whipped cream properties

Commercially available UHT cream was tempered at 4 degrees C for 24 h and whipped for different times: 3. 6. 9 and 12 nun. The following properties of cream were measured: rheological and interfacial properties. overrun and size distribution of air bubbles. The whipping process changes the properties of cream, which exhibits viscoelastic behaviour with a high influence of elastic component. The air bubbles incorporated during the process result in forming stronger foam containing smaller bubbles. and also give a higher overrun. These changes are observed around 9 min of whipping. when the amount of air is sufficient to create a stable structure. Further whipping reduces the overrun and the foam partly collapses: this may be caused by aggregation of fat droplets. (c) 2005 Elsevier Ltd. All rights reserved.

Assessment of the quality of cottage cheese produced from standard and protein-fortified skim milk

The effects of milk protein fortification on the texture and microstructure of cottage cheese curd were evaluated. Protein powder (92.6% protein) was added to the skim milk at a level of 0.4% (w/w) to produce curds. Control curds with no protein powder addition were also produced. These curds were analysed for differences in yield, total solids, curd size, texture and structure. It was found that the addition of protein powder contributed to a significant yield increase, which can be attributed to increased water retention, with better curd size distribution. Control curds were firmer than the fortified curds and the structure showed less open-pore structure as revealed by electron microscopy. However, the addition of dressing masked the textural differences, and a sensory panel was unable to distinguish between cheeses produced from fortified milk and controls.