Advances in sequential data assimilation and numerical weather forecasting: an Ensemble Transform Kalman-Bucy Filter, a study on clustering in deterministic Ensemble Square Root Filters, and a test of a new time stepping scheme in an atmospheric model

This dissertation deals with aspects of sequential data assimilation (in particular ensemble Kalman filtering) and numerical weather forecasting. In the first part, the recently formulated Ensemble Kalman-Bucy (EnKBF) filter is revisited. It is shown that the previously used numerical integration scheme fails when the magnitude of the background error covariance grows beyond that of the observational error covariance in the forecast window. Therefore, we present a suitable integration scheme that handles the stiffening of the differential equations involved and doesn’t represent further computational expense. Moreover, a transform-based alternative to the EnKBF is developed: under this scheme, the operations are performed in the ensemble space instead of in the state space. Advantages of this formulation are explained. For the first time, the EnKBF is implemented in an atmospheric model. The second part of this work deals with ensemble clustering, a phenomenon that arises when performing data assimilation using of deterministic ensemble square root filters in highly nonlinear forecast models. Namely, an M-member ensemble detaches into an outlier and a cluster of M-1 members. Previous works may suggest that this issue represents a failure of EnSRFs; this work dispels that notion. It is shown that ensemble clustering can be reverted also due to nonlinear processes, in particular the alternation between nonlinear expansion and compression of the ensemble for different regions of the attractor. Some EnSRFs that use random rotations have been developed to overcome this issue; these formulations are analyzed and their advantages and disadvantages with respect to common EnSRFs are discussed. The third and last part contains the implementation of the Robert-Asselin-Williams (RAW) filter in an atmospheric model. The RAW filter is an improvement to the widely popular Robert-Asselin filter that successfully suppresses spurious computational waves while avoiding any distortion in the mean value of the function. Using statistical significance tests both at the local and field level, it is shown that the climatology of the SPEEDY model is not modified by the changed time stepping scheme; hence, no retuning of the parameterizations is required. It is found the accuracy of the medium-term forecasts is increased by using the RAW filter.

Ensemble clustering in deterministic ensemble Kalman filters

Ensemble clustering (EC) can arise in data assimilation with ensemble square root filters (EnSRFs) using non-linear models: an M-member ensemble splits into a single outlier and a cluster of M−1 members. The stochastic Ensemble Kalman Filter does not present this problem. Modifications to the EnSRFs by a periodic resampling of the ensemble through random rotations have been proposed to address it. We introduce a metric to quantify the presence of EC and present evidence to dispel the notion that EC leads to filter failure. Starting from a univariate model, we show that EC is not a permanent but transient phenomenon; it occurs intermittently in non-linear models. We perform a series of data assimilation experiments using a standard EnSRF and a modified EnSRF by a resampling though random rotations. The modified EnSRF thus alleviates issues associated with EC at the cost of traceability of individual ensemble trajectories and cannot use some of algorithms that enhance performance of standard EnSRF. In the non-linear regimes of low-dimensional models, the analysis root mean square error of the standard EnSRF slowly grows with ensemble size if the size is larger than the dimension of the model state. However, we do not observe this problem in a more complex model that uses an ensemble size much smaller than the dimension of the model state, along with inflation and localisation. Overall, we find that transient EC does not handicap the performance of the standard EnSRF.

Structural and morphological studies of the dipeptide based L-Pro-L-Val organocatalytic gels and their rheological behaviour

Organocatalytic gels based on the dipeptide sequence L-Pro-L-Val have been studied by two different FTIR techniques. This suggests a different arrangement of the gelator molecules in the self-assembled fibers depending on the organic solvent employed. In acetonitrile and nitromethane the structure of the supramolecular aggregates is similar and provides similar catalytic properties (supramolecularenhancement of basicity). In contrast, the self-assembled fibers obtained in toluene clearly presented a different molecular arrangement consistent with its different catalytic behaviour (enamine-based catalysis). In addition these gels have been studied by microscopy and rheology.

Faunal evidence of a Holocene pluvial phase in Southern Arabia with remarks on the morphological variability of Helenina anderseni

Although foraminifera have been found living in inland saline lakes isolated from the sea, this phenomenon has rarely been recognized in the fossil record. This study documents the occurrence of benthic foraminifera in Holocene lake sediments located nearly 500 km inland from the Red Sea, in the Al-Mundafan region of southern Saudi Arabia. The lake formed during a regional pluvial period, 10,500–6000 yr BP. The presence of foraminifera and brackish charophytes in the studied section represent an interval when the lake was slightly brackish due to high evaporation. The studied sediments yielded a bispecific benthic foraminiferal fauna comprised of Helenina anderseni and Trichohyalus aguayoi, as well as the brackish charophyte genus Lamprothamnium. The benthic foraminifera are species characteristic of mangrove swamps, salt marshes, and lagoons, which are environments currently widespread along the Red Sea coasts. Because the Al Mundafan area was never connected to the sea during the Quaternary, wading birds must have been the vector that transported the foraminifera to the paleolake.

