Impact of long-range correlations on trend detection in total ozone

Total ozone trends are typically studied using linear regression models that assume a first-order autoregression of the residuals [so-called AR(1) models]. We consider total ozone time series over 60°S–60°N from 1979 to 2005 and show that most latitude bands exhibit long-range correlated (LRC) behavior, meaning that ozone autocorrelation functions decay by a power law rather than exponentially as in AR(1). At such latitudes the uncertainties of total ozone trends are greater than those obtained from AR(1) models and the expected time required to detect ozone recovery correspondingly longer. We find no evidence of LRC behavior in southern middle-and high-subpolar latitudes (45°–60°S), where the long-term ozone decline attributable to anthropogenic chlorine is the greatest. We thus confirm an earlier prediction based on an AR(1) analysis that this region (especially the highest latitudes, and especially the South Atlantic) is the optimal location for the detection of ozone recovery, with a statistically significant ozone increase attributable to chlorine likely to be detectable by the end of the next decade. In northern middle and high latitudes, on the other hand, there is clear evidence of LRC behavior. This increases the uncertainties on the long-term trend attributable to anthropogenic chlorine by about a factor of 1.5 and lengthens the expected time to detect ozone recovery by a similar amount (from ∼2030 to ∼2045). If the long-term changes in ozone are instead fit by a piecewise-linear trend rather than by stratospheric chlorine loading, then the strong decrease of northern middle- and high-latitude ozone during the first half of the 1990s and its subsequent increase in the second half of the 1990s projects more strongly on the trend and makes a smaller contribution to the noise. This both increases the trend and weakens the LRC behavior at these latitudes, to the extent that ozone recovery (according to this model, and in the sense of a statistically significant ozone increase) is already on the verge of being detected. The implications of this rather controversial interpretation are discussed.

O_3-N_2O correlations from the Atmospheric Chemistry Experiment: revisiting a diagnostic of transport and chemistry in the stratosphere

Our knowledge of stratospheric O3-N2O correlations is extended, and their potential for model-measurement comparison assessed, using data from the Atmospheric Chemistry Experiment (ACE) satellite and the Canadian Middle Atmosphere Model (CMAM). ACE provides the first comprehensive data set for the investigation of interhemispheric, interseasonal, and height-resolved differences of the O_3-N_2O correlation structure. By subsampling the CMAM data, the representativeness of the ACE data is evaluated. In the middle stratosphere, where the correlations are not compact and therefore mainly reflect the data sampling, joint probability density functions provide a detailed picture of key aspects of transport and mixing, but also trace polar ozone loss. CMAM captures these important features, but exhibits a displacement of the tropical pipe into the Southern Hemisphere (SH). Below about 21 km, the ACE data generally confirm the compactness of the correlations, although chemical ozone loss tends to destroy the compactness during late winter/spring, especially in the SH. This allows a quantitative comparison of the correlation slopes in the lower and lowermost stratosphere (LMS), which exhibit distinct seasonal cycles that reveal the different balances between diabatic descent and horizontal mixing in these two regions in the Northern Hemisphere (NH), reconciling differences found in aircraft measurements, and the strong role of chemical ozone loss in the SH. The seasonal cycles are qualitatively well reproduced by CMAM, although their amplitude is too weak in the NH LMS. The correlation slopes allow a "chemical" definition of the LMS, which is found to vary substantially in vertical extent with season.

Correlations of long-lived chemical species in a middle atmosphere general circulation model

Correlations between various chemical species simulated by the Canadian Middle Atmosphere Model, a general circulation model with fully interactive chemistry, are considered in order to investigate the general conditions under which compact correlations can be expected to form. At the same time, the analysis serves to validate the model. The results are compared to previous work on this subject, both from theoretical studies and from atmospheric measurements made from space and from aircraft. The results highlight the importance of having a data set with good spatial coverage when working with correlations and provide a background against which the compactness of correlations obtained from atmospheric measurements can be confirmed. It is shown that for long-lived species, distinct correlations are found in the model in the tropics, the extratropics, and the Antarctic winter vortex. Under these conditions, sparse sampling such as arises from occultation instruments is nevertheless suitable to define a chemical correlation within each region even from a single day of measurements, provided a sufficient range of mixing ratio values is sampled. In practice, this means a large vertical extent, though the requirements are less stringent at more poleward latitudes.

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.

Cross-correlations and cross-bicorrelations in Sterling exchange rates

This paper proposes two new tests for linear and nonlinear lead/lag relationships between time series based on the concepts of cross-correlations and cross-bicorrelations, respectively. The tests are then applied to a set of Sterling-denominated exchange rates. Our analysis indicates that there existed periods during the post-Bretton Woods era where the temporal relationship between different exchange rates was strong, although these periods have become less frequent over the past 20 years. In particular, our results demonstrate the episodic nature of the nonlinearity, and have implications for the speed of flow of information between financial series. The method generalises recently proposed tests for nonlinearity to the multivariate context.

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.

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.

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.

Drizzle and Fog Analysis in the Sao Paulo Metropolitan Area: Changes 1933-2005 and Correlations with other Climate Factors

In order to improve our understanding of climate change, the aim of this research project was to study the climatology and the time trends of drizzle and fog events in the Sao Paulo Metropolitan Area, and the possible connections of this variability with the sea surface temperature (SST) of the Atlantic and Pacific Oceans. The climatology of both phenomena presents differences and similarities. Fog shows a marked maximum frequency in winter and a minimum frequency in summer, while the seasonal differences of drizzle occurrence are less pronounced, there is a maximum in spring, whereas the other seasons present smaller and similar numbers of events. Both phenomena present a negative trend from 1933 to 2005 which is stronger for fog events. A multivariate statistical analysis indicates that the South Atlantic SST could increase warm temperature advection to the continent. This could be one of the responsible factors for the negative tendency in the number of both fog and drizzle events.

