scholarly journals Search for evidence of trend slow-down in the long-term TOMS/SBUV total ozone data record: the importance of instrument drift uncertainty

2006 ◽  
Vol 6 (12) ◽  
pp. 4057-4065 ◽  
Author(s):  
R. S. Stolarski ◽  
S. M. Frith

Abstract. We have developed a merged ozone data set (MOD) for the period October 1978 through June 2006 combining total ozone measurements (Version 8 retrieval) from the TOMS (Nimbus 7, Earth Probe) and SBUV/SBUV2 (Nimbus 7, NOAA 9/11/16) series of satellite instruments. We use the MOD data set to search for evidence of ozone recovery in response to the observed leveling off of chlorine and bromine compounds in the stratosphere. A crucial step in any time series analysis is the evaluation of uncertainties. In addition to the standard statistical time series uncertainties, we evaluate the possible instrument drift uncertainty for the MOD data set. We combine these two sources of uncertainty and apply them to a cumulative sum of residuals (CUSUM) analysis for trend slow-down. For the extra-polar mean between 60° S and 60° N, the apparent slow-down in trend is found to be clearly significant if instrument uncertainties are ignored. When instrument uncertainties are added, the slow-down becomes marginally significant at the 2σ level. For the mid-latitudes of the northern hemisphere (30° to 60° N) the trend slow-down is highly significant at the 2σ level, while in the southern hemisphere the trend slow-down has yet to meet the 2σ significance criterion. The rate of change of chlorine/bromine compounds is similar in both hemispheres, and we expect the ozone response to be similar in both hemispheres as well. The asymmetry in the trend slow-down between hemispheres likely reflects the influence of dynamical variability, and thus a clearly statistically significant response of total ozone to the leveling off of chlorine and bromine in the stratosphere is not yet indicated.

2006 ◽  
Vol 6 (3) ◽  
pp. 3883-3912 ◽  
Author(s):  
R. S. Stolarski ◽  
S. Frith

Abstract. We have developed a merged ozone data (MOD) data set for the period October 1978 through October 2005 combining total ozone measurements (version 8 retrieval) from the TOMS (Nimbus 7, Meteor 3, and Earth Probe) and SBUV/SBUV2 (Nimbus 7, NOAA 9/11/16) series of satellite instruments. We use MOD to search for evidence of ozone recovery in response to the observed leveling off of chlorine compounds in the stratosphere. A crucial step in any time series analysis is the evaluation of uncertainties. In addition to the standard statistical time-series uncertainties, we evaluate the possible instrumental drift uncertainty for the MOD data set. We combine these two sources of uncertainty and apply them to a cumulative sum of residuals (CUSUM) analysis for trend slow-down. For the quasi-global mean between 60° S and 60° N, the apparent slow-down in trend is found to be clearly significant if instrument uncertainties are ignored. When instrument uncertainties are added, the slow-down becomes marginally significant at the 2σ level. For the mid-latitudes of the northern hemisphere (30° to 60° N) the trend slow-down is significant. For the mid-latitudes of the southern hemisphere (30° to 60° S) it is not significant. The fingerprint of ozone recovery expected from model calculations suggests both northern and southern mid-latitude total ozone levels should recover together. Our result fails this fingerprint test and is therefore not a demonstration of the response of total ozone to the leveling off of chlorine.


2010 ◽  
Vol 10 (20) ◽  
pp. 10021-10031 ◽  
Author(s):  
H. E. Rieder ◽  
J. Staehelin ◽  
J. A. Maeder ◽  
T. Peter ◽  
M. Ribatet ◽  
...  

