Parameter estimation using a particle method: inferring mixing coefficients from sea level observations

This paper presents a first attempt to estimate mixing parameters from sea level observations using a particle method based on importance sampling. The method is applied to an ensemble of 128 members of model simulations with a global ocean general circulation model of high complexity. Idealized twin experiments demonstrate that the method is able to accurately reconstruct mixing parameters from an observed mean sea level field when mixing is assumed to be spatially homogeneous. An experiment with inhomogeneous eddy coefficients fails because of the limited ensemble size. This is overcome by the introduction of local weighting, which is able to capture spatial variations in mixing qualitatively. As the sensitivity of sea level for variations in mixing is higher for low values of mixing coefficients, the method works relatively well in regions of low eddy activity.

Decadal-scale rainfall variability in Ethiopia recorded in an annually laminated, Holocene-age, stalagmite

An annually laminated, uranium-series dated, Holocene stalagmite from southeast Ethiopia has been analysed for growth rate and δ13C and δ18O variations at annual to biennial resolution, in order to provide the first long duration proxy record of decadal-scale rainfall variability in this climatically sensitive region. Our study site (10°N) is climatically influenced by both summer (June—August) and spring (March—May) rainfall caused by the annual movement of the Inter-Tropical Convergence Zone (ITCZ) and modulated by large-scale anomalies in the atmospheric circulation and in ocean temperatures. Here we show that stalagmite growth, episodic throughout the last 7800 years, demonstrates decadal-scale (8—25 yr) variability in both growth rate and δ 18O. A hydrological model was employed and indicates that this decadal variability is due to variations in the relative amounts of rainfall in the two rain seasons. Our record, unique in its combination of length (a total of ~1000 years), annual chronology and high resolution δ18O, shows for the first time that such decadal-scale variability in rainfall in this region has occurred through the Holocene, which implies persistent decadal-scale variability for the large-scale atmospheric and oceanic driving factors.

Determining the influence of rainfall patterns and carbendazim on the surface activity of the earthworm Lumbricus terrestris

Carbendazim is highly toxic to earthworms and is used as a standard control substance when running field-based trials of pesticides, but results using carbendazim are highly variable. In the present study, impacts of timing of rainfall events following carbendazim application on earthworms were investigated. Lumbricus terrestris were maintained in soil columns to which carbendazim and then deionized water (a rainfall substitute) were applied. Carbendazim was applied at 4 kg/ha, the rate recommended in pesticide field trials. Three rainfall regimes were investigated: initial and delayed heavy rainfall 24 h and 6 d after carbendazim application, and frequent rainfall every 48 h. Earthworm mortality and movement of carbendazim through the soil was assessed 14 d after carbendazim application. No detectable movement of carbendazim occurred through the soil in any of the treatments or controls. Mortality in the initial heavy and frequent rainfall was significantly higher (approximately 55%) than in the delayed rainfall treatment (approximately 25%). This was due to reduced bioavailability of carbendazim in the latter treatment due to a prolonged period of sorption of carbendazim to soil particles before rainfall events. The impact of carbendazim application on earthworm surface activity was assessed using video cameras. Carbendazim applications significantly reduced surface activity due to avoidance behavior of the earthworms. Surface activity reductions were least in the delayed rainfall treatment due to the reduced bioavailability of the carbendazim. The nature of rainfall events' impacts on the response of earthworms to carbendazim applications, and details of rainfall events preceding and following applications during field trials should be made at a higher level of resolution than is currently practiced according to standard International Organization for Standardization protocols.

Maize yield and rainfall on different spatial and temporal scales in Southern Brazil

This study aimed to establish relationships between maize yield and rainfall on different temporal and spatial scales, in order to provide a basis for crop monitoring and modelling. A 16-year series of maize yield and daily rainfall from 11 municipalities and micro-regions of Rio Grande do Sul State was used. Correlation and regression analyses were used to determine associations between crop yield and rainfall for the entire crop cycle, from tasseling to 30 days after, and from 5 days before tasseling to 40 days after. Close relationships between maize yield and rainfall were found, particularly during the reproductive period (45-day period comprising the flowering and grain filling). Relationships were closer on a regional scale than at smaller scales. Implications of the crop-rainfall relationships for crop modelling are discussed.

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.

Effects of summer rainfall manipulations on the abundance and vertical distribution of herbivorous soil macro-invertebrates

Soil invertebrate communities are likely to be highly vulnerable to low soil moisture, caused by a reduction in summer rainfall which is predicted for some regions under current climate change scenarios. However, the effects of changes in summer rainfall on soil invertebrate assemblages have rarely been tested experimentally. In this study, samples were taken in 2003 and 2004 from a long-running field experiment, to investigate the impact of 10 years of experimental summer drought and increased summer rainfall manipulations on the soil fauna of a calcareous grassland. Summer drought altered the soil invertebrate assemblage in the autumn, immediately following treatment application, but by the following spring treatment effects were no longer apparent. The two most common root herbivore species responded differently to the summer rainfall manipulations. Larvae of the dominant root-chewing species, Agriotes lineatus, were more numerous under enhanced rainfall in both the spring and autumn. In contrast, abundance of the Coccoidea Lecanopsis formicarum was unaffected by the rainfall manipulations. The responses of root herbivores to an increased incidence of summer droughts are therefore likely to vary, depending on their feeding strategy and life history. (c) 2007 Elsevier Masson SAS. All rights reserved.

