On the evolution of the solar wind between 1 and 5 AU at the time of the Cassini Jupiter flyby: Multispacecraft observations of interplanetary coronal mass ejections including the formation of a merged interaction region

The Cassini flyby of Jupiter occurred at a time near solar maximum. Consequently, the pre-Jupiter data set reveals clear and numerous transient perturbations to the Parker Spiral solar wind structure. Limited plasma data are available at Cassini for this period due to pointing restrictions imposed on the instrument. This renders the identification of the nature of such structures ambiguous, as determinations based on the magnetic field data alone are unreliable. However, a fortuitous alignment of the planets during this encounter allowed us to trace these structures back to those observed previously by the Wind spacecraft near the Earth. Of the phenomena that we are satisfactorily able to trace back to their manifestation at 1 AU, two are identified as being due to interplanetary coronal mass ejections. One event at Cassini is shown to be a merged interaction region, which is formed from the compression of a magnetic cloud by two anomalously fast solar wind streams. The flux-rope structure associated with this magnetic cloud is not as apparent at Cassini and has most likely been compressed and deformed. Confirmation of the validity of the ballistic projections used here is provided by results obtained from a one-dimensional magnetohydrodynamic projection of solar wind parameters measured upstream near the Earth. It is found that when the Earth and Cassini are within a few tens of degrees in heliospheric longitude, the results of this one-dimensional model predict the actual conditions measured at 5 AU to an impressive degree. Finally, the validity of the use of such one-dimensional projections in obtaining quasi-solar wind parameters at the outer planets is discussed.

Global scale water stages from SAR imagery (GlosSARi) to support global flood forecasting

Despite the success of studies attempting to integrate remotely sensed data and flood modelling and the need to provide near-real time data routinely on a global scale as well as setting up online data archives, there is to date a lack of spatially and temporally distributed hydraulic parameters to support ongoing efforts in modelling. Therefore, the objective of this project is to provide a global evaluation and benchmark data set of floodplain water stages with uncertainties and assimilation in a large scale flood model using space-borne radar imagery. An algorithm is developed for automated retrieval of water stages with uncertainties from a sequence of radar imagery and data are assimilated in a flood model using the Tewkesbury 2007 flood event as a feasibility study. The retrieval method that we employ is based on possibility theory which is an extension of fuzzy sets and that encompasses probability theory. In our case we first attempt to identify main sources of uncertainty in the retrieval of water stages from radar imagery for which we define physically meaningful ranges of parameter values. Possibilities of values are then computed for each parameter using a triangular ‘membership’ function. This procedure allows the computation of possible values of water stages at maximum flood extents along a river at many different locations. At a later stage in the project these data are then used in assimilation, calibration or validation of a flood model. The application is subsequently extended to a global scale using wide swath radar imagery and a simple global flood forecasting model thereby providing improved river discharge estimates to update the latter.

Retrieval of effective ice crystal size in the infrared: sensitivity study and global measurements from TIROS-N Operational Vertical Sounder

The difference between cirrus emissivities at 8 and 11 μm is sensitive to the mean effective ice crystal size of the cirrus cloud, De. By using single scattering properties of ice crystals shaped as planar polycrystals, diameters of up to about 70 μm can be retrieved, instead of up to 45 μm assuming spheres or hexagonal columns. The method described in this article is used for a global determination of mean effective ice crystal sizes of cirrus clouds from TOVS satellite observations. A sensitivity study of the De retrieval to uncertainties in hypotheses on ice crystal shape, size distributions, and temperature profiles, as well as in vertical and horizontal cloud heterogeneities shows that uncertainties can be as large as 30%. However, the TOVS data set is one of few data sets which provides global and long-term coverage. Having analyzed the years 1987–1991, it was found that measured effective ice crystal diameters De are stable from year to year. For 1990 a global median De of 53.5 μm was determined. Averages distinguishing ocean/land, season, and latitude lie between 23 μm in winter over Northern Hemisphere midlatitude land and 64 μm in the tropics. In general, larger Des are found in regions with higher atmospheric water vapor and for cirrus with a smaller effective emissivity.

Lagrangian validation of numerical drifter trajectories using drifting buoys: Application to the Agulhas system

The skill of numerical Lagrangian drifter trajectories in three numerical models is assessed by comparing these numerically obtained paths to the trajectories of drifting buoys in the real ocean. The skill assessment is performed using the two-sample Kolmogorov–Smirnov statistical test. To demonstrate the assessment procedure, it is applied to three different models of the Agulhas region. The test can either be performed using crossing positions of one-dimensional sections in order to test model performance in specific locations, or using the total two-dimensional data set of trajectories. The test yields four quantities: a binary decision of model skill, a confidence level which can be used as a measure of goodness-of-fit of the model, a test statistic which can be used to determine the sensitivity of the confidence level, and cumulative distribution functions that aid in the qualitative analysis. The ordering of models by their confidence levels is the same as the ordering based on the qualitative analysis, which suggests that the method is suited for model validation. Only one of the three models, a 1/10° two-way nested regional ocean model, might have skill in the Agulhas region. The other two models, a 1/2° global model and a 1/8° assimilative model, might have skill only on some sections in the region