Abstract. In this study ideas from extreme value theory are for the first time applied in the field of stratospheric ozone research, because statistical analysis showed that previously used concepts assuming a Gaussian distribution (e.g. fixed deviations from mean values) of total ozone data do not adequately address the structure of the extremes. We show that statistical extreme value methods are appropriate to identify ozone extremes and to describe the tails of the Arosa (Switzerland) total ozone time series. In order to accommodate the seasonal cycle in total ozone, a daily moving threshold was determined and used, with tools from extreme value theory, to analyse the frequency of days with extreme low (termed ELOs) and high (termed EHOs) total ozone at Arosa. The analysis shows that the Generalized Pareto Distribution (GPD) provides an appropriate model for the frequency distribution of total ozone above or below a mathematically well-defined threshold, thus providing a statistical description of ELOs and EHOs. The results show an increase in ELOs and a decrease in EHOs during the last decades. The fitted model represents the tails of the total ozone data set with high accuracy over the entire range (including absolute monthly minima and maxima), and enables a precise computation of the frequency distribution of ozone mini-holes (using constant thresholds). Analyzing the tails instead of a small fraction of days below constant thresholds provides deeper insight into the time series properties. Fingerprints of dynamical (e.g. ENSO, NAO) and chemical features (e.g. strong polar vortex ozone loss), and major volcanic eruptions, can be identified in the observed frequency of extreme events throughout the time series. Overall the new approach to analysis of extremes provides more information on time series properties and variability than previous approaches that use only monthly averages and/or mini-holes and mini-highs.


2015 ◽  
Vol 8 (11) ◽  
pp. 4845-4850 ◽  
Author(s):  
R. D. McPeters ◽  
S. Frith ◽  
G. J. Labow

Abstract. The ozone data record from the Ozone Monitoring Instrument (OMI) onboard the NASA Earth Observing System (EOS) Aura satellite has proven to be very stable over the 10-plus years of operation. The OMI total column ozone processed through the Total Ozone Mapping Spectrometer (TOMS) ozone retrieval algorithm (version 8.5) has been compared with ground-based measurements and with ozone from a series of SBUV/2 (Solar Backscatter Ultraviolet) instruments. Comparison with an ensemble of Brewer–Dobson sites shows an absolute offset of about 1.5 % and almost no relative trend. Comparison with a merged ozone data set (MOD) created by combining data from a series of SBUV/2 instruments again shows an offset, of about 1 %, and a relative trend of less than 0.5 % over 10 years. The offset is mostly due to the use of the old Bass–Paur ozone cross sections in the OMI retrievals rather than the Brion–Daumont–Malicet cross sections that are now recommended. The bias in the Southern Hemisphere is smaller than that in the Northern Hemisphere, 0.9 % vs. 1.5 %, for reasons that are not completely understood. When OMI was compared with the European realization of a multi-instrument ozone time series, the GTO (GOME type Total Ozone) data set, there was a small trend of about −0.85 % decade−1. Since all the comparisons of OMI relative to other ozone measuring systems show relative trends that are less than 1 % decade−1, we conclude that the OMI total column ozone data are sufficiently stable that they can be used in studies of ozone trends.


2010 ◽  
Vol 10 (5) ◽  
pp. 12765-12794 ◽  
Author(s):  
H. E. Rieder ◽  
J. Staehelin ◽  
J. A. Maeder ◽  
T. Peter ◽  
M. Ribatet ◽  
...  

Abstract. In this study ideas from extreme value theory are for the first time applied in the field of stratospheric ozone research, because statistical analysis showed that previously used concepts assuming a Gaussian distribution (e.g. fixed deviations from mean values) of total ozone data do not adequately address the structure of the extremes. We show that statistical extreme value methods are appropriate to identify ozone extremes and to describe the tails of the Arosa (Switzerland) total ozone time series. In order to accommodate the seasonal cycle in total ozone, a daily moving threshold was determined and used, with tools from extreme value theory, to analyse the frequency of days with extreme low (termed ELOs) and high (termed EHOs) total ozone at Arosa. The analysis shows that the Generalized Pareto Distribution (GPD) provides an appropriate model for the frequency distribution of total ozone above or below a mathematically well-defined threshold, thus providing a statistical description of ELOs and EHOs. The results show an increase in ELOs and a decrease in EHOs during the last decades. The fitted model represents the tails of the total ozone data set with high accuracy over the entire range (including absolute monthly minima and maxima), and enables a precise computation of the frequency distribution of ozone mini-holes (using constant thresholds). Analyzing the tails instead of a small fraction of days below constant thresholds provides deeper insight into the time series properties. Fingerprints of dynamical (e.g. ENSO, NAO) and chemical features (e.g. strong polar vortex ozone loss), and major volcanic eruptions, can be identified in the observed frequency of extreme events throughout the time series. Overall the new approach to analysis of extremes provides more information on time series properties and variability than previous approaches that use only monthly averages and/or mini-holes and mini-highs.