Cooled infrared filters and dichroics for the sea and land surface temperature radiometer (SLSTR)

The Sea and Land Surface Temperature Radiometer (SLSTR) is a nine channel visible and infrared high precision radiometer designed to provide climate data of global sea and land surface temperatures. The SLSTR payload is destined to fly on the Ocean and Medium-Resolution Land Mission for the ESA/EU Global Monitoring for Environment and Security (GMES) Programme Sentinel-3 mission to measure the sea and land temperature and topography for near real-time environmental and atmospheric climate monitoring of the Earth. In this paper we describe the optical layout of infrared optics in the instrument, spectral thin-film multilayer design, and system channel throughput analysis for the combined interference filter and dichroic beamsplitter coatings to discriminate wavelengths at 3.74, 10.85 & 12.0 μm. The rationale for selection of thin-film materials, deposition technique, and environmental testing, inclusive of humidity, thermal cycling and ionizing radiation testing are also described.

Viruses in soils: morphological diversity and abundance in the rhizosphere

Soil viruses are potentially of great importance as they may influence the ecology and evolution of soil biological communities through both an ability to transfer genes from host to host and as a potential cause of microbial mortality. Despite this importance, the area of soil virology is understudied. Here, we report the isolation and preliminary characterisation of viruses from soils in the Dundee area of Scotland. Different virus morphotypes including tailed, polyhedral (spherical), rod shaped, filamentous and bacilliform particles were detected in the soil samples. An apparent predominance of small spherical and filamentous bacteriophages was observed, whereas tailed bacteriophages were significantly less abundant. In this report, we also present observations and characterisation of viruses from different soil functional domains surrounding wheat roots: rhizosheath, rhizosphere and bulk soil. In spite of the differences in abundance of bacterial communities in these domains, no significant variations in viral population structure in terms of morphology and abundance were found. Typically, there were approximately 1.1–1.2 × 109 virions g−1 dry weight, implicating remarkable differences in virus-to-bacteria ratios in domains close to roots, rhizosphere and rhizosheath (approximately 0.27) and in bulk soil (approximately 4.68).

Sometimes they use it, sometimes they don’t: an epistemological discussion of L2 morphological production and its use as a competence measurement

The main goal of all approaches to adult second language acquisition (SLA) is to accurately describe and explain the overall acquisition process. To accomplish this, SLA researchers must come to agree on some key issues. In this commentary, I defend the necessity of the competence/performance distinction and how this relates to why an examination of morphological production presents challenges for SLA research. I suggest that such a methodology is meaningful only when it is dovetailed with procedures that test for related syntactic/semantic knowledge.

Ensemble transform Kalman-Bucy filters

Two recent works have adapted the Kalman–Bucy filter into an ensemble setting. In the first formulation, the ensemble of perturbations is updated by the solution of an ordinary differential equation (ODE) in pseudo-time, while the mean is updated as in the standard Kalman filter. In the second formulation, the full ensemble is updated in the analysis step as the solution of single set of ODEs in pseudo-time. Neither requires matrix inversions except for the frequently diagonal observation error covariance. We analyse the behaviour of the ODEs involved in these formulations. We demonstrate that they stiffen for large magnitudes of the ratio of background error to observational error variance, and that using the integration scheme proposed in both formulations can lead to failure. A numerical integration scheme that is both stable and is not computationally expensive is proposed. We develop transform-based alternatives for these Bucy-type approaches so that the integrations are computed in ensemble space where the variables are weights (of dimension equal to the ensemble size) rather than model variables. Finally, the performance of our ensemble transform Kalman–Bucy implementations is evaluated using three models: the 3-variable Lorenz 1963 model, the 40-variable Lorenz 1996 model, and a medium complexity atmospheric general circulation model known as SPEEDY. The results from all three models are encouraging and warrant further exploration of these assimilation techniques.

A molecular, morphological and cytological investigation of the identity of non-native Ludwigia (Onagraceae) populations in Britain

A combined molecular, morphological and cytological analysis was used to study the identity and number of species of Ludwigia section Oligospermum occurring in British waterways. Only one taxon was identified for which the name L. grandiflora subsp. hexapetala (Hook. & Arn.) G.L. Nesom & Kartesz is preferred. A chromosome count of 2n = 80 was made for all plants tested and DNA evidence demonstrates that at least two clones are present in Britain. Morphological characters to differentiate L.grandiflora subsp. hexapetala and L. peploides subsp. montevidensis (Spreng.) P.H. Raven are provided. Though the production of fruit in Britain by apparently isolated populations is noted, repeated introduction into the wild from gardens is judged to be primarily responsible for the British distribution of the taxon. Legislative implications are considered and an amendment to Schedule 9 of the UK Wildlife and Countryside Act (1981) is strongly advocated.