Adaptation to climate change and variability: farmer responses to intra-seasonal precipitation trends in South Africa

We describe the nature of recent (50 year) rainfall variability in the summer rainfall zone, South Africa, and how variability is recognised and responded to on the ground by farmers. Using daily rainfall data and self-organising mapping (SOM) we identify 12 internally homogeneous rainfall regions displaying differing parameters of precipitation change. Three regions, characterised by changing onset and timing of rains, rainfall frequencies and intensities, in Limpopo, North West and KwaZulu Natal provinces, were selected to investigate farmer perceptions of, and responses to, rainfall parameter changes. Village and household level analyses demonstrate that the trends and variabilities in precipitation parameters differentiated by the SOM analysis were clearly recognised by people living in the areas in which they occurred. A range of specific coping and adaptation strategies are employed by farmers to respond to climate shifts, some generic across regions and some facilitated by specific local factors. The study has begun to understand the complexity of coping and adaptation, and the factors that influence the decisions that are taken.

The effects of environmental perturbation and measurement error on estimates of the shape parameter in the theta-logistic model of population regulation

The theta-logistic is a widely used generalisation of the logistic model of regulated biological processes which is used in particular to model population regulation. Then the parameter theta gives the shape of the relationship between per-capita population growth rate and population size. Estimation of theta from population counts is however subject to bias, particularly when there are measurement errors. Here we identify factors disposing towards accurate estimation of theta by simulation of populations regulated according to the theta-logistic model. Factors investigated were measurement error, environmental perturbation and length of time series. Large measurement errors bias estimates of theta towards zero. Where estimated theta is close to zero, the estimated annual return rate may help resolve whether this is due to bias. Environmental perturbations help yield unbiased estimates of theta. Where environmental perturbations are large, estimates of theta are likely to be reliable even when measurement errors are also large. By contrast where the environment is relatively constant, unbiased estimates of theta can only be obtained if populations are counted precisely Our results have practical conclusions for the design of long-term population surveys. Estimation of the precision of population counts would be valuable, and could be achieved in practice by repeating counts in at least some years. Increasing the length of time series beyond ten or 20 years yields only small benefits. if populations are measured with appropriate accuracy, given the level of environmental perturbation, unbiased estimates can be obtained from relatively short censuses. These conclusions are optimistic for estimation of theta. (C) 2008 Elsevier B.V All rights reserved.

A swarm pattern transformation method based on macroscopic parameter operations

The work reported in this paper is motivated by biomimetic inspiration - the transformation of patterns. The major issue addressed is the development of feasible methods for transformation based on a macroscopic tool. The general requirement for the feasibility of the transformation method is determined by classifying pattern formation approaches an their characteristics. A formal definition for pattern transformation is provided and four special cases namely, elementary and geometric transformation based on repositioning all and some robotic agents are introduced. A feasible method for transforming patterns geometrically, based on the macroscopic parameter operation of a swarm is considered. The transformation method is applied to a swarm model which lends itself to the transformation technique. Simulation studies are developed to validate the feasibility of the approach, and do indeed confirm the approach.

Riemann solvers with primitive parameter vectors for two-dimensional compressible flows of a real gas

An analysis of averaging procedures is presented for an approximate Riemann solver for the equations governing the compressible flow of a real gas. This study extends earlier work for the Euler equations with ideal gases.

Geostatistical analysis of rainfall

Rainfall can be modeled as a spatially correlated random field superimposed on a background mean value; therefore, geostatistical methods are appropriate for the analysis of rain gauge data. Nevertheless, there are certain typical features of these data that must be taken into account to produce useful results, including the generally non-Gaussian mixed distribution, the inhomogeneity and low density of observations, and the temporal and spatial variability of spatial correlation patterns. Many studies show that rigorous geostatistical analysis performs better than other available interpolation techniques for rain gauge data. Important elements are the use of climatological variograms and the appropriate treatment of rainy and nonrainy areas. Benefits of geostatistical analysis for rainfall include ease of estimating areal averages, estimation of uncertainties, and the possibility of using secondary information (e.g., topography). Geostatistical analysis also facilitates the generation of ensembles of rainfall fields that are consistent with a given set of observations, allowing for a more realistic exploration of errors and their propagation in downstream models, such as those used for agricultural or hydrological forecasting. This article provides a review of geostatistical methods used for kriging, exemplified where appropriate by daily rain gauge data from Ethiopia.