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.

Evaluation of the sea ice simulation in a new coupled atmosphere-ocean climate model (HadGEM1)

[ 1] A rapid increase in the variety, quality, and quantity of observations in polar regions is leading to a significant improvement in the understanding of sea ice dynamic and thermodynamic processes and their representation in global climate models. We assess the simulation of sea ice in the new Hadley Centre Global Environmental Model (HadGEM1) against the latest available observations. The HadGEM1 sea ice component uses elastic-viscous-plastic dynamics, multiple ice thickness categories, and zero-layer thermodynamics. The model evaluation is focused on the mean state of the key variables of ice concentration, thickness, velocity, and albedo. The model shows good agreement with observational data sets. The variability of the ice forced by the North Atlantic Oscillation is also found to agree with observations.

Evaluation of the S(T) assimilation method with the Argo dataset

Recent observations from the Argo dataset of temperature and salinity profiles are used to evaluate a series of 3-year data assimilation experiments in a global ice–ocean general circulation model. The experiments are designed to evaluate a new data assimilation system whereby salinity is assimilated along isotherms, S(T ). In addition, the role of a balancing salinity increment to maintain water mass properties is investigated. This balancing increment is found to effectively prevent spurious mixing in tropical regions induced by univariate temperature assimilation, allowing the correction of isotherm geometries without adversely influencing temperature–salinity relationships. In addition, the balancing increment is able to correct a fresh bias associated with a weak subtropical gyre in the North Atlantic using only temperature observations. The S(T ) assimilation method is found to provide an important improvement over conventional depth level assimilation, with lower root-mean-squared forecast errors over the upper 500 m in the tropical Atlantic and Pacific Oceans. An additional set of experiments is performed whereby Argo data are withheld and used for independent evaluation. The most significant improvements from Argo assimilation are found in less well-observed regions (Indian, South Atlantic and South Pacific Oceans). When Argo salinity data are assimilated in addition to temperature, improvements to modelled temperature fields are obtained due to corrections to model density gradients and the resulting circulation. It is found that observations from the Argo array provide an invaluable tool for both correcting modelled water mass properties through data assimilation and for evaluating the assimilation methods themselves.

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.

An ex-ante evaluation of the economic impact of Bt cotton adoption by spanish farmers facing the EU cotton sector reform

Cotton production in the European Union (EU) is limited to areas of Greece and Southern Spain (Andalusia). The 2004 reform of the EU cotton policy severely affected the profitability of the crop. In this article we analyze how the introduction of genetically modified (GM), insect-resistant cotton varieties (Bt cotton) might help EU cotton farmers to increase profitability and therefore face the cotton policy reform. We first study farmers’ attitudes toward adoption of Bt cotton varieties through a survey conducted in Andalusia (Southern Spain). The results show a positive attitude of Andalusian cotton farmers toward the Bt cotton varieties. Second, we perform an ex-ante analysis of the effects of introducing Bt cotton in Andalusia. Finally, we integrate the analysis of the effects of Bt cotton with the analysis of the EU cotton reform. Our results show that despite the significant economic benefits of Bt cotton, the current policy reform is likely to jeopardize the profitability of cotton production in the EU.

An evaluation of the use of Cottonseed cake in the diet of growing pigs

An experiment was conducted to determine the effects of including cottonseed cake in rations for weaned growing pigs. Thirty-two Landrace x Large White pigs, weighing 20-24 kg, were included in four blocks formed on the basis of initial weight within sex in an otherwise completely randomized block design. The pigs were killed when they reached a live weight of 75.0 +/- 2.0 kg and the half careases were analysed into cuts and the weights of the organs were recorded. An estimate of the productivity of the pigs on each diet was calculated. Cottonseed cake reduced the voluntary feed intake (p < 0.001) and live weight gains (p < 0.001) and increased the heart, kidney and liver weights (p < 0.01). The pigs on the soya bean-based control diet took the shortest time to reach slaughter weight. The result was probably in part due to lysine deficiency and in part to the effect of free gossypol. It was found that it is at present cost-effective to include cottonseed cake in pig weaner grower diets up to 300 g/kg in Cameroon.

Statistical evaluation of the fixed concentration procedure for acute inhalation toxicity assessment

The conventional method for the assessment of acute inhalation toxicity (OECD Test Guideline 403, 1981) uses death of animals as an endpoint to identify the median lethal concentration (LC50). A new OECD Testing Guideline called the Fixed Concentration Procedure (FCP) is being prepared to provide an alternative to Test Guideline 403. Unlike Test Guideline 403, the FCP does not provide a point estimate of the LC50, but aims to identify an airborne exposure level that causes clear signs of nonlethal toxicity. This is then used to assign classification according to the new Globally Harmonized System of Classification and Labelling scheme (GHS). The FCP has been validated using statistical simulation rather than byin vivo testing. The statistical simulation approach predicts the GHS classification outcome and the numbers of deaths and animals used in the test for imaginary substances with a range of LC50 values and dose response curve slopes. This paper describes the FCP and reports the results from the statistical simulation study assessing its properties. It is shown that the procedure will be completed with considerably less death and suffering than Test Guideline 403, and will classify substances either in the same or a more stringent GHS class than that assigned on the basis of the LC50 value.

