Measures of observation impact in non-Gaussian data assimilation

ABSTRACT Non-Gaussian/non-linear data assimilation is becoming an increasingly important area of research in the Geosciences as the resolution and non-linearity of models are increased and more and more non-linear observation operators are being used. In this study, we look at the effect of relaxing the assumption of a Gaussian prior on the impact of observations within the data assimilation system. Three different measures of observation impact are studied: the sensitivity of the posterior mean to the observations, mutual information and relative entropy. The sensitivity of the posterior mean is derived analytically when the prior is modelled by a simplified Gaussian mixture and the observation errors are Gaussian. It is found that the sensitivity is a strong function of the value of the observation and proportional to the posterior variance. Similarly, relative entropy is found to be a strong function of the value of the observation. However, the errors in estimating these two measures using a Gaussian approximation to the prior can differ significantly. This hampers conclusions about the effect of the non-Gaussian prior on observation impact. Mutual information does not depend on the value of the observation and is seen to be close to its Gaussian approximation. These findings are illustrated with the particle filter applied to the Lorenz ’63 system. This article is concluded with a discussion of the appropriateness of these measures of observation impact for different situations.

Dietary supplementation with Bifidobacterium lactis NCC2818 from weaning reduces local immunoglobulin production in lymphoid-associated tissues but increases systemic antibodies in healthy neonates

Weaning is associated with a major shift in the microbial community of the intestine, and this instability may make it more acquiescent than the adult microbiota to long-term changes. Modulation achieved through dietary interventions may have potentially beneficial effects on the developing immune system, which is driven primarily by the microbiota. The specific aim of the present study was to determine whether immune development could be modified by dietary supplementation with the human probiotic Bifidobacterium lactis NCC2818 in a tractable model of weaning in infants. Piglets were reared by their mothers before being weaned onto a solid diet supplemented with B. lactis NCC2818, while sibling controls did not receive supplementation. Probiotic supplementation resulted in a reduction in IgA (P,0·0005) and IgM (P,0·009) production by mucosal tissues but had no effect on IgG production (P.0·05). Probiotic- supplemented pigs had more mast cells than unsupplemented littermates (P,0·0001), although numbers in both groups were low. In addition, the supplemented piglets made stronger serum IgG responses to fed and injected antigens (P,0·05). The present findings are consistent with B. lactis NCC2818 reducing intestinal permeability induced by weaning, and suggest that the piglet is a valuable intermediate between rodent models and human infants. The results also strongly suggest that measures of the effect of probiotic supplementation on the immune system need to be interpreted carefully as proxy measures of health benefit. However, they are useful in developing an understanding of the mechanism of action of probiotic strains, an important factor in predicting favourable health outcomes of nutritional intervention.

The Marlborough mound, Wiltshire. A further Neolithic monumental mound by the River Kennet.

Recent radiocarbon dates obtained from two soil cores taken through the Marlborough Castle mound, Wiltshire, show the main body of it to be a contemporaneous monument to Silbury Hill, dating to the second half of the 3rd millennium cal bc. In light of these dates, this paper considers the sequence identified within the cores, which includes two possible flood events early in the construction of the mound. It also describes four cores taken through the surrounding ditch, as well as small-scale work to the north-east of the mound. The topographic location of the mound in a low-lying area and close to rivers and springs is discussed, and the potential for Late Neolithic sites nearby is set out, with the land to the south of the mound identified as an area for future research. The paper ends with the prospect that other apparent mottes in Wiltshire and beyond may well also have prehistoric origins

Priority research questions for the UK food system

The rise of food security up international political, societal and academic agendas has led to increasing interest in novel means of improving primary food production and reducing waste. There are however, also many ‘post-farm gate’ activities that are critical to food security, including processing, packaging, distributing, retailing, cooking and consuming. These activities all affect a range of important food security elements, notably availability, affordability and other aspects of access, nutrition and safety. Addressing the challenge of universal food security, in the context of a number of other policy goals (e.g. social, economic and environmental sustainability), is of keen interest to a range of UK stakeholders but requires an up-to-date evidence base and continuous innovation. An exercise was therefore conducted, under the auspices of the UK Global Food Security Programme, to identify priority research questions with a focus on the UK food system (though the outcomes may be broadly applicable to other developed nations). Emphasis was placed on incorporating a wide range of perspectives (‘world views’) from different stakeholder groups: policy, private sector, non-governmental organisations, advocacy groups and academia. A total of 456 individuals submitted 820 questions from which 100 were selected by a process of online voting and a three-stage workshop voting exercise. These 100 final questions were sorted into 10 themes and the ‘top’ question for each theme identified by a further voting exercise. This step also allowed four different stakeholder groups to select the top 7–8 questions from their perspectives. Results of these voting exercises are presented. It is clear from the wide range of questions prioritised in this exercise that the different stakeholder groups identified specific research needs on a range of post-farm gate activities and food security outcomes. Evidence needs related to food affordability, nutrition and food safety (all key elements of food security) featured highly in the exercise. While there were some questions relating to climate impacts on production, other important topics for food security (e.g. trade, transport, preference and cultural needs) were not viewed as strongly by the participants.

