Two methods for modelling the propagation of terahertz radiation in a layered structure

Modelling the interaction of terahertz(THz) radiation with biological tissueposes many interesting problems. THzradiation is neither obviously described byan electric field distribution or anensemble of photons and biological tissueis an inhomogeneous medium with anelectronic permittivity that is bothspatially and frequency dependent making ita complex system to model.A three-layer system of parallel-sidedslabs has been used as the system throughwhich the passage of THz radiation has beensimulated. Two modelling approaches havebeen developed a thin film matrix model anda Monte Carlo model. The source data foreach of these methods, taken at the sametime as the data recorded to experimentallyverify them, was a THz spectrum that hadpassed though air only.Experimental verification of these twomodels was carried out using athree-layered in vitro phantom. Simulatedtransmission spectrum data was compared toexperimental transmission spectrum datafirst to determine and then to compare theaccuracy of the two methods. Goodagreement was found, with typical resultshaving a correlation coefficient of 0.90for the thin film matrix model and 0.78 forthe Monte Carlo model over the full THzspectrum. Further work is underway toimprove the models above 1 THz.

The inhibitory effect of natural microflora of food on growth of Listeria monocytogenes in enrichment broths

The aims of this study were to (i) compare the inhibitory effects of the natural microflora of different foods on the growth of Listeria monocytogenes during enrichment in selective and non-selective broths; (ii) to isolate and identify components of the microflora of the most inhibitory food; and (iii) to determine which of these components was most inhibitory to growth of L. monocytogenes in co-culture studies. Growth of an antibioticresistant marker strain of L. monocytogenes was examined during enrichment of a range of different foods in Tryptone Soya Broth (TSB), Half Fraser Broth (HFB) and Oxoid Novel Enrichment (ONE) Broth. Inhibition of L. monocytogenes was greatest in the presence of minced beef, salami and soft cheese and least with prepared fresh salad and chicken pâté. For any particular food the numbers of L. monocytogenes present after 24 h enrichment in different broths increased in the order: TSB, HFB and ONE Broth. Numbers of L. monocytogenes recovered after enrichment in TSB were inversely related to the initial aerobic plate count (APC) in the food but with only a moderate coefficient of determination (R2) of 0.51 implying that microbial numbers and the composition of the microflora both influenced the degree of inhibition of L. monocytogenes. In HFB and ONE Broth the relationship between APC and final L. monocytogenes counts was weaker. The microflora of TSB after 24 h enrichment of minced beef consisted of lactic acid bacteria, Brochothrix thermosphacta, Pseudomonas spp., Enterobacteriaceae, and enterococci. In co-culture studies of L. monocytogenes with different components of the microflora in TSB, the lactic acid bacteria were the most inhibitory followed by the Enterobacteriaceae. The least inhibitory organisms were Pseudomonas sp., enterococci and B. thermosphacta. In HFB and ONE Broth the growth of Gram-negative organisms was inhibited but lactic acid bacteria still reached high numbers after 24 h. A more detailed study of the growth of low numbers of L. monocytogenes during enrichment of minced beef in TSB revealed that growth of L. monocytogenes ceased at a cell concentration of about 102 cfu/ml when lactic acid bacteria entered stationary phase. However in ONE Broth growth of lactic acid bacteria was slower than in TSB with a longer lag time allowing L. monocytogenes to achieve much higher numbers before lactic acid bacteria reached stationary phase. This work has identified the relative inhibitory effects of different components of a natural food microflora and shown that the ability of low numbers of L. monocytogenes to achieve high cell concentrations is highly dependent on the extent to which enrichment media are able to inhibit or delay growth of the more effective competitors.

Effect of milk fat concentration and gel firmness on syneresis during curd stirring in cheese-making

An experiment was undertaken to investigate the effect of milk fat level (0%, 2.5% and 5.0% w/w) and gel firmness level at cutting (5, 35 and 65 Pa) on indices of syneresis, while curd was undergoing stirring. The curd moisture content, yield of whey, fat in whey and casein fines in whey were measured at fixed intervals between 5 and 75 min after cutting the gel. The casein level in milk and clotting conditions was kept constant in all trials. The trials were carried out using recombined whole milk in an 11 L cheese vat. The fat level in milk had a large negative effect on the yield of whey. A clear effect of gel firmness on casein fines was observed. The best overall prediction, in terms of coefficient of determination, was for curd moisture content using milk fat concentration, time after gel cutting and set-to-cut time (R2 = 0.95).

