A data assimilation method for using low-resolution Earth observation data in heterogeneous ecosystems

We present an approach for dealing with coarse-resolution Earth observations (EO) in terrestrial ecosystem data assimilation schemes. The use of coarse-scale observations in ecological data assimilation schemes is complicated by spatial heterogeneity and nonlinear processes in natural ecosystems. If these complications are not appropriately dealt with, then the data assimilation will produce biased results. The “disaggregation” approach that we describe in this paper combines frequent coarse-resolution observations with temporally sparse fine-resolution measurements. We demonstrate the approach using a demonstration data set based on measurements of an Arctic ecosystem. In this example, normalized difference vegetation index observations are assimilated into a “zero-order” model of leaf area index and carbon uptake. The disaggregation approach conserves key ecosystem characteristics regardless of the observation resolution and estimates the carbon uptake to within 1% of the demonstration data set “truth.” Assimilating the same data in the normal manner, but without the disaggregation approach, results in carbon uptake being underestimated by 58% at an observation resolution of 250 m. The disaggregation method allows the combination of multiresolution EO and improves in spatial resolution if observations are located on a grid that shifts from one observation time to the next. Additionally, the approach is not tied to a particular data assimilation scheme, model, or EO product and can cope with complex observation distributions, as it makes no implicit assumptions of normality.

European plant phenology and climate as seen in a 20-year AVHRR land-surface parameter dataset

Vegetation distribution and state have been measured since 1981 by the AVHRR (Advanced Very High Resolution Radiometer) instrument through satellite remote sensing. In this study a correction method is applied to the Pathfinder NDVI (Normalized Difference Vegetation Index) data to create a continuous European vegetation phenology dataset of a 10-day temporal and 0.1° spatial resolution; additionally, land surface parameters for use in biosphere–atmosphere modelling are derived. The analysis of time-series from this dataset reveals, for the years 1982–2001, strong seasonal and interannual variability in European land surface vegetation state. Phenological metrics indicate a late and short growing season for the years 1985–1987, in addition to early and prolonged activity in the years 1989, 1990, 1994 and 1995. These variations are in close agreement with findings from phenological measurements at the surface; spring phenology is also shown to correlate particularly well with anomalies in winter temperature and winter North Atlantic Oscillation (NAO) index. Nevertheless, phenological metrics, which display considerable regional differences, could only be determined for vegetation with a seasonal behaviour. Trends in the phenological phases reveal a general shift to earlier (−0.54 days year−1) and prolonged (0.96 days year−1) growing periods which are statistically significant, especially for central Europe.

Some numerical experiments on the effect of the variation of static stability in two-layer quasi-geostrophic models

A method to solve a quasi-geostrophic two-layer model including the variation of static stability is presented. The divergent part of the wind is incorporated by means of an iterative procedure. The procedure is rather fast and the time of computation is only 60–70% longer than for the usual two-layer model. The method of solution is justified by the conservation of the difference between the gross static stability and the kinetic energy. To eliminate the side-boundary conditions the experiments have been performed on a zonal channel model. The investigation falls mainly into three parts: The first part (section 5) contains a discussion of the significance of some physically inconsistent approximations. It is shown that physical inconsistencies are rather serious and for these inconsistent models which were studied the total kinetic energy increased faster than the gross static stability. In the next part (section 6) we are studying the effect of a Jacobian difference operator which conserves the total kinetic energy. The use of this operator in two-layer models will give a slight improvement but probably does not have any practical use in short periodic forecasts. It is also shown that the energy-conservative operator will change the wave-speed in an erroneous way if the wave-number or the grid-length is large in the meridional direction. In the final part (section 7) we investigate the behaviour of baroclinic waves for some different initial states and for two energy-consistent models, one with constant and one with variable static stability. According to the linear theory the waves adjust rather rapidly in such a way that the temperature wave will lag behind the pressure wave independent of the initial configuration. Thus, both models give rise to a baroclinic development even if the initial state is quasi-barotropic. The effect of the variation of static stability is very small, qualitative differences in the development are only observed during the first 12 hours. For an amplifying wave we will get a stabilization over the troughs and an instabilization over the ridges.