The time-course of morphological constraints: Evidence from eye-movements during reading

Lexical compounds in English are constrained in that the non-head noun can be an irregular but not a regular plural (e.g. mice eater vs. *rats eater), a contrast that has been argued to derive from a morphological constraint on modifiers inside compounds. In addition, bare nouns are preferred over plural forms inside compounds (e.g. mouse eater vs. mice eater), a contrast that has been ascribed to the semantics of compounds. Measuring eyemovements during reading, this study examined how morphological and semantic information become available over time during the processing of a compound. We found that the morphological constraint affected both early and late eye-movement measures, whereas the semantic constraint for singular non-heads only affected late measures of processing. These results indicate that morphological information becomes available earlier than semantic information during the processing of compounds.

The time-course of morphological constraints: a study of plurals inside derived words

The avoidance of regular but not irregular plurals inside compounds (e.g. *rats eater vs. mice eater) has been one of the most widely studied morphological phenomena in the psycholinguistics literature. To examine whether the constraints that are responsible for this contrast have any general significance beyond compounding, we investigated derived word forms containing regular and irregular plurals in two experiments. Experiment 1 was an offline acceptability judgment task, and experiment 2 measured eye movements during reading derived words containing regular and irregular plurals and uninflected base nouns. The results from both experiments show that the constraint against regular plurals inside compounds generalizes to derived words. We argue that this constraint cannot be reduced to phonological properties, but is instead morphological in nature. The eye-movement data provide detailed information on the time-course of processing derived word forms indicating that early stages of processing are affected by a general constraint that disallows inflected words from feeding derivational processes, and that the more specific constraint against regular plurals comes in at a subsequent later stage of processing. We argue that these results are consistent with stage-based models of language processing.

The time-course of morphological processing in a second language

We report findings from psycholinguistic experiments investigating the detailed timing of processing morphologically complex words by proficient adult second (L2) language learners of English in comparison to adult native (L1) speakers of English. The first study employed the masked priming technique to investigate -ed forms with a group of advanced Arabic-speaking learners of English. The results replicate previously found L1/L2 differences in morphological priming, even though in the present experiment an extra temporal delay was offered after the presentation of the prime words. The second study examined the timing of constraints against inflected forms inside derived words in English using the eye-movement monitoring technique and an additional acceptability judgment task with highly advanced Dutch L2 learners of English in comparison to adult L1 English controls. Whilst offline the L2 learners performed native-like, the eye-movement data showed that their online processing was not affected by the morphological constraint against regular plurals inside derived words in the same way as in native speakers. Taken together, these findings indicate that L2 learners are not just slower than native speakers in processing morphologically complex words, but that the L2 comprehension system employs real-time grammatical analysis (in this case, morphological information) less than the L1 system.

Flow-induced morphological instability of a mushy layer

A morphological instability of a mushy layer due to a forced flow in the melt is analysed. The instability is caused by flow induced in the mushy layer by Bernoulli suction at the crests of a sinusoidally perturbed mush–melt interface. The flow in the mushy layer advects heat away from crests which promotes solidification. Two linear stability analyses are presented: the fundamental mechanism for instability is elucidated by considering the case of uniform flow of an inviscid melt; a more complete analysis is then presented for the case of a parallel shear flow of a viscous melt. The novel instability mechanism we analyse here is contrasted with that investigated by Gilpin et al. (1980) and is found to be more potent for the case of newly forming sea ice.

Palaeobiology of the medieval population of Albano (Rome, Italy): a combined morphological and biomolecular approach

Recent years have seen increased interest in skeletal populations from the Imperial Roman Age in Italy, but much less is known about diet and standards of living in the subsequent medieval period. To fill this gap, we conducted a morphological analysis of human remains from Albano, an Italian town near Rome, as well as a stable isotope analysis of bone collagen to reconstruct diet. The sample was recovered from a Medieval cemetery (1040–1220 cal. yr. BP) located in the gardens of the historical Palazzo Doria Pamphili in Albano. A minimum number of 40 individuals (31 adults and 9 sub-adults) were examined using standard methods. Though the general health status of the population was good, signs of cribra orbitalia and diffuse enthesopathies were noted during the morphological examination. Stable carbon and nitrogen isotope analyses of the bone collagen from 24 adult humans and three faunal bones indicate that the diet of the population may be described as predominantly terrestrial and C3-plant based although the data for some of the individuals suggest a modest consumption of C4-(millet) based or aquatic proteins. No evidence of significant dietary differences between the sexes was found. The comparison of the isotope data from Albano with those from populations recovered in the same region is consistent with a shift from a terrestrial, possibly plant foods-dominated subsistence in the Early Middle Ages to a diet with a higher contribution from animal proteins, both terrestrial and aquatic, in the Later Middle Ages.