Efficient non-linear data assimilation in geophysical fluid dynamics

New ways of combining observations with numerical models are discussed in which the size of the state space can be very large, and the model can be highly nonlinear. Also the observations of the system can be related to the model variables in highly nonlinear ways, making this data-assimilation (or inverse) problem highly nonlinear. First we discuss the connection between data assimilation and inverse problems, including regularization. We explore the choice of proposal density in a Particle Filter and show how the ’curse of dimensionality’ might be beaten. In the standard Particle Filter ensembles of model runs are propagated forward in time until observations are encountered, rendering it a pure Monte-Carlo method. In large-dimensional systems this is very inefficient and very large numbers of model runs are needed to solve the data-assimilation problem realistically. In our approach we steer all model runs towards the observations resulting in a much more efficient method. By further ’ensuring almost equal weight’ we avoid performing model runs that are useless in the end. Results are shown for the 40 and 1000 dimensional Lorenz 1995 model.

Continental shift? An analysis of convergence trends in European real estate equities

European economic and political integration have been recognised as having implications for patterns of performance in national real estate and capital markets and have generated a wide body of research and commentary. In 1999, progress towards monetary integration within the European Union culminated in the introduction of a common currency and monetary policy. This paper investigates the effects of this ‘event’ on the behaviour of stock returns in European real estate companies. A range of statistical tests is applied to the performance of European property companies to test for changes in segmentation, co-movement and causality. The results suggest that, relative to the wider equity markets, the dispersion of performance is higher, correlations are lower, a common contemporaneous factor has much lower explanatory power whilst lead-lag relationships are stronger. Consequently, the evidence of transmission of monetary integration to real estate securities is less noticeable than to general securities. Less and slower integration is attributed to the relatively small size of the real estate securities market and the local and national nature of the majority of the companies’ portfolios.

Can carbon surface oxidation shift the pore size distribution curve calculated from Ar, N(2) and CO(2) adsorption isotherms? Simulation results for a realistic carbon model

Using the virtual porous carbon model proposed by Harris et al, we study the effect of carbon surface oxidation on the pore size distribution (PSD) curve determined from simulated Ar, N(2) and CO(2) isotherms. It is assumed that surface oxidation is not destructive for the carbon skeleton, and that all pores are accessible for studied molecules (i.e., only the effect of the change of surface chemical composition is studied). The results obtained show two important things, i.e., oxidation of the carbon surface very slightly changes the absolute porosity (calculated from the geometric method of Bhattacharya and Gubbins (BG)); however, PSD curves calculated from simulated isotherms are to a greater or lesser extent affected by the presence of surface oxides. The most reliable results are obtained from Ar adsorption data. Not only is adsorption of this adsorbate practically independent from the presence of surface oxides, but, more importantly, for this molecule one can apply the slit-like model of pores as the first approach to recover the average pore diameter of a real carbon structure. For nitrogen, the effect of carbon surface chemical composition is observed due to the quadrupole moment of this molecule, and this effect shifts the PSD curves compared to Ar. The largest differences are seen for CO2, and it is clearly demonstrated that the PSD curves obtained from adsorption isotherms of this molecule contain artificial peaks and the average pore diameter is strongly influenced by the presence of electrostatic adsorbate-adsorbate as well as adsorbate-adsorbent interactions.

Viscous coupling of shear-free turbulence across nearly flat fluid interfaces

