On the statistical modeling of persistence in total ozone anomalies

Geophysical time series sometimes exhibit serial correlations that are stronger than can be captured by the commonly used first‐order autoregressive model. In this study we demonstrate that a power law statistical model serves as a useful upper bound for the persistence of total ozone anomalies on monthly to interannual timescales. Such a model is usually characterized by the Hurst exponent. We show that the estimation of the Hurst exponent in time series of total ozone is sensitive to various choices made in the statistical analysis, especially whether and how the deterministic (including periodic) signals are filtered from the time series, and the frequency range over which the estimation is made. In particular, care must be taken to ensure that the estimate of the Hurst exponent accurately represents the low‐frequency limit of the spectrum, which is the part that is relevant to long‐term correlations and the uncertainty of estimated trends. Otherwise, spurious results can be obtained. Based on this analysis, and using an updated equivalent effective stratospheric chlorine (EESC) function, we predict that an increase in total ozone attributable to EESC should be detectable at the 95% confidence level by 2015 at the latest in southern midlatitudes, and by 2020–2025 at the latest over 30°–45°N, with the time to detection increasing rapidly with latitude north of this range.

Decline and recovery of total column ozone using a multimodel time series analysis

Simulations of 15 coupled chemistry climate models, for the period 1960–2100, are presented. The models include a detailed stratosphere, as well as including a realistic representation of the tropospheric climate. The simulations assume a consistent set of changing greenhouse gas concentrations, as well as temporally varying chlorofluorocarbon concentrations in accordance with observations for the past and expectations for the future. The ozone results are analyzed using a nonparametric additive statistical model. Comparisons are made with observations for the recent past, and the recovery of ozone, indicated by a return to 1960 and 1980 values, is investigated as a function of latitude. Although chlorine amounts are simulated to return to 1980 values by about 2050, with only weak latitudinal variations, column ozone amounts recover at different rates due to the influence of greenhouse gas changes. In the tropics, simulated peak ozone amounts occur by about 2050 and thereafter total ozone column declines. Consequently, simulated ozone does not recover to values which existed prior to the early 1980s. The results also show a distinct hemispheric asymmetry, with recovery to 1980 values in the Northern Hemisphere extratropics ahead of the chlorine return by about 20 years. In the Southern Hemisphere midlatitudes, ozone is simulated to return to 1980 levels only 10 years ahead of chlorine. In the Antarctic, annually averaged ozone recovers at about the same rate as chlorine in high latitudes and hence does not return to 1960s values until the last decade of the simulations.

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.

Summertime total ozone variations over middle and polar latitudes

The statistical relationship between springtime and summertime ozone over middle and polar latitudes is analyzed using zonally averaged total ozone data. Shortterm variations in springtime midlatitude ozone demonstrate only a modest correlation with springtime polar ozone variations. However by early summer, ozone variations throughout the extratropics are highly correlated. Analysis of correlation functions indicates that springtime midlatitude ozone, not polar ozone, is the best predictor for summertime polar ozone. Long-term total ozone trends at middle and high latitudes are also different for spring and nearly identical for summer. About 39% of the observed southern midlatitude ozone decline in December can be attributed to the polar ozone depletion up to November. In the Northern Hemisphere, the corresponding contribution is about 15%, but the error bars are too large to make an accurate estimate.

Persistence and photochemical decay of springtime total ozone anomalies in the Canadian Middle Atmosphere Model