ASYM20: a program for force constant and normal coordinate calculations, with a critical review of the theory involved

The theory of harmonic force constant refinement calculations is reviewed, and a general-purpose program for force constant and normal coordinate calculations is described. The program, called ASYM20. is available through Quantum Chemistry Program Exchange. It will work on molecules of any symmetry containing up to 20 atoms and will produce results on a series of isotopomers as desired. The vibrational secular equations are solved in either nonredundant valence internal coordinates or symmetry coordinates. As well as calculating the (harmonic) vibrational wavenumbers and normal coordinates, the program will calculate centrifugal distortion constants, Coriolis zeta constants, harmonic contributions to the α′s. root-mean-square amplitudes of vibration, and other quantities related to gas electron-diffraction studies and thermodynamic properties. The program will work in either a predict mode, in which it calculates results from an input force field, or in a refine mode, in which it refines an input force field by least squares to fit observed data on the quantities mentioned above. Predicate values of the force constants may be included in the data set for a least-squares refinement. The program is written in FORTRAN for use on a PC or a mainframe computer. Operation is mainly controlled by steering indices in the input data file, but some interactive control is also implemented.

Bayesian estimation of willingness-to-pay where respondents mis-report their preferences

We introduce a modified conditional logit model that takes account of uncertainty associated with mis-reporting in revealed preference experiments estimating willingness-to-pay (WTP). Like Hausman et al. [Journal of Econometrics (1988) Vol. 87, pp. 239-269], our model captures the extent and direction of uncertainty by respondents. Using a Bayesian methodology, we apply our model to a choice modelling (CM) data set examining UK consumer preferences for non-pesticide food. We compare the results of our model with the Hausman model. WTP estimates are produced for different groups of consumers and we find that modified estimates of WTP, that take account of mis-reporting, are substantially revised downwards. We find a significant proportion of respondents mis-reporting in favour of the non-pesticide option. Finally, with this data set, Bayes factors suggest that our model is preferred to the Hausman model.

Modifying willingness to pay estimates where respondents mis-report their preferences

The likelihood for the Logit model is modified, so as to take account of uncertainty associated with mis-reporting in stated preference experiments estimating willingness to pay (WTP). Monte Carlo results demonstrate the bias imparted to estimates where there is mis-reporting. The approach is applied to a data set examining consumer preferences for food produced employing a nonpesticide technology. Our modified approach leads to WTP that are substantially downwardly revised.

Design and optimisation of a large-area process-based model for annual crops

The formulation of a new process-based crop model, the general large-area model (GLAM) for annual crops is presented. The model has been designed to operate on spatial scales commensurate with those of global and regional climate models. It aims to simulate the impact of climate on crop yield. Procedures for model parameter determination and optimisation are described, and demonstrated for the prediction of groundnut (i.e. peanut; Arachis hypogaea L.) yields across India for the period 1966-1989. Optimal parameters (e.g. extinction coefficient, transpiration efficiency, rate of change of harvest index) were stable over space and time, provided the estimate of the yield technology trend was based on the full 24-year period. The model has two location-specific parameters, the planting date, and the yield gap parameter. The latter varies spatially and is determined by calibration. The optimal value varies slightly when different input data are used. The model was tested using a historical data set on a 2.5degrees x 2.5degrees grid to simulate yields. Three sites are examined in detail-grid cells from Gujarat in the west, Andhra Pradesh towards the south, and Uttar Pradesh in the north. Agreement between observed and modelled yield was variable, with correlation coefficients of 0.74, 0.42 and 0, respectively. Skill was highest where the climate signal was greatest, and correlations were comparable to or greater than correlations with seasonal mean rainfall. Yields from all 35 cells were aggregated to simulate all-India yield. The correlation coefficient between observed and simulated yields was 0.76, and the root mean square error was 8.4% of the mean yield. The model can be easily extended to any annual crop for the investigation of the impacts of climate variability (or change) on crop yield over large areas. (C) 2004 Elsevier B.V. All rights reserved.

Climate change and the flowering time of annual crops

