The effect of nonlinearity on the variational assimilation of satellite observations using a simple column model

Cloud imagery is not currently used in numerical weather prediction (NWP) to extract the type of dynamical information that experienced forecasters have extracted subjectively for many years. For example, rapidly developing mid-latitude cyclones have characteristic signatures in the cloud imagery that are most fully appreciated from a sequence of images rather than from a single image. The Met Office is currently developing a technique to extract dynamical development information from satellite imagery using their full incremental 4D-Var (four-dimensional variational data assimilation) system. We investigate a simplified form of this technique in a fully nonlinear framework. We convert information on the vertical wind field, w(z), and profiles of temperature, T(z, t), and total water content, qt (z, t), as functions of height, z, and time, t, to a single brightness temperature by defining a 2D (vertical and time) variational assimilation testbed. The profiles of w, T and qt are updated using a simple vertical advection scheme. We define a basic cloud scheme to obtain the fractional cloud amount and, when combined with the temperature field, we convert this information into a brightness temperature, having developed a simple radiative transfer scheme. With the exception of some matrix inversion routines, all our code is developed from scratch. Throughout the development process we test all aspects of our 2D assimilation system, and then run identical twin experiments to try and recover information on the vertical velocity, from a sequence of observations of brightness temperature. This thesis contains a comprehensive description of our nonlinear models and assimilation system, and the first experimental results.

On the 'downward control' of extratropical diabatic circulations by eddy-induced mean zonal forces

The situation considered is that of a zonally symmetric model of the middle atmosphere subject to a given quasi-steady zonal force F̄, conceived to be the result of irreversible angular momentum transfer due to the upward propagation and breaking of Rossby and gravity waves together with any other dissipative eddy effects that may be relevant. The model's diabatic heating is assumed to have the qualitative character of a relaxation toward some radiatively determined temperature field. To the extent that the force F̄ may be regarded as given, and the extratropical angular momentum distribution is realistic, the extratropical diabatic mass flow across a given isentropic surface may be regarded as controlled exclusively by the F̄ distribution above that surface (implying control by the eddy dissipation above that surface and not, for instance, by the frequency of tropopause folding below). This “downward control” principle expresses a critical part of the dynamical chain of cause and effect governing the average rate at which photochemical products like ozone become available for folding into, or otherwise descending into, the extratropical troposphere. The dynamical facts expressed by the principle are also relevant, for instance, to understanding the seasonal-mean rate of upwelling of water vapor to the summer mesopause, and the interhemispheric differences in stratospheric tracer transport. The robustness of the principle is examined when F̄ is time-dependent. For a global-scale, zonally symmetric diabatic circulation with a Brewer-Dobson-like horizontal structure given by the second zonally symmetric Hough mode, with Rossby height HR = 13 km in an isothermal atmosphere with density scale height H = 7 km, the vertical partitioning of the unsteady part of the mass circulation caused by fluctuations in F̄ confined to a shallow layer LF̄ is always at least 84% downward. It is 90% downward when the force fluctuates sinusoidally on twice the radiative relaxation timescale and 95% if five times slower. The time-dependent adjustment when F̄ is changed suddenly is elucidated, extending the work of Dickinson (1968), when the atmosphere is unbounded above and below. Above the forcing, the adjustment is characterized by decay of the meridional mass circulation cell at a rate proportional to the radiative relaxation rate τr−1 divided by {1 + (4H2/HR2)}. This decay is related to the boundedness of the angular momentum that can be taken up by the finite mass of air above LF̄ without causing an ever-increasing departure from thermal wind balance. Below the forcing, the meridional mass circulation cell penetrates downward at a speed τr−1 HR2/H. For the second Hough mode, the time for downward penetration through one density scale height is about 6 days if the radiative relaxation time is 20 days, the latter being representative of the lower stratosphere. At any given altitude, a steady state is approached. The effect of a rigid lower boundary on the time-dependent adjustment is also considered. If a frictional planetary boundary layer is present then a steady state is ultimately approached everywhere, with the mass circulation extending downward from LF̄ and closing via the boundary layer. Satellite observations of temperature and ozone are used in conjunction with a radiative transfer scheme to estimate the altitudes from which the lower stratospheric diabatic vertical velocity is controlled by the effective F̄ in the real atmosphere. The data appear to indicate that about 80% of the effective control is usually exerted from below 40 km but with significant exceptions up to 70 km (in the high latitude southern hemispheric winter). The implications for numerical modelling of chemical transport are noted.

