Determining the effect of asymmetric data on the variogram. I. Underlying asymmetry

Matheron's usual variogram estimator can result in unreliable variograms when data are strongly asymmetric or skewed. Asymmetry in a distribution can arise from a long tail of values in the underlying process or from outliers that belong to another population that contaminate the primary process. This paper examines the effects of underlying asymmetry on the variogram and on the accuracy of prediction, and the second one examines the effects arising from outliers. Standard geostatistical texts suggest ways of dealing with underlying asymmetry; however, this is based on informed intuition rather than detailed investigation. To determine whether the methods generally used to deal with underlying asymmetry are appropriate, the effects of different coefficients of skewness on the shape of the experimental variogram and on the model parameters were investigated. Simulated annealing was used to create normally distributed random fields of different size from variograms with different nugget:sill ratios. These data were then modified to give different degrees of asymmetry and the experimental variogram was computed in each case. The effects of standard data transformations on the form of the variogram were also investigated. Cross-validation was used to assess quantitatively the performance of the different variogram models for kriging. The results showed that the shape of the variogram was affected by the degree of asymmetry, and that the effect increased as the size of data set decreased. Transformations of the data were more effective in reducing the skewness coefficient in the larger sets of data. Cross-validation confirmed that variogram models from transformed data were more suitable for kriging than were those from the raw asymmetric data. The results of this study have implications for the 'standard best practice' in dealing with asymmetry in data for geostatistical analyses. (C) 2007 Elsevier Ltd. All rights reserved.

Determining the effect of asymmetric data on the variogram. II. Outliers

Asymmetry in a distribution can arise from a long tail of values in the underlying process or from outliers that belong to another population that contaminate the primary process. The first paper of this series examined the effects of the former on the variogram and this paper examines the effects of asymmetry arising from outliers. Simulated annealing was used to create normally distributed random fields of different size that are realizations of known processes described by variograms with different nugget:sill ratios. These primary data sets were then contaminated with randomly located and spatially aggregated outliers from a secondary process to produce different degrees of asymmetry. Experimental variograms were computed from these data by Matheron's estimator and by three robust estimators. The effects of standard data transformations on the coefficient of skewness and on the variogram were also investigated. Cross-validation was used to assess the performance of models fitted to experimental variograms computed from a range of data contaminated by outliers for kriging. The results showed that where skewness was caused by outliers the variograms retained their general shape, but showed an increase in the nugget and sill variances and nugget:sill ratios. This effect was only slightly more for the smallest data set than for the two larger data sets and there was little difference between the results for the latter. Overall, the effect of size of data set was small for all analyses. The nugget:sill ratio showed a consistent decrease after transformation to both square roots and logarithms; the decrease was generally larger for the latter, however. Aggregated outliers had different effects on the variogram shape from those that were randomly located, and this also depended on whether they were aggregated near to the edge or the centre of the field. The results of cross-validation showed that the robust estimators and the removal of outliers were the most effective ways of dealing with outliers for variogram estimation and kriging. (C) 2007 Elsevier Ltd. All rights reserved.

Effect of time of day and rabbit strain on patterns of aortic wall permeability

Lipid deposits occur more frequently downstream of branch points than upstream in immature rabbit and human aortas but the opposite pattern is seen in mature vessels. These distributions correlate spatially with age-related patterns of aortic permeability, observed in rabbits, and may be determined by them. The mature but not the immature pattern of permeability is dependent on endogenous nitric oxide synthesis. Although the transport patterns have hitherto seemed robust, recent studies have given the upstream pattern in some mature rabbits but the downstream pattern in others. Here we show that transport in mature rabbits is significantly skewed to the downstream pattern in the afternoon compared with the morning (P < 0.05), and switches from a downstream to an upstream pattern at around 21 months in rabbits of the Murex strain, but at twice this age in Highgate rabbits (P < 0.001). The effect of time of day was not explained by changes in nitric oxide production, assessed from plasma levels of nitrate and nitrate, nor did it correlate with conduit artery tone, assessed from the shape of the peripheral pulse wave. The effect of strain could not be explained by variation in nitric oxide production nor by differences in wall structure. The effects of time of day and rabbit strain on permeability patterns explain recent discrepancies, provide a useful tool for investigating underlying mechanisms and may have implications for human disease.

Effect of psychological treatment on attentional bias in eating disorders

Objective: The aims of these studies were (a) to investigate the relationship between attentional bias and eating disorders and (b) examine the impact of psychological treatment on attentional bias. Method: The first study compared performance on a pictorial dot probe of 82 female patients with clinical eating disorders and 44 healthy female controls. The second study compared the performance of 31 patients with eating disorder on the same task before and after receiving 20 weeks of standardized cognitive behavior therapy. Twenty-four patients with eating disorder served as wait-list controls. Results: With the exception of neutral shape stimuli, attentional biases for eating, shape, and weight stimuli were greater in the patient sample than the healthy controls. The second study found that attentional biases significantly reduced after active treatment only. Conclusion: Attentional biases may be an expression of the eating disorder. The question of whether such biases warrant specific intervention requires further investigation. (C) 2008 by Wiley Periodicals, Inc.

