Estimating surface CO2 fluxes from space-borne dry air mole fraction observations using an Ensemble Kalman filter

We have developed an ensemble Kalman Filter (EnKF) to estimate 8-day regional surface fluxes of CO2 from space-borne CO2 dry-air mole fraction observations (XCO2) and evaluate the approach using a series of synthetic experiments, in preparation for data from the NASA Orbiting Carbon Observatory (OCO). The 32-day duty cycle of OCO alternates every 16 days between nadir and glint measurements of backscattered solar radiation at short-wave infrared wavelengths. The EnKF uses an ensemble of states to represent the error covariances to estimate 8-day CO2 surface fluxes over 144 geographical regions. We use a 12×8-day lag window, recognising that XCO2 measurements include surface flux information from prior time windows. The observation operator that relates surface CO2 fluxes to atmospheric distributions of XCO2 includes: a) the GEOS-Chem transport model that relates surface fluxes to global 3-D distributions of CO2 concentrations, which are sampled at the time and location of OCO measurements that are cloud-free and have aerosol optical depths <0.3; and b) scene-dependent averaging kernels that relate the CO2 profiles to XCO2, accounting for differences between nadir and glint measurements, and the associated scene-dependent observation errors. We show that OCO XCO2 measurements significantly reduce the uncertainties of surface CO2 flux estimates. Glint measurements are generally better at constraining ocean CO2 flux estimates. Nadir XCO2 measurements over the terrestrial tropics are sparse throughout the year because of either clouds or smoke. Glint measurements provide the most effective constraint for estimating tropical terrestrial CO2 fluxes by accurately sampling fresh continental outflow over neighbouring oceans. We also present results from sensitivity experiments that investigate how flux estimates change with 1) bias and unbiased errors, 2) alternative duty cycles, 3) measurement density and correlations, 4) the spatial resolution of estimated flux estimates, and 5) reducing the length of the lag window and the size of the ensemble. At the revision stage of this manuscript, the OCO instrument failed to reach its orbit after it was launched on 24 February 2009. The EnKF formulation presented here is also applicable to GOSAT measurements of CO2 and CH4.

On the energetics of stratified turbulent mixing, irreversible thermodynamics, Boussinesq models and the ocean heat engine controversy

In this paper, the available potential energy (APE) framework of Winters et al. (J. Fluid Mech., vol. 289, 1995, p. 115) is extended to the fully compressible Navier– Stokes equations, with the aims of clarifying (i) the nature of the energy conversions taking place in turbulent thermally stratified fluids; and (ii) the role of surface buoyancy fluxes in the Munk & Wunsch (Deep-Sea Res., vol. 45, 1998, p. 1977) constraint on the mechanical energy sources of stirring required to maintain diapycnal mixing in the oceans. The new framework reveals that the observed turbulent rate of increase in the background gravitational potential energy GPEr , commonly thought to occur at the expense of the diffusively dissipated APE, actually occurs at the expense of internal energy, as in the laminar case. The APE dissipated by molecular diffusion, on the other hand, is found to be converted into internal energy (IE), similar to the viscously dissipated kinetic energy KE. Turbulent stirring, therefore, does not introduce a new APE/GPEr mechanical-to-mechanical energy conversion, but simply enhances the existing IE/GPEr conversion rate, in addition to enhancing the viscous dissipation and the entropy production rates. This, in turn, implies that molecular diffusion contributes to the dissipation of the available mechanical energy ME =APE +KE, along with viscous dissipation. This result has important implications for the interpretation of the concepts of mixing efficiency γmixing and flux Richardson number Rf , for which new physically based definitions are proposed and contrasted with previous definitions. The new framework allows for a more rigorous and general re-derivation from the first principles of Munk & Wunsch (1998, hereafter MW98)’s constraint, also valid for a non-Boussinesq ocean: G(KE) ≈ 1 − ξ Rf ξ Rf Wr, forcing = 1 + (1 − ξ )γmixing ξ γmixing Wr, forcing , where G(KE) is the work rate done by the mechanical forcing, Wr, forcing is the rate of loss of GPEr due to high-latitude cooling and ξ is a nonlinearity parameter such that ξ =1 for a linear equation of state (as considered by MW98), but ξ <1 otherwise. The most important result is that G(APE), the work rate done by the surface buoyancy fluxes, must be numerically as large as Wr, forcing and, therefore, as important as the mechanical forcing in stirring and driving the oceans. As a consequence, the overall mixing efficiency of the oceans is likely to be larger than the value γmixing =0.2 presently used, thereby possibly eliminating the apparent shortfall in mechanical stirring energy that results from using γmixing =0.2 in the above formula.

Identifying management zones in agricultural fields using spatially constrained classification of soil and ancillary data

