Conditioning model output statistics of regional climate model precipitation on circulation patterns

Dynamical downscaling of Global Climate Models (GCMs) through regional climate models (RCMs) potentially improves the usability of the output for hydrological impact studies. However, a further downscaling or interpolation of precipitation from RCMs is often needed to match the precipitation characteristics at the local scale. This study analysed three Model Output Statistics (MOS) techniques to adjust RCM precipitation; (1) a simple direct method (DM), (2) quantile-quantile mapping (QM) and (3) a distribution-based scaling (DBS) approach. The modelled precipitation was daily means from 16 RCMs driven by ERA40 reanalysis data over the 1961–2000 provided by the ENSEMBLES (ENSEMBLE-based Predictions of Climate Changes and their Impacts) project over a small catchment located in the Midlands, UK. All methods were conditioned on the entire time series, separate months and using an objective classification of Lamb's weather types. The performance of the MOS techniques were assessed regarding temporal and spatial characteristics of the precipitation fields, as well as modelled runoff using the HBV rainfall-runoff model. The results indicate that the DBS conditioned on classification patterns performed better than the other methods, however an ensemble approach in terms of both climate models and downscaling methods is recommended to account for uncertainties in the MOS methods.

Pre-conditioning the epigenetic response to high vapor pressure deficit increases the drought tolerance of Arabidopsis thaliana

Epigenetic modification of the genome via cytosine methylation is a dynamic process that responds to changes in the growing environment. This modification can also be heritable. The combination of both properties means that there is the potential for the life experiences of the parental generation to modify the methylation profiles of their offspring and so potentially to ‘pre-condition’ them to better accommodate abiotic conditions encountered by their parents. We recently identified high vapor pressure deficit (vpd)-induced DNA methylation at two gene loci in the stomatal development pathway and an associated reduction in leaf stomatal frequency.1 Here, we test whether this epigenetic modification pre-conditioned parents and their offspring to the more severe water stress of periodic drought. We found that three generations of high vpd-grown plants were better able to withstand periodic drought stress over two generations. This resistance was not directly associated with de novo methylation of the target stomata genes, but was associated with the cmt3 mutant’s inability to maintain asymmetric sequence context methylation. If our finding applies widely, it could have significant implications for evolutionary biology and breeding for stressful environments.

Long-term conditioning of soil by plantation eucalypts and pines does not affect growth of the native jarrah tree

Plant species can condition the physico-chemical and biological properties of soil in ways that modify plant growth via plant–soil feedback (PSF). Plant growth can be positively affected, negatively affected or neutrally affected by soil conditioning by the same or other plant species. Soil conditioning by other plant species has particular relevance to ecological restoration of historic ecosystems because sites set aside for restoration are often conditioned by other, potentially non-native, plant species. We investigated changes in properties of jarrah forest soils after long-term (35 years) conditioning by pines (Pinus radiata), Sydney blue gums (Eucalyptus saligna), both non-native, plantation trees, and jarrah (Eucalyptus marginata; dominant native tree). Then, we tested the influence of the conditioned soils on the growth of jarrah seedlings. Blue gums and pines similarly conditioned the physico-chemical properties of soils, which differed from soil conditioning caused by jarrah. Especially important were the differences in conditioning of the properties C:N ratio, pH, and available K. The two eucalypt species similarly conditioned the biological properties of soil (i.e. community level physiological profile, numbers of fungal-feeding nematodes, omnivorous nematodes, and nematode channel ratio), and these differed from conditioning caused by pines. Species-specific conditioning of soil did not translate into differences in the amounts of biomass produced by jarrah seedlings and a neutral PSF was observed. In summary, we found that decades of soil conditioning by non-native plantation trees did not influence the growth of jarrah seedlings and will therefore not limit restoration of jarrah following the removal of the plantation trees.

Soil conditioning and plant-soil feedbacks in a modified forest ecosystem are soil-context dependent

