Adding constrained discontinuities to Gaussian process models of wind fields

Gaussian Processes provide good prior models for spatial data, but can be too smooth. In many physical situations there are discontinuities along bounding surfaces, for example fronts in near-surface wind fields. We describe a modelling method for such a constrained discontinuity and demonstrate how to infer the model parameters in wind fields with MCMC sampling.

Adding constrained discontinuities to Gaussian process models of wind fields

Gaussian Processes provide good prior models for spatial data, but can be too smooth. In many physical situations there are discontinuities along bounding surfaces, for example fronts in near-surface wind fields. We describe a modelling method for such a constrained discontinuity and demonstrate how to infer the model parameters in wind fields with MCMC sampling.

A brief review of magnetic fields from the human visual system

Both the eye and brain generate magnetic fields when stimulated with a variety of visual cues. These magnetic fields can be measured with a magnetometer; a device which uses superconducting technology. The application of this technique to measuring the magnetooculogram, magnetoretinogram and visually evoked fields from the brain is described. So far the main use of this technique has been in pure research. Its potential for diagnosing ocular and neurological diseases is discussed.

Automated assessment of visual fields and their inter-relation to evoked potentials in visual disorders

The Octopus Automated Perimeter was validated in a comparative study and found to offer many advantages in the assessment of the visual field. The visual evoked potential was investigated in an extensive study using a variety of stimulus parameters to simulate hemianopia and central visual field defects. The scalp topography was recorded topographically and a technique to compute the source derivation of the scalp potential was developed. This enabled clarification of the expected scalp distribution to half field stimulation using different electrode montages. The visual evoked potential following full field stimulation was found to be asymmetrical around the midline with a bias over the left occiput particularly when the foveal polar projections of the occipital cortex were preferentially stimulated. The half field response reflected the distribution asymmetry. Masking of the central 3° resulted in a response which was approximately symmetrical around the midline but there was no evidence of the PNP-complex. A method for visual field quantification was developed based on the neural representation of visual space (Drasdo and Peaston 1982) in an attempt to relate visual field depravation with the resultant visual evoked potentials. There was no form of simple, diffuse summation between the scalp potential and the cortical generators. It was, however, possible to quantify the degree of scalp potential attenuation for M-scaled full field stimuli. The results obtained from patients exhibiting pre-chiasmal lesions suggested that the PNP-complex is not scotomatous in nature but confirmed that it is most likely to be related to specific diseases (Harding and Crews 1982). There was a strong correlation between the percentage information loss of the visual field and the diagnostic value of the visual evoked potential in patients exhibiting chiasmal lesions.

Relaxation procedures for an iterative MFS algorithm for the stable reconstruction of elastic fields from Cauchy data in two-dimensional isotropic linear elasticity

We investigate two numerical procedures for the Cauchy problem in linear elasticity, involving the relaxation of either the given boundary displacements (Dirichlet data) or the prescribed boundary tractions (Neumann data) on the over-specified boundary, in the alternating iterative algorithm of Kozlov et al. (1991). The two mixed direct (well-posed) problems associated with each iteration are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method, while the optimal value of the regularization parameter is chosen via the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases the iterative procedure at the point where the accumulation of noise becomes dominant and the errors in predicting the exact solutions increase, is also presented. The MFS-based iterative algorithms with relaxation are tested for Cauchy problems for isotropic linear elastic materials in various geometries to confirm the numerical convergence, stability, accuracy and computational efficiency of the proposed method.

Polarisation dynamics of vector soliton molecules in mode locked fibre laser

Two fundamental laser physics phenomena - dissipative soliton and polarisation of light are recently merged to the concept of vector dissipative soliton (VDS), viz. train of short pulses with specific state of polarisation (SOP) and shape defined by an interplay between anisotropy, gain/loss, dispersion, and nonlinearity. Emergence of VDSs is both of the fundamental scientific interest and is also a promising technique for control of dynamic SOPs important for numerous applications from nano-optics to high capacity fibre optic communications. Using specially designed and developed fast polarimeter, we present here the first experimental results on SOP evolution of vector soliton molecules with periodic polarisation switching between two and three SOPs and superposition of polarisation switching with SOP precessing. The underlying physics presents an interplay between linear and circular birefringence of a laser cavity along with light induced anisotropy caused by polarisation hole burning.

Fiber optic sensing of magnetic fields utilizing femtosecond laser sculpted microslots and long period gratings coated with Terfenol-D

Fiber optic sensors are fabricated for detecting static magnetic fields. The sensors consist of a UV inscribed long period grating with two 50 micron long microslots. The microslots are fabricated using the femtosecond laser based inscribe and etch technique. The microslots and the fiber surface are coated with a magnetostrictive material Terfenol-D. A spectral sensitivity of 1.15 pm/mT was measured in transmission with a working resolution of ±0.2 mT for a static magnetic field strength below 10 mT. These devices also present a different response when the spatial orientation of the fiber was adjusted relative to the magnetic field lines.

Dynamics of vector rogue waves in a fiber laser with a ring cavity