The effect of non-uniform radiative damping on the zonal-mean dynamics of the extratropical middle atmosphere

The effect of spatial and temporal variations in the radiative damping rate on the response to an imposed forcing or diabatic heating is examined in a zonal-mean model of the middle atmosphere. Attention is restricted to the extratropics, where a linear approach is viable. It is found that regions with weak radiative damping rates are more sensitive in terms of temperature to the remote influence of the diabatic circulation. The delay in the response in such regions can mean that ‘downward’ control is not achieved on seasonal time-scales. A seasonal variation in the radiative damping rate modulates the evolution of the response and leaves a transient-like signature in the annual mean temperature field. Several idealized examples are considered, motivated by topical questions. It is found that wave drag outside the polar vortex can significantly affect the temperatures in its interior, so that high-latitude, high-altitude gravity-wave drag is not the only mechanism for warming the southern hemisphere polar vortex. Diabatic mass transport through the 100 hPa surface is found to lag the seasonal evolution of the wave drag that drives the transport, and thus cannot be considered to be in the downward control regime. On the other hand, the seasonal variation of the radiative damping rate is found to make only a weak contribution to the annual mean temperature increase that has been observed above the ozone hole. Copyright © 2002 Royal Meteorological Society.

Effect of high temperature and water stress on groundnuts under field conditions

Groundnuts cultivated in the semiarid tropics are often exposed to water stress (mid-season and end season) and high temperature (> 34 °C) during the critical stages of flowering and pod development. This study evaluated the effects of both water stress and high temperature under field conditions at ICRISAT, India. Treatments included two irrigations (full irrigation, 100 % of crop evapotranspiration; and water stress, 40 % of crop evapotranspiration), four temperature treatments from a combination of two sowing dates and heat tunnels with mean temperatures from sowing to maturity of 26.3° (T1), 27.3° (T2), 29.0° (T3) and 29.7 °C (T4) and two genotypes TMV2 and ICGS 11. The heat tunnels were capable of raising the day temperature by > 10 °C compared to ambient. During the 20-day high-temperature treatment at flowering, mean temperatures were 33.8° (T1), 41.6° (T2), 38.7° (T3) and 43.5°C (T4). The effects of water stress and high temperature were additive and temporary for both vegetative and pod yield, and disappeared as soon as high-temperature stress was removed. Water use efficiency was significantly affected by the main effects of temperature and cultivar and not by water stress treatments. Genotypic differences for tolerance to high temperature can be attributed to differences in flowering pattern, flower number, peg-set and harvest index. It can be inferred from this study that genotypes that are tolerant to water stress are also tolerant to high temperature under field conditions. In addition, genotypes with an ability to establish greater biomass and with a significantly greater partitioning of biomass to pod yield would be suitable for sustaining higher yields in semiarid tropics with high temperature and water stress.

Using Distributed Temperature Sensing to monitor field scale dynamics of ground surface temperature and related substrate heat flux

We present one of the first studies of the use of Distributed Temperature Sensing (DTS) along fibre-optic cables to purposely monitor spatial and temporal variations in ground surface temperature (GST) and soil temperature, and provide an estimate of the heat flux at the base of the canopy layer and in the soil. Our field site was at a groundwater-fed wet meadow in the Netherlands covered by a canopy layer (between 0-0.5 m thickness) consisting of grass and sedges. At this site, we ran a single cable across the surface in parallel 40 m sections spaced by 2 m, to create a 40×40 m monitoring field for GST. We also buried a short length (≈10 m) of cable to depth of 0.1±0.02 m to measure soil temperature. We monitored the temperature along the entire cable continuously over a two-day period and captured the diurnal course of GST, and how it was affected by rainfall and canopy structure. The diurnal GST range, as observed by the DTS system, varied between 20.94 and 35.08◦C; precipitation events acted to suppress the range of GST. The spatial distribution of GST correlated with canopy vegetation height during both day and night. Using estimates of thermal inertia, combined with a harmonic analysis of GST and soil temperature, substrate and soil-heat fluxes were determined. Our observations demonstrate how the use of DTS shows great promise in better characterising area-average substrate/soil heat flux, their spatiotemporal variability, and how this variability is affected by canopy structure. The DTS system is able to provide a much richer data set than could be obtained from point temperature sensors. Furthermore, substrate heat fluxes derived from GST measurements may be able to provide improved closure of the land surface energy balance in micrometeorological field studies. This will enhance our understanding of how hydrometeorological processes interact with near-surface heat fluxes.