Crop production is inherently sensitive to variability in climate. Temperature is a major determinant of the rate of plant development and, under climate change, warmer temperatures that shorten development stages of determinate crops will most probably reduce the yield of a given variety. Earlier crop flowering and maturity have been observed and documented in recent decades, and these are often associated with warmer (spring) temperatures. However, farm management practices have also changed and the attribution of observed changes in phenology to climate change per se is difficult. Increases in atmospheric [CO2] often advance the time of flowering by a few days, but measurements in FACE (free air CO2 enrichment) field-based experiments suggest that elevated [CO2] has little or no effect on the rate of development other than small advances in development associated with a warmer canopy temperature. The rate of development (inverse of the duration from sowing to flowering) is largely determined by responses to temperature and photoperiod, and the effects of temperature and of photoperiod at optimum and suboptimum temperatures can be quantified and predicted. However, responses to temperature, and more particularly photoperiod, at supraoptimal temperature are not well understood. Analysis of a comprehensive data set of time to tassel initiation in maize (Zea mays) with a wide range of photoperiods above and below the optimum suggests that photoperiod modulates the negative effects of temperature above the optimum. A simulation analysis of the effects of prescribed increases in temperature (0-6 degrees C in + 1 degrees C steps) and temperature variability (0% and + 50%) on days to tassel initiation showed that tassel initiation occurs later, and variability was increased, as the temperature exceeds the optimum in models both with and without photoperiod sensitivity. However, the inclusion of photoperiod sensitivity above the optimum temperature resulted in a higher apparent optimum temperature and less variability in the time of tassel initiation. Given the importance of changes in plant development for crop yield under climate change, the effects of photoperiod and temperature on development rates above the optimum temperature clearly merit further research, and some of the knowledge gaps are identified herein.

Emergence of Chenopodium album and Stellaria media of different origins under different climatic conditions

The emergence behaviour of weed species in relation to cultural and meteorological events was studied. Dissimilarities between populations in dormancy and germination ecology, between-year maturation conditions and seed quality and burial site climate all contribute to potentially unpredictable variability. Therefore, a weed emergence data set was produced for weed seeds of Stellaria media and Chenopodium album matured and collected from three populations (Italy, Sweden and UK). The seeds were collected in two consecutive seasons (1999 and 2000) and subsequently buried in the autumn of the same year of maturation in eight contrasting climatic locations throughout Europe and the USA. The experiment sought to explore and explain differences between the three populations in their emergence behaviour. Evidence was demonstrated of synchrony in the timing of the emergence of different populations of a species at a given burial site. The relative magnitudes of emergence from the three populations at a given burial site in a given year were generally similar across all the burial sites in the study. The resulting data set was also used to construct a simple weed emergence model, which was tested for its application to the range of different burial environments and populations. The study demonstrated the possibility of using a simple thermal time-based model to describe part of the emergence behaviour across different burial sites, seed populations and seasons, and a simple winter chilling relationship to adjust for the magnitude of the flush of emergence at a given burial site. This study demonstrates the possibility of developing robust generic models for simple predictions of emergence timing across populations.

Genetic parameters and evaluations from single- and multiple-trait analysis of dairy cow fertility and milk production

Genetic parameters and breeding values for dairy cow fertility were estimated from 62 443 lactation records. Two-trait analysis of fertility and milk yield was investigated as a method to estimate fertility breeding values when culling or selection based on milk yield in early lactation determines presence or absence of fertility observations in later lactations. Fertility traits were calving interval, intervals from calving to first service, calving to conception and first to last service, conception success to first service and number of services per conception. Milk production traits were 305-day milk, fat and protein yield. For fertility traits, range of estimates of heritability (h(2)) was 0.012 to 0.028 and of permanent environmental variance (c(2)) was 0.016 to 0.032. Genetic correlations (r(g)) among fertility traits were generally high ( > 0.70). Genetic correlations of fertility with milk production traits were unfavourable (range -0.11 to 0.46). Single and two-trait analyses of fertility were compared using the same data set. The estimates of h(2) and c(2) were similar for two types of analyses. However, there were differences between estimated breeding values and rankings for the same trait from single versus multi-trait analyses. The range for rank correlation was 0.69-0.83 for all animals in the pedigree and 0.89-0.96 for sires with more than 25 daughters. As single-trait method is biased due to selection on milk yield, a multi-trait evaluation of fertility with milk yield is recommended. (C) 2002 Elsevier Science B.V. All rights reserved.

An evaluation of different growth functions for describing the profile of live weight with time (age) in meat and egg strains of chicken

A total of 86 profiles from meat and egg strains of chickens (male and female) were used in this study. Different flexible growth functions were evaluated with regard to their ability to describe the relationship between live weight and age and were compared with the Gompertz and logistic equations, which have a fixed point of inflection. Six growth functions were used: Gompertz, logistic, Lopez, Richards, France, and von Bertalanffy. A comparative analysis was carried out based on model behavior and statistical performance. The results of this study confirmed the initial concern about the limitation of a fixed point of inflection, such as in the Gompertz equation. Therefore, consideration of flexible growth functions as an alternatives to the simpler equations (with a fixed point of inflection) for describing the relationship between live weight and age are recommended for the following reasons: they are easy to fit, they very often give a closer fit to data points because of their flexibility and therefore a smaller RSS value, than the simpler models, and they encompasses simpler models for the addition of an extra parameter, which is especially important when the behavior of a particular data set is not defined previously.

