A Bayesian approach for multi-model downscaling: Seasonal forecasting of regional rainfall and river flows in South America

This study addresses three issues: spatial downscaling, calibration, and combination of seasonal predictions produced by different coupled ocean-atmosphere climate models. It examines the feasibility Of using a Bayesian procedure for producing combined, well-calibrated downscaled seasonal rainfall forecasts for two regions in South America and river flow forecasts for the Parana river in the south of Brazil and the Tocantins river in the north of Brazil. These forecasts are important for national electricity generation management and planning. A Bayesian procedure, referred to here as forecast assimilation, is used to combine and calibrate the rainfall predictions produced by three climate models. Forecast assimilation is able to improve the skill of 3-month lead November-December-January multi-model rainfall predictions over the two South American regions. Improvements are noted in forecast seasonal mean values and uncertainty estimates. River flow forecasts are less skilful than rainfall forecasts. This is partially because natural river flow is a derived quantity that is sensitive to hydrological as well as meteorological processes, and to human intervention in the form of reservoir management.

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.

General concepts in meteorology and dynamics

The aim of this chapter is to give a general overview of the atmospheric circulation, highlighting the main concepts that are important for a basic understanding of meteorology and atmospheric dynamics relevant to atmospheric data assimilation.

Characterising the background errors for the boundary-layer capping inversion

The one-dimensional variational assimilation of vertical temperature information in the presence of a boundary-layer capping inversion is studied. For an optimal analysis of the vertical temperature profile, an accurate representation of the background error covariances is essential. The background error covariances are highly flow-dependent due to the variability in the presence, structure and height of the boundary-layer capping inversion. Flow-dependent estimates of the background error covariances are shown by studying the spread in an ensemble of forecasts. A forecast of the temperature profile (used as a background state) may have a significant error in the position of the capping inversion with respect to observations. It is shown that the assimilation of observations may weaken the inversion structure in the analysis if only magnitude errors are accounted for as is the case for traditional data assimilation methods used for operational weather prediction. The positional error is treated explicitly here in a new data assimilation scheme to reduce positional error, in addition to the traditional framework to reduce magnitude error. The distribution of the positional error of the background inversion is estimated for use with the new scheme.

Northern hemisphere snow: measurement, modelling and predictability