An integrated model of soil-canopy spectral radiances, photosynthesis, fluorescence, temperature and energy balance

This paper presents the model SCOPE (Soil Canopy Observation, Photochemistry and Energy fluxes), which is a vertical (1-D) integrated radiative transfer and energy balance model. The model links visible to thermal infrared radiance spectra (0.4 to 50 μm) as observed above the canopy to the fluxes of water, heat and carbon dioxide, as a function of vegetation structure, and the vertical profiles of temperature. Output of the model is the spectrum of outgoing radiation in the viewing direction and the turbulent heat fluxes, photosynthesis and chlorophyll fluorescence. A special routine is dedicated to the calculation of photosynthesis rate and chlorophyll fluorescence at the leaf level as a function of net radiation and leaf temperature. The fluorescence contributions from individual leaves are integrated over the canopy layer to calculate top-of-canopy fluorescence. The calculation of radiative transfer and the energy balance is fully integrated, allowing for feedback between leaf temperatures, leaf chlorophyll fluorescence and radiative fluxes. Leaf temperatures are calculated on the basis of energy balance closure. Model simulations were evaluated against observations reported in the literature and against data collected during field campaigns. These evaluations showed that SCOPE is able to reproduce realistic radiance spectra, directional radiance and energy balance fluxes. The model may be applied for the design of algorithms for the retrieval of evapotranspiration from optical and thermal earth observation data, for validation of existing methods to monitor vegetation functioning, to help interpret canopy fluorescence measurements, and to study the relationships between synoptic observations with diurnally integrated quantities. The model has been implemented in Matlab and has a modular design, thus allowing for great flexibility and scalability.

Assessment of the effect of time and temperature on the availability of residual phosphate in a glasshouse study of four soils using the Olsen method

The Olsen method is an indicator of plant-available phosphorus (P). The effect of time and temperature on residual phosphate in soils was measured using the Olsen method in a pot experiment. Four soils were investigated: two from Pakistan and one each from England (calcareous) and Colombia (acidic). Two levels of residual phosphate were developed in each soil after addition of phosphate by incubation at either 10degreesC or 45degreesC. The amount of phosphate added was based on the P maximum of each soil, calculated using the Langmuir equation. Rvegrass was used as the test crop. The pooled data for the four soils incubated at 10degreesC showed good correlation between Olsen P and dry matter yield or P uptake (r(2) = 0.85 and 0.77, respectively), whereas at 45 degreesC, each soil had its own relationship and pooled data did not show correlation of Olsen P with dry matter yield or P uptake. When the data at both temperatures were pooled, Olsen P was a good indicator of yield and uptake for the English soil. For the Pakistani soils, Olsen P after 45 degreesC treatment was an underestimate relative to the 10 degreesC data and for the Colombian soil it was an overestimate. The reasons for these differences need to be explored further before high temperature incubation can be used to simulate long-term changes in the field.

The diurnal course of soil moisture as measured by various dielectric sensors: Effects of soil temperature and the implications for evaporation estimates

Soil moisture content, theta, of a bare and vegetated UK gravelly sandy loam soil (in situ and repacked in small lysimeters) was measured using various dielectric instruments (single-sensor ThetaProbes, multi-sensor Profile Probes, and Aquaflex Sensors), at depths ranging between 0.03 and I m, during the summers of 2001 (in situ soil) and 2002 (mini-lysimeters). Half-hourly values of evaporation, E, were calculated from diurnal changes in total soil profile water content, using the soil water balance equation. For the bare soil field, Profile Probes and ML2x ThetaProbes indicated a diurnal course of theta that did not concur with typical soil physical observations: surface layer soil moisture content increased from early morning until about midday, after which theta declined, generally until the early evening. The unexpected course of theta was positively correlated to soil temperature, T-s, also at deeper depths. Aquaflex and ML1 ThetaProbe (older models) outputs, however, reflected common observations: 0 increased slightly during the night (capillary rise) and decreased from the morning until late afternoon (as a result of evaporation). For the vegetated plot, the spurious diurnal theta fluctuations were less obvious, because canopy shading resulted in lower amplitudes of T-s. The unrealistic theta profiles measured for the bare and vegetated field sites caused diurnal estimates of E to attain downward daytime and upward night-time values. In the mini-lysimeters, at medium to high moisture contents, theta values measured by (ML2x) ThetaProbes followed a relatively realistic course, and predictions of E from diurnal changes in vertically integrated theta generally compared well with lysimeter estimates of E. However, time courses of theta and E became comparable to those observed for the field plots when the soil in the lysimeters reached relatively low values of theta. Attempts to correct measured theta for fluctuations in T, revealed that no generally applicable formula could be derived. (c) 2005 Elsevier B.V. All rights reserved.