A statistical evaluation of the fixed dose procedure

The fixed-dose procedure (FDP) was introduced as OECD Test Guideline 420 in 1992, as an alternative to the conventional median lethal dose (LD50) test for the assessment of acute oral toxicity (OECD Test Guideline 401). The FDP uses fewer animals and causes less suffering than the conventional test, while providing information on the acute toxicity to allow substances to be ranked according to the EU hazard classification system. Recently the FDP has been revised, with the aim of providing further reductions and refinements, and classification according to the criteria of the Globally Harmonized Hazard Classification and Labelling scheme (GHS). This paper describes the revised FDP and analyses its properties, as determined by a statistical modelling approach. The analysis shows that the revised FDP classifies substances for acute oral toxicity generally in the same, or a more stringent, hazard class as that based on the LD50 value, according to either the GHS or the EU classification scheme. The likelihood of achieving the same classification is greatest for substances with a steep dose-response curve and median toxic dose (TD50) close to the LD50. The revised FDP usually requires five or six animals with two or fewer dying as a result of treatment in most cases.

An evaluation of the costs of making specific secondary metabolites: does the yield penalty incurred by host plant resistance to insects due to competition for resources?

The presumption that the synthesis of 'defence' compounds in plants must incur some 'trade-off' or penalty in terms of annual crop yields has been used to explain observed inverse correlations between resistance to herbivores and rates of growth or photosynthesis. An analysis of the cost of making secondary compounds suggests that this accounts for only a small part of the overall carbon budget of annual crop plants. Even the highest reported amounts of secondary metabolites found in different crop species (flavonoids, allylisothiocyanates, hydroxamic acids, 2-tridecanone) represent a carbon demand that can be satisfied by less than an hour's photosynthesis. Similar considerations apply to secondary compounds containing nitrogen or sulphur, which are unlikely to represent a major investment compared to the cost of making proteins, the major demand for these elements. Decreases in growth and photosynthesis in response to stress are more likely the result of programmed down-regulation. Observed correlations between yield and low contents of unpalatable or toxic compounds may be the result of parallel selection during the refinement of crop species by humans.

Evaluation of the cumulative evidence for freedom from BSE in birth cohorts

Substantial resources are used for surveillance of bovine spongiform encephalopathy (BSE) despite an extremely low detection rate, especially in healthy slaughtered cattle. We have developed a method based on the geometric waiting time distribution to establish and update the statistical evidence for BSE-freedom for defined birth cohorts using continued surveillance data. The results suggest that currently (data included till September 2004) a birth cohort of Danish cattle born after March 1999 is free from BSE with probability (power) of 0.8746 or 0.8509, depending on the choice of a model for the diagnostic sensitivity. These results apply to an assumed design prevalence of 1 in 10,000 and account for prevalence heterogeneity. The age-dependent, diagnostic sensitivity for the detection of BSE has been identified as major determinant of the power. The incorporation of heterogeneity was deemed adequate on scientific grounds and led to improved power values. We propose our model as a decision tool for possible future modification of the BSE surveillance and discuss public health and international trade implications.

Evaluation of the binding ability of a novel dioxatetraazamacrocyclic receptor that contains two phenanthroline units: selective uptake of carboxylate anions

The novel dioxatetraaza macrocycle [26]phen(2)N(4)O(2), which incorporates two phenanthroline units, has been synthesized, and its acid-base behavior has been evaluated by potentiometric and H-1 NMR methods. Six protonation constants were determined, and the protonation sequence was established by NMR. The location of the fifth proton on the phen nitrogen was confirmed by X-ray determinations of the crystal structures of the receptor as bromide and chloride salts. The two compounds have the general molecular formula {(H-5[26]phen(2)N(4)O(2))X-n(H2O)(5-n)}X(n-1)(.)mH(2)O, where X = Cl, n = 3, and m = 6 or X = Br, n = 4, and m = 5.5. In the solid state, the (H-5[26]phen(2)N(4)O(2))(5+) cation adopts a "horseshoe" topology with sufficient room to encapsulate three or four halogen anions through the several N-(HX)-X-... hydrogen-bonding interactions. Two supermolecules {(H-5[26]phen(2)N(4)O(2))X-n(H2O)(5-n)}((5-n)+) form an interpenetrating dimeric species, which was also found by ESI mass spectrum. Binding studies of the protonated macrocycle with aliphatic (ox(2-), mal(2-), suc(2-), cit(3-), cta(3-)) and aromatic (bzc(-), naphc(-), anthc(-), pyrc(-), ph(2-), iph(2-), tph(2-), btc(3-)) anions were determined in water by potentiometric methods. These studies were complemented by H-1 NMR titrations in D2O of the receptor with selected anions. The H-i[26]phen(2)N(4)O(2)(i+) receptor can selectively uptake highly charged or extended aromatic carboxylate anions, such as btc(3-) and pyrc(-), in the pH ranges of 4.0-8.5 and < 4.0, respectively, from aqueous solution that contain the remaining anions as pollutants or contaminants. To obtain further insight into these structural and experimental findings, molecular dynamics (MD) simulations were carried out in water solution.