Aims There is potential for altered plant-soil feedback (PSF) to develop in human-modified ecosystems but empirical data to test this idea are limited. Here, we compared the PSF operating in jarrah forest soil restored after bauxite mining in Western Australia with that operating in unmined soil. Methods Native seedlings of jarrah (Eucalyptus marginata), acacia (Acacia pulchella), and bossiaea (Bossiaea ornata) were grown in unmined and restored soils to measure conditioning of chemical and biological properties as compared with unplanted control soils. Subsequently, acacia and bossiaea were grown in soils conditioned by their own or by jarrah seedlings to determine the net PSF. Results In unmined soil, the three plant species conditioned the chemical properties but had little effect on the biological properties. In comparison, jarrah and bossiaea conditioned different properties of restored soil while acacia did not condition this soil. In unmined soil, neutral PSF was observed, whereas in restored soil, negative PSF was associated with acacia and bossiaea. Conclusions Soil conditioning was influenced by soil context and plant species. The net PSF was influenced by soil context, not by plant species and it was different in restored and unmined soils. The results have practical implications for ecosystem restoration after human activities.

Conditioning of the weak-constraint variational data assimilation problem for numerical weather prediction

4-Dimensional Variational Data Assimilation (4DVAR) assimilates observations through the minimisation of a least-squares objective function, which is constrained by the model flow. We refer to 4DVAR as strong-constraint 4DVAR (sc4DVAR) in this thesis as it assumes the model is perfect. Relaxing this assumption gives rise to weak-constraint 4DVAR (wc4DVAR), leading to a different minimisation problem with more degrees of freedom. We consider two wc4DVAR formulations in this thesis, the model error formulation and state estimation formulation. The 4DVAR objective function is traditionally solved using gradient-based iterative methods. The principle method used in Numerical Weather Prediction today is the Gauss-Newton approach. This method introduces a linearised `inner-loop' objective function, which upon convergence, updates the solution of the non-linear `outer-loop' objective function. This requires many evaluations of the objective function and its gradient, which emphasises the importance of the Hessian. The eigenvalues and eigenvectors of the Hessian provide insight into the degree of convexity of the objective function, while also indicating the difficulty one may encounter while iterative solving 4DVAR. The condition number of the Hessian is an appropriate measure for the sensitivity of the problem to input data. The condition number can also indicate the rate of convergence and solution accuracy of the minimisation algorithm. This thesis investigates the sensitivity of the solution process minimising both wc4DVAR objective functions to the internal assimilation parameters composing the problem. We gain insight into these sensitivities by bounding the condition number of the Hessians of both objective functions. We also precondition the model error objective function and show improved convergence. We show that both formulations' sensitivities are related to error variance balance, assimilation window length and correlation length-scales using the bounds. We further demonstrate this through numerical experiments on the condition number and data assimilation experiments using linear and non-linear chaotic toy models.

Conditioning and preconditioning of the optimal state estimation problem

Optimal state estimation is a method that requires minimising a weighted, nonlinear, least-squares objective function in order to obtain the best estimate of the current state of a dynamical system. Often the minimisation is non-trivial due to the large scale of the problem, the relative sparsity of the observations and the nonlinearity of the objective function. To simplify the problem the solution is often found via a sequence of linearised objective functions. The condition number of the Hessian of the linearised problem is an important indicator of the convergence rate of the minimisation and the expected accuracy of the solution. In the standard formulation the convergence is slow, indicating an ill-conditioned objective function. A transformation to different variables is often used to ameliorate the conditioning of the Hessian by changing, or preconditioning, the Hessian. There is only sparse information in the literature for describing the causes of ill-conditioning of the optimal state estimation problem and explaining the effect of preconditioning on the condition number. This paper derives descriptive theoretical bounds on the condition number of both the unpreconditioned and preconditioned system in order to better understand the conditioning of the problem. We use these bounds to explain why the standard objective function is often ill-conditioned and why a standard preconditioning reduces the condition number. We also use the bounds on the preconditioned Hessian to understand the main factors that affect the conditioning of the system. We illustrate the results with simple numerical experiments.