Site-specific management requires accurate knowledge of the spatial variation in a range of soil properties within fields. This involves considerable sampling effort, which is costly. Ancillary data, such as crop yield, elevation and apparent electrical conductivity (ECa) of the soil, can provide insight into the spatial variation of some soil properties. A multivariate classification with spatial constraint imposed by the variogram was used to classify data from two arable crop fields. The yield data comprised 5 years of crop yield, and the ancillary data 3 years of yield data, elevation and ECa. Information on soil chemical and physical properties was provided by intensive surveys of the soil. Multivariate variograms computed from these data were used to constrain sites spatially within classes to increase their contiguity. The constrained classifications resulted in coherent classes, and those based on the ancillary data were similar to those from the soil properties. The ancillary data seemed to identify areas in the field where the soil is reasonably homogeneous. The results of targeted sampling showed that these classes could be used as a basis for management and to guide future sampling of the soil.

Two Independent Triggers for the Indian Ocean Dipole/Zonal Mode in a Coupled GCM

The question of whether and how tropical Indian Ocean dipole or zonal mode (IOZM) interannual variability is independent of El Nino-Southern Oscillation (ENSO) variability in the Pacific is addressed in a comparison of twin 200-yr runs of a coupled climate model. The first is a reference simulation, and the second has ENSO-scale variability suppressed with a constraint on the tropical Pacific wind stress. The IOZM can exist in the model without ENSO, and the composite evolution of the main anomalies in the Indian Ocean in the two simulations is virtually identical. Its growth depends on a positive feedback between anomalous equatorial easterly winds, upwelling equatorial and coastal Kelvin waves reducing the thermocline depth and sea surface temperature off the coast of Sumatra, and the atmospheric dynamical response to the subsequently reduced convection. Two IOZM triggers in the boreal spring are found. The first is an anomalous Hadley circulation over the eastern tropical Indian Ocean and Maritime Continent, with an early northward penetration of the Southern Hemisphere southeasterly trades. This situation grows out of cooler sea surface temperatures in the southeastern tropical Indian Ocean left behind by a reinforcement of the late austral summer winds. The second trigger is a consequence of a zonal shift in the center of convection associated with a developing El Nino, a Walker cell anomaly. The first trigger is the only one present in the constrained simulation and is similar to the evolution of anomalies in 1994, when the IOZM occurred in the absence of a Pacific El Nino state. The presence of these two triggers-the first independent of ENSO and the second phase locking the IOZM to El Nino-allows an understanding of both the existence of IOZM events when Pacific conditions are neutral and the significant correlation between the IOZM and El Nino.

Projected gradient approach to the numerical solution of the SCoTLASS

The SCoTLASS problem-principal component analysis modified so that the components satisfy the Least Absolute Shrinkage and Selection Operator (LASSO) constraint-is reformulated as a dynamical system on the unit sphere. The LASSO inequality constraint is tackled by exterior penalty function. A globally convergent algorithm is developed based on the projected gradient approach. The algorithm is illustrated numerically and discussed on a well-known data set. (c) 2004 Elsevier B.V. All rights reserved.

In search of simple structures in climate: Simplifying EOFs

Empirical orthogonal functions (EOFs) are widely used in climate research to identify dominant patterns of variability and to reduce the dimensionality of climate data. EOFs, however, can be difficult to interpret. Rotated empirical orthogonal functions (REOFs) have been proposed as more physical entities with simpler patterns than EOFs. This study presents a new approach for finding climate patterns with simple structures that overcomes the problems encountered with rotation. The method achieves simplicity of the patterns by using the main properties of EOFs and REOFs simultaneously. Orthogonal patterns that maximise variance subject to a constraint that induces a form of simplicity are found. The simplified empirical orthogonal function (SEOF) patterns, being more 'local'. are constrained to have zero loadings outside the main centre of action. The method is applied to winter Northern Hemisphere (NH) monthly mean sea level pressure (SLP) reanalyses over the period 1948-2000. The 'simplified' leading patterns of variability are identified and compared to the leading patterns obtained from EOFs and REOFs. Copyright (C) 2005 Royal Meteorological Society.

The elliptic sine-Gordon equation in a half plane

We consider boundary value problems for the elliptic sine-Gordon equation posed in the half plane y > 0. This problem was considered in Gutshabash and Lipovskii (1994 J. Math. Sci. 68 197–201) using the classical inverse scattering transform approach. Given the limitations of this approach, the results obtained rely on a nonlinear constraint on the spectral data derived heuristically by analogy with the linearized case. We revisit the analysis of such problems using a recent generalization of the inverse scattering transform known as the Fokas method, and show that the nonlinear constraint of Gutshabash and Lipovskii (1994 J. Math. Sci. 68 197–201) is a consequence of the so-called global relation. We also show that this relation implies a stronger constraint on the spectral data, and in particular that no choice of boundary conditions can be associated with a decaying (possibly mod 2π) solution analogous to the pure soliton solutions of the usual, time-dependent sine-Gordon equation. We also briefly indicate how, in contrast to the evolutionary case, the elliptic sine-Gordon equation posed in the half plane does not admit linearisable boundary conditions.