The persistence and decay of springtime total ozone anomalies over the entire extratropics (midlatitudes plus polar regions) is analysed using results from the Canadian Middle Atmosphere Model (CMAM), a comprehensive chemistry-climate model. As in the observations, interannual anomalies established through winter and spring persist with very high correlation coefficients (above 0.8) through summer until early autumn, while decaying in amplitude as a result of photochemical relaxation in the quiescent summertime stratosphere. The persistence and decay of the ozone anomalies in CMAM agrees extremely well with observations, even in the southern hemisphere when the model is run without heterogeneous chemistry (in which case there is no ozone hole and the seasonal cycle of ozone is quite different from observations). However in a version of CMAM with strong vertical diffusion, the northern hemisphere anomalies decay far too rapidly compared to observations. This shows that ozone anomaly persistence and decay does not depend on how the springtime anomalies are created or on their magnitude, but reflects the transport and photochemical decay in the model. The seasonality of the long-term trends over the entire extratropics is found to be explained by the persistence of the interannual anomalies, as in the observations, demonstrating that summertime ozone trends reflect winter/spring trends rather than any change in summertime ozone chemistry. However this mechanism fails in the northern hemisphere midlatitudes because of the relatively large impact, compared to observations, of the CMAM polar anomalies. As in the southern hemisphere, the influence of polar ozone loss in CMAM increases the midlatitude summertime loss, leading to a relatively weak seasonal dependence of ozone loss in the Northern Hemisphere compared to the observations.

Seasonal persistence of midlatitude total ozone anomalies

Temporal autocorrelations of monthly mean total ozone anomalies over the 35–60°S and 35–60°N latitude bands reveal that anomalies established in the wintertime midlatitude ozone buildup persist (with photochemical decay) until the end of the following autumn, and then are rapidly erased once the next winter's buildup begins. The photochemical decay rate is found to be identical between the two hemispheres. High predictability of ozone through late summer exists based on the late-spring values. In the northern hemisphere, extending the 1979–2001 springtime ozone trend to other months through regression based on the seasonal persistence of anomalies captures the seasonality of the ozone trends remarkably well. In the southern hemisphere, the springtime trend only accounts for part of the summertime trends. There is a strong correlation between the ozone anomalies in northern hemisphere spring and those in the subsequent southern hemisphere spring, but not the converse.

Geographical variation in total and inorganic arsenic content of polished (white) rice

An extensive data set of total arsenic analysis for 901 polished (white) grain samples, originating from 10 countries from 4 continents, was compiled. The samples represented the baseline (i.e., not specifically collected from arsenic contaminated areas), and all were for market sale in major conurbations. Median total arsenic contents of rice varied 7-fold, with Egypt (0.04 mg/kg) and India (0.07 mg/kg) having the lowest arsenic content while the U.S. (0.25 mg/kg) and France (0.28 mg/kg) had the highest content. Global distribution of total arsenic in rice was modeled by weighting each country’s arsenic distribution by that country’s contribution to global production. A subset of 63 samples from Bangladesh, China, India, Italy, and the U.S. was analyzed for arsenic species. The relationship between inorganic arsenic content versus total arsenic content significantly differed among countries, with Bangladesh and India having the steepest slope in linear regression, and the U.S. having the shallowest slope. Using country-specific rice consumption data, daily intake of inorganic arsenic was estimated and the associated internal cancer risk was calculated using the U.S. Environmental Protection Agency (EPA) cancer slope. Median excess internal cancer risks posed by inorganic arsenic ranged 30-fold for the 5 countries examined, being 0.7 per 10,000 for Italians to 22 per 10,000 for Bangladeshis, when a 60 kg person was considered.

Sensitivity of global sulphate aerosol production to changes in oxidant concentrations and climate

The oxidation of SO2 to sulphate aerosol is an important process to include in climate models, and uncertainties caused by ignoring feedback mechanisms affecting the oxidants concerned need to be investigated. Here we present the results of an investigation into the sensitivity of sulphate concentrations to oxidant changes (from changes in climate and in emissions of oxidant precursors) and to changes in climate, in a version of HadGAM1 (the atmosphere-only version of HadGEM1) with an improved sulphur cycle scheme. We find that, when oxidants alone are changed, the global total sulphate burden decreases by approximately 3%, due mainly to a reduction in the OH burden. When climate alone is changed, our results show that the global total sulphate burden increases by approximately 9%; we conclude that this is probably attributable to reduced precipitation in regions of high sulphate abundance. When both oxidants and climate are changed simultaneously, we find that the effects of the two changes combine approximately linearly.