The spatial and temporal shape of oculomotor inhibition

Selecting a stimulus as the target for a goal-directed movement involves inhibiting other competing possible responses. Inhibition has generally proved hard to study behaviorally, because it results in no measurable output. The effect of distractors on the shape of oculomotor and manual trajectories provide evidence of such inhibition. Individual saccades may deviate initially either towards, or away from, a competing distractor - the direction and extent of this deviation depends upon saccade latency, target predictability and the target to distractor separation. The experiment reported here used these effects to show how inhibition of distractor locations develops over time. Distractors could be presented at various distances from unpredictable and predictable targets in two separate experiments. The deviation of saccade trajectories was compared between trials with and without distractors. Inhibition was measured by saccade trajectory deviation. Inhibition was found to increase as the distractor distance from target decreased but was found to increase with saccade latency at all distractor distances (albeit to different peaks). Surprisingly, no differences were found between unpredictable and predictable targets perhaps because our saccade latencies were generally long (similar to 260-280 ms.). We conclude that oculomotor inhibition of saccades to possible target objects involves the same mechanisms for all distractor distances and target types. (C) 2009 Elsevier Ltd. All rights reserved.

Rheology of magmas with bimodal crystal size and shape distributions: insights from analog experiments

Magmas in volcanic conduits commonly contain microlites in association with preexisting phenocrysts, as often indicated by volcanic rock textures. In this study, we present two different experiments that inves- tigate the flow behavior of these bidisperse systems. In the first experiments, rotational rheometric methods are used to determine the rheology of monodisperse and polydisperse suspensions consisting of smaller, prolate particles (microlites) and larger, equant particles (phenocrysts) in a bubble‐free Newtonian liquid (silicate melt). Our data show that increasing the relative proportion of prolate microlites to equant pheno- crysts in a magma at constant total particle content can increase the relative viscosity by up to three orders of magnitude. Consequently, the rheological effect of particles in magmas cannot be modeled by assuming a monodisperse population of particles. We propose a new model that uses interpolated parameters based on the relative proportions of small and large particles and produces a considerably improved fit to the data than earlier models. In a second series of experiments we investigate the textures produced by shearing bimodal suspensions in gradually solidifying epoxy resin in a concentric cylinder setup. The resulting textures show the prolate particles are aligned with the flow lines and spherical particles are found in well‐organized strings, with sphere‐depleted shear bands in high‐shear regions. These observations may explain the measured variation in the shear thinning and yield stress behavior with increasing solid fraction and particle aspect ratio. The implications for magma flow are discussed, and rheological results and tex- tural observations are compared with observations on natural samples.

The effect of annealing on the local structure of glassy polycarbonate

Annealing of polycarbonate glasses at temperatures belowTg leads to an increase in yield stress and a drop in the impact strength. Although such behaviour may be related to the corresponding reduction in free volume upon annealing, variations in the wide-angle X-ray scattering curves indicate some modification to the local structure. The area of an intrachain peak at s ∼ 0.7 Å−1 is monitored with respect to annealing temperature and time. It is proposed that the variations may be described by an increasing level of interlocking or (nesting) between neighbouring chain segments, a process which is a natural consequence of the molecular shape of polycarbonate.

The effect of wind shear and curvature on the gravity wave drag produced by a ridge

The analytical model proposed by Teixeira, Miranda, and Valente is modified to calculate the gravity wave drag exerted by a stratified flow over a 2D mountain ridge. The drag is found to be more strongly affected by the vertical variation of the background velocity than for an axisymmetric mountain. In the hydrostatic approximation, the corrections to the drag due to this effect do not depend on the detailed shape of the ridge as long as this is exactly 2D. Besides the drag, all the perturbed quantities of the flow at the surface, including the pressure, may be calculated analytically.

Effect of shear rupture on aggregate scale formation in sea ice

A discrete element model is used to study shear rupture of sea ice under convergent wind stresses. The model includes compressive, tensile, and shear rupture of viscous elastic joints connecting floes that move under the action of the wind stresses. The adopted shear rupture is governed by Coulomb’s criterion. The ice pack is a 400 km long square domain consisting of 4 km size floes. In the standard case with tensile strength 10 times smaller than the compressive strength, under uniaxial compression the failure regime is mainly shear rupture with the most probable scenario corresponding to that with the minimum failure work. The orientation of cracks delineating formed aggregates is bimodal with the peaks around the angles given by the wing crack theory determining diamond-shaped blocks. The ice block (floe aggregate) size decreases as the wind stress gradient increases since the elastic strain energy grows faster leading to a higher speed of crack propagation. As the tensile strength grows, shear rupture becomes harder to attain and compressive failure becomes equally important leading to elongation of blocks perpendicular to the compression direction and the blocks grow larger. In the standard case, as the wind stress confinement ratio increases the failure mode changes at a confinement ratio within 0.2–0.4, which corresponds to the analytical critical confinement ratio of 0.32. Below this value, the cracks are bimodal delineating diamond shape aggregates, while above this value failure becomes isotropic and is determined by small-scale stress anomalies due to irregularities in floe shape.