Vibrational energy levels of water

Several quartic force fields and a full sextic anharmonic force field for H,O have been determined from high-quality ab initio calculations, the highest at the aug-cc-pVQZ CCSD(T) level of theory. These force fields have been used to determine vibrational excited state band origins up to 15 000 cm - ’ above the zero-point level, using both a perturbation-resonancea pproach and a variational approach. An optimisedq uartic force field hasb eeno btained by least squares refinement of our best ab initio results to fit the observed overtone levels of 5 symmetrically substituted isotopomers of water (Hi60, Hi70, HisO, D,O, and T,O) with an rms error of less than 10 cm-‘, using the perturbation-resonancem odel for the vibrational calculation. Predicatel east squaresr efinement was usedt o provide a loose constraint of the refined force field to the ab initio results. The results obtained prove the viability of the perturbation-resonancem odel for usei n larger molecular systemsa nd also highlight someo f its weaknesse

The calculation of force constants and normal coordinates III: constrained force fields for the methyl halides

A method is discussed for imposing any desired constraint on the force field obtained in a force constant refinement calculation. The application of this method to force constant refinement calculations for the methyl halide molecules is reported. All available data on the vibration frequencies, Coriolis interaction constants and centrifugal stretching constants of CH3X and CD3X molecules were used in the refinements, but despite this apparent abundance of data it was found that constraints were necessary in order to obtain a unique solution to the force field. The results of unconstrained calculations, and of three different constrained calculations, are reported in this paper. The constrained models reported are a Urey—Bradley force field, a modified valence force field, and a constraint based on orbital-following bond-hybridization arguments developed in the following paper. The results are discussed, and compared with previous results for these molecules. The third of the above models is found to reproduce the observed data better than either of the first two, and additional reasons are given for preferring this solution to the force field for the methyl halide molecules.

The calculation of force constants and normal co-ordinates--VII Normal co-ordinates in the E species of allene

Force constant and normal co-ordinate calculations are reported for the E species vibrations of the allene molecule. Data on the fundamental vibration frequencies of allene-h4, allene-d4 and allene-1.1-d2 and on the five experimentally determined Coriolis zeta constants of C3H4 and C3D4, were used in a force constant refinement procedure. Allowing for product and sum rules this gives 21 independent data which were used to refine to the most general harmonic force field (10 parameters) with one constraint (in the absence of any constraints the refinement was not satisfactory). The results have been used to calculate the complete ζz Coriolis interaction matrix for the allene-1.1-d2 molecule, and hence to calculate the expected rotational structure of the perpendicular bending vibrations of this molecule; the good agreement obtained with the observed spectra is a check on our results.

Farmers' attitudes towards techniques for improving oestrus detection in dairy herds in South West England

Unidentified heats contribute to declining fertility rates in English dairy herds. Several techniques have been advocated to improve heat detection rates. Despite demonstrable technical efficacy and cost-effectiveness, uptake is low. A study in South West England used the Theory of Reasoned Action (TORA) to explore dairy farmers' attitudes and beliefs towards heat detection techniques. Few farmers were convinced that following prescribed observation times, milk progesterone testing and using pedometers would fit their system or improve on their current heat detection practices. Perceived difficulty of using a technique was not a constraint on adoption. Without promotion that addresses identified barriers and drivers to adoption, little change in current practice can be expected. (c) 2006 Elsevier B.V. All rights reserved.

Measuring and comparing constraints to improved biosecurity amongst GB farmers, veterinarians and the auxiliary industries

Constraints to the introduction of enhanced biosecurity systems are rarely considered in sufficient detail when population medicine specialists initiate new control schemes. The main objective of our research was to investigate and compare the different attitudes constraining improvement in biosecurity for cattle and sheep farmers, practising veterinary surgeons and the auxiliary industries in Great Britain (GB). This study was carried out utilizing farmer focus groups, a questionnaire survey of veterinary practitioners and a telephone survey of auxiliary industry representatives. It appears that farmers and veterinarians have their own relatively clear definitions for biosecurity in relation to some major diseases threatening GB agriculture. Overall, farmers believe that other stakeholders, such as the government, should make a greater contribution towards biosecurity within GB. Conversely, veterinary practitioners saw their clients' ability or willingness to invest in biosecurity measures as a major constraint. Veterinary practitioners also felt that there was need for additional proof of efficacy and/or the potential economic benefits of proposed farm biosecurity practices better demonstrated. Auxiliary industries, in general, were not certain of their role in biosecurity although study participants highlighted zoonoses as part of the issue and offered that most of the constraints operated at farm level. (C) 2008 Elsevier B.V. All rights reserved.