Clinical correlations of human cytomegalovirus strains and viral load in kidney transplant recipients

Little is known about clinical differences associated with cytomegalovirus (CMV) infection by distinct strains in renal transplant patients. Different clinical pictures may be associated with specific viral genotypes. viral load, as well as host factors. The objective of this study was to identify CMV strains to determine viral load (antigenemia), and their correlation with clinical data in renal transplant recipients. Seventy-one patients were enrolled, comprising 91 samples. After selection, polymorphonuclear cells were used to amplify and sequence the gB region of CMV DNA. The sequences were analyzed to ascertain the frequency of different genotypes. Additionally, the results of this Study showed that the gB coding gene presents a great variability, revealing a variety of patterns: classical gB (1.4%), gB1V (46.4%), classical gB2 (35.2%), gB2V (2.8%), gB3 (1.4%), classical gB4 (4.9%) and gB4V (4.9%). The mean viral load in kidney transplant patient was 75.1 positive cells (1-1000). A higher viral load was observed in patients with genotype 4 infection. Statistically significant differences were detected between gB1 and gB4 (p=0.010), and between gB2 and gB4 (p=0.021). The average numbers of positive cells in relation to clinical presentation were: 34.5 in asymptomatic, 49.5 in CMV associated syndrome and 120.7 in patients with invasive disease (p=0.048). As a group, gB1 was the most frequent strain and revealed a potential risk for developing invasive disease. Viral load also seemed to be important as a marker associated with clinical presentation of the disease. (C) 2008 Elsevier B.V. All rights reserved.

A closer look at the influence of tubular initial conditions on two-particle correlations

In a recent paper, the hydrodynamic code NEXSPheRIO was used in conjunction with STAR analysis methods to study two-particle correlations as a function of Delta(eta) and Delta phi. The various structures observed in the data were reproduced. In this work, we discuss the origin of these structures as well as present new results.

Hydrodynamics: Fluctuating initial conditions and two-particle correlations

Event-by-event hydrodynamics (or hydrodynamics with fluctuating initial conditions) has been developed in the past few years. Here we discuss how it may help to understand the various structures observed in two-particle correlations. (C) 2010 Elsevier B.V. All rights reserved.

Protecting clean critical points by local disorder correlations

We show that a broad class of quantum critical points can be stable against locally correlated disorder even if they are unstable against uncorrelated disorder. Although this result seemingly contradicts the Harris criterion, it follows naturally from the absence of a random-mass term in the associated order parameter field theory. We illustrate the general concept with explicit calculations for quantum spin-chain models. Instead of the infinite-randomness physics induced by uncorrelated disorder, we find that weak locally correlated disorder is irrelevant. For larger disorder, we find a line of critical points with unusual properties such as an increase of the entanglement entropy with the disorder strength. We also propose experimental realizations in the context of quantum magnetism and cold-atom physics. Copyright (C) EPLA, 2011

Electronic Structure and Spectro-Structural Correlations of Fe(III)Zn(II) Biomimetics for Purple Acid Phosphatases: Relevance to DNA Cleavage and Cytotoxic Activity

Purple acid phosphatases (PAPs) are a group of metallohydrolases that contain a dinuclear Fe(II)M(II) center (M(II) = Fe, Mn, Zn) in the active site and are able to catalyze the hydrolysis of a variety of phosphoric acid esters. The dinuclear complex [(H(2)O)Fe(III)(mu-OH)Zn(II)(L-H)](CIO(4))(2) (2) with the ligand 2-[N-bis(2-pyridylmethyl)aminomethyl]-4-methyl-6-[N-(2-pyridylmethyl)(2-hydroxybenzyl) aminomethyl]phenol (H(2)L-H) has recently been prepared and is found to closely mimic the coordination environment of the Fe(III)Zn(II) active site found in red kidney bean PAP (Neves et al. J. Am. Chem. Soc. 2007, 129, 7486). The biomimetic shows significant catalytic activity in hydrolytic reactions. By using a variety of structural, spectroscopic, and computational techniques the electronic structure of the Fe(III) center of this biomimetic complex was determined. In the solid state the electronic ground state reflects the rhombically distorted Fe(III)N(2)O(4) octahedron with a dominant tetragonal compression align ad along the mu-OH-Fe-O(phenolate) direction. To probe the role of the Fe-O(phenolate) bond, the phenolate moiety was modified to contain electron-donating or -withdrawing groups (-CH(3), -H, -Br, -NO(2)) in the 5-position. Tie effects of the substituents on the electronic properties of the biomimetic complexes were studied with a range of experimental and computational techniques. This study establishes benchmarks against accurate crystallographic struck ral information using spectroscopic techniques that are not restricted to single crystals. Kinetic studies on the hydrolysis reaction revealed that the phosphodiesterase activity increases in the order -NO(2)<- Br <- H <- CH(3) when 2,4-bis(dinitrophenyl)phosphate (2,4-bdnpp) was used as substrate, and a linear free energy relationship is found when log(k(cat)/k(0)) is plotted against the Hammett parameter a. However, nuclease activity measurements in the cleavage of double stranded DNA showed that the complexes containing the electron-withdrawing -NO(2) and electron-donating CH3 groups are the most active while the cytotoxic activity of the biomimetics on leukemia and lung tumoral cells is highest for complexes with electron-donating groups.