Estimating interchannel observation-error correlations for IASI radiance data in the Met Office system

The optimal utilisation of hyper-spectral satellite observations in numerical weather prediction is often inhibited by incorrectly assuming independent interchannel observation errors. However, in order to represent these observation-error covariance structures, an accurate knowledge of the true variances and correlations is needed. This structure is likely to vary with observation type and assimilation system. The work in this article presents the initial results for the estimation of IASI interchannel observation-error correlations when the data are processed in the Met Office one-dimensional (1D-Var) and four-dimensional (4D-Var) variational assimilation systems. The method used to calculate the observation errors is a post-analysis diagnostic which utilises the background and analysis departures from the two systems. The results show significant differences in the source and structure of the observation errors when processed in the two different assimilation systems, but also highlight some common features. When the observations are processed in 1D-Var, the diagnosed error variances are approximately half the size of the error variances used in the current operational system and are very close in size to the instrument noise, suggesting that this is the main source of error. The errors contain no consistent correlations, with the exception of a handful of spectrally close channels. When the observations are processed in 4D-Var, we again find that the observation errors are being overestimated operationally, but the overestimation is significantly larger for many channels. In contrast to 1D-Var, the diagnosed error variances are often larger than the instrument noise in 4D-Var. It is postulated that horizontal errors of representation, not seen in 1D-Var, are a significant contributor to the overall error here. Finally, observation errors diagnosed from 4D-Var are found to contain strong, consistent correlation structures for channels sensitive to water vapour and surface properties.

Estimating correlated observation error statistics using an ensemble transform Kalman filter

For certain observing types, such as those that are remotely sensed, the observation errors are correlated and these correlations are state- and time-dependent. In this work, we develop a method for diagnosing and incorporating spatially correlated and time-dependent observation error in an ensemble data assimilation system. The method combines an ensemble transform Kalman filter with a method that uses statistical averages of background and analysis innovations to provide an estimate of the observation error covariance matrix. To evaluate the performance of the method, we perform identical twin experiments using the Lorenz ’96 and Kuramoto-Sivashinsky models. Using our approach, a good approximation to the true observation error covariance can be recovered in cases where the initial estimate of the error covariance is incorrect. Spatial observation error covariances where the length scale of the true covariance changes slowly in time can also be captured. We find that using the estimated correlated observation error in the assimilation improves the analysis.

Effects of action observation on corticospinal excitability: muscle specificity, direction, and timing of the mirror response

Many human behaviours and pathologies have been attributed to the putative mirror neuron system, a neural system that is active during both the observation and execution of actions. While there are now a very large number of papers on the mirror neuron system, variations in the methods and analyses employed by researchers mean that the basic characteristics of the mirror response are not clear. This review focuses on three important aspects of the mirror response, as measured by modulations in corticospinal excitability: (1) muscle specificity, (2) direction, and (3) timing of modulation. We focus mainly on electromyographic (EMG) data gathered following single-pulse transcranial magnetic stimulation (TMS), because this method provides precise information regarding these three aspects of the response. Data from paired-pulse TMS paradigms and peripheral nerve stimulation (PNS) are also considered when we discuss the possible mechanisms underlying the mirror response. In this systematic review of the literature, we examine the findings of 85 TMS and PNS studies of the human mirror response, and consider the limitations and advantages of the different methodological approaches these have adopted in relation to discrepancies between their findings. We conclude by proposing a testable model of how action observation modulates corticospinal excitability in humans. Specifically, we propose that action observation elicits an early, non-specific facilitation of corticospinal excitability (at around 90 ms from action onset), followed by a later modulation of activity specific to the muscles involved in the observed action (from around 200 ms). Testing this model will greatly advance our understanding of the mirror mechanism and provide a more stable grounding on which to base inferences about its role in human behaviour.

A novel ensemble method for retrieving properties of warm cloud in 3-D using ground-based scanning radar and zenith radiances