Impact of serving method on the consumption of nutritional supplement drinks: randomised trial in older adults with cognitive impairment

Aim: To analyse the influence of serving method on compliance and consumption of nutritional supplement drinks in older adults with cognitive impairment. Background: Oral nutritional supplement drinks have positive benefits on increasing nutritional status within undernourished elderly people leading to weight gain. However, consumption of these drinks is low and therefore limits their effectiveness. Design: This study was a non blind randomised control trial where participants either consumed nutritional supplement drinks in a glass/beaker or consumed them through a straw inserted directly into the container. Method: Participants with longstanding cognitive impairment were recruited from nursing homes (n=31) and hospitals (n=14). Participants were randomised to serving method. Nursing and care staff were instructed to give the supplement drinks three times per day on alternate days over a week by the allocated serving method. The researcher weighed the amount of supplement drink remaining after consumption. Data were collected over 12 months in 2011-2012. Results: 45 people participated in this study mean age 86.7 (SD 7.5 ) years. After randomisation there was no significant difference between the baseline characteristics of the two groups. Participants randomised to consume nutritional drinks from a glass / beaker drank significantly more than those who consumed them via a straw inserted directly into the container. However, supplements allocated to be given in a glass/beaker were more frequently omitted. Conclusion: Nutritional supplement drinks should be given to people with dementia who are able to feed themselves in a glass or a beaker if staffing resources allow (NIHR CSP ref 31101).

A method for under-sampled ecological network data analysis: plant-pollination as case study

In this paper, we develop a method, termed the Interaction Distribution (ID) method, for analysis of quantitative ecological network data. In many cases, quantitative network data sets are under-sampled, i.e. many interactions are poorly sampled or remain unobserved. Hence, the output of statistical analyses may fail to differentiate between patterns that are statistical artefacts and those which are real characteristics of ecological networks. The ID method can support assessment and inference of under-sampled ecological network data. In the current paper, we illustrate and discuss the ID method based on the properties of plant-animal pollination data sets of flower visitation frequencies. However, the ID method may be applied to other types of ecological networks. The method can supplement existing network analyses based on two definitions of the underlying probabilities for each combination of pollinator and plant species: (1), pi,j: the probability for a visit made by the i’th pollinator species to take place on the j’th plant species; (2), qi,j: the probability for a visit received by the j’th plant species to be made by the i’th pollinator. The method applies the Dirichlet distribution to estimate these two probabilities, based on a given empirical data set. The estimated mean values for pi,j and qi,j reflect the relative differences between recorded numbers of visits for different pollinator and plant species, and the estimated uncertainty of pi,j and qi,j decreases with higher numbers of recorded visits.

Assessment of inter-examiner agreement and variability in the manual classification of auditory brainstem response

Abstract Background: The analysis of the Auditory Brainstem Response (ABR) is of fundamental importance to the investigation of the auditory system behaviour, though its interpretation has a subjective nature because of the manual process employed in its study and the clinical experience required for its analysis. When analysing the ABR, clinicians are often interested in the identification of ABR signal components referred to as Jewett waves. In particular, the detection and study of the time when these waves occur (i.e., the wave latency) is a practical tool for the diagnosis of disorders affecting the auditory system. Significant differences in inter-examiner results may lead to completely distinct clinical interpretations of the state of the auditory system. In this context, the aim of this research was to evaluate the inter-examiner agreement and variability in the manual classification of ABR. Methods: A total of 160 ABR data samples were collected, for four different stimulus intensity (80dBHL, 60dBHL, 40dBHL and 20dBHL), from 10 normal-hearing subjects (5 men and 5 women, from 20 to 52 years). Four examiners with expertise in the manual classification of ABR components participated in the study. The Bland-Altman statistical method was employed for the assessment of inter-examiner agreement and variability. The mean, standard deviation and error for the bias, which is the difference between examiners’ annotations, were estimated for each pair of examiners. Scatter plots and histograms were employed for data visualization and analysis. Results: In most comparisons the differences between examiner’s annotations were below 0.1 ms, which is clinically acceptable. In four cases, it was found a large error and standard deviation (>0.1 ms) that indicate the presence of outliers and thus, discrepancies between examiners. Conclusions: Our results quantify the inter-examiner agreement and variability of the manual analysis of ABR data, and they also allows for the determination of different patterns of manual ABR analysis.