Assessment of the dietary intake of total flavan-3-ols, monomeric flavan-3-ols, proanthocyanidins and theaflavins in the European Union

Dietary interventions with flavan-3-ols have shown beneficial effects on vascular function. The translation of these findings into the context of the health of the general public requires detailed information on habitual dietary intake. However, only limited data are currently available for European populations. Therefore, in the present study, we assessed the habitual intake of flavan-3-ol monomers, proanthocyanidins (PA) and theaflavins in the European Union (EU) and determined their main food sources using the EFSA (European Food Safety Authority) Comprehensive European Food Consumption Database. Data for adults aged 18–64 years were available from fourteen European countries, and intake was determined using the FLAVIOLA Flavanol Food Composition Database, developed for the present study and based on the latest US Department of Agriculture and Phenol-Explorer databases. The mean habitual intake of flavan-3-ol monomers, theaflavins and PA ranged from 181 mg/d (Czech Republic) to 793 mg/d (Ireland). The highest intakes of flavan-3-ol monomers and theaflavins were observed in Ireland (191/505 mg/d) and the lowest intakes in Spain (24/9 mg/d). In contrast, the daily intake of PA was highest in Spain (175 mg/d) and lowest in The Netherlands (96 mg/d). Main sources were tea (62 %), pome fruits (11 %), berries (3 %) and cocoa products (3 %). Tea was the major single contributor to monomer intake (75 %), followed by pome fruits (6 %). Pome fruits were also the main source of PA (28 %). The present study provides important data on the population-based intake of flavanols in the EU and demonstrates that dietary intake amounts for flavan-3-ol monomers, PA and theaflavins vary significantly across European countries. The average habitual intake of flavan-3-ols is considerably below the amounts used in most dietary intervention studies.

Contributions to twentieth century total column ozone change from halocarbons, tropospheric ozone precursors, and climate change

We investigate ozone changes from preindustrial times to the present using a chemistry-climate model. The influence of changes in physical climate, ozone-depleting substances, N2O, and tropospheric ozone precursors is estimated using equilibrium simulations with these different factors set at either preindustrial or present-day values. When these effects are combined, the entire decrease in total column ozone from preindustrial to present day is very small (–1.8 DU) in the global annual average, though with significant decreases in total column ozone over large parts of the Southern Hemisphere during austral spring and widespread increases in column ozone over the Northern Hemisphere during boreal summer. A significant contribution to the total ozone column change is the increase in lower stratospheric ozone associated with the increase in ozone precursors (5.9 DU). Also noteworthy is the near cancellation of the global average climate change effect on ozone (3.5 DU) by the increase in N2O (–3.9 DU).

The effect of adding stinging nettle (Urtica dioica) haylage to a total mixed ration on performance and rumen function of lactating dairy cows