We present a novel method for retrieving high-resolution, three-dimensional (3-D) nonprecipitating cloud fields in both overcast and broken-cloud situations. The method uses scanning cloud radar and multiwavelength zenith radiances to obtain gridded 3-D liquid water content (LWC) and effective radius (re) and 2-D column mean droplet number concentration (Nd). By using an adaption of the ensemble Kalman filter, radiances are used to constrain the optical properties of the clouds using a forward model that employs full 3-D radiative transfer while also providing full error statistics given the uncertainty in the observations. To evaluate the new method, we first perform retrievals using synthetic measurements from a challenging cumulus cloud field produced by a large-eddy simulation snapshot. Uncertainty due to measurement error in overhead clouds is estimated at 20% in LWC and 6% in re, but the true error can be greater due to uncertainties in the assumed droplet size distribution and radiative transfer. Over the entire domain, LWC and re are retrieved with average error 0.05–0.08 g m-3 and ~2 μm, respectively, depending on the number of radiance channels used. The method is then evaluated using real data from the Atmospheric Radiation Measurement program Mobile Facility at the Azores. Two case studies are considered, one stratocumulus and one cumulus. Where available, the liquid water path retrieved directly above the observation site was found to be in good agreement with independent values obtained from microwave radiometer measurements, with an error of 20 g m-2.

Identification of domestication-related loci associated with flowering time and seed size in soybean with the RAD-seq genotyping method

Flowering time and seed size are traits related to domestication. However, identification of domestication-related loci/genes of controlling the traits in soybean is rarely reported. In this study, we identified a total of 48 domestication-related loci based on RAD-seq genotyping of a natural population comprising 286 accessions. Among these, four on chromosome 12 and additional two on chromosomes 11 and 15 were associated with flowering time, and four on chromosomes 11 and 16 were associated with seed size. Of the five genes associated with flowering time and the three genes associated with seed size, three genes Glyma11g18720, Glyma11g15480 and Glyma15g35080 were homologous to Arabidopsis genes, additional five genes were found for the first time to be associated with these two traits. Glyma11g18720 and Glyma05g28130 were co-expressed with five genes homologous to flowering time genes in Arabidopsis, and Glyma11g15480 was co-expressed with 24 genes homologous to seed development genes in Arabidopsis. This study indicates that integration of population divergence analysis, genome-wide association study and expression analysis is an efficient approach to identify candidate domestication-related genes.

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.

Selecting algorithms for Earth observation of climate within the European Space Agency Climate Change initiative: introduction to a special issue

This special issue is focused on the assessment of algorithms for the observation of Earth’s climate from environ- mental satellites. Climate data records derived by remote sensing are increasingly a key source of insight into the workings of and changes in Earth’s climate system. Producers of data sets must devote considerable effort and expertise to maximise the true climate signals in their products and minimise effects of data processing choices and changing sensors. A key choice is the selection of algorithm(s) for classification and/or retrieval of the climate variable. Within the European Space Agency Climate Change Initiative, science teams undertook systematic assessment of algorithms for a range of essential climate variables. The papers in the special issue report some of these exercises (for ocean colour, aerosol, ozone, greenhouse gases, clouds, soil moisture, sea surface temper- ature and glaciers). The contributions show that assessment exercises must be designed with care, considering issues such as the relative importance of different aspects of data quality (accuracy, precision, stability, sensitivity, coverage, etc.), the availability and degree of independence of validation data and the limitations of validation in characterising some important aspects of data (such as long-term stability or spatial coherence). As well as re- quiring a significant investment of expertise and effort, systematic comparisons are found to be highly valuable. They reveal the relative strengths and weaknesses of different algorithmic approaches under different observa- tional contexts, and help ensure that scientific conclusions drawn from climate data records are not influenced by observational artifacts, but are robust.

A systematic method of parameterisation estimation using data assimilation

In numerical weather prediction, parameterisations are used to simulate missing physics in the model. These can be due to a lack of scientific understanding or a lack of computing power available to address all the known physical processes. Parameterisations are sources of large uncertainty in a model as parameter values used in these parameterisations cannot be measured directly and hence are often not well known; and the parameterisations themselves are also approximations of the processes present in the true atmosphere. Whilst there are many efficient and effective methods for combined state/parameter estimation in data assimilation (DA), such as state augmentation, these are not effective at estimating the structure of parameterisations. A new method of parameterisation estimation is proposed that uses sequential DA methods to estimate errors in the numerical models at each space-time point for each model equation. These errors are then fitted to pre-determined functional forms of missing physics or parameterisations that are based upon prior information. We applied the method to a one-dimensional advection model with additive model error, and it is shown that the method can accurately estimate parameterisations, with consistent error estimates. Furthermore, it is shown how the method depends on the quality of the DA results. The results indicate that this new method is a powerful tool in systematic model improvement.