Effects of phosphorus intake on phosphorus flow in growing pigs: Application and comparison of two models

A comparison of the models of Vitti et al. (2000, J. Anim. Sci. 78, 2706-2712) and Fernandez (1995c, Livest. Prod. Sci. 41, 255-261) was carried out using two data sets on growing pigs as input. The two models compared were based on similar basic principles, although their aims and calculations differed. The Vitti model employs the rate:state formalism and describes phosphorus (P) flow between four pools representing P content in gut, blood, bone and soft tissue in growing goats. The Fernandez model describes flow and fractional recirculation between P pools in gut, blood and bone in growing pigs. The results from both models showed similar trends for P absorption from gut to blood and net retention in bone with increasing P intake, with the exception of the 65 kg results from Date Set 2 calculated using the FernAndez model. Endogenous loss from blood back to gut increased faster with increasing P intake in the FernAndez than in the Vitti model for Data Set 1. However, for Data Set 2, endogenous loss increased with increasing P intake using the Vitti model, but decreased when calculated using the FernAndez model. Incorporation of P into bone was not influenced by intake in the FernAndez model, while in the Vitti model there was an increasing trend. The FernAndez model produced a pattern of decreasing resorption in bone with increasing P intake, with one of the data sets, which was not observed when using the Vitti model. The pigs maintained their P homeostasis in blood by regulation of P excretion in urine. (c) 2005 Elsevier Ltd. All rights reserved.

A predictive model for early-warning of Septoria leaf blotch on winter wheat

Disease-weather relationships influencing Septoria leaf blotch (SLB) preceding growth stage (GS) 31 were identified using data from 12 sites in the UK covering 8 years. Based on these relationships, an early-warning predictive model for SLB on winter wheat was formulated to predict the occurrence of a damaging epidemic (defined as disease severity of 5% or > 5% on the top three leaf layers). The final model was based on accumulated rain > 3 mm in the 80-day period preceding GS 31 (roughly from early-February to the end of April) and accumulated minimum temperature with a 0A degrees C base in the 50-day period starting from 120 days preceding GS 31 (approximately January and February). The model was validated on an independent data set on which the prediction accuracy was influenced by cultivar resistance. Over all observations, the model had a true positive proportion of 0.61, a true negative proportion of 0.73, a sensitivity of 0.83, and a specificity of 0.18. True negative proportion increased to 0.85 for resistant cultivars and decreased to 0.50 for susceptible cultivars. Potential fungicide savings are most likely to be made with resistant cultivars, but such benefits would need to be identified with an in-depth evaluation.

Sperm competition and the evolution of testes size in terrestrial mammalian carnivores

Summary Understanding the factors influencing variation in the degree of sperm competition is a key question underlying the mechanisms driving sexual conflict. Previous behavioural and comparative studies have indicated that carnivores appear to have evolved under sperm competition but an analysis of the predictors of the level of sperm competition is missing. In this study, we use phylogenetic comparative methods to investigate life-history parameters predicted to affect the degree of sperm competition in terrestrial carnivores using variation in relative testes size (RTS, after controlling for body size allometry) as a measure of the level of sperm competition. Due to a paucity of consistent data across taxa, we used three measures of RTS: testes mass (n = 40 species), testes and epididymes mass combined (n = 38), and testes volume (n = 48). We also created a derived data set (n = 79) with testes mass estimated from regression analyses on the other measures of testes size. Carnivores with shorter mating seasons had relatively larger testes, consistent with the hypothesis that sperm competition is greater when the degree of female oestrous synchrony is high. This relationship was stronger in spontaneous versus induced ovulators, suggesting higher sperm competition levels in spontaneous ovulators. This is the first comparative study to show this within mammalian taxa. Neither social mating system nor reproductive lifespan were significantly associated with variation in RTS and hence are poor predictors of sperm competition levels. None of the above relationships were found to be significant for the testes and epididymes mass combined data set, but our understanding of the role of the epididymis in sperm competition is too limited to draw any conclusions. Finally, we consistently found a significant phylogenetic signal in all analyses, indicating that phylogeny has played a significant role in the evolution of carnivore testes size and, therefore, in shaping levels of sperm competition. Our results shed new light into the factors affecting levels of sperm competition in terrestrial carnivores by showing that the degree of oestrous synchrony and ovulation type interact to predict variation in RTS.