Observation impact in data assimilation: the effect of non-Gaussian observation error.

Data assimilation methods which avoid the assumption of Gaussian error statistics are being developed for geoscience applications. We investigate how the relaxation of the Gaussian assumption affects the impact observations have within the assimilation process. The effect of non-Gaussian observation error (described by the likelihood) is compared to previously published work studying the effect of a non-Gaussian prior. The observation impact is measured in three ways: the sensitivity of the analysis to the observations, the mutual information, and the relative entropy. These three measures have all been studied in the case of Gaussian data assimilation and, in this case, have a known analytical form. It is shown that the analysis sensitivity can also be derived analytically when at least one of the prior or likelihood is Gaussian. This derivation shows an interesting asymmetry in the relationship between analysis sensitivity and analysis error covariance when the two different sources of non-Gaussian structure are considered (likelihood vs. prior). This is illustrated for a simple scalar case and used to infer the effect of the non-Gaussian structure on mutual information and relative entropy, which are more natural choices of metric in non-Gaussian data assimilation. It is concluded that approximating non-Gaussian error distributions as Gaussian can give significantly erroneous estimates of observation impact. The degree of the error depends not only on the nature of the non-Gaussian structure, but also on the metric used to measure the observation impact and the source of the non-Gaussian structure.

Challenges in quantifying changes in the global water cycle

Understanding observed changes to the global water cycle is key to predicting future climate changes and their impacts. While many datasets document crucial variables such as precipitation, ocean salinity, runoff, and humidity, most are uncertain for determining long-term changes. In situ networks provide long time-series over land but are sparse in many regions, particularly the tropics. Satellite and reanalysis datasets provide global coverage, but their long-term stability is lacking. However, comparisons of changes among related variables can give insights into the robustness of observed changes. For example, ocean salinity, interpreted with an understanding of ocean processes, can help cross-validate precipitation. Observational evidence for human influences on the water cycle is emerging, but uncertainties resulting from internal variability and observational errors are too large to determine whether the observed and simulated changes are consistent. Improvements to the in situ and satellite observing networks that monitor the changing water cycle are required, yet continued data coverage is threatened by funding reductions. Uncertainty both in the role of anthropogenic aerosols, and due to large climate variability presently limits confidence in attribution of observed changes.

The science case for an orbital mission to Uranus: Exploring the origins and evolution of ice giant planets

Giant planets helped to shape the conditions we see in the Solar System today and they account for more than 99% of the mass of the Sun’s planetary system. They can be subdivided into the Ice Giants (Uranus and Neptune) and the Gas Giants (Jupiter and Saturn), which differ from each other in a number of fundamental ways. Uranus, in particular is the most challenging to our understanding of planetary formation and evolution, with its large obliquity, low self-luminosity, highly asymmetrical internal field, and puzzling internal structure. Uranus also has a rich planetary system consisting of a system of inner natural satellites and complex ring system, five major natural icy satellites, a system of irregular moons with varied dynamical histories, and a highly asymmetrical magnetosphere. Voyager 2 is the only spacecraft to have explored Uranus, with a flyby in 1986, and no mission is currently planned to this enigmatic system. However, a mission to the uranian system would open a new window on the origin and evolution of the Solar System and would provide crucial information on a wide variety of physicochemical processes in our Solar System. These have clear implications for understanding exoplanetary systems. In this paper we describe the science case for an orbital mission to Uranus with an atmospheric entry probe to sample the composition and atmospheric physics in Uranus’ atmosphere. The characteristics of such an orbiter and a strawman scientific payload are described and we discuss the technical challenges for such a mission. This paper is based on a white paper submitted to the European Space Agency’s call for science themes for its large-class mission programme in 2013.

Exploring strategies for coupled 4D-Var data assimilation using an idealised atmosphere-ocean model

Operational forecasting centres are currently developing data assimilation systems for coupled atmosphere-ocean models. Strongly coupled assimilation, in which a single assimilation system is applied to a coupled model, presents significant technical and scientific challenges. Hence weakly coupled assimilation systems are being developed as a first step, in which the coupled model is used to compare the current state estimate with observations, but corrections to the atmosphere and ocean initial conditions are then calculated independently. In this paper we provide a comprehensive description of the different coupled assimilation methodologies in the context of four dimensional variational assimilation (4D-Var) and use an idealised framework to assess the expected benefits of moving towards coupled data assimilation. We implement an incremental 4D-Var system within an idealised single column atmosphere-ocean model. The system has the capability to run both strongly and weakly coupled assimilations as well as uncoupled atmosphere or ocean only assimilations, thus allowing a systematic comparison of the different strategies for treating the coupled data assimilation problem. We present results from a series of identical twin experiments devised to investigate the behaviour and sensitivities of the different approaches. Overall, our study demonstrates the potential benefits that may be expected from coupled data assimilation. When compared to uncoupled initialisation, coupled assimilation is able to produce more balanced initial analysis fields, thus reducing initialisation shock and its impact on the subsequent forecast. Single observation experiments demonstrate how coupled assimilation systems are able to pass information between the atmosphere and ocean and therefore use near-surface data to greater effect. We show that much of this benefit may also be gained from a weakly coupled assimilation system, but that this can be sensitive to the parameters used in the assimilation.