The interactions between shear-free turbulence in two regions (denoted as + and − on either side of a nearly flat horizontal interface are shown here to be controlled by several mechanisms, which depend on the magnitudes of the ratios of the densities, ρ+/ρ−, and kinematic viscosities of the fluids, μ+/μ−, and the root mean square (r.m.s.) velocities of the turbulence, u0+/u0−, above and below the interface. This study focuses on gas–liquid interfaces so that ρ+/ρ− ≪ 1 and also on where turbulence is generated either above or below the interface so that u0+/u0− is either very large or very small. It is assumed that vertical buoyancy forces across the interface are much larger than internal forces so that the interface is nearly flat, and coupling between turbulence on either side of the interface is determined by viscous stresses. A formal linearized rapid-distortion analysis with viscous effects is developed by extending the previous study by Hunt & Graham (J. Fluid Mech., vol. 84, 1978, pp. 209–235) of shear-free turbulence near rigid plane boundaries. The physical processes accounted for in our model include both the blocking effect of the interface on normal components of the turbulence and the viscous coupling of the horizontal field across thin interfacial viscous boundary layers. The horizontal divergence in the perturbation velocity field in the viscous layer drives weak inviscid irrotational velocity fluctuations outside the viscous boundary layers in a mechanism analogous to Ekman pumping. The analysis shows the following. (i) The blocking effects are similar to those near rigid boundaries on each side of the interface, but through the action of the thin viscous layers above and below the interface, the horizontal and vertical velocity components differ from those near a rigid surface and are correlated or anti-correlated respectively. (ii) Because of the growth of the viscous layers on either side of the interface, the ratio uI/u0, where uI is the r.m.s. of the interfacial velocity fluctuations and u0 the r.m.s. of the homogeneous turbulence far from the interface, does not vary with time. If the turbulence is driven in the lower layer with ρ+/ρ− ≪ 1 and u0+/u0− ≪ 1, then uI/u0− ~ 1 when Re (=u0−L−/ν−) ≫ 1 and R = (ρ−/ρ+)(v−/v+)1/2 ≫ 1. If the turbulence is driven in the upper layer with ρ+/ρ− ≪ 1 and u0+/u0− ≫ 1, then uI/u0+ ~ 1/(1 + R). (iii) Nonlinear effects become significant over periods greater than Lagrangian time scales. When turbulence is generated in the lower layer, and the Reynolds number is high enough, motions in the upper viscous layer are turbulent. The horizontal vorticity tends to decrease, and the vertical vorticity of the eddies dominates their asymptotic structure. When turbulence is generated in the upper layer, and the Reynolds number is less than about 106–107, the fluctuations in the viscous layer do not become turbulent. Nonlinear processes at the interface increase the ratio uI/u0+ for sheared or shear-free turbulence in the gas above its linear value of uI/u0+ ~ 1/(1 + R) to (ρ+/ρ−)1/2 ~ 1/30 for air–water interfaces. This estimate agrees with the direct numerical simulation results from Lombardi, De Angelis & Bannerjee (Phys. Fluids, vol. 8, no. 6, 1996, pp. 1643–1665). Because the linear viscous–inertial coupling mechanism is still significant, the eddy motions on either side of the interface have a similar horizontal structure, although their vertical structure differs.

Weaning diet induces sustained metabolic phenotype shift in the pig and influences host response to Bifidobacterium lactis NCC2818

Background The process of weaning causes a major shift in intestinal microbiota and is a critical period for developing appropriate immune responses in young mammals.Objective To use a new systems approach to provide an overview of host metabolism and the developing immune system in response to nutritional intervention around the weaning period.Design Piglets (n=14) were weaned onto either an egg-based or soya-based diet at 3 weeks until 7 weeks, when all piglets were switched onto a fish-based diet. Half the animals on each weaning diet received Bifidobacterium lactis NCC2818 supplementation from weaning onwards. Immunoglobulin production from immunologically relevant intestinal sites was quantified and the urinary (1)H NMR metabolic profile was obtained from each animal at post mortem (11 weeks).Results Different weaning diets induced divergent and sustained shifts in the metabolic phenotype, which resulted in the alteration of urinary gut microbial co-metabolites, even after 4 weeks of dietary standardisation. B lactis NCC2818 supplementation affected the systemic metabolism of the different weaning diet groups over and above the effects of diet. Additionally, production of gut mucosa-associated IgA and IgM was found to depend upon the weaning diet and on B lactis NCC2818 supplementation.ConclusionThe correlation of urinary (1)H NMR metabolic profile with mucosal immunoglobulin production was demonstrated, thus confirming the value of this multi-platform approach in uncovering non-invasive biomarkers of immunity. This has clear potential for translation into human healthcare with the development of urine testing as a means of assessing mucosal immune status. This might lead to early diagnosis of intestinal dysbiosis and with subsequent intervention, arrest disease development. This system enhances our overall understanding of pathologies under supra-organismal control.

Elevated levels of corticotrophin-releasing factor binding protein in the blood of patients suffering from arthritis and septicaemia and the presence of novel ligands in synovial fluid.

In view of the reported inflammatory effects of corticotrophin-releasing factor (CRF) and the associated regulatory elements in the gene of its binding protein (BP), we postulate that both BP as well as novel BP-ligands other than CRF may be involved in inflammatory disease. We have investigated BP in the blood of patients with arthritis and septicaemia and have attempted to identify CRF and other BP-ligands in synovial fluid. The BP was found to be significantly elevated in the blood of patients with rheumatoid arthritis and septicaemia. There was less BP-ligand and CRF in synovial fluid from patients with rheumatoid arthritis that from those with osteo- or psoriatic arthritis. There was at least 10-fold more BP-ligand than CRF in the fluid of all three groups of patients. A small amount of immunoreactive human (h)CRF, eluting in the expected position of CRF-41, was detected after high-pressure liquid chromatography of arthritic synovial fluid; however, the bulk of material with BP-ligand binding activity eluted earlier, suggesting that synovial fluid contained novel peptides that interacted with the BP. These results would suggest that the BP and its ligands could play an endocrine immunomodulatory role in inflammatory disease.

Liquid-vapour homogenisation of fluid inclusions in stalagmites: Evaluation of a new thermometer for paleoclimate research

We present a new approach to determine palaeotemperatures (mean annual surface temperatures) based on measurements of the liquid–vapour homogenisation temperature of fluid inclusions in stalagmites. The aim of this study is to explore the potential and the limitations of this new palaeothermometer and to develop a reliable methodology for routine applications in palaeoclimate research. Therefore, we have investigated recent fluid inclusions from the top part of actively growing stalagmites that have formed at temperatures close to the present-day cave air temperature. A precondition for measuring homogenisation temperatures of originally monophase inclusions is the nucleation of a vapour bubble by means of single ultra-short laser pulses. Based on the observed homogenisation temperatures (Th(obs)) and measurements of the vapour bubble diameter at a known temperature, we calculated stalagmite formation temperatures (Tf) by applying a thermodynamic model that takes into account the effect of surface tension on liquid–vapour homogenisation. Results from recent stalagmite samples demonstrate that calculated stalagmite formation temperatures match the present-day cave air temperature within ± 0.2 °C. To avoid artificially induced changes of the fluid density we defined specific demands on the selection, handling and preparation of the stalagmite samples. Application of the method is restricted to stalagmites that formed at cave temperatures greater than ~ 9–11 °C.

Determination of Holocene cave temperatures from Kr and Xe concentrations in stalagmite fluid inclusions

From the concentrations of dissolved atmospheric noble gases in water, a so-called “noble gas temperature” (NGT) can be determined that corresponds to the temperature of the water when it was last in contact with the atmosphere. Here we demonstrate that the NGT concept is applicable to water inclusions in cave stalagmites, and yields NGTs that are in good agreement with the ambient air temperatures in the caves. We analysed samples from two Holocene and one undated stalagmite. The three stalagmites originate from three caves located in different climatic regions having modern mean annual air temperatures of 27 °C, 12 °C and 8 °C, respectively. In about half of the samples analysed Kr and Xe concentrations originated entirely from the two well-defined noble gas components air-saturated water and atmospheric air, which allowed NGTs to be determined successfully from Kr and Xe concentrations. One stalagmite seems to be particularly suitable for NGT determination, as almost all of its samples yielded the modern cave temperature. Notably, this stalagmite contains a high proportion of primary water inclusions, which seem to preserve the temperature-dependent signature well in their Kr and Xe concentrations. In future work on stalagmites detailed microscopic inspection of the fluid inclusions prior to noble gas analysis is therefore likely to be crucial in increasing the number of successful NGT determinations.

Borrowing and shift-induced interference: contrasting patterns in French-Germanic contact in Brussels and Strasbourg

The main aim of the present article is to test hypotheses derived from the model for contact- induced language change as formulated in Thomason and Kaufman (1988 et seq.). As the model correctly predicts the asymmetries between the mutual influences of the Germanic and the Romance varieties in Brussels and Strasbourg it is a very powerful tool for describing the contact patterns in these cities. The analysis shows that the contact patterns are very similar, both from a quantitative and from a qualitative point of view, despite important differences in the sociolinguistic situation of both cities. The striking similarities in the outcome of language contact seem to find a plausible explanation in the fact that the language contact situations in both cities are similar from a typological point of view: in each city a variety of French is in contact with a Germanic variety (Alsatian and Brussels Dutch). Thus, the claim of the present article is that the structure of the languages plays a more prominent role in the outcome of language contact than the sociolinguistic history of the speakers.

Atlantic Ocean influence on a shift in European climate in the 1990s

European climate exhibits variability on a wide range of timescales. Understanding the nature and drivers of this variability is an essential step in developing robust climate predictions and risk assessments. The Atlantic Ocean has been suggested as an important driver of variability in European climate on decadal timescales1, but the importance of this influence in recent decades has been unclear, partly because of difficulties in separating the influence of the Atlantic Ocean from other contributions, for example, from the tropical Pacific Ocean and the stratosphere. Here we analyse four data sets derived from observations to show that, during the 1990s, there was a substantial shift in European climate towards a pattern characterized by anomalously wet summers in northern Europe, and hot, dry, summers in southern Europe, with related shifts in spring and autumn. These changes in climate coincided with a substantial warming of the North Atlantic Ocean, towards a state last seen in the 1950s. The patterns of European climate change in the 1990s are consistent with earlier changes attributed to the influence of the North Atlantic Ocean, and provide compelling evidence that the Atlantic Ocean was the key driver. Our results suggest that the recent pattern of anomalies in European climate will persist as long as the North Atlantic Ocean remains anomalously warm.

Accurate analysis of noble gas concentrations in small water samples and its application to fluid inclusions in stalagmites

The concentrations of dissolved noble gases in water are widely used as a climate proxy to determine noble gas temperatures (NGTs); i.e., the temperature of the water when gas exchange last occurred. In this paper we make a step forward to apply this principle to fluid inclusions in stalagmites in order to reconstruct the cave temperature prevailing at the time when the inclusion was formed. We present an analytical protocol that allows us accurately to determine noble gas concentrations and isotope ratios in stalagmites, and which includes a precise manometrical determination of the mass of water liberated from fluid inclusions. Most important for NGT determination is to reduce the amount of noble gases liberated from air inclusions, as they mask the temperature-dependent noble gas signal from the water inclusions. We demonstrate that offline pre-crushing in air to subsequently extract noble gases and water from the samples by heating is appropriate to separate gases released from air and water inclusions. Although a large fraction of recent samples analysed by this technique yields NGTs close to present-day cave temperatures, the interpretation of measured noble gas concentrations in terms of NGTs is not yet feasible using the available least squares fitting models. This is because the noble gas concentrations in stalagmites are not only composed of the two components air and air saturated water (ASW), which these models are able to account for. The observed enrichments in heavy noble gases are interpreted as being due to adsorption during sample preparation in air, whereas the excess in He and Ne is interpreted as an additional noble gas component that is bound in voids in the crystallographic structure of the calcite crystals. As a consequence of our study's findings, NGTs will have to be determined in the future using the concentrations of Ar, Kr and Xe only. This needs to be achieved by further optimizing the sample preparation to minimize atmospheric contamination and to further reduce the amount of noble gases released from air inclusions.

The unsteady flow of a weakly compressible fluid in a thin porous layer II: three-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a three-dimensional layer, composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Numerical solution of this three-dimensional evolution problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l, a situation which occurs frequently in the application to oil and gas reservoir recovery and which leads to significant stiffness in the numerical problem. Under the assumption that $\epsilon\propto h/l\ll 1$, we show that, to leading order in $\epsilon$, the pressure field varies only in the horizontal directions away from the wells (the outer region). We construct asymptotic expansions in $\epsilon$ in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive expressions for all significant process quantities. The only computations required are for the solution of non-stiff linear, elliptic, two-dimensional boundary-value, and eigenvalue problems. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the layer, $\epsilon$, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighbourhood of wells and away from wells.

The unsteady flow of a weakly compressible fluid in a thin porous layer III: three-dimensional computations

We describe a novel method for determining the pressure and velocity fields for a weakly compressible fluid flowing in a thin three-dimensional layer composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Our approach uses the method of matched asymptotic expansions to derive expressions for all significant process quantities, the computation of which requires only the solution of linear, elliptic, two-dimensional boundary value and eigenvalue problems. In this article, we provide full implementation details and present numerical results demonstrating the efficiency and accuracy of our scheme.

Hamiltonian geophysical fluid dynamics

Hamiltonian dynamics describes the evolution of conservative physical systems. Originally developed as a generalization of Newtonian mechanics, describing gravitationally driven motion from the simple pendulum to celestial mechanics, it also applies to such diverse areas of physics as quantum mechanics, quantum field theory, statistical mechanics, electromagnetism, and optics – in short, to any physical system for which dissipation is negligible. Dynamical meteorology consists of the fundamental laws of physics, including Newton’s second law. For many purposes, diabatic and viscous processes can be neglected and the equations are then conservative. (For example, in idealized modeling studies, dissipation is often only present for numerical reasons and is kept as small as possible.) In such cases dynamical meteorology obeys Hamiltonian dynamics. Even when nonconservative processes are not negligible, it often turns out that separate analysis of the conservative dynamics, which fully describes the nonlinear interactions, is essential for an understanding of the complete system, and the Hamiltonian description can play a useful role in this respect. Energy budgets and momentum transfer by waves are but two examples.