In vitro studies found that inclusion of dried stinging nettle (Urtica dioica) at 100 mg/g dry matter (DM) increased the pH of a rumen fluid inoculated fermentation buffer by 30% and the effect was persistent for 7 days. Our objective was to evaluate the effects of adding stinging nettle haylage to a total mixed ration on feed intake, eating and rumination activity, rumen pH, milk yield, and milk composition of lactating dairy cows. Six lactating Holstein-Friesian cows were used in a replicated 3 × 3 Latin Square design experiment with 3 treatments and 3 week periods. Treatments were a control (C) high-starch (311 g/kg DM) total mixed ration diet and two treatment diets containing 50 (N5) and 100 (N10) g nettle haylage (DM/kg) as a replacement for ryegrass silage (Lolium perenne). There was an increase (linear, P < 0.010) in the proportion of large particles and a reduction in medium (linear, P = 0.045) and fine particles (linear, P = 0.026) in the diet offered with increasing nettle inclusion. A numerical decrease (linear, P = 0.106) in DM intake (DMI) was observed as nettle inclusion in the diet increased. Milk yield averaged 20.3 kg/day and was not affected by diet. There was a decrease (quadratic, P = 0.01) in the time animals spent ruminating as nettle inclusion in the diet increased, in spite of an increase in the number of boli produced daily for the N5 diet (quadratic, P = 0.031). Animals fed the N10 diet spent less time with a rumen pH below 5.5 (P < 0.05) than cows fed the N5 diet. Averaged over an 8.5 h sampling period, there were no changes in the concentration or proportions of acetate or propionate in the rumen, but feeding nettle haylage reduced the concentrations of n-butyrate (quadratic, P < 0.001), i-butyrate (linear, P < 0.009) and n-caproate (linear, P < 0.003). Milk and fat and protein corrected milk yield were not affected when nettles replaced ryegrass silage in the diet of lactating dairy cows, despite a numerical reduction in feed intake. Rumination activity was reduced by the addition of nettle haylage to the diet, which may reflect differences in fibre structure between the nettle haylage and ryegrass silage fed. Changes observed in rumen pH suggest potential benefits of feeding nettle haylage for reducing rumen acidosis. However, the extent to which these effects were due to the fermentability and structure of the nettle haylage compared to the ryegrass silage fed, or a bioactive component of the nettles, is not certain

Reconciliation of halogen-induced ozone loss with the total-column ozone record

The observed depletion of the ozone layer from the 1980s onwards is attributed to halogen source gases emitted by human activities. However, the precision of this attribution is complicated by year-to-year variations in meteorology, that is, dynamical variability, and by changes in tropospheric ozone concentrations. As such, key aspects of the total-column ozone record, which combines changes in both tropospheric and stratospheric ozone, remain unexplained, such as the apparent absence of a decline in total-column ozone levels before 1980, and of any long-term decline in total-column ozone levels in the tropics. Here we use a chemistry–climate model to estimate changes in halogen-induced ozone loss between 1960 and 2010; the model is constrained by observed meteorology to remove the effects of dynamical variability, and driven by emissions of tropospheric ozone precursors to separate out changes in tropospheric ozone. We show that halogen-induced ozone loss closely followed stratospheric halogen loading over the studied period. Pronounced enhancements in ozone loss were apparent in both hemispheres following the volcanic eruptions of El Chichon and, in particular, Mount Pinatubo, which significantly enhanced stratospheric aerosol loads. We further show that approximately 40% of the long-term non-volcanic ozone loss occurred before 1980, and that long-term ozone loss also occurred in the tropical stratosphere. Finally, we show that halogen-induced ozone loss has declined by over 10% since stratospheric halogen loading peaked in the late 1990s, indicating that the recovery of the ozone layer is well underway.

The impacts of European and Asian anthropogenic sulphur dioxide emissions on Sahel rainfall

In this study, the atmospheric component of a state-of-the-art climate model (HadGEM2-ES) has been used to investigate the impacts of regional anthropogenic sulphur dioxide emissions on boreal summer Sahel rainfall. The study focuses on the transient response of the West African monsoon (WAM) to a sudden change in regional anthropogenic sulphur dioxide emissions, including land surface feedbacks, but without sea surface temperature (SST) feedbacks. The response occurs in two distinct phases: 1) fast adjustment of the atmosphere on a time scale of days to weeks (up to 3 weeks) through aerosol-radiation and aerosol-cloud interactions with weak hydrological cycle changes and surface feedbacks. 2) adjustment of the atmosphere and land surface with significant local hydrological cycle changes and changes in atmospheric circulation (beyond 3 weeks). European emissions lead to an increase in shortwave (SW) scattering by increased sulphate burden, leading to a decrease in surface downward SW radiation which causes surface cooling over North Africa, a weakening of the Saharan heat low and WAM, and a decrease in Sahel precipitation. In contrast, Asian emissions lead to very little change in sulphate burden over North Africa, but they induce an adjustment of the Walker Circulation which leads again to a weakening of the WAM and a decrease in Sahel precipitation. The responses to European and Asian emissions during the second phase exhibit similar large scale patterns of anomalous atmospheric circulation and hydrological variables, suggesting a preferred response. The results support the idea that sulphate aerosol emissions contributed to the observed decline in Sahel precipitation in the second half of the twentieth century.