On integral equation and least squares methods for scattering by diffraction gratings

In this paper we consider the scattering of a plane acoustic or electromagnetic wave by a one-dimensional, periodic rough surface. We restrict the discussion to the case when the boundary is sound soft in the acoustic case, perfectly reflecting with TE polarization in the EM case, so that the total field vanishes on the boundary. We propose a uniquely solvable first kind integral equation formulation of the problem, which amounts to a requirement that the normal derivative of the Green's representation formula for the total field vanish on a horizontal line below the scattering surface. We then discuss the numerical solution by Galerkin's method of this (ill-posed) integral equation. We point out that, with two particular choices of the trial and test spaces, we recover the so-called SC (spectral-coordinate) and SS (spectral-spectral) numerical schemes of DeSanto et al., Waves Random Media, 8, 315-414 1998. We next propose a new Galerkin scheme, a modification of the SS method that we term the SS* method, which is an instance of the well-known dual least squares Galerkin method. We show that the SS* method is always well-defined and is optimally convergent as the size of the approximation space increases. Moreover, we make a connection with the classical least squares method, in which the coefficients in the Rayleigh expansion of the solution are determined by enforcing the boundary condition in a least squares sense, pointing out that the linear system to be solved in the SS* method is identical to that in the least squares method. Using this connection we show that (reflecting the ill-posed nature of the integral equation solved) the condition number of the linear system in the SS* and least squares methods approaches infinity as the approximation space increases in size. We also provide theoretical error bounds on the condition number and on the errors induced in the numerical solution computed as a result of ill-conditioning. Numerical results confirm the convergence of the SS* method and illustrate the ill-conditioning that arises.

The seasonal forecast of electricity demand: A hierarchical Bayesian model with climatological weather generator

In this paper we focus on the one year ahead prediction of the electricity peak-demand daily trajectory during the winter season in Central England and Wales. We define a Bayesian hierarchical model for predicting the winter trajectories and present results based on the past observed weather. Thanks to the flexibility of the Bayesian approach, we are able to produce the marginal posterior distributions of all the predictands of interest. This is a fundamental progress with respect to the classical methods. The results are encouraging in both skill and representation of uncertainty. Further extensions are straightforward at least in principle. The main two of those consist in conditioning the weather generator model with respect to additional information like the knowledge of the first part of the winter and/or the seasonal weather forecast. Copyright (C) 2006 John Wiley & Sons, Ltd.

Resolution of the uncertainties in the radiative forcing of HFC-134a

HFC-134a (CF3CH2F) is the most rapidly growing hydrofluorocarbon in terms of atmospheric abundance. It is currently used in a large number of household refrigerators and air-conditioning systems and its concentration in the atmosphere is forecast to increase substantially over the next 50–100 years. Previous estimates of its radiative forcing per unit concentration have differed significantly 25%. This paper uses a two-step approach to resolve this discrepancy. In the first step six independent absorption cross section datasets are analysed. We find that, for the integrated cross section in the spectral bands that contribute most to the radiative forcing, the differences between the various datasets are typically smaller than 5% and that the dependence on pressure and temperature is not significant. A “recommended'' HFC-134a infrared absorption spectrum was obtained based on the average band intensities of the strongest bands. In the second step, the “recommended'' HFC-134a spectrum was used in six different radiative transfer models to calculate the HFC-134a radiative forcing efficiency. The clear-sky instantaneous radiative forcing, using a single global and annual mean profile, differed by 8%, between the 6 models, and the latitudinally-resolved adjusted cloudy sky radiative forcing estimates differed by a similar amount.

The social impacts of heat waves