Selective environmental stress caused by magmatic sulfur emissions from continental flood basalt eruptions

Several biotic crises during the past 300 million years have been linked to episodes of continental flood basalt volcanism, and in particular to the release of massive quantities of magmatic sulphur gas species. Flood basalt provinces were typically formed by numerous individual eruptions, each lasting years to decades. However, the environmental impact of these eruptions may have been limited by the occurrence of quiescent periods that lasted hundreds to thousands of years. Here we use a global aerosol model to quantify the sulphur-induced environmental effects of individual, decade-long flood basalt eruptions representative of the Columbia River Basalt Group, 16.5–14.5 million years ago, and the Deccan Traps, 65 million years ago. For a decade-long eruption of Deccan scale, we calculate a decadal-mean reduction in global surface temperature of 4.5 K, which would recover within 50 years after an eruption ceased unless climate feedbacks were very different in deep-time climates. Acid mists and fogs could have caused immediate damage to vegetation in some regions, but acid-sensitive land and marine ecosystems were well-buffered against volcanic sulphur deposition effects even during century-long eruptions. We conclude that magmatic sulphur from flood basalt eruptions would have caused a biotic crisis only if eruption frequencies and lava discharge rates had been high and sustained for several centuries at a time.

An idealized study of coupled atmosphere-ocean 4D-Var in the presence of model error

Atmosphere only and ocean only variational data assimilation (DA) schemes are able to use window lengths that are optimal for the error growth rate, non-linearity and observation density of the respective systems. Typical window lengths are 6-12 hours for the atmosphere and 2-10 days for the ocean. However, in the implementation of coupled DA schemes it has been necessary to match the window length of the ocean to that of the atmosphere, which may potentially sacrifice the accuracy of the ocean analysis in order to provide a more balanced coupled state. This paper investigates how extending the window length in the presence of model error affects both the analysis of the coupled state and the initialized forecast when using coupled DA with differing degrees of coupling. Results are illustrated using an idealized single column model of the coupled atmosphere-ocean system. It is found that the analysis error from an uncoupled DA scheme can be smaller than that from a coupled analysis at the initial time, due to faster error growth in the coupled system. However, this does not necessarily lead to a more accurate forecast due to imbalances in the coupled state. Instead coupled DA is more able to update the initial state to reduce the impact of the model error on the accuracy of the forecast. The effect of model error is potentially most detrimental in the weakly coupled formulation due to the inconsistency between the coupled model used in the outer loop and uncoupled models used in the inner loop.

Elevation-dependent warming in mountain regions of the world

There is growing evidence that the rate of warming is amplified with elevation, such that high-mountain environments experience more rapid changes in temperature than environments at lower elevations. Elevation-dependent warming (EDW) can accelerate the rate of change in mountain ecosystems, cryospheric systems, hydrological regimes and biodiversity. Here we review important mechanisms that contribute towards EDW: snow albedo and surface-based feedbacks; water vapour changes and latent heat release; surface water vapour and radiative flux changes; surface heat loss and temperature change; and aerosols. All lead to enhanced warming with elevation (or at a critical elevation), and it is believed that combinations of these mechanisms may account for contrasting regional patterns of EDW. We discuss future needs to increase knowledge of mountain temperature trends and their controlling mechanisms through improved observations, satellite-based remote sensing and model simulations.

Comparação dos achados histológicos do testículo intra-abdominal antes e após a divisão dos vasos espermáticos

Introdução: a videolaparoscopia proporciona adequado diagnóstico do testículo impalpável e difundiu a orquidopexia estagiada pela técnica de Fowler- Stephens; entretanto, não há relatos de que a divisão dos vasos espermáticos poderia ocasionar alterações histológicas no testículo intra-abdominal. O objetivo do presente trabalho foi avaliar a viabilidade e comparar a histologia e o volume dos testículos intraabdominais antes e após a divisão dos vasos espermáticos. Métodos: foram avaliados 44 testículos de 35 pacientes com idades variando de 4 a 168 meses com testículos intra-abdominais, que foram submetidos a videolaparoscopia para diagnóstico e tratamento. Foi realizada biópsia e medição do volume do testículo antes e após a ligadura e divisão dos vasos espermáticos, com intervalo de seis meses, pela técnica de Fowler-Stephens. A análise volumétrica e a histologia dos testículos foram comparadas. Resultados: 40% dos testículos impalpáveis se localizaram na região intraabdominal; 97,7% dos testículos permaneceram viáveis. Os achados histológicos e o volume testicular antes e após a divisão dos vasos espermáticos não demonstraram diferença estatisticamente significativa (p > 0,05). Conclusões: o estudo foi eficaz para avaliar a viabilidade dos testículos. A divisão dos vasos espermáticos não ocasionou alterações no volume e nos achados histológicos dos testículos intra-abdominais. Unitermos: criptorquia, testículo intra-abdominal, videolaparoscopia, histologia testicular.