A standardized approach for the definition and reporting of vertical resolution of the ozone and temperature lidar profiles contributing to the Network for the Detection for Atmospheric Composition ...Change (NDACC) database is proposed. Two standardized definitions homogeneously and unequivocally describing the impact of vertical filtering are recommended. The first proposed definition is based on the width of the response to a finite-impulse-type perturbation. The response is computed by convolving the filter coefficients with an impulse function, namely, a Kronecker delta function for smoothing filters, and a Heaviside step function for derivative filters. Once the response has been computed, the proposed standardized definition of vertical resolution is given by Δz = δz × HFWHM, where δz is the lidar's sampling resolution and HFWHM is the full width at half maximum (FWHM) of the response, measured in sampling intervals. The second proposed definition relates to digital filtering theory. After applying a Laplace transform to a set of filter coefficients, the filter's gain characterizing the effect of the filter on the signal in the frequency domain is computed, from which the cut-off frequency fC, defined as the frequency at which the gain equals 0.5, is computed. Vertical resolution is then defined by Δz = δz∕(2fC). Unlike common practice in the field of spectral analysis, a factor 2fC instead of fC is used here to yield vertical resolution values nearly equal to the values obtained with the impulse response definition using the same filter coefficients. When using either of the proposed definitions, unsmoothed signals yield the best possible vertical resolution Δz = δz (one sampling bin). Numerical tools were developed to support the implementation of these definitions across all NDACC lidar groups. The tools consist of ready-to-use “plug-in” routines written in several programming languages that can be inserted into any lidar data processing software and called each time a filtering operation occurs in the data processing chain. When data processing implies multiple smoothing operations, the filtering information is analytically propagated through the multiple calls to the routines in order for the standardized values of vertical resolution to remain theoretically and numerically exact at the very end of data processing.
We present a comprehensive analysis of the trends of stratospheric ozone in the midlatitudes and subtropics. The analysis is performed using ground‐based and space‐based measurements over the light ...detection and ranging stations for the period 1985–2012. Also, trends are estimated for the zonal mean data made from a merged satellite data set, Global OZone Chemistry And Related trace gas Data records for the Stratosphere, over 1979–2012. The linear trends in stratospheric ozone are estimated using piecewise linear trend (PWLT) functions. The ozone trends during the increasing phase of halogens (before 1997) range from −0.2 ± 0.08 to −1 ± 0.07% yr−1 in the midlatitudes and −0.2 ± 0.06 to −0.7 ± 0.05 % yr−1 in the subtropics at 15–45 km, depending on altitude. In 1997–2012, the PWLT analyses show a positive trend, significantly different from zero at the 95% confidence intervals, toward ozone recovery in the middle‐ and low‐latitude upper stratosphere (35–45 km), and the trends are about +0.5 ± 0.07% yr−1 at midlatitudes and about +0.3 ± 0.05% yr−1 at subtropical latitudes. However, negative and insignificant trends are estimated in the lower stratosphere (15–20 km) over 1997–2012 in the midlatitudes, mainly due to the dynamics, as demonstrated by the large (50–60%) contributions from the quasi‐biennial oscillation, El Niño–Southern Oscillation, and planetary wave activity to recent ozone changes. This suggests that the ozone changes are governed by the interannual variations in meteorology and dynamics of the regions; these factors will influence the recovery detection time and the behavior of the recovery path to pre‐1980 levels.
Key Points
Presents the midlatitude and subtropical ozone trends
Shows clear recovery signal in the upper stratosphere
Lower stratospheric recovery is now masked by the dynamics
A standardized approach for the definition, propagation, and reporting of uncertainty in the ozone differential absorption lidar data products contributing to the Network for the Detection for ...Atmospheric Composition Change (NDACC) database is proposed. One essential aspect of the proposed approach is the propagation in parallel of all independent uncertainty components through the data processing chain before they are combined together to form the ozone combined standard uncertainty. The independent uncertainty components contributing to the overall budget include random noise associated with signal detection, uncertainty due to saturation correction, background noise extraction, the absorption cross sections of O3, NO2, SO2, and O2, the molecular extinction cross sections, and the number densities of the air, NO2, and SO2. The expression of the individual uncertainty components and their step-by-step propagation through the ozone differential absorption lidar (DIAL) processing chain are thoroughly estimated. All sources of uncertainty except detection noise imply correlated terms in the vertical dimension, which requires knowledge of the covariance matrix when the lidar signal is vertically filtered. In addition, the covariance terms must be taken into account if the same detection hardware is shared by the lidar receiver channels at the absorbed and non-absorbed wavelengths. The ozone uncertainty budget is presented as much as possible in a generic form (i.e., as a function of instrument performance and wavelength) so that all NDACC ozone DIAL investigators across the network can estimate, for their own instrument and in a straightforward manner, the expected impact of each reviewed uncertainty component. In addition, two actual examples of full uncertainty budget are provided, using nighttime measurements from the tropospheric ozone DIAL located at the Jet Propulsion Laboratory (JPL) Table Mountain Facility, California, and nighttime measurements from the JPL stratospheric ozone DIAL located at Mauna Loa Observatory, Hawai'i.
The ozone profile records of a large number of limb and occultation satellite instruments are widely used to address several key questions in ozone research. Further progress in some domains depends ...on a more detailed understanding of these data sets, especially of their long-term stability and their mutual consistency. To this end, we made a systematic assessment of 14 limb and occultation sounders that, together, provide more than three decades of global ozone profile measurements. In particular, we considered the latest operational Level-2 records by SAGE II, SAGE III, HALOE, UARS MLS, Aura MLS, POAM II, POAM III, OSIRIS, SMR, GOMOS, MIPAS, SCIAMACHY, ACE-FTS and MAESTRO. Central to our work is a consistent and robust analysis of the comparisons against the ground-based ozonesonde and stratospheric ozone lidar networks. It allowed us to investigate, from the troposphere up to the stratopause, the following main aspects of satellite data quality: long-term stability, overall bias and short-term variability, together with their dependence on geophysical parameters and profile representation. In addition, it permitted us to quantify the overall consistency between the ozone profilers. Generally, we found that between 20 and 40 km the satellite ozone measurement biases are smaller than ±5 %, the short-term variabilities are less than 5–12 % and the drifts are at most ±5 % decade−1 (or even ±3 % decade−1 for a few records). The agreement with ground-based data degrades somewhat towards the stratopause and especially towards the tropopause where natural variability and low ozone abundances impede a more precise analysis. In part of the stratosphere a few records deviate from the preceding general conclusions; we identified biases of 10 % and more (POAM II and SCIAMACHY), markedly higher single-profile variability (SMR and SCIAMACHY) and significant long-term drifts (SCIAMACHY, OSIRIS, HALOE and possibly GOMOS and SMR as well). Furthermore, we reflected on the repercussions of our findings for the construction, analysis and interpretation of merged data records. Most notably, the discrepancies between several recent ozone profile trend assessments can be mostly explained by instrumental drift. This clearly demonstrates the need for systematic comprehensive multi-instrument comparison analyses.
A standardized approach for the definition, propagation, and reporting of uncertainty in the temperature lidar data products contributing to the Network for the Detection for Atmospheric Composition ...Change (NDACC) database is proposed. One important aspect of the proposed approach is the ability to propagate all independent uncertainty components in parallel through the data processing chain. The individual uncertainty components are then combined together at the very last stage of processing to form the temperature combined standard uncertainty. The identified uncertainty sources comprise major components such as signal detection, saturation correction, background noise extraction, temperature tie-on at the top of the profile, and absorption by ozone if working in the visible spectrum, as well as other components such as molecular extinction, the acceleration of gravity, and the molecular mass of air, whose magnitudes depend on the instrument, data processing algorithm, and altitude range of interest. The expression of the individual uncertainty components and their step-by-step propagation through the temperature data processing chain are thoroughly estimated, taking into account the effect of vertical filtering and the merging of multiple channels. All sources of uncertainty except detection noise imply correlated terms in the vertical dimension, which means that covariance terms must be taken into account when vertical filtering is applied and when temperature is integrated from the top of the profile. Quantitatively, the uncertainty budget is presented in a generic form (i.e., as a function of instrument performance and wavelength), so that any NDACC temperature lidar investigator can easily estimate the expected impact of individual uncertainty components in the case of their own instrument. Using this standardized approach, an example of uncertainty budget is provided for the Jet Propulsion Laboratory (JPL) lidar at Mauna Loa Observatory, Hawai'i, which is typical of the NDACC temperature lidars transmitting at 355 nm. The combined temperature uncertainty ranges between 0.1 and 1 K below 60 km, with detection noise, saturation correction, and molecular extinction correction being the three dominant sources of uncertainty. Above 60 km and up to 10 km below the top of the profile, the total uncertainty increases exponentially from 1 to 10 K due to the combined effect of random noise and temperature tie-on. In the top 10 km of the profile, the accuracy of the profile mainly depends on that of the tie-on temperature. All other uncertainty components remain below 0.1 K throughout the entire profile (15–90 km), except the background noise correction uncertainty, which peaks around 0.3–0.5 K. It should be kept in mind that these quantitative estimates may be very different for other lidar instruments, depending on their altitude range and the wavelengths used.
The validation of ozone profiles retrieved by satellite instruments through comparison with data from ground-based instruments is important to monitor the evolution of the satellite instrument, to ...assist algorithm development and to allow multi-mission trend analyses. In this study we compare ozone profiles derived from GOMOS night-time observations with measurements from lidar, microwave radiometer and balloon sonde. Collocated pairs are analysed for dependence on several geophysical and instrument observational parameters. Validation results are presented for the operational ESA level 2 data (GOMOS version 5.00) obtained during nearly seven years of observations and a comparison using a smaller dataset from the previous processor (version 4.02) is also included. The profiles obtained from dark limb measurements (solar zenith angle >107°) when the provided processing flag is properly considered match the ground-based measurements within ±2 percent over the altitude range 20 to 40 km. Outside this range, the pairs start to deviate more and there is a latitudinal dependence: in the polar region where there is a higher amount of straylight contamination, differences start to occur lower in the mesosphere than in the tropics, whereas for the lower part of the stratosphere the opposite happens: the profiles in the tropics reach less far down as the signal reduces faster because of the higher altitude at which the maximum ozone concentration is found compared to the mid and polar latitudes. Also the bias is shifting from mostly negative in the polar region to more positive in the tropics Profiles measured under "twilight" conditions are often matching the ground-based measurements very well, but care has to be taken in all cases when dealing with "straylight" contaminated profiles. For the selection criteria applied here (data within 800 km, 3 degrees in equivalent latitude, 20 h (5 h above 50 km) and a relative ozone error in the GOMOS data of 20% or less), no dependence was found on stellar magnitude, star temperature, nor the azimuth angle of the line of sight. No evidence of a temporal trend was seen either in the bias or frequency of outliers, but a comparison applying less strict data selection criteria might show differently.
Upper stratospheric ozone anomalies from the satellite-borne Solar Backscatter Ultra-Violet (SBUV), Stratospheric Aerosol and Gas Experiment II (SAGE II), Halogen Occultation Experiment (HALOE), ...Global Ozone Monitoring by Occultation of Stars (GOMOS), and Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) instruments agree within 5% or better with ground-based data from lidars and microwave radiometers at five stations of the Network for the Detection of Atmospheric Composition Change (NDACC), from 45°S to 48°N. From 1979 until the late 1990s, all available data show a clear decline of ozone near 40 km, by 10%-15%. This decline has not continued in the last 10 years. At some sites, ozone at 40 km appears to have increased since 2000, consistent with the beginning decline of stratospheric chlorine. The phaseout of chlorofluorocarbons after the International Montreal Protocol in 1987 has been successful, and is now showing positive effects on ozone in the upper stratosphere. Temperature anomalies near 40 km altitude from European Centre for Medium Range Weather Forecast reanalyses (ERA-40), from National Centers for Environmental Prediction (NCEP) operational analyses, and from HALOE and lidar measurements show good consistency at the five stations, within about 3 K. Since about 1985, upper stratospheric temperatures have been fluctuating around a constant level at all five NDACC stations. This non-decline of upper stratospheric temperatures is a significant change from the more or less linear cooling of the upper stratosphere up until the mid-1990s, reported in previous trend assessments. It is also at odds with the almost linear 1 K per decade cooling simulated over the entire 1979-2010 period by chemistry-climate models (CCMs). The same CCM simulations, however, track the historical ozone anomalies quite well, including the change of ozone tendency in the late 1990s.
Traditional validation of atmospheric profiles is based on the intercomparison of two or more data sets in predefined ranges or classes of a given observational characteristic such as latitude or ...solar zenith angle. In this study we trained a self-organising map (SOM) with a full time series of relative difference profiles of SCIAMACHY limb v5.02 and lidar ozone profiles from seven observation sites. Each individual observation characteristic was then mapped to the obtained SOM to investigate to which degree variation in this characteristic is explanatory for the variation seen in the SOM map. For the studied data sets, altitude-dependent relations for the global data set were found between the difference profiles and studied variables. From the lowest altitude studied (18 km) ascending, the most influencing factors were found to be longitude, followed by solar zenith angle and latitude, sensor age and again solar zenith angle together with the day of the year at the highest altitudes studied here (up to 45 km). After accounting for both latitude and longitude, residual partial correlations with a reduced magnitude are seen for various factors. However, (partial) correlations cannot point out which (combination) of the factors drives the observed differences between the ground-based and satellite ozone profiles as most of the factors are inter-related. Clustering into three classes showed that there are also some local dependencies, with for instance one cluster having a much stronger correlation with the sensor age (days since launch) between 36 and 42 km. The proposed SOM-based approach provides a powerful tool for the exploration of differences between data sets without being limited to a priori defined data subsets.
The long-term evolution of stratospheric ozone at different stations in the low and mid-latitudes is investigated. The analysis is performed by comparing the collocated profiles of ozone lidars, at ...the northern mid-latitudes (Meteorological Observatory Hohenpeißenberg, Haute-Provence Observatory, Tsukuba and Table Mountain Facility), tropics (Mauna Loa Observatory) and southern mid-latitudes (Lauder), with ozonesondes and space-borne sensors (SBUV(/2), SAGE II, HALOE, UARS MLS and Aura MLS), extracted around the stations. Relative differences are calculated to find biases and temporal drifts in the measurements. All measurement techniques show their best agreement with respect to the lidar at 20–40 km, where the differences and drifts are generally within ±5% and ±0.5% yr−1, respectively, at most stations. In addition, the stability of the long-term ozone observations (lidar, SBUV(/2), SAGE II and HALOE) is evaluated by the cross-comparison of each data set. In general, all lidars and SBUV(/2) exhibit near-zero drifts and the comparison between SAGE II and HALOE shows larger, but insignificant drifts. The RMS of the drifts of lidar and SBUV(/2) is 0.22 and 0.27% yr−1, respectively at 20–40 km. The average drifts of the long-term data sets, derived from various comparisons, are less than ±0.3% yr−1 in the 20–40 km altitude at all stations. A combined time series of the relative differences between SAGE II, HALOE and Aura MLS with respect to lidar data at six sites is constructed, to obtain long-term data sets lasting up to 27 years. The relative drifts derived from these combined data are very small, within ±0.2% yr−1.
In this paper we evaluate the potential of ENVISAT-Medium Resolution Imaging Spectrometer (MERIS) fused images for land-cover mapping and vegetation status assessment in heterogeneous landscapes. A ...series of MERIS fused images (15 spectral bands; 25 m pixel size) is created using the linear mixing model and a Landsat Thematic Mapper (TM) image acquired over the Netherlands. First, the fused images are classified to produce a map of the eight main land-cover types of the Netherlands. Subsequently, the maps are validated using the Dutch land-cover/land-use database as a reference. Then, the fused image with the highest overall classification accuracy is selected as the best fused image. Finally, the best fused image is used to compute three vegetation indices: the normalized difference vegetation index (NDVI) and two indices specifically designed to monitor vegetation status using MERIS data: the MERIS terrestrial chlorophyll index (MTCI) and the MERIS global vegetation index (MGVI).
Results indicate that the selected data fusion approach is able to downscale MERIS data to a Landsat-like spatial resolution. The spectral information in the fused images originates fully from MERIS and is not influenced by the TM data. Classification results for the TM and for the best fused image are similar and, when comparing spectrally similar images (i.e. TM with no short-wave infrared bands), the results of the fused image outperform those of TM. With respect to the vegetation indices, a good correlation was found between the NDVI computed from TM and from the best fused image (in spite of the spectral differences between these two sensors). In addition, results show the potential of using MERIS vegetation indices computed from fused images to monitor individual fields. This is not possible using the original MERIS full resolution image. Therefore, we conclude that MERIS-TM fused images are very useful to map heterogeneous landscapes.