2014 ◽  
Vol 14 (13) ◽  
pp. 7059-7074 ◽  
Author(s):  
W. Chehade ◽  
M. Weber ◽  
J. P. Burrows

Abstract. The study presents a long-term statistical trend analysis of total ozone data sets obtained from various satellites. A multi-variate linear regression was applied to annual mean zonal mean data using various natural and anthropogenic explanatory variables that represent dynamical and chemical processes which modify global ozone distributions in a changing climate. The study investigated the magnitude and zonal distribution of the different atmospheric chemical and dynamical factors contributing to long-term total ozone changes. The regression model included the equivalent effective stratospheric chlorine (EESC), the 11-year solar cycle, the quasi-biennial oscillation (QBO), stratospheric aerosol loading describing the effects from major volcanic eruptions, the El Niño–Southern Oscillation (ENSO), the Arctic and Antarctic oscillation (AO/AAO), and accumulated eddy heat flux (EHF), the latter representing changes due to the Brewer–Dobson circulation. The total ozone column data set used here comprises the Solar Backscater Ultraviolet SBUV/SBUV-2 merged ozone data set (MOD) V8.6, the merged data set of the Solar Backscaterr Ultraviolet, the Total Ozone Mapping Spectrometer and the Ozone Monitoring Instrument SBUV/TOMS/OMI (1979–2012) MOD V8.0 and the merged data set of the Global Ozone Monitoring Experiment, the Scanning Imaging Absorption spectroMeter for Atmospheric ChartograpHY and the Global Ozone Monitoring Experiment 2 GOME/SCIAMACHY/GOME-2 (GSG) (1995–2012). The trend analysis was performed for twenty-six 5° wide latitude bands from 65° S to 65° N, and the analysis explained most of the ozone variability to within 70 to 90%. The results show that QBO dominates the ozone variability in the tropics (±7 DU) while at higher latitudes, the dynamical indices, AO/AAO and eddy heat flux, have substantial influence on total ozone variations by up to ±10 DU. The contribution from volcanic aerosols is only prominent during the major eruption periods (El Chichón and Mt. Pinatubo), and together with the ENSO signal, is more evident in the Northern Hemisphere. The signature of the solar cycle covers all latitudes and contributes about 10 DU from solar maximum to solar minimum. EESC is found to be a main contributor to the long-term ozone decline and the trend changes after the end of the 1990s. From the EESC fits, statistically significant upward trends after 1997 were found in the extratropics, which points at the slowing of ozone decline and the onset of ozone recovery. The EESC based trends are compared with the trends obtained from the statistical piecewise linear trend (PWLT) model (known as hockey stick) with a turnaround in 1997 to examine the differences between both approaches. In case of the SBUV merged V8.6 data the EESC and PWLT trends before and after 1997 are in good agreement (within 2 σ), however, the positive post-1997 linear trends from the PWLT regression are not significant within 2 σ. A sensitivity study is carried out by comparing the regression results, using SBUV/SBUV-2 MOD V8.6 merged time series (1979–2012) and a merged data set combining SBUV/SBUV-2 (1979–June 1995) and GOME/SCIAMACHY/GOME-2 ("GSG") WFDOAS (Weighting Function DOAS) (July 1995–2012) as well as SBUV/TOMS/OMI MOD V8.0 (1979–2012) in the regression analysis in order to investigate the uncertainty in the long-term trends due to different ozone data sets and data versions. Replacing the late SBUV/SBUV-2 merged data record with GSG data (unscaled and adjusted) leads to very similar results demonstrating the high consistency between satellite data sets. However, the comparison of the new SBUV/SBUV-2 MOD V8.6 with the MOD V8.0 and MOD8.6/GSG data showed somewhat smaller sensitivities with regard to several proxies as well as the linear EESC trends. On the other hand, the PWLT trends after 1997 show some differences, however, within the 2 σ error bars the PWLT trends agree with each other for all three data sets.