Measuring and comparing constraints to improved biosecurity amongst GB farmers, veterinarians and the auxiliary industries

Constraints to the introduction of enhanced biosecurity systems are rarely considered in sufficient detail when population medicine specialists initiate new control schemes. The main objective of our research was to investigate and compare the different attitudes constraining improvement in biosecurity for cattle and sheep farmers, practising veterinary surgeons and the auxiliary industries in Great Britain (GB). This study was carried out utilizing farmer focus groups, a questionnaire survey of veterinary practitioners and a telephone survey of auxiliary industry representatives. It appears that farmers and veterinarians have their own relatively clear definitions for biosecurity in relation to some major diseases threatening GB agriculture. Overall, farmers believe that other stakeholders, such as the government, should make a greater contribution towards biosecurity within GB. Conversely, veterinary practitioners saw their clients' ability or willingness to invest in biosecurity measures as a major constraint. Veterinary practitioners also felt that there was need for additional proof of efficacy and/or the potential economic benefits of proposed farm biosecurity practices better demonstrated. Auxiliary industries, in general, were not certain of their role in biosecurity although study participants highlighted zoonoses as part of the issue and offered that most of the constraints operated at farm level. (C) 2008 Elsevier B.V. All rights reserved.

The state of veterinary epidemiology and economics

Recent developments in the fields of veterinary epidemiology and economics are critically reviewed and assessed. The impacts of recent technological developments in diagnosis, genetic characterisation, data processing and statistical analysis are evaluated. It is concluded that the acquisition and availability of data remains the principal constraint to the application of available techniques in veterinary epidemiology and economics, especially at population level. As more commercial producers use computerised management systems, the availability of data for analysis within herds is improving. However, consistency of recording and diagnosis remains problematic. Recent trends to the development of national livestock databases intended to provide reassurance to consumers of the safety and traceability of livestock products are potentially valuable sources of data that could lead to much more effective application of veterinary epidemiology and economics. These opportunities will be greatly enhanced if data from different sources, such as movement recording, official animal health programmes, quality assurance schemes, production recording and breed societies can be integrated. However, in order to realise such integrated databases, it will be necessary to provide absolute control of user access to guarantee data security and confidentiality. The potential applications of integrated livestock databases in analysis, modelling, decision-support, and providing management information for veterinary services and livestock producers are discussed. (c) 2004 Elsevier B.V. All rights reserved.

The effect of wilting and soaking Eupatorium adenophorum on its digestibility in vitro and voluntary intake by goats

An experiment was conducted to determine what effect simple treatments might have on the voluntary intake by goats in Nepal of Eupatorium adenophorum, an invasive weed that is usually only consumed by goats to a very limited extent. Samples of E. adenophorum were collected and either untreated, soaked for 2 h or wilted for 2 h before being oven dried (60 degrees C) and ground. Soaking and wilting had little effect on the chemical composition of E. adenophorum, but did increase (P=0.036) its in vitro organic matter degradability, by approximately 8%. The short-term intake rate (STIR) of treated and untreated E. adenophorum was then estimated with eight goats. Soaking time (from 2 to 24 h) was not related to STIR (r = -0.111, P=0.198), but the time E. adenophorum was left to wilt (from 2 to 48h), was positively related to STIR (r=0.521, P<0.001), with values of STIR (g dry matter/min kg goat liveweight(0.75)) being 0.405, 0.649,1.058, S.E.M. 0.088 for E. adenophorum, that had been wilted for 0, 24 and 48 h respectively (P<0.001). Liveweight change of goats and voluntary intake of E. adenophorum by goats was then estimated with 24 goats. E. adenophorum was fed either unwilted, or wilted for 24 or 48 h. It was fed as the sole forage or as a 3:1 mixture (dry matter basis) with Ficus cunia. There was a linear (P<0.001) and quadratic (P<0.01) increase in the intake of total forage and E. adenophorum with wilting time of E. adenophorum. Offering Ficus cunia increased total forage intake, but decreased E. adenophorum intake (P<0.05). After four weeks, there was virtually no change in goat liveweight and no significant difference between treatments. The results suggest that wilting E adenophorum for 24 h could increase its intake by goats, and thereby increase its usefulness, as a potential source of forage in the dry season of Nepal, when forage scarcity is a common constraint to livestock production. (C) 2007 Elsevier B.V. All rights reserved.