A fast two-grid and finite section method for a class of integral equations on the real line with application to an acoustic scattering problem in the half-plane

We consider the numerical treatment of second kind integral equations on the real line of the form ∅(s) = ∫_(-∞)^(+∞)▒〖κ(s-t)z(t)ϕ(t)dt,s=R〗 (abbreviated ϕ= ψ+K_z ϕ) in which K ϵ L_1 (R), z ϵ L_∞ (R) and ψ ϵ BC(R), the space of bounded continuous functions on R, are assumed known and ϕ ϵ BC(R) is to be determined. We first derive sharp error estimates for the finite section approximation (reducing the range of integration to [-A, A]) via bounds on (1-K_z )^(-1)as an operator on spaces of weighted continuous functions. Numerical solution by a simple discrete collocation method on a uniform grid on R is then analysed: in the case when z is compactly supported this leads to a coefficient matrix which allows a rapid matrix-vector multiply via the FFT. To utilise this possibility we propose a modified two-grid iteration, a feature of which is that the coarse grid matrix is approximated by a banded matrix, and analyse convergence and computational cost. In cases where z is not compactly supported a combined finite section and two-grid algorithm can be applied and we extend the analysis to this case. As an application we consider acoustic scattering in the half-plane with a Robin or impedance boundary condition which we formulate as a boundary integral equation of the class studied. Our final result is that if z (related to the boundary impedance in the application) takes values in an appropriate compact subset Q of the complex plane, then the difference between ϕ(s)and its finite section approximation computed numerically using the iterative scheme proposed is ≤C_1 [kh log⁡〖(1⁄kh)+(1-Θ)^((-1)⁄2) (kA)^((-1)⁄2) 〗 ] in the interval [-ΘA,ΘA](Θ<1) for kh sufficiently small, where k is the wavenumber and h the grid spacing. Moreover this numerical approximation can be computed in ≤C_2 N log⁡N operations, where N = 2A/h is the number of degrees of freedom. The values of the constants C1 and C2 depend only on the set Q and not on the wavenumber k or the support of z.

Restricting microbial exposure in early life negates the immune benefits associated with gut colonization in environments of high microbial diversity

Background: Acquisition of the intestinal microbiota in early life corresponds with the development of the mucosal immune system. Recent work on caesarean-delivered infants revealed that early microbial composition is influenced by birthing method and environment. Furthermore, we have confirmed that early-life environment strongly influences both the adult gut microbiota and development of the gut immune system. Here, we address the impact of limiting microbial exposure after initial colonization on the development of adult gut immunity. Methodology/Principal Findings: Piglets were born in indoor or outdoor rearing units, allowing natural colonization in the immediate period after birth, prior to transfer to high-health status isolators. Strikingly, gut closure and morphological development were strongly affected by isolator-rearing, independent of indoor or outdoor origins of piglets. Isolator-reared animals showed extensive vacuolation and disorganization of the gut epithelium, inferring that normal gut closure requires maturation factors present in maternal milk. Although morphological maturation and gut closure were delayed in isolatorreared animals, these hard-wired events occurred later in development. Type I IFN, IL-22, IL-23 and Th17 pathways were increased in indoor-isolator compared to outdoor-isolator animals during early life, indicating greater immune activation in pigs originating from indoor environments reflecting differences in the early microbiota. This difference was less apparent later in development due to enhanced immune activation and convergence of the microbiota in all isolator-reared animals. This correlated with elevation of Type I IFN pathways in both groups, although T cell pathways were still more affected in indoor-reared animals. Conclusions/Significance: Environmental factors, in particular microbial exposure, influence expression of a large number of immune-related genes. However, the homeostatic effects of microbial colonization in outdoor environments require sustained microbial exposure throughout development. Gut development in high-hygiene environments negatively impacts on normal succession of the gut microbiota and promotes innate immune activation which may impair immune homeostasis.