2018 ◽  
Vol 10 (3) ◽  
pp. 1451-1456 ◽  
Author(s):  
Marion Maturilli ◽  
Kerstin Ebell

Abstract. Clouds are a key factor for the Arctic amplification of global warming, but their actual appearance and distribution are still afflicted by large uncertainty. On the Arctic-wide scale, large discrepancies are found between the various reanalyses and satellite products, respectively. Although ground-based observations by remote sensing are limited to point measurements, they have the advantage of obtaining extended time series of vertically resolved cloud properties. Here, we present a 25-year data record of cloud base height measured by ceilometer at the Ny-Ålesund, Svalbard, Arctic site. We explain the composition of the three sub-periods with different instrumentation contributing to the data set, and show examples of potential application areas. Linked to cyclonic activity, the cloud base height provides essential information for the interpretation of the surface radiation balance and contributes to the understanding of meteorological processes. Furthermore, it is a useful auxiliary component for the analysis of advanced technologies that provide insight into cloud microphysical properties, like the cloud radar. The long-term time series also allows derivation of an annual cycle of the cloud occurrence frequency, revealing the more frequent cloud cover in summer and the lowest cloud cover amount in April. However, as the use of different ceilometer instruments over the years potentially imposed inhomogeneities onto the data record, any long-term trend analysis should be avoided. The Ny-Ålesund cloud base height data from August 1992 to July 2017 are provided in a high temporal resolution of 5 min (1 min) before (after) July 1998, respectively, at the PANGAEA repository (https://doi.org/10.1594/PANGAEA.880300).


2008 ◽  
Vol 8 (11) ◽  
pp. 2847-2857 ◽  
Author(s):  
J. W. Krzyścin ◽  
J. L. Borkowski

Abstract. The total ozone data over Europe are available for only few ground-based stations in the pre-satellite era disallowing examination of the spatial trend variability over the whole continent. A need of having gridded ozone data for a trend analysis and input to radiative transfer models stimulated a reconstruction of the daily ozone values since January 1950. Description of the reconstruction model and its validation were a subject of our previous paper. The data base used was built within the objectives of the COST action 726 "Long-term changes and climatology of UV radiation over Europe". Here we focus on trend analyses. The long-term variability of total ozone is discussed using results of a flexible trend model applied to the reconstructed total ozone data for the period 1950–2004. The trend pattern, which comprises both anthropogenic and "natural" component, is not a priori assumed but it comes from a smooth curve fit to the zonal monthly means and monthly grid values. The ozone long-term changes are calculated separately for cold (October–next year April) and warm (May–September) seasons. The confidence intervals for the estimated ozone changes are derived by the block bootstrapping. The statistically significant negative trends are found almost over the whole Europe only in the period 1985–1994. Negative trends up to −3% per decade appeared over small areas in earlier periods when the anthropogenic forcing on the ozone layer was weak . The statistically positive trends are found only during warm seasons 1995–2004 over Svalbard archipelago. The reduction of ozone level in 2004 relative to that before the satellite era is not dramatic, i.e., up to ~−5% and ~−3.5% in the cold and warm subperiod, respectively. Present ozone level is still depleted over many popular resorts in southern Europe and northern Africa. For high latitude regions the trend overturning could be inferred in last decade (1995–2004) as the ozone depleted areas are not found there in 2004 in spite of substantial ozone depletion in the period 1985–1994.


2021 ◽  
Author(s):  
Annette Dietmaier ◽  
Thomas Baumann

<p>The European Water Framework Directive (WFD) commits EU member states to achieve a good qualitative and quantitative status of all their water bodies.  WFD provides a list of actions to be taken to achieve the goal of good status.  However, this list disregards the specific conditions under which deep (> 400 m b.g.l.) groundwater aquifers form and exist.  In particular, deep groundwater fluid composition is influenced by interaction with the rock matrix and other geofluids, and may assume a bad status without anthropogenic influences. Thus, a new concept with directions of monitoring and modelling this specific kind of aquifers is needed. Their status evaluation must be based on the effects induced by their exploitation. Here, we analyze long-term real-life production data series to detect changes in the hydrochemical deep groundwater characteristics which might be triggered by balneological and geothermal exploitation. We aim to use these insights to design a set of criteria with which the status of deep groundwater aquifers can be quantitatively and qualitatively determined. Our analysis is based on a unique long-term hydrochemical data set, taken from 8 balneological and geothermal sites in the molasse basin of Lower Bavaria, Germany, and Upper Austria. It is focused on a predefined set of annual hydrochemical concentration values. The data range dates back to 1937. Our methods include developing threshold corridors, within which a good status can be assumed, and developing cluster analyses, correlation, and piper diagram analyses. We observed strong fluctuations in the hydrochemical characteristics of the molasse basin deep groundwater during the last decades. Special interest is put on fluctuations that seem to have a clear start and end date, and to be correlated with other exploitation activities in the region. For example, during the period between 1990 and 2020, bicarbonate and sodium values displayed a clear increase, followed by a distinct dip to below-average values and a subsequent return to average values at site F. During the same time, these values showed striking irregularities at site B. Furthermore, we observed fluctuations in several locations, which come close to disqualifying quality thresholds, commonly used in German balneology. Our preliminary results prove the importance of using long-term (multiple decades) time series analysis to better inform quality and quantity assessments for deep groundwater bodies: most fluctuations would stay undetected within a < 5 year time series window, but become a distinct irregularity when viewed in the context of multiple decades. In the next steps, a quality assessment matrix and threshold corridors will be developed, which take into account methods to identify these fluctuations. This will ultimately aid in assessing the sustainability of deep groundwater exploitation and reservoir management for balneological and geothermal uses.</p>


2018 ◽  
Vol 611 ◽  
pp. A85 ◽  
Author(s):  
R. Silvotti ◽  
S. Schuh ◽  
S.-L. Kim ◽  
R. Lutz ◽  
M. Reed ◽  
...  

V391 Peg (alias HS 2201+2610) is a subdwarf B (sdB) pulsating star that shows both p- and g-modes. By studying the arrival times of the p-mode maxima and minima through the O–C method, in a previous article the presence of a planet was inferred with an orbital period of 3.2 years and a minimum mass of 3.2 MJup. Here we present an updated O–C analysis using a larger data set of 1066 h of photometric time series (~2.5× larger in terms of the number of data points), which covers the period between 1999 and 2012 (compared with 1999–2006 of the previous analysis). Up to the end of 2008, the new O–C diagram of the main pulsation frequency (f1) is compatible with (and improves) the previous two-component solution representing the long-term variation of the pulsation period (parabolic component) and the giant planet (sine wave component). Since 2009, the O–C trend of f1 changes, and the time derivative of the pulsation period (p.) passes from positive to negative; the reason of this change of regime is not clear and could be related to nonlinear interactions between different pulsation modes. With the new data, the O–C diagram of the secondary pulsation frequency (f2) continues to show two components (parabola and sine wave), like in the previous analysis. Various solutions are proposed to fit the O–C diagrams of f1 and f2, but in all of them, the sinusoidal components of f1 and f2 differ or at least agree less well than before. The nice agreement found previously was a coincidence due to various small effects that are carefully analyzed. Now, with a larger dataset, the presence of a planet is more uncertain and would require confirmation with an independent method. The new data allow us to improve the measurement of p. for f1 and f2: using only the data up to the end of 2008, we obtain p.1 = (1.34 ± 0.04) × 10−12 and p.2 = (1.62 ± 0.22) × 10−12. The long-term variation of the two main pulsation periods (and the change of sign of p.1) is visible also in direct measurements made over several years. The absence of peaks near f1 in the Fourier transform and the secondary peak close to f2 confirm a previous identification as l = 0 and l = 1, respectively, and suggest a stellar rotation period of about 40 days. The new data allow constraining the main g-mode pulsation periods of the star.


Sign in / Sign up

Export Citation Format

Share Document