The inclusion of forage mixtures in the diet of growing dairy heifers: impacts on digestion, energy utilisation, and methane emissions

Intensive farming focusing on monoculture grass species to maximise forage production has led to a reduction in the extent and diversity of species-rich grasslands. However, plant communities with higher species number (richness) are a potential strategy for more sustainable production and mitigation of greenhouse gas (GHG) emissions. Research has indicated the need to understand opportunities that forage mixtures can offer sustainable ruminant production systems. The objective of the two experiments reported here were to evaluate multiple species forage mixtures in comparison to ryegrass-dominant pasture, when conserved or grazed, on digestion, energy utilisation, N excretion, and methane emissions by growing 10–15 month old heifers. Experiment 1 was a 4 × 4 Latin square design with five week periods. Four forage treatments of: (1) ryegrass (control); permanent pasture with perennial ryegrass (Lolium perenne); (2) clover; a ryegrass:red clover (Trifolium pratense) mixture; (3) trefoil; a ryegrass:birdsfoot trefoil (Lotus corniculatus) mixture; and (4) flowers; a ryegrass:wild flower mixture of predominately sorrel (Rumex acetosa), ox-eye daisy (Leucanthemum vulgare), yarrow (Achillea millefolium), knapweed (Centaurea nigra) and ribwort plantain (Plantago lanceolata), were fed as haylages to four dairy heifers. Measurements included digestibility, N excretion, and energy utilisation (including methane emissions measured in respiration chambers). Experiment 2 used 12 different dairy heifers grazing three of the same forage treatments used to make haylage in experiment 1 (ryegrass, clover and flowers) and methane emissions were estimated using the sulphur hexafluoride (SF6) tracer technique. Distribution of ryegrass to other species (dry matter (DM) basis) was approximately 70:30 (clover), 80:20 (trefoil), and 40:60 (flowers) for experiment 1. During the first and second grazing rotations (respectively) in experiment 2, perennial ryegrass accounted for 95 and 98% of DM in ryegrass, and 84 and 52% of DM in clover, with red clover accounting for almost all of the remainder. In the flowers mixture, perennial ryegrass was 52% of the DM in the first grazing rotation and only 30% in the second, with a variety of other flower species occupying the remainder. Across both experiments, compared to the forage mixtures (clover, trefoil and flowers), ryegrass had a higher crude protein (CP) content (P < 0.001, 187 vs. 115 g kg −1 DM) and DM intake (P < 0.05, 9.0 vs. 8.1 kg day −1). Heifers in experiment 1 fed ryegrass, compared to the forage mixtures, had greater total tract digestibility (g kg −1) of DM (DMD; P < 0.008, 713 vs. 641) and CP (CPD, P < 0.001, 699 vs. 475), and used more intake energy (%) for body tissue deposition (P < 0.05, 2.6 vs. −4.9). For both experiments, heifers fed flowers differed the most compared to the ryegrass control for a number of measurements. Compared to ryegrass, flowers had 40% lower CP content (P < 0.001, 113 vs. 187 g kg −1), 18% lower DMD (P < 0.01, 585 vs. 713 g kg −1), 42% lower CPD (P < 0.001, 407 vs. 699 g kg −1), and 10% lower methane yield (P < 0.05, 22.6 vs. 25.1 g kg −1 DM intake). This study has shown inclusion of flowers in forage mixtures resulted in a lower CP concentration, digestibility and intake. These differences were due in part to sward management and maturity at harvest. Further research is needed to determine how best to exploit the potential environmental benefits of forage mixtures in sustainable ruminant production systems.