Theoretical insight into diagnosing observation error correlations using observation-minus-background and observation-minus-analysis statistics

To improve the quantity and impact of observations used in data assimilation it is necessary to take into account the full, potentially correlated, observation error statistics. A number of methods for estimating correlated observation errors exist, but a popular method is a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. The accuracy of the results it yields is unknown as the diagnostic is sensitive to the difference between the exact background and exact observation error covariances and those that are chosen for use within the assimilation. It has often been stated in the literature that the results using this diagnostic are only valid when the background and observation error correlation length scales are well separated. Here we develop new theory relating to the diagnostic. For observations on a 1D periodic domain we are able to the show the effect of changes in the assumed error statistics used in the assimilation on the estimated observation error covariance matrix. We also provide bounds for the estimated observation error variance and eigenvalues of the estimated observation error correlation matrix. We demonstrate that it is still possible to obtain useful results from the diagnostic when the background and observation error length scales are similar. In general, our results suggest that when correlated observation errors are treated as uncorrelated in the assimilation, the diagnostic will underestimate the correlation length scale. We support our theoretical results with simple illustrative examples. These results have potential use for interpreting the derived covariances estimated using an operational system.

Structure and electrical transport properties of the ordered skutterudites MGe1.5S1.5 (M = Co, Rh, Ir)

High-resolution powder neutron diffraction data collected for the skutterudites MGe1.5S1.5 (M=Co, Rh, Ir) reveal that these materials adopt an ordered skutterudite structure (space group R3¯), in which the anions are ordered in layers perpendicular to the [111] direction. In this ordered structure, the anions form two-crystallographically distinct four-membered rings, with stoichiometry Ge2S2, in which the Ge and S atoms are trans to each other. The transport properties of these materials, which are p-type semiconductors, are discussed in the light of the structural results. The effect of iron substitution in CoGe1.5S1.5 has been investigated. While doping of CoGe1.5S1.5 has a marked effect on both the electrical resistivity and the Seebeck coefficient, these ternary skutterudites exhibit significantly higher electrical resistivities than their binary counterparts.

Designing a sampling scheme to reveal correlations between weeds and soil properties at multiple spatial scales

Weeds tend to aggregate in patches within fields and there is evidence that this is partly owing to variation in soil properties. Because the processes driving soil heterogeneity operate at different scales, the strength of the relationships between soil properties and weed density would also be expected to be scale-dependent. Quantifying these effects of scale on weed patch dynamics is essential to guide the design of discrete sampling protocols for mapping weed distribution. We have developed a general method that uses novel within-field nested sampling and residual maximum likelihood (REML) estimation to explore scale-dependent relationships between weeds and soil properties. We have validated the method using a case study of Alopecurus myosuroides in winter wheat. Using REML, we partitioned the variance and covariance into scale-specific components and estimated the correlations between the weed counts and soil properties at each scale. We used variograms to quantify the spatial structure in the data and to map variables by kriging. Our methodology successfully captured the effect of scale on a number of edaphic drivers of weed patchiness. The overall Pearson correlations between A. myosuroides and soil organic matter and clay content were weak and masked the stronger correlations at >50 m. Knowing how the variance was partitioned across the spatial scales we optimized the sampling design to focus sampling effort at those scales that contributed most to the total variance. The methods have the potential to guide patch spraying of weeds by identifying areas of the field that are vulnerable to weed establishment.

Predicting sizes of hexagonal and gyroid metal nanostructures from liquid crystal templating

We describe a method to predict and control the lattice parameters of hexagonal and gyroid mesoporous materials formed by liquid crystal templating. In the first part, we describe a geometric model with which the lattice parameters of different liquid crystal mesophases can be predicted as a function of their water/surfactant/oil volume fractions, based on certain geometric parameters relating to the constituent surfactant molecules. We demonstrate the application of this model to the lamellar (LR), hexagonal (H1), and gyroid bicontinuous cubic (V1) mesophases formed by the binary Brij-56 (C16EO10)/water system and the ternary Brij-56/hexadecane/water system. In this way, we demonstrate predictable and independent control over the size of the cylinders (with hexadecane) and their spacing (with water). In the second part, we produce mesoporous platinum using as templates hexagonal and gyroid phases with different compositions and show that in each case the symmetry and lattice parameter of the metal nanostructure faithfully replicate those of the liquid crystal template, which is itself in agreement with the model. This demonstrates a rational control over the geometry, size, and spacing of pores in a mesoporous metal.

Diagnosing observation error correlations for Doppler radar radial winds in the Met Office UKV model using observation-minus-background and observation-minus-analysis statistics

With the development of convection-permitting numerical weather prediction the efficient use of high resolution observations in data assimilation is becoming increasingly important. The operational assimilation of these observations, such as Dopplerradar radial winds, is now common, though to avoid violating the assumption of un- correlated observation errors the observation density is severely reduced. To improve the quantity of observations used and the impact that they have on the forecast will require the introduction of the full, potentially correlated, error statistics. In this work, observation error statistics are calculated for the Doppler radar radial winds that are assimilated into the Met Office high resolution UK model using a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. This is the first in-depth study using the diagnostic to estimate both horizontal and along-beam correlated observation errors. By considering the new results obtained it is found that the Doppler radar radial wind error standard deviations are similar to those used operationally and increase as the observation height increases. Surprisingly the estimated observation error correlation length scales are longer than the operational thinning distance. They are dependent on both the height of the observation and on the distance of the observation away from the radar. Further tests show that the long correlations cannot be attributed to the use of superobservations or the background error covariance matrix used in the assimilation. The large horizontal correlation length scales are, however, in part, a result of using a simplified observation operator.