Effect of forage conservation method on plasma lipids, mammary lipogenesis, and milk fatty acid composition in lactating cows fed diets containing a 60:40 forage-to-concentrate ratio

The effects of forage conservation method on plasma lipids, mammary lipogenesis, and milk fat were examined in 2 complementary experiments. Treatments comprised fresh grass, hay, or untreated (UTS) or formic acid treated silage (FAS) prepared from the same grass sward. Preparation of conserved forages coincided with the collection of samples from cows fed fresh grass. In the first experiment, 5 multiparous Finnish Ayrshire cows (229 d in milk) were used to compare a diet based on fresh grass followed by hay during 2 consecutive 14-d periods, separated by a 5-d transition during which extensively wilted grass was fed. In the second experiment, 5 multiparous Finnish Ayrshire cows (53 d in milk) were assigned to 1 of 2 blocks and allocated treatments according to a replicated 3 × 3 Latin square design, with 14-d periods to compare hay, UTS, and FAS. Cows received 7 or 9 kg/d of the same concentrate in experiments 1 and 2, respectively. Arterial concentrations of triacylglycerol (TAG) and phospholipid were higher in cows fed fresh grass, UTS, and FAS compared with hay. Nonesterified fatty acid (NEFA) concentrations and the relative abundance of 18:2n-6 and 18:3n-3 in TAG of arterial blood were also higher in cows fed fresh grass than conserved forages. On all diets, TAG was the principle source of fatty acids (FA) for milk fat synthesis, whereas mammary extraction of NEFA was negligible, except during zero-grazing, which was associated with a lower, albeit positive calculated energy balance. Mammary FA uptake was higher and the synthesis of 16:0 lower in cows fed fresh grass than hay. Conservation of grass by drying or ensiling had no influence on mammary extraction of TAG and NEFA, despite an increase in milk fat secretion for silages compared with hay and for FAS than UTS. Relative to hay, milk fat from fresh grass contained lower 12:0, 14:0, and 16:0 and higher S3,R7,R11,15-tetramethyl-16:0, cis-9 18:1, trans-11 18:1, cis-9,trans-11 18:2, 18:2n-6, and 18:3n-3 concentrations. Even though conserved forages altered mammary lipogenesis, differences in milk FA composition were relatively minor, other than a higher enrichment of S3,R7,R11,15-tetramethyl-16:0 in milk from silages compared with hay. In conclusion, differences in milk fat composition on fresh grass relative to conserved forages were associated with a lower energy balance, increased uptake of preformed FA, and decreased synthesis of 16:0 de novo in the mammary glands, in the absence of alterations in stearoyl-coenzyme A desaturase activity.

Equator-to-pole temperature differences and the extra-tropical storm track responses of the CMIP5 climate models