Validation of the GRAPE single view aerosol retrieval for ATSR-2 and insights into the long term global AOD trend over the ocean

The Global Retrieval of ATSR Cloud Parameters and Evaluation (GRAPE) project has produced a global data-set of cloud and aerosol properties from the Along Track Scanning Radiometer-2 (ATSR-2) instrument, covering the time period 1995�2001. This paper presents the validation of aerosol optical depths (AODs) over the ocean from this product against AERONET sun-photometer measurements, as well as a comparison to the Advanced Very High Resolution Radiometer (AVHRR) optical depth product produced by the Global Aerosol Climatology Project (GACP). The GRAPE AOD over ocean is found to be in good agreement with AERONET measurements, with a Pearson's correlation coefficient of 0.79 and a best-fit slope of 1.0±0.1, but with a positive bias of 0.08±0.04. Although the GRAPE and GACP datasets show reasonable agreement, there are significant differences. These discrepancies are explored, and suggest that the downward trend in AOD reported by GACP may arise from changes in sampling due to the orbital drift of the AVHRR instruments.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain

Scintillometry is an established technique for determining large areal average sensible heat fluxes. The scintillometer measurement is related to sensible heat flux via Monin–Obukhov similarity theory, which was developed for ideal homogeneous land surfaces. In this study it is shown that judicious application of scintillometry over heterogeneous mixed agriculture on undulating topography yields valid results when compared to eddy covariance (EC). A large aperture scintillometer (LAS) over a 2.4 km path was compared with four EC stations measuring sensible (H) and latent (LvE) heat fluxes over different vegetation (cereals and grass) which when aggregated were representative of the LAS source area. The partitioning of available energy into H and LvE varied strongly for different vegetation types, with H varying by a factor of three between senesced winter wheat and grass pasture. The LAS derived H agrees (one-to-one within the experimental uncertainty) with H aggregated from EC with a high coefficient of determination of 0.94. Chronological analysis shows individual fields may have a varying contribution to the areal average sensible heat flux on short (weekly) time scales due to phenological development and changing soil moisture conditions. Using spatially aggregated measurements of net radiation and soil heat flux with H from the LAS, the areal averaged latent heat flux (LvELAS) was calculated as the residual of the surface energy balance. The regression of LvELAS against aggregated LvE from the EC stations has a slope of 0.94, close to ideal, and demonstrates that this is an accurate method for the landscape-scale estimation of evaporation over heterogeneous complex topography.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KMGAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KMGAP is based on the PRA model framework (P¨oschl-Rudich- Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modeled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system and the computational constraints, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmospheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270K is close to unity (Winkler et al., 2006). Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for efficient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Standardised ileal digestibility of crude protein and amino acids of UK-grown peas and faba beans by broilers

In view of the increasing interest in home-grown legumes as components of diets for non-ruminant livestock and in an attempt to reduce the reliance on imported soya bean meal (SBM), two experiments were conducted to evaluate samples of peas and faba beans for their standardised ileal digestibility (SID) of amino acids determined with young broiler chicks. Experiment 1 evaluated six faba bean and seven pea cultivars and Experiment 2 evaluated two faba bean and three pea cultivars as well as a sample of soya bean meal provided as a reference material. Peas and beans were added at 750g/kg as the only source of protein/amino acids in a semi-synthetic diet containing the inert marker titanium dioxide; SBM was added, in a control diet, at 500g/kg. Each diet was fed to six replicates of a cage containing two Ross-type broilers for 96h at which point birds were culled allowing removal of ileal digesta. Chemical analyses allowed the calculation of the coefficient of SID of amino acids. There were no differences between samples of the same pulse species (P>0.05) but peas had higher values (P<0.05), similar to SBM, than beans. Trypsin inhibitor content (expressed as g trypsin inhibitor units/mg sample) of all pea samples was low and in the range 0.83–1.77mg/kg. There was relatively little variation in bean tannin content and composition amongst the coloured-flowered varieties; however, the white-flowered cultivar had no tannins. There was no correlation between tannin content and coefficient of SID. The content of SID of amino acids (g/kg legume) was higher in SBM when compared with peas and beans by virtue of having higher total concentrations.

Canopy characteristics of contrasting clones of cacao (Theobroma cacao)

Canopy characteristics (leaf area index, fractional light interception, extinction coefficient) of mature trees of ten clonally propagated cacao cultivars were measured over a period of 14 months at an experimental site in Bahia, Brazil. Differences in leaf area index between clones became more pronounced over time. When an approximately constant leaf area index was reached (after about nine months), LAI varied between clones from 2.8 to 4.5. Clonal differences in the relationship between leaf area index and fractional light interception implied differences in canopy architecture, as reflected by the range of extinction coefficients (mean values ranged from 0.63 for the clone TSH-565 to 0.82 for CC-10). The results demonstrate the potential for breeding more photosynthetically efficient cacao canopies.

Extending Murty interferometry to the Terahertz part of the spectrum

We discuss some novel technologies that enable the implementation of shearing interferometry at the terahertz part of the spectrum. Possible applications include the direct measurement of lens parameters, the measurement of refractive index of materials that are transparent to terahertz frequencies, determination of homogeneity of samples, measurement of optical distortions and the non-contact evaluation of thermal expansion coefficient of materials buried inside media that are opaque to optical or infrared frequencies but transparent to THz frequencies. The introduction of a shear to a Gaussian free-space propagating terahertz beam in a controlled manner also makes possible a range of new encoding and optical signal processing modalities.

The global distribution of phytoplankton size spectrum and size classes from their light-absorption spectra derived from satellite data

The absorption spectra of phytoplankton in the visible domain hold implicit information on the phytoplankton community structure. Here we use this information to retrieve quantitative information on phytoplankton size structure by developing a novel method to compute the exponent of an assumed power-law for their particle-size spectrum. This quantity, in combination with total chlorophyll-a concentration, can be used to estimate the fractional concentration of chlorophyll in any arbitrarily-defined size class of phytoplankton. We further define and derive expressions for two distinct measures of cell size of mixed populations, namely, the average spherical diameter of a bio-optically equivalent homogeneous population of cells of equal size, and the average equivalent spherical diameter of a population of cells that follow a power-law particle-size distribution. The method relies on measurements of two quantities of a phytoplankton sample: the concentration of chlorophyll-a, which is an operational index of phytoplankton biomass, and the total absorption coefficient of phytoplankton in the red peak of visible spectrum at 676 nm. A sensitivity analysis confirms that the relative errors in the estimates of the exponent of particle size spectra are reasonably low. The exponents of phytoplankton size spectra, estimated for a large set of in situ data from a variety of oceanic environments (~ 2400 samples), are within a reasonable range; and the estimated fractions of chlorophyll in pico-, nano- and micro-phytoplankton are generally consistent with those obtained by an independent, indirect method based on diagnostic pigments determined using high-performance liquid chromatography. The estimates of cell size for in situ samples dominated by different phytoplankton types (diatoms, prymnesiophytes, Prochlorococcus, other cyanobacteria and green algae) yield nominal sizes consistent with the taxonomic classification. To estimate the same quantities from satellite-derived ocean-colour data, we combine our method with algorithms for obtaining inherent optical properties from remote sensing. The spatial distribution of the size-spectrum exponent and the chlorophyll fractions of pico-, nano- and micro-phytoplankton estimated from satellite remote sensing are in agreement with the current understanding of the biogeography of phytoplankton functional types in the global oceans. This study contributes to our understanding of the distribution and time evolution of phytoplankton size structure in the global oceans.

Retrieval of phytoplankton size from bio-optical measurements: theory and applications

The absorption coefficient of a substance distributed as discrete particles in suspension is less than that of the same material dissolved uniformly in a medium—a phenomenon commonly referred to as the flattening effect. The decrease in the absorption coefficient owing to flattening effect depends on the concentration of the absorbing pigment inside the particle, the specific absorption coefficient of the pigment within the particle, and on the diameter of the particle, if the particles are assumed to be spherical. For phytoplankton cells in the ocean, with diameters ranging from less than 1 µm to more than 100 µm, the flattening effect is variable, and sometimes pronounced, as has been well documented in the literature. Here, we demonstrate how the in vivo absorption coefficient of phytoplankton cells per unit concentration of its major pigment, chlorophyll a, can be used to determine the average cell size of the phytoplankton population. Sensitivity analyses are carried out to evaluate the errors in the estimated diameter owing to potential errors in the model assumptions. Cell sizes computed for field samples using the model are compared qualitatively with indirect estimates of size classes derived from high performance liquid chromatography data. Also, the results are compared quantitatively against measurements of cell size in laboratory cultures. The method developed is easy-to-apply as an operational tool for in situ observations, and has the potential for application to remote sensing of ocean colour data.

Estimates of global multicomponent aerosol optical depth and direct radiative perturbation in the Laboratoire de Météorologie Dynamique general circulation model

The global cycle of multicomponent aerosols including sulfate, black carbon (BC),organic matter (OM), mineral dust, and sea salt is simulated in the Laboratoire de Me´te´orologie Dynamique general circulation model (LMDZT GCM). The seasonal open biomass burning emissions for simulation years 2000–2001 are scaled from climatological emissions in proportion to satellite detected fire counts. The emissions of dust and sea salt are parameterized online in the model. The comparison of model-predicted monthly mean aerosol optical depth (AOD) at 500 nm with Aerosol Robotic Network (AERONET) shows good agreement with a correlation coefficient of 0.57(N = 1324) and 76% of data points falling within a factor of 2 deviation. The correlation coefficient for daily mean values drops to 0.49 (N = 23,680). The absorption AOD (ta at 670 nm) estimated in the model is poorly correlated with measurements (r = 0.27, N = 349). It is biased low by 24% as compared to AERONET. The model reproduces the prominent features in the monthly mean AOD retrievals from Moderate Resolution Imaging Spectroradiometer (MODIS). The agreement between the model and MODIS is better over source and outflow regions (i.e., within a factor of 2).There is an underestimation of the model by up to a factor of 3 to 5 over some remote oceans. The largest contribution to global annual average AOD (0.12 at 550 nm) is from sulfate (0.043 or 35%), followed by sea salt (0.027 or 23%), dust (0.026 or 22%),OM (0.021 or 17%), and BC (0.004 or 3%). The atmospheric aerosol absorption is predominantly contributed by BC and is about 3% of the total AOD. The globally and annually averaged shortwave (SW) direct aerosol radiative perturbation (DARP) in clear-sky conditions is �2.17 Wm�2 and is about a factor of 2 larger than in all-sky conditions (�1.04 Wm�2). The net DARP (SW + LW) by all aerosols is �1.46 and �0.59 Wm�2 in clear- and all-sky conditions, respectively. Use of realistic, less absorbing in SW, optical properties for dust results in negative forcing over the dust-dominated regions.

Method of determining a geophysical-scale sea ice rheology from laboratory measurements

We present a methodology that allows a sea ice rheology, suitable for use in a General Circulation Model (GCM), to be determined from laboratory and tank experiments on sea ice when combined with a kinematic model of deformation. The laboratory experiments determine a material rheology for sea ice, and would investigate a nonlinear friction law of the form τ ∝ σ n⅔, instead of the more familiar Amonton's law, τ = μσn (τ is the shear stress, μ is the coefficient of friction and σ n is the normal stress). The modelling approach considers a representative region R containing ice floes (or floe aggregates), separated by flaws. The deformation of R is imposed and the motion of the floes determined using a kinematic model, which will be motivated from SAR observations. Deformation of the flaws is inferred from the floe motion and stress determined from the material rheology. The stress over R is then determined from the area-weighted contribution from flaws and floes

Ozonolysis of methyl oleate monolayers at the air–water interface: oxidation kinetics, reaction products and atmospheric implications

Ozonolysis of methyl oleate monolayers at the air–water interface results in surprisingly rapid loss of material through cleavage of the C[double bond, length as m-dash]C bond and evaporation/dissolution of reaction products. We determine using neutron reflectometry a rate coefficient of (5.7 ± 0.9) × 10−10 cm2 molecule−1 s−1 and an uptake coefficient of [similar]3 × 10−5 for the oxidation of a methyl ester monolayer: the atmospheric lifetime is [similar]10 min. We obtained direct experimental evidence that <2% of organic material remains at the surface on atmospheric timescales. Therefore known long atmospheric residence times of unsaturated fatty acids suggest that these molecules cannot be present at the interface throughout their ageing cycle, i.e. the reported atmospheric longevity is likely to be attributed to presence in the bulk and viscosity-limited reactive loss. Possible reaction products were characterized by ellipsometry and uncertainties in the atmospheric fate of organic surfactants such as oleic acid and its methyl ester are discussed. Our results suggest that a minor change to the structure of the molecule (fatty acid vs. its methyl ester) considerably impacts on reactivity and fate of the organic film.