Stratospheric temperature trends: impact of ozone variability and the QBO

In most climate simulations used by the Intergovernmental Panel on Climate Change 2007 fourth assessment report, stratospheric processes are only poorly represented. For example, climatological or simple specifications of time-varying ozone concentrations are imposed and the quasi-biennial oscillation (QBO) of equatorial stratospheric zonal wind is absent. Here we investigate the impact of an improved stratospheric representation using two sets of perturbed simulations with the Hadley Centre coupled ocean atmosphere model HadGEM1 with natural and anthropogenic forcings for the 1979–2003 period. In the first set of simulations, the usual zonal mean ozone climatology with superimposed trends is replaced with a time series of observed zonal mean ozone distributions that includes interannual variability associated with the solar cycle, QBO and volcanic eruptions. In addition to this, the second set of perturbed simulations includes a scheme in which the stratospheric zonal wind in the tropics is relaxed to appropriate zonal mean values obtained from the ERA-40 re-analysis, thus forcing a QBO. Both of these changes are applied strictly to the stratosphere only. The improved ozone field results in an improved simulation of the stepwise temperature transitions observed in the lower stratosphere in the aftermath of the two major recent volcanic eruptions. The contribution of the solar cycle signal in the ozone field to this improved representation of the stepwise cooling is discussed. The improved ozone field and also the QBO result in an improved simulation of observed trends, both globally and at tropical latitudes. The Eulerian upwelling in the lower stratosphere in the equatorial region is enhanced by the improved ozone field and is affected by the QBO relaxation, yet neither induces a significant change in the upwelling trend.

Simulation of the impact of high temperature stress on annual crop yields

Brief periods of high temperature which occur near flowering can severely reduce the yield of annual crops such as wheat and groundnut. A parameterisation of this well-documented effect is presented for groundnut (i.e. peanut; Arachis hypogaeaL.). This parameterisation was combined with an existing crop model, allowing the impact of season-mean temperature, and of brief high-temperature episodes at various times near flowering, to be both independently and jointly examined. The extended crop model was tested with independent data from controlled environment experiments and field experiments. The impact of total crop duration was captured, with simulated duration being within 5% of observations for the range of season-mean temperatures used (20-28 degrees C). In simulations across nine differently timed high temperature events, eight of the absolute differences between observed and simulated yield were less than 10% of the control (no-stress) yield. The parameterisation of high temperature stress also allows the simulation of heat tolerance across different genotypes. Three parameter sets, representing tolerant, moderately sensitive and sensitive genotypes were developed and assessed. The new parameterisation can be used in climate change studies to estimate the impact of heat stress on yield. It can also be used to assess the potential for adaptation of cropping systems to increased temperature threshold exceedance via the choice of genotype characteristics. (c) 2005 Elsevier B.V. All rights reserved.

Effects of restricted water availability and increased temperature on the grain filling, drying and quality of winter wheat

Experiments in controlled environments examined the effects of the timing and severity of drought, and increased temperature, on grain development of Hereward winter wheat. Environmental effects on grain specific weight, protein content, Hagberg Falling Number, SDS-sedimentation volume, and sulphur content were also studied. Drought and increased temperature applied before the end of grain filling shortened the grain filling period and reduced grain yield, mean grain weight and specific weight. Grain filling was most severely affected by drought between days 1-14 after anthesis. Protein content was increased by stresses before the end of grain growth, because nitrogen harvest index was less severely affected than was dry matter harvest index. Hagberg Falling Number was increased to the greatest extent by stresses applied 15-28 days after anthesis. Treatment effects on grain sulphur content were similar to those on protein content, such that N:S ratio was not significantly affected by drought nor temperature stresses. The effects of restricted water on grain yield and quality were linearly related to soil moisture between 44 and about 73% field capacity (FC) from days 15-28. Drought stress (but not temperature stress) before the end of grain filling decreased SDS-sedimentation volume relative to drought applied later. (C) 2003 Elsevier Science Ltd. All rights reserved.