The effects of the 2003 European heat wave have highlighted the need for society to prepare itself for and cope more effectively with heat waves. This is particularly important in the context of predicted climate change and the likelihood of more frequent extreme climate events; to date, heat as a natural hazard has been largely ignored. In order to develop better coping strategies, this report explores the factors that shape the social impacts of heat waves, and sets out a programme of research to address the considerable knowledge gaps in this area. Heat waves, or periods of anomalous warmth, do not affect everyone; it is the vulnerable individuals or sectors of society who will most experience their effects. The main factors of vulnerability are being elderly, living alone, having a pre-existing disease, being immobile or suffering from mental illness and being economically disadvantaged. The synergistic effects of such factors may prove fatal for some. Heat waves have discernible impacts on society including a rise in mortality, an increased strain on infrastructure (power, water and transport) and a possible rise in social disturbance. Wider impacts may include effects on the retail industry, ecosystem services and tourism. Adapting to more frequent heat waves should include soft engineering options and, where possible, avoid the widespread use of air conditioning which could prove unsustainable in energy terms. Strategies for coping with heat include changing the way in which urban areas are developed or re-developed, and setting up heat watch warning systems based around weather and seasonal climate forecasting and intervention strategies. Although heat waves have discernible effects on society, much remains unknown about their wider social impacts, diffuse health issues and how to manage them.

Estimators in capture–recapture studies with two sources

This paper investigates the applications of capture–recapture methods to human populations. Capture–recapture methods are commonly used in estimating the size of wildlife populations but can also be used in epidemiology and social sciences, for estimating prevalence of a particular disease or the size of the homeless population in a certain area. Here we focus on estimating the prevalence of infectious diseases. Several estimators of population size are considered: the Lincoln–Petersen estimator and its modified version, the Chapman estimator, Chao’s lower bound estimator, the Zelterman’s estimator, McKendrick’s moment estimator and the maximum likelihood estimator. In order to evaluate these estimators, they are applied to real, three-source, capture-recapture data. By conditioning on each of the sources of three source data, we have been able to compare the estimators with the true value that they are estimating. The Chapman and Chao estimators were compared in terms of their relative bias. A variance formula derived through conditioning is suggested for Chao’s estimator, and normal 95% confidence intervals are calculated for this and the Chapman estimator. We then compare the coverage of the respective confidence intervals. Furthermore, a simulation study is included to compare Chao’s and Chapman’s estimator. Results indicate that Chao’s estimator is less biased than Chapman’s estimator unless both sources are independent. Chao’s estimator has also the smaller mean squared error. Finally, the implications and limitations of the above methods are discussed, with suggestions for further development.

Inheritance of time to flowering in cowpea (Vigna unguiculata (L.) Walp.)

Time to flowering and maturity is an important adaptive feature in annual crops, including cowpeas (Vigna unguiculata (L.) Walp.). In West and Central Africa, photoperiod is the most important environmental variable affecting time to flowering in cowpea. The inheritance of time from sowing to flowering (f) in cowpeas was studied by crossing a photoperiod-sensitive genotype Kanannnado to a photoperiod-insensitive variety IT97D-941-1. Sufficient seed of F-1, F-2, F-3 and backcross populations were generated. The parental, F-1, F-2, F-3 and the backcross populations were screened for f under long natural days (mean daylength 13.4 h per day) in the field and the parents, F-1, F-2 and backcross populations under short day (10 h per day) conditions. The result of the screening showed that photoperiod in the field was long enough to delay flowering of photoperiod-sensitive genotypes. Photoperiod-sensitivity was found to be partially dominant to insensitivity. Frequency distribution of the trait in the various populations indicated quantitative inheritance. Additive (d) and additive x dominance (j) interactions were the most important gene actions conditioning time to flowering. A narrow sense heritability of 86% was estimated for this trait. This will result in 26 days gain in time to flowering with 5% selection intensity from the F-2 to F-3 generation. At least seven major gene pairs, with an average delay of 6 days each, were estimated to control time to flowering in this cross.

Estimators in capture-recapture studies with two sources