The effect of climate policy on the impacts of climate change on river flows in the UK

This paper compares the effects of two indicative climate mitigation policies on river flows in six catchments in the UK with two scenarios representing un-mitigated emissions. It considers the consequences of uncertainty in both the pattern of catchment climate change as represented by different climate models and hydrological model parameterisation on the effects of mitigation policy. Mitigation policy has little effect on estimated flow magnitudes in 2030. By 2050 a mitigation policy which achieves a 2oC temperature rise target reduces impacts on low flows by 20-25% compared to a business-as-usual emissions scenario which increases temperatures by 4oC by the end of the 21st century, but this is small compared to the range in impacts between different climate model scenarios. However, the analysis also demonstrates that an early peak in emissions would reduce impacts by 40-60% by 2080 (compared with the 4oC pathway), easing the adaptation challenge over the long term, and can delay by several decades the impacts that would be experienced from around 2050 in the absence of policy. The estimated proportion of impacts avoided varies between climate model patterns and, to a lesser extent, hydrological model parameterisations, due to variations in the projected shape of the relationship between climate forcing and hydrological response.

The effect of chosen extraterrestrial solar spectrum on clear-sky atmospheric absorption and heating rates in the near infrared

The extraterrestrial solar spectrum (ESS) is an important component in near infrared (near-IR) radiative transfer calculations. However, the impact of a particular choice of the ESS in these regions has been given very little attention. A line-by-line (LBL) transfer model has been used to calculate the absorbed solar irradiance and solar heating rates in the near-IR from 2000-10000 cm−1(1-5 μm) using different ESS. For overhead sun conditions in a mid-latitude summer atmosphere, the absorbed irradiances could differ by up to about 11 Wm−2 (8.2%) while the tropospheric and stratospheric heating rates could differ by up to about 0.13 K day−1 (8.1%) and 0.19 K day−1 (7.6%). The spectral shape of the ESS also has a small but non-negligible impact on these factors in the near-IR.

Effect of uncertainty in surface mass balance--elevation feedback on projections of the future sea level contribution of the Greenland ice sheet

We apply a new parameterisation of the Greenland ice sheet (GrIS) feedback between surface mass balance (SMB: the sum of surface accumulation and surface ablation) and surface elevation in the MAR regional climate model (Edwards et al., 2014) to projections of future climate change using five ice sheet models (ISMs). The MAR (Modèle Atmosphérique Régional: Fettweis, 2007) climate projections are for 2000–2199, forced by the ECHAM5 and HadCM3 global climate models (GCMs) under the SRES A1B emissions scenario. The additional sea level contribution due to the SMB– elevation feedback averaged over five ISM projections for ECHAM5 and three for HadCM3 is 4.3% (best estimate; 95% credibility interval 1.8–6.9 %) at 2100, and 9.6% (best estimate; 95% credibility interval 3.6–16.0 %) at 2200. In all results the elevation feedback is significantly positive, amplifying the GrIS sea level contribution relative to the MAR projections in which the ice sheet topography is fixed: the lower bounds of our 95% credibility intervals (CIs) for sea level contributions are larger than the “no feedback” case for all ISMs and GCMs. Our method is novel in sea level projections because we propagate three types of modelling uncertainty – GCM and ISM structural uncertainties, and elevation feedback parameterisation uncertainty – along the causal chain, from SRES scenario to sea level, within a coherent experimental design and statistical framework. The relative contributions to uncertainty depend on the timescale of interest. At 2100, the GCM uncertainty is largest, but by 2200 both the ISM and parameterisation uncertainties are larger. We also perform a perturbed parameter ensemble with one ISM to estimate the shape of the projected sea level probability distribution; our results indicate that the probability density is slightly skewed towards higher sea level contributions.

Effect of temperature and strain rate on the compressive behaviour of supramolecular polyurethane

Supramolecular polyurethanes (SPUs) possess thermoresponsive and thermoreversible properties, and those characteristics are highly desirable in both bulk commodity and value-added applications such as adhesives, shape-memory materials, healable coatings and lightweight, impact-resistant structures (e.g. protection for mobile electronics). A better understanding of the mechanical properties, especially the rate and temperature sensitivity, of these materials are required to assess their suitability for different applications. In this paper, a newly developed SPU with tuneable thermal properties was studied, and the response of this SPU to compressive loading over strain rates from 10−3 to 104 s−1 was presented. Furthermore, the effect of temperature on the mechanical response was also demonstrated. The sample was tested using an Instron mechanical testing machine for quasi-static loading, a home-made hydraulic system for moderate rates and a traditional split Hopkinson pressure bars (SHPBs) for high strain rates. Results showed that the compression stress-strain behaviour was affected significantly by the thermoresponsive nature of SPU, but that, as expected for polymeric materials, the general trends of the temperature and the rate dependence mirror each other. However, this behaviour is more complicated than observed for many other polymeric materials, as a result of the richer range of transitions that influence the behaviour over the range of temperatures and strain rates tested.