Northern hemisphere snow water equivalent (SWE) distribution from remote sensing (SSM/I), the ERA40 reanalysis product and the HadCM3 general circulation model are compared. Large differences are seen in the February climatologies, particularly over Siberia. The SSM/I retrieval algorithm may be overestimating SWE in this region, while comparison with independent runoff estimates suggest that HadCM3 is underestimating SWE. Treatment of snow grain size and vegetation parameterizations are concerns with the remotely sensed data. For this reason, ERA40 is used as `truth' for the following experiments. Despite the climatology differences, HadCM3 is able to reproduce the distribution of ERA40 SWE anomalies when assimilating ERA40 anomaly fields of temperature, sea level pressure, atmospheric winds and ocean temperature and salinity. However when forecasts are released from these assimilated initial states, the SWE anomaly distribution diverges rapidly from that of ERA40. No predictability is seen from one season to another. Strong links between European SWE distribution and the North Atlantic Oscillation (NAO) are seen, but forecasts of this index by the assimilation scheme are poor. Longer term relationships between SWE and the NAO, and SWE and the El Ni\~no-Southern Oscillation (ENSO) are also investigated in a multi-century run of HadCM3. SWE is impacted by ENSO in the Himalayas and North America, while the NAO affects SWE in North America and Europe. While significant connections with the NAO index were only present in DJF (and to an extent SON), the link between ENSO and February SWE distribution was seen to exist from the previous JJA ENSO index onwards. This represents a long lead time for SWE prediction for hydrological applications such as flood and wildfire forecasting. Further work is required to develop reliable large scale observation-based SWE datasets with which to test these model-derived connections.

Can 4D-Var use dynamical information from targeted observations of a baroclinic structure?

Targeted observations are generally taken in regions of high baroclinicity, but often show little impact. One plausible explanation is that important dynamical information, such as upshear tilt, is not extracted from the targeted observations by the data assimilation scheme and used to correct initial condition error. This is investigated by generating pseudo targeted observations which contain a singular vector (SV) structure that is not present in the background field or routine observations, i.e. assuming that the background has an initial condition error with tilted growing structure. Experiments were performed for a single case-study with varying numbers of pseudo targeted observations. These were assimilated by the Met Office four-dimensional variational (4D-Var) data assimilation scheme, which uses a 6 h window for observations and background-error covariances calculated using the National Meteorological Centre (NMC) method. The forecasts were run using the operational Met Office Unified Model on a 24 km grid. The results presented clearly demonstrate that a 6 h window 4D-Var system is capable of extracting baroclinic information from a limited set of observations and using it to correct initial condition error. To capture the SV structure well (projection of 0.72 in total energy), 50 sondes over an area of 1×106 km2 were required. When the SV was represented by only eight sondes along an example targeting flight track covering a smaller area, the projection onto the SV structure was lower; the resulting forecast perturbations showed an SV structure with increased tilt and reduced initial energy. The total energy contained in the perturbations decreased as the SV structure was less well described by the set of observations (i.e. as fewer pseudo observations were assimilated). The assimilated perturbation had lower energy than the SV unless the pseudo observations were assimilated with the dropsonde observation errors halved from operational values. Copyright © 2010 Royal Meteorological Society

Estimating discharge from inundation imagery

Remote sensing from space-borne platforms is often seen as an appealing method of monitoring components of the hydrological cycle, including river discharge, due to its spatial coverage. However, data from these platforms is often less than ideal because the geophysical properties of interest are rarely measured directly and the measurements that are taken can be subject to significant errors. This study assimilated water levels derived from a TerraSAR-X synthetic aperture radar image and digital aerial photography with simulations from a two dimensional hydraulic model to estimate discharge, inundation extent, depths and velocities at the confluence of the rivers Severn and Avon, UK. An ensemble Kalman filter was used to assimilate spot heights water levels derived by intersecting shorelines from the imagery with a digital elevation model. Discharge was estimated from the ensemble of simulations using state augmentation and then compared with gauge data. Assimilating the real data reduced the error between analyzed mean water levels and levels from three gauging stations to less than 0.3 m, which is less than typically found in post event water marks data from the field at these scales. Measurement bias was evident, but the method still provided a means of improving estimates of discharge for high flows where gauge data are unavailable or of poor quality. Posterior estimates of discharge had standard deviations between 63.3 m3s-1 and 52.7 m3s-1, which were below 15% of the gauged flows along the reach. Therefore, assuming a roughness uncertainty of 0.03-0.05 and no model structural errors discharge could be estimated by the EnKF with accuracy similar to that arguably expected from gauging stations during flood events. Quality control prior to assimilation, where measurements were rejected for being in areas of high topographic slope or close to tall vegetation and trees, was found to be essential. The study demonstrates the potential, but also the significant limitations of currently available imagery to reduce discharge uncertainty in un-gauged or poorly gauged basins when combined with model simulations in a data assimilation framework.

Use of a crop model ensemble to quantify CO2 stimulation of water-stressed and well-watered crops

Increased atmospheric concentrations of carbon dioxide (CO2) will benefit the yield of most crops. Two free air CO2 enrichment (FACE) meta-analyses have shown increases in yield of between 0 and 73% for C3 crops. Despite this large range, few crop modelling studies quantify the uncertainty inherent in the parameterisation of crop growth and development. We present a novel perturbed-parameter method of crop model simulation, which uses some constraints from observations, that does this. The model used is the groundnut (i.e. peanut; Arachis hypogaea L.) version of the general large-area model for annual crops (GLAM). The conclusions are of relevance to C3 crops in general. The increases in yield simulated by GLAM for doubled CO2 were between 16 and 62%. The difference in mean percentage increase between well-watered and water-stressed simulations was 6.8. These results were compared to FACE and controlled environment studies, and to sensitivity tests on two other crop models of differing levels of complexity: CROPGRO, and the groundnut model of Hammer et al. [Hammer, G.L., Sinclair, T.R., Boote, K.J., Wright, G.C., Meinke, H., Bell, M.J., 1995. A peanut simulation model. I. Model development and testing. Agron. J. 87, 1085-1093]. The relationship between CO2 and water stress in the experiments and in the models was examined. From a physiological perspective, water-stressed crops are expected to show greater CO2 stimulation than well-watered crops. This expectation has been cited in literature. However, this result is not seen consistently in either the FACE studies or in the crop models. In contrast, leaf-level models of assimilation do consistently show this result. An analysis of the evidence from these models and from the data suggests that scale (canopy versus leaf), model calibration, and model complexity are factors in determining the sign and magnitude of the interaction between CO2 and water stress. We conclude from our study that the statement that 'water-stressed crops show greater CO2 stimulation than well-watered crops' cannot be held to be universally true. We also conclude, preliminarily, that the relationship between water stress and assimilation varies with scale. Accordingly, we provide some suggestions on how studies of a similar nature, using crop models of a range of complexity, could contribute further to understanding the roles of model calibration, model complexity and scale. (C) 2008 Elsevier B.V. All rights reserved.

Major advances in fundamental dairy cattle nutrition

Fundamental nutrition seeks to describe the complex biochemical reactions involved in assimilation and processing of nutrients by various tissues and organs, and to quantify nutrient movement (flux) through those processes. Over the last 25 yr, considerable progress has been made in increasing our understanding of metabolism in dairy cattle. Major advances have been made at all levels of biological organization, including the whole animal, organ systems, tissues, cells, and molecules. At the whole-animal level, progress has been made in delineating metabolism during late pregnancy and the transition to lactation, as well as in whole-body use of energy-yielding substrates and amino acids for growth in young calves. An explosion of research using multicatheterization techniques has led to better quantitative descriptions of nutrient use by tissues of the portal-drained viscera (digestive tract, pancreas, and associated adipose tissues) and liver. Isolated tissue preparations have provided important information on the interrelationships among glucose, fatty acid, and amino acid metabolism in liver, adipose tissue, and mammary gland, as well as the regulation of these pathways during different physiological states. Finally, the last 25 yr has witnessed the birth of "molecular biology" approaches to understanding fundamental nutrition. Although measurements of mRNA abundance for proteins of interest already have provided new insights into regulation of metabolism, the next 25 yr will likely see remarkable advances as these techniques continue to be applied to problems of dairy cattle biology. Integration of the "omics" technologies (functional genomics, proteomics, and metabolomics) with measurements of tissue metabolism obtained by other methods is a particularly exciting prospect for the future. The result should be improved animal health and well being, more efficient dairy production, and better models to predict nutritional requirements and provide rations to meet those requirements.

Modelling simultaneously water content and dry matter dynamics of wheat grains

A model was devised to describe simultaneously the grain masses of water and dry matter against thermal time during grain filling and maturation of winter wheat. The model accounted for a linear increase in water mass of duration anthesis-m(1) (end of rapid water assimilation phase) and rate a, followed by a more stable water mass until in,, after which water mass declined rapidly at rate e. Grain dry matter was described as a linear increase of rate bgf until a maximum size (maxgf) was attained at m(2).The model was fitted to plot data from weekly samples of grains taken from replicated field experiments investigating effects of grain position (apical or medial), fungicide (five contrasting treatments), sowing date (early or late), cultivar (Malacca or Shamrock) and season (2001/2002 and 2002/2003) on grain filling. The model accounted for between 83 and 99% of the variation ( 2) when fitted to data from individual plots, and between 97 and 99% when fitted to treatment means. Endosperm cell number of grains from early-sown plots in the first season were also counted. Differences in maxgf between grain positions and also between cultivars were mostly the result of effects on bgf and were empirically associated with water mass at nil. Fungicide application controlled S. tritici and powdery mildew infection, delayed flag leaf senescence, increased water mass at m(1) (wm(1)), and also increased m(2), bgf and maxgf. Fungicide effects on water mass were detected before fungicide effects on dry matter, but comparison of the effects of individual fungicide treatments showed no evidence that effects on wm(1), nor on endosperm cell numbers at about m(1), were required for fungicide effects on maxgf, (c) 2005 Elsevier B.V. All rights reserved.

Differential effects of temperature on fruit development and bean quality of contrasting genotypes of cacao (Theobroma cacao)

The effects of temperature and light integral on fruit growth and development of five cacao genotypes (Amelonado, AMAZ 15/15, SCA 6, SPEC 54/1 and UF 676) were studied in semi-controlled environment glasshouses in which the thermal regimes of cacao-growing regions of Brazil, Ghana and Malaysia were simulated. Fruit losses because of physiological will (cherelle will) were greater at higher temperatures and also differed significantly between genotypes, reflecting genetic differences in competition for assimilates between vegetative and reproductive components. Short-term measurements of fruit growth indicated faster growth rates at higher temperatures. In addition, a significant negative linear relationship between temperature and development time was observed. There was an effect of genotype on this relationship, such that time to fruit maturation at a given temperature was greatest for the clone UF 676 and least for AMAZ 15/15. Analysis of base temperatures, derived from these relationships indicated genetic variability in sensitivity of cacao fruit growth to temperature (base temperatures ranged from 7.5 degrees C for Amelonado and AMAZ 15/15 to 12.9 for SPEC 54/1). Final fruit size was a positive function of beam number for all genotypes and a positive function of light integral for Amelonado in the Malaysia simulated environment (where the temperature was almost constant). In simulated environments where temperature was the main variable (Brazil and Ghana) increases in temperature resulted in a significant decrease in final pod size for one genotype (Amelonado) in Brazil and for two genotypes (SPEC 54/1 and UF 676) in Ghana. It was hypothesised that pod growth duration (mediated by temperature), assimilation and beam number are all determinants of final pod size but that under specific conditions one of these factors may override the others. There was variability between genotypes in the response of beam size and beam lipid content to temperature. Negative relationships between temperature and bean size were found for Amelonado and UF 676. Lipid concentration was a curvilinear function of temperature for Amelonado and UF 676, with optimal temperatures of 23 degrees C and 24 degrees C, respectively. The variability observed here of different cacao genotypes to temperature highlights the need and opportunities for appropriate matching of planting material with local environments.

Shoot dieback in clipped young Golden Leyland (Cupressocyparis leylandii) trees - a physiological mechanism?

Shoot dieback is a problem in frequently trimmed Leyland hedges and is increasingly affecting gardeners’ choice of hedge trees, having a negative effect on a conifer nursery industry. Some damage can be attributed to the feeding by aphids, but it is unclear if there are also underlying physiological causes. In this study, we tested the hypothesis that shoot-clipping of conifer trees during adverse growing conditions (i.e. high air temperature and low soil moisture) could be leading to shoot ‘dieback’. Three-year-old Golden Leyland Cypress (x Cupressocyparis leylandii ‘Excalibur Gold’) plants were subjected to either a well-watered or droughted irrigation regime and placed in either a ‘hot’ (average day temperature = 40°C) or a ‘cool’ (average day temperature = 27°C) glasshouse compartment. Half of the plants from each glasshouse were clipped on Day 14 and again on Day 50. Measurements of soil moisture content (SMC), net CO2 assimilation rate (A), stomatal conductance (gs), branchlet xylem water potential (XWP), plant height and foliage colour were made. Within the clipped and unclipped treatments of both glasshouse compartments, plants from the droughted regime had significantly lower values for A, gs and XWP than those from the well-watered regime. However, there was no difference in these parameters between the hot and cool glasshouse compartments. The trends seen for A, gs and XWP of all treatments generally mirrored changes in SMC indicating a direct effect of water supply on these parameters. By the end of the experiment the overall foliage colour of plants from the hot glasshouse was darker than that of plants from the cool glasshouse and the overall foliage colour was also darker following shoot clipping. In general, shoot clipping led to increases in A, gs XWP and SMC. This may be due to the reduction in total leaf area leading to a greater supply of water for the remaining leaves. No shoot ‘dieback’ was observed in any treatment in response to drought stress or shoot-clipping.

Root pruning as a means to encourage root growth in two ornamental shrubs, Buddleja davidii ‘Summer Beauty’ and Cistus ‘Snow Fire'

The inability of a plant to grow roots rapidly upon transplanting is one of the main factors contributing to poor establishment. In bare-rooted trees, treatments such as root pruning or application of the plant hormone auxin [e.g., indole butyric acid (IBA)] can promote root growth and aid long-term establishment. There is little information on ornamental containerised plants, however, other than the anecdotal notion that 'teasing' the roots out of the rootsoil mass before transplanting can be beneficial. In the present study we tested the ability of various root-pruning treatments and application of IBA to encourage new root and shoot growth in two shrub species, commonly produced in containers - Buddleja davidii 'Summer Beauty' and Cistus 'Snow Fire'. In a number of experiments, young plants were exposed to root manipulation (teasing, light pruning, or two types of heavy pruning) and/or treatment with IBA (at 500 or 1,000 mg l-1) before being transplanted into larger containers containing a medium of 1:1:1 (v/v/v) fine bark, sand and loam. Leaf stomatal conductance (gl) was measured 20 min, and 1, 2, 4 and 6 h after root manipulation. Net leaf CO2 assimilation (A) was measured frequently during the first week after transplanting, then at regular intervals up to 8 weeks after transplanting. Plants were harvested 8 weeks after transplanting, and root and shoot weights were measured. In both species, light root pruning alone, or in combination with 500 mg l-1 IBA, was most effective in stimulating root growth. In contrast, teasing, which is commonly used, showed no positive effect on root growth in Buddleja, and decreased new root growth in Cistus. The requirement for exogenous auxin to encourage new root growth varied between experiments and appeared to be influenced by the age and developmental stage of the plants. There were no consistent responses between root treatments and net CO2 assimilation rates, and changes in root weight were not closely correlated with changes in assimilation. The mechanisms whereby new root growth is sustained are discussed.