This paper aims to understand the physical processes causing the large spread in the storm track projections of the CMIP5 climate models. In particular, the relationship between the climate change responses of the storm tracks, as measured by the 2–6 day mean sea level pressure variance, and the equator-to-pole temperature differences at upper- and lower-tropospheric levels is investigated. In the southern hemisphere the responses of the upper- and lower-tropospheric temperature differences are correlated across the models and as a result they share similar associations with the storm track responses. There are large regions in which the storm track responses are correlated with the temperature difference responses, and a simple linear regression model based on the temperature differences at either level captures the spatial pattern of the mean storm track response as well explaining between 30 and 60 % of the inter-model variance of the storm track responses. In the northern hemisphere the responses of the two temperature differences are not significantly correlated and their associations with the storm track responses are more complicated. In summer, the responses of the lower-tropospheric temperature differences dominate the inter-model spread of the storm track responses. In winter, the responses of the upper- and lower-temperature differences both play a role. The results suggest that there is potential to reduce the spread in storm track responses by constraining the relative magnitudes of the warming in the tropical and polar regions.

Modelling of urban sensible heat flux at multiple spatial scales : a demonstration using airborne hyperspectral imagery of Shanghai and a temperature–emissivity separation approach

Urbanization related alterations to the surface energy balance impact urban warming (‘heat islands’), the growth of the boundary layer, and many other biophysical processes. Traditionally, in situ heat flux measures have been used to quantify such processes, but these typically represent only a small local-scale area within the heterogeneous urban environment. For this reason, remote sensing approaches are very attractive for elucidating more spatially representative information. Here we use hyperspectral imagery from a new airborne sensor, the Operative Modular Imaging Spectrometer (OMIS), along with a survey map and meteorological data, to derive the land cover information and surface parameters required to map spatial variations in turbulent sensible heat flux (QH). The results from two spatially-explicit flux retrieval methods which use contrasting approaches and, to a large degree, different input data are compared for a central urban area of Shanghai, China: (1) the Local-scale Urban Meteorological Parameterization Scheme (LUMPS) and (2) an Aerodynamic Resistance Method (ARM). Sensible heat fluxes are determined at the full 6 m spatial resolution of the OMIS sensor, and at lower resolutions via pixel aggregation and spatial averaging. At the 6 m spatial resolution, the sensible heat flux of rooftop dominated pixels exceeds that of roads, water and vegetated areas, with values peaking at ∼ 350 W m− 2, whilst the storage heat flux is greatest for road dominated pixels (peaking at around 420 W m− 2). We investigate the use of both OMIS-derived land surface temperatures made using a Temperature–Emissivity Separation (TES) approach, and land surface temperatures estimated from air temperature measures. Sensible heat flux differences from the two approaches over the entire 2 × 2 km study area are less than 30 W m− 2, suggesting that methods employing either strategy maybe practica1 when operated using low spatial resolution (e.g. 1 km) data. Due to the differing methodologies, direct comparisons between results obtained with the LUMPS and ARM methods are most sensibly made at reduced spatial scales. At 30 m spatial resolution, both approaches produce similar results, with the smallest difference being less than 15 W m− 2 in mean QH averaged over the entire study area. This is encouraging given the differing architecture and data requirements of the LUMPS and ARM methods. Furthermore, in terms of mean study QH, the results obtained by averaging the original 6 m spatial resolution LUMPS-derived QH values to 30 and 90 m spatial resolution are within ∼ 5 W m− 2 of those derived from averaging the original surface parameter maps prior to input into LUMPS, suggesting that that use of much lower spatial resolution spaceborne imagery data, for example from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is likely to be a practical solution for heat flux determination in urban areas.

Global snow mass measurements and the effect of stratigraphic detail on inversion of microwave brightness temperatures