A pulse–pulse interaction that leads to rogue wave (RW) generation in lasers was previously attributed either to soliton–soliton or soliton–dispersive-wave interaction. The beating between polarization modes in the absence of a saturable absorber causes similar effects. Accounting for these polarization modes in a laser resonator is the purpose of the distributed vector model of laser resonators. Furthermore, high pump power, high amplitude, and short pulse duration are not necessary conditions to observe pulse attraction, repulsion, and collisions and the resonance exchange of energy between among them. The regimes of interest can be tuned just by changing the birefringence in the cavity with the pump power slightly higher than the laser threshold. This allows the observation of a wide range of RW patterns in the same experiment, as well as to classify them. The dynamics of the interaction between pulses leads us to the conclusion that all of these effects occur due to nonlinearity induced by the inverse population in the active fiber as well as an intrinsic nonlinearity in the passive part of the cavity. Most of the mechanisms of pulse–pulse interaction were found to be mutually exclusive. This means that all the observed RW patterns, namely, the “lonely,” “twins,” “three sisters,” and “cross,” are probably different cases of the same process.

Bayesian retrieval of scatterometer wind fields

The retrieval of wind fields from scatterometer observations has traditionally been separated into two phases; local wind vector retrieval and ambiguity removal. Operationally, a forward model relating wind vector to backscatter is inverted, typically using look up tables, to retrieve up to four local wind vector solutions. A heuristic procedure, using numerical weather prediction forecast wind vectors and, often, some neighbourhood comparison is then used to select the correct solution. In this paper we develop a Bayesian method for wind field retrieval, and show how a direct local inverse model, relating backscatter to wind vector, improves the wind vector retrieval accuracy. We compare these results with the operational U.K. Meteorological Office retrievals, our own CMOD4 retrievals and a neural network based local forward model retrieval. We suggest that the neural network based inverse model, which is extremely fast to use, improves upon current forward models when used in a variational data assimilation scheme.

Splines and vector splines

The thrust of this report concerns spline theory and some of the background to spline theory and follows the development in (Wahba, 1991). We also review methods for determining hyper-parameters, such as the smoothing parameter, by Generalised Cross Validation. Splines have an advantage over Gaussian Process based procedures in that we can readily impose atmospherically sensible smoothness constraints and maintain computational efficiency. Vector splines enable us to penalise gradients of vorticity and divergence in wind fields. Two similar techniques are summarised and improvements based on robust error functions and restricted numbers of basis functions given. A final, brief discussion of the application of vector splines to the problem of scatterometer data assimilation highlights the problems of ambiguous solutions.

A three dimensional numerical investigation of hillslope flow processes

Physically based distributed models of catchment hydrology are likely to be made available as engineering tools in the near future. Although these models are based on theoretically acceptable equations of continuity, there are still limitations in the present modelling strategy. Of interest to this thesis are the current modelling assumptions made concerning the effects of soil spatial variability, including formations producing distinct zones of preferential flow. The thesis contains a review of current physically based modelling strategies and a field based assessment of soil spatial variability. In order to investigate the effects of soil nonuniformity a fully three dimensional model of variability saturated flow in porous media is developed. The model is based on a Galerkin finite element approximation to Richards equation. Accessibility to a vector processor permits numerical solutions on grids containing several thousand node points. The model is applied to a single hillslope segment under various degrees of soil spatial variability. Such variability is introduced by generating random fields of saturated hydraulic conductivity using the turning bands method. Similar experiments are performed under conditions of preferred soil moisture movement. The results show that the influence of soil variability on subsurface flow may be less significant than suggested in the literature, due to the integrating effects of three dimensional flow. Under conditions of widespread infiltration excess runoff, the results indicate a greater significance of soil nonuniformity. The recognition of zones of preferential flow is also shown to be an important factor in accurate rainfall-runoff modelling. Using the results of various fields of soil variability, experiments are carried out to assess the validity of the commonly used concept of `effective parameters'. The results of these experiments suggest that such a concept may be valid in modelling subsurface flow. However, the effective parameter is observed to be event dependent when the dominating mechanism is infiltration excess runoff.

Characterisation and performance of a Terfenol-D coated femtosecond laser inscribed optical fibre Bragg sensor with a laser ablated microslot for the detection of static magnetic fields

We present a novel device for the characterisation of static magnetic fields through monitoring wavelength shifts of femtosecond inscribed fibre Bragg grating and micromachined slot, coated with Terfenol-D. The device was sensitive to static magnetic fields and can be used as a vectoral sensor for the detection of magnetic fields as low as 0.046 mT with a resolution of ± 0.3mT in transmission and ± 0.7mT in reflection. The use of a femtosecond laser to both inscribe the FBGs and micromachine the slot in a single stage prior to coating the device significantly simplifies the fabrication.

Vector solitons with locked and precessing states of polarization

We demonstrate experimentally new families of vector solitons with locked and precessing states of polarization for fundamental and multipulse soliton operations in a carbon nanotube mode-locked fiber laser with anomalous dispersion laser cavity.

Femtosecond laser inscribed Bragg sensor in Terfenol-D coated optical fibre with ablated microslot for the detection of static magnetic fields

A novel device for the detection and characterisation of static magnetic fields is presented. It consists of a femtosecond laser inscribed fibre Bragg grating (FBG) that is incorporated into an optical fibre with a femtosecond laser micromachined slot. The symmetry of the fibre is broken by the micro-slot, producing non-uniform strain across the fibre cross section. The sensing region is coated with Terfenol-D making the device sensitive to static magnetic fields, whereas the symmetry breaking results in a vectorial sensor for the detection of magnetic fields as low as 0.046 mT with a resolution of ±0.3mT in transmission and ±0.7mT in reflection. The sensor output is directly wavelength encoded from the FBG filtering, leading to simple demodulation through the monitoring of wavelength shifts that result as the fibre structure changes shape in response to the external magnetic field. The use of a femtosecond laser to both inscribe the FBG and micro-machine the slot in a single stage, prior to coating the device, significantly simplifies the sensor fabrication.