Role of temperature, moisture and Trichoderma species on the survival of Fusarium oxysporum ciceri in the rainfed areas of Pakistan

It has been observed in the present study that when spores of Trichoderma harzianum (Th-2) isolate were applied in the sandy clay loam soil and continuously incubated for 4 months at 25 degrees C and 35 degrees C and at three water potentials, -0.03 MPa, -0.3 MPa and <-50 MPa, it has resulted in significantly reduced (P<0.05), growth of Fusarium oxysporum ciceri (Foc) on branches of chickpea plant. The pathogen population was greatly reduced in the moist soil (43 MPa) when compared with the wet soil (-0.03 MPa) at both temperatures which was indicated by greater colonization and growth of T. harzanum-2 on the branch pieces of chickpea plants. The pathogen was completely eradicated from the chickpea branch pieces, after 6 months at 35 degrees C in the moist soil. In air-dry soil (<-50 MPa), Foc survived in 100% of the branch pieces even after 6 months at both temperatures. When chickpea plant branch pieces having pathogen was sprayed with Th-2 antagonistic isolates of Trichoderma spp., the Th-2 isolate killed the pathogen up to minimum level (10-12%) after 5 months at 35 degrees C in the sandy clay loam soil. It can be concluded that in chickpea growing rainfed areas of Pakistan having sandy clay loam soil, Foc can be controlled by using specific Trichoderma spp., especially in the summer season as after harvest of the crop the temperature increased up and there is rainfall during this period which makes the soil moist. This practice will be able to reduce the inoculum of Foc during this hot period as field remain fallow till next crop is sown in most of the chickpea growing rainfed areas of Pakistan.

The effects of infection by Tetracapsuloides bryosalmonae (Myxozoa) and temperature on Fredericella sultana (Bryozoa)

The myxozoan, Tetracapsuloides bryosalmonae, exploits freshwater bryozoans as definitive hosts, occurring as cryptic stages in bryozoan colonies during covert infections and as spore-forming sacs during overt infections. Spores released from sacs are infective to salmonid fish, causing the devastating Proliferative Kidney Disease (PKD). We undertook laboratory studies using mesocosm systems running at 10, 14 and 20 degrees C to determine how infection by T bryosalmonae and water temperature influence fitness of one of its most important bryozoan hosts, Fredericella sultana, over a period of 4 weeks. The effects of infection were context-dependent and often undetectable. Covert infections appear to pose very low energetic costs. Thus, we found that growth of covertly infected F. sultana colonies was similar to that of uninfected colonies regardless of temperature, as was the propensity to produce dormant resting stages (statoblasts). Production of statoblasts, however, was associated with decreased growth. Overt infections imposed greater effects on correlates of host fitness by: (i) reducing growth rates at the two higher temperatures: (ii) increasing mortality rates at the highest temperature: (iii) inhibiting statoblast production. Our results indicate that parasitism should have a relatively small effect on host fitness in the field as the negative effects of infection were mainly expressed in environmentally extreme conditions (20 degrees C for 4 weeks). The generally low virulence of T. bryosalmonae is similar to that recently demonstrated for another myxozoan endoparasite of freshwater bryozoans. The unique opportunity for extensive vertical transmission in these colonial invertebrate hosts couples the reproductive interests of host and parasite and may well give rise to the low virulence that characterises these systems. Our study implies that climate change can be expected to exacerbate PKD outbreaks and increase the geographic range of PKD as a result of the combined responses of T. bryosalmonae and its bryozoan hosts to higher temperatures. Crown Copyright (C) 2009 Published by Elsevier Ltd. All rights reserved.

Mutations in the NS2B and NS3 genes affect mouse neuroinvasiveness of a Western European field strain of tick-borne encephalitis virus

An attenuated strain (263) of the tick-borne encephalitis virus, isolated from field ticks, was either serially subcultured, 5 times in mice, or at 40 degrees C in PS cells, producing 2 independent strains, 263-m5 and 263-TR with identical genomes; both strains exhibited increased plaque size, neuroinvasiveness and temperature-resistance. Sequencing revealed two unique amino acid substitutions, one mapping close to the catalytic site of the viral protease. These observations imply that virus adaptation from ticks to mammals occurs by selection of pre-existing virulent variants from the quasispecies population rather than by the emergence of new random mutations. The significance of these observations is discussed. (c) 2008 Elsevier Inc. All rights reserved.