Snow provides large seasonal storage of freshwater, and information about the distribution of snow mass as Snow Water Equivalent (SWE) is important for hydrological planning and detecting climate change impacts. Large regional disagreements remain between estimates from reanalyses, remote sensing and modelling. Assimilating passive microwave information improves SWE estimates in many regions but the assimilation must account for how microwave scattering depends on snow stratigraphy. Physical snow models can estimate snow stratigraphy, but users must consider the computational expense of model complexity versus acceptable errors. Using data from the National Aeronautics and Space Administration Cold Land Processes Experiment (NASA CLPX) and the Helsinki University of Technology (HUT) microwave emission model of layered snowpacks, it is shown that simulations of the brightness temperature difference between 19 GHz and 37 GHz vertically polarised microwaves are consistent with Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and Special Sensor Microwave Imager (SSM/I) retrievals once known stratigraphic information is used. Simulated brightness temperature differences for an individual snow profile depend on the provided stratigraphic detail. Relative to a profile defined at the 10 cm resolution of density and temperature measurements, the error introduced by simplification to a single layer of average properties increases approximately linearly with snow mass. If this brightness temperature error is converted into SWE using a traditional retrieval method then it is equivalent to ±13 mm SWE (7% of total) at a depth of 100 cm. This error is reduced to ±5.6 mm SWE (3 % of total) for a two-layer model.

A new rapid detection threshold method for use with older adults: reducing fatigue whilst maintaining accuracy

Taste and smell detection threshold measurements are frequently time consuming especially when the method involves reversing the concentrations presented to replicate and improve accuracy of results. These multiple replications are likely to cause sensory and cognitive fatigue which may be more pronounced in elderly populations. A new rapid detection threshold methodology was developed that quickly located the likely position of each individuals sensory detection threshold then refined this by providing multiple concentrations around this point to determine their threshold. This study evaluates the reliability and validity of this method. Findings indicate that this new rapid detection threshold methodology was appropriate to identify differences in sensory detection thresholds between different populations and has positive benefits in providing a shorter assessment of detection thresholds. The results indicated that this method is appropriate at determining individual as well as group detection thresholds.

Survival analysis statistics applied to threshold data obtained from the ascending forced-choice method of limits

Sensory thresholds are often collected through ascending forced-choice methods. Group thresholds are important for comparing stimuli or populations; yet, the method has two problems. An individual may correctly guess the correct answer at any concentration step and might detect correctly at low concentrations but become adapted or fatigued at higher concentrations. The survival analysis method deals with both issues. Individual sequences of incorrect and correct answers are adjusted, taking into account the group performance at each concentration. The technique reduces the chance probability where there are consecutive correct answers. Adjusted sequences are submitted to survival analysis to determine group thresholds. The technique was applied to an aroma threshold and a taste threshold study. It resulted in group thresholds similar to ASTM or logarithmic regression procedures. Significant differences in taste thresholds between younger and older adults were determined. The approach provides a more robust technique over previous estimation methods.

Sex/gender differences and autism: setting the scene for future research

Objective The relationship between sex/gender differences and autism has attracted a variety of research ranging from clinical, neurobiological to etiological, stimulated by the male bias in autism prevalence. Findings are complex and do not always relate to each other in a straightforward manner. Distinct but interlinked questions on the relationship between sex/gender differences and autism remain under addressed. To better understand the implications from existing research and to help design future studies, we propose a four-level conceptual framework to clarify the embedded themes. Method We searched PubMed for publications before September 2014 using search terms “‘sex OR gender OR females’ AND autism.” 1,906 citations were screened for relevance, along with publications identified via additional literature reviews, resulting in 329 reports that were reviewed. Results Level 1 “Nosological and diagnostic challenges” concerns the question “How should autism be defined and diagnosed in males and females?” Level 2 “Sex/gender-independent and sex/gender-dependent characteristics” addresses the question “What are the similarities and differences between males and females with autism?” Level 3 “General models of etiology: liability and threshold” asks the question “How is the liability for developing autism linked to sex/gender?” Level 4 “Specific etiological-developmental mechanisms” focuses on the question “What etiological-developmental mechanisms of autism are implicated by sex/gender and/or sexual/gender differentiation?” Conclusions Using this conceptual framework, findings can be more clearly summarized, and the implications of the links between findings from different levels can become clearer. Based on this four-level framework, we suggest future research directions, methodology, and specific topics in sex/gender differences and autism.