Reliability of a method for evaluating porosity in denture base resins

Background: The method of porosity analysis by water absorption has been carried out by the storage of the specimens in pure water, but it does not exclude the potential plasticising effect of the water generating unreal values of porosity. Objective: The present study evaluated the reliability of this method of porosity analysis in polymethylmethacrylate denture base resins by the determination of the most satisfactory solution for storage (S), where the plasticising effect was excluded. Materials and methods: Two specimen shapes (rectangular and maxillary denture base) and two denture base resins, water bath-polymerised (Classico) and microwave-polymerised (Acron MC) were used. Saturated anhydrous calcium chloride solutions (25%, 50%, 75%) and distilled water were used for specimen storage. Sorption isotherms were used to determine S. Porosity factor (PF) and diffusion coefficient (D) were calculated within S and for the groups stored in distilled water. anova and Tukey tests were performed to identify significant differences in PF results and Kruskal-Wallis test and Dunn multiple comparison post hoc test, for D results (alpha = 0.05). Results: For Acron MC denture base shape, FP results were 0.24% (S 50%) and 1.37% (distilled water); for rectangular shape FP was 0.35% (S 75%) and 0.19% (distilled water). For Classico denture base shape, FP results were 0.54% (S 75%) and 1.21% (distilled water); for rectangular shape FP was 0.7% (S 50%) and 1.32% (distilled water). FP results were similar in S and distilled water only for Acron MC rectangular shape (p > 0.05). D results in distilled water were statistically higher than S for all groups. Conclusions: The results of the study suggest that an adequate solution for storing specimens must be used to measure porosity by water absorption, based on excluding the plasticising effect.

Paleointensity data from Early Cretaceous Ponta Grossa dikes (Brazil) using a multisample method

Definition of the long-term variation of the geomagnetic virtual dipole moment requires more reliable paleointensity results. Here, we applied a multisample protocol to the study of the 130.5 Ma Ponta Grossa basaltic dikes (southern Brazil) that carry a very stable dual-polarity magnetic component. The magnetic stability of the samples wits checked using thermomagnetic curves and by monitoring the magnetic Susceptibility evolution through the paleointensity experiments. Twelve sites containing the least alterable samples were chosen for the paleointensity measurements. Although these rocks failed stepwise double-heating experiments, they yielded coherent results in the multisample method for all sites but one. The coherent sites show low to moderate field intensities between 5.7 +/- 0.2 and 26.4 +/- 0.7 mu T (average 13.4 +/- 1.9 mu T). Virtual dipole moments for these sites range from 1.3 +/- 0.04 to 6.0 +/- 0.2 x 10(22) A m(2) (average 2.9 +/- 0.5 x 10(22) A m(2)). Our results agree with the tendency for low dipole moments during the Early Cretaceous, immediately prior to the Cretaceous Normal Superchron (CNS). The available paleointensity database shows a strong variability of the field between 80 and 160 Ma. There seems to be no firm evidence for a Mesozoic Dipole Low, but a long-term tendency does emerge from the data with the highest dipole moments Occurring at the middle of the CNS.

First observation of the Cabibbo-suppressed decays Xi(+)(c)->Sigma(+)pi(-)pi(+) and Xi(+)(c)->Sigma(-)pi(+)pi(+) and measurement of their branching ratios

We report the first observation of two Cabibbo-suppressed ecay modes, Xi(+)(c) -> Sigma(+)pi(-)pi(+) and Xi c+ -> Sigma(-)pi(+)pi(+). We observe 59 +/- 14 over a background of 87, and 22 +/- 8 over a background of 13 events, respectively, for the signals. The data were accumulated using the SELEX spectrometer during the 1996-1997 fixed target run at Fermilab, chiefly from a 600 Gev/c Sigma(-) beam. The branching ratios of the decays relative to the Cabibbo-favored Xi c+ -> Xi(-)pi(+)pi(+) are measured to be B(Xi(+)(c) -> Sigma(+)pi(-)pi(+))/B(Xi(+)(c) -> Xi(-)pi(+)pi(+)) = 0.48 +/- 0.20, and B(Xi(+)(c) -> Sigma(-)pi(+)pi(+))/B(Xi(+)(c) -> Sigma(-)pi(+)pi(+)) = 0.18 +/- 0.09, respectively. We also report branching ratios for the same decay modes of the Delta(+)(c) relative to Delta(+)(c) -> pK(-)pi(+.) (C) 2008 Elsevier B.V. All rights reserved.