This paper investigates the applications of capture-recapture methods to human populations. Capture-recapture methods are commonly used in estimating the size of wildlife populations but can also be used in epidemiology and social sciences, for estimating prevalence of a particular disease or the size of the homeless population in a certain area. Here we focus on estimating the prevalence of infectious diseases. Several estimators of population size are considered: the Lincoln-Petersen estimator and its modified version, the Chapman estimator, Chao's lower bound estimator, the Zelterman's estimator, McKendrick's moment estimator and the maximum likelihood estimator. In order to evaluate these estimators, they are applied to real, three-source, capture-recapture data. By conditioning on each of the sources of three source data, we have been able to compare the estimators with the true value that they are estimating. The Chapman and Chao estimators were compared in terms of their relative bias. A variance formula derived through conditioning is suggested for Chao's estimator, and normal 95% confidence intervals are calculated for this and the Chapman estimator. We then compare the coverage of the respective confidence intervals. Furthermore, a simulation study is included to compare Chao's and Chapman's estimator. Results indicate that Chao's estimator is less biased than Chapman's estimator unless both sources are independent. Chao's estimator has also the smaller mean squared error. Finally, the implications and limitations of the above methods are discussed, with suggestions for further development.

The seasonal forecast of electricity demand: A hierarchical Bayesian model with climatological weather generator

In this paper we focus on the one year ahead prediction of the electricity peak-demand daily trajectory during the winter season in Central England and Wales. We define a Bayesian hierarchical model for predicting the winter trajectories and present results based on the past observed weather. Thanks to the flexibility of the Bayesian approach, we are able to produce the marginal posterior distributions of all the predictands of interest. This is a fundamental progress with respect to the classical methods. The results are encouraging in both skill and representation of uncertainty. Further extensions are straightforward at least in principle. The main two of those consist in conditioning the weather generator model with respect to additional information like the knowledge of the first part of the winter and/or the seasonal weather forecast. Copyright (C) 2006 John Wiley & Sons, Ltd.

Bahrain's Formula-1 racing circuit: energy and environmental considerations

The Bahrain International Circuit (BIC) and complex, at latitude 26.00N and longitude 51.54E, was built in 483 days and cost 150 million US$. The circuit consists of six different individual tracks with a 3.66 km outer track (involving 10 turns) and a 2.55 km inner track (having six turns). The complex has been designed to host a variety of other sporting activities. Fifty thousand spectators, including 10,500 in the main grandstand, can be accommodated simultaneously. State-of-the art on-site media and broadcast facilities are available. The noise level emitted from vehicles on the circuit during the Formula-1 event, on April 4th 2004, was acceptable and caused no physical disturbance to the fans in the VIP lounges or to scholars studying at the University of Bahrain's Shakeir Campus, which is only 1.5 km away from the circuit. The sound-intensity level (SIL) recorded on the balcony of the VIP lounge was 128 dB(A) and was 80 dB(A) inside the lounge. The calculated SIL immediately outside the lecture halls of the University of Bahrain was 70 dB(A) and 65 dB(A) within them. Thus racing at BIC can proceed without significantly disturbing the academic-learning process. The purchased electricity demand by the BIC complex peaked (at 4.5 MW) during the first Formula-1 event on April 4th 2004. The reverse-osmosis (RO) plant at the BIC provides 1000 m(3) of desalinated water per day for landscape irrigation. Renewable-energy inputs, (i.e., via solar and wind power), at the BIC could be harnessed to generate electricity for water desalination, air conditioning, lighting as well as for irrigation. If the covering of the BIC complex was covered by adhesively fixed modern photovoltaic cells, then similar to 1.2 MW of solar electricity could be generated. If two horizontal-axis, at 150 m height above the ground, three 75m bladed, wind turbines were to be installed at the BIC, then the output could reach 4 MW. Furthermore, if 10,000 Jojoba trees (a species renowned for having a low demand for water, needing only five irrigations per year in Bahrain and which remain green throughout the year) are planted near the circuit, then the local micro-climate would be improved with respect to human comfort as well as the local environment becoming cleaner.