Observation of an Antimatter Hypernucleus

Nuclear collisions recreate conditions in the universe microseconds after the Big Bang. Only a very small fraction of the emitted fragments are light nuclei, but these states are of fundamental interest. We report the observation of antihypertritons-comprising an antiproton, an antineutron, and an antilambda hyperon-produced by colliding gold nuclei at high energy. Our analysis yields 70 +/- 17 antihypertritons (3/Lambda(H) over bar) and 157 +/- 30 hypertritons ((3)(Lambda)H). The measured yields of (3)(Lambda)H (3/Lambda(H) over bar) and (3)He ((3)(He) over bar) are similar, suggesting an equilibrium in coordinate and momentum space populations of up, down, and strange quarks and antiquarks, unlike the pattern observed at lower collision energies. The production and properties of antinuclei, and of nuclei containing strange quarks, have implications spanning nuclear and particle physics, astrophysics, and cosmology.

Observation of the antimatter helium-4 nucleus

High-energy nuclear collisions create an energy density similar to that of the Universe microseconds after the Big Bang(1); in both cases, matter and antimatter are formed with comparable abundance. However, the relatively short-lived expansion in nuclear collisions allows antimatter to decouple quickly from matter, and avoid annihilation. Thus, a high-energy accelerator of heavy nuclei provides an efficient means of producing and studying antimatter. The antimatter helium-4 nucleus ((4)(He) over bar), also known as the anti-alpha ((alpha) over bar), consists of two antiprotons and two antineutrons (baryon number B = -4). It has not been observed previously, although the alpha-particle was identified a century ago by Rutherford and is present in cosmic radiation at the ten per cent level(2). Antimatter nuclei with B -1 have been observed only as rare products of interactions at particle accelerators, where the rate of antinucleus production in high-energy collisions decreases by a factor of about 1,000 with each additional antinucleon(3-5). Here we report the observation of (4)<(He) over bar, the heaviest observed antinucleus to date. In total, 18 (4)(He) over bar counts were detected at the STAR experiment at the Relativistic Heavy Ion Collider (RHIC; ref. 6) in 10(9) recorded gold-on-gold (Au+Au) collisions at centre-of-mass energies of 200 GeV and 62 GeV per nucleon-nucleon pair. The yield is consistent with expectations from thermodynamic(7) and coalescent nucleosynthesis(8) models, providing an indication of the production rate of even heavier antimatter nuclei and a benchmark for possible future observations of (4)(He) over bar in cosmic radiation.

Mice and larvae tracking using a particle filter with an auto-adjustable observation model

This paper proposes a novel way to combine different observation models in a particle filter framework. This, so called, auto-adjustable observation model, enhance the particle filter accuracy when the tracked objects overlap without infringing a great runtime penalty to the whole tracking system. The approach has been tested under two important real world situations related to animal behavior: mice and larvae tracking. The proposal was compared to some state-of-art approaches and the results show, under the datasets tested, that a good trade-off between accuracy and runtime can be achieved using an auto-adjustable observation model. (C) 2009 Elsevier B.V. All rights reserved.

Vad sker i elevers parskrivande vid datorn? : Observation av fyra elevpar i årskurs ett

The aim with this study is to explore if and how linguistic learning takes places when four pairs ofstudents in grade one writes on the computer. Observations of eight students took place duringtwelve lessons to investigate what happened during their writing processes. This study isimportant because students’ computer writing is a relatively new phenomenon in the Swedishprimary school.It is clear that the peer collaboration creates a complex interaction that becomes the platform, andalso the prerequisite, for successful language training. The method can clearly give students profitsin their linguistic development. When working in pairs, children support each other in writing andreading. They also get to train letter knowledge, phonemic awareness and the correspondencebetween phonemes and graphemes. Although the method creates good possibilities for readingand writing development, it cannot be the only way to learn reading. There are also risks thatstudents get stuck in different positions that can affect their progresses, which is important thatthe teachers consider while creating writing pairs. Teachers need to identify student’s reading andwriting development processes and also be able to create reading and writing instruction thatprovides different learning strategies.