An observing system simulation experiment for climate monitoring with GNSS radio occulation data: setup and testbed study

The long-term stability, high accuracy, all-weather capability, high vertical resolution, and global coverage of Global Navigation Satellite System (GNSS) radio occultation (RO) suggests it as a promising tool for global monitoring of atmospheric temperature change. With the aim to investigate and quantify how well a GNSS RO observing system is able to detect climate trends, we are currently performing an (climate) observing system simulation experiment over the 25-year period 2001 to 2025, which involves quasi-realistic modeling of the neutral atmosphere and the ionosphere. We carried out two climate simulations with the general circulation model MAECHAM5 (Middle Atmosphere European Centre/Hamburg Model Version 5) of the MPI-M Hamburg, covering the period 2001–2025: One control run with natural variability only and one run also including anthropogenic forcings due to greenhouse gases, sulfate aerosols, and tropospheric ozone. On the basis of this, we perform quasi-realistic simulations of RO observables for a small GNSS receiver constellation (six satellites), state-of-the-art data processing for atmospheric profiles retrieval, and a statistical analysis of temperature trends in both the “observed” climatology and the “true” climatology. Here we describe the setup of the experiment and results from a test bed study conducted to obtain a basic set of realistic estimates of observational errors (instrument- and retrieval processing-related errors) and sampling errors (due to spatial-temporal undersampling). The test bed results, obtained for a typical summer season and compared to the climatic 2001–2025 trends from the MAECHAM5 simulation including anthropogenic forcing, were found encouraging for performing the full 25-year experiment. They indicated that observational and sampling errors (both contributing about 0.2 K) are consistent with recent estimates of these errors from real RO data and that they should be sufficiently small for monitoring expected temperature trends in the global atmosphere over the next 10 to 20 years in most regions of the upper troposphere and lower stratosphere (UTLS). Inspection of the MAECHAM5 trends in different RO-accessible atmospheric parameters (microwave refractivity and pressure/geopotential height in addition to temperature) indicates complementary climate change sensitivity in different regions of the UTLS so that optimized climate monitoring shall combine information from all climatic key variables retrievable from GNSS RO data.

Jupiter's and Saturn's convectively driven banded jets in the laboratory

The banded patterns of cloud and wind are among the most striking features of the atmospheres of Jupiter and Saturn, but their dynamical origin remains poorly understood. Most approaches towards understanding zonation so far (also in the terrestrial oceans) have used highly idealized models to show that it might originate from dynamical anisotropy in a shallow turbulent fluid layer due to the planetary β-effect. Here we report the results of laboratory experiments, conducted on a 14-m diameter turntable, which quantitatively confirm that multiple zonal jets may indeed be generated and maintained by this mechanism in the presence of deep convection and a topographic β-effect. At the very small values of Ekman number (≤2 × 10−5) and large local Reynolds numbers (≥2000, based on jet scales) achieved, the kinetic energy spectra suggest the presence of both energy-cascading and enstrophy-cascading inertial ranges in addition to the zonation near twice the Rhines wave number.

Diurnal cycle of land surface temperature in a desert encroachment zone as observed from satellites

Climate variability in the African Soudano-Sahel savanna zone has attracted much attention because of the persistence of anomalously low rainfall. Past efforts to monitor the climate of this region have focused on rainfall and vegetation conditions, while land surface temperature (LST) has received less attention. Remote sensing of LST is feasible and possible at global scale. Most remotely sensed estimates of LST are based on the National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR) that are limited in their ability to capture the full diurnal cycle. Although more frequent observations are available from past geostationary satellites, their spatial resolution is coarser than that of polar orbiting satellites. In this study, the improved capabilities of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on the METEOSAT Second Generation (MSG) instrument are used to remotely sense the LST in the African Soudano-Sahel savanna zone at a resolution of 3 km and 15 minutes. In support of the Radiative Atmospheric Divergence using the ARM Mobile Facility (AMF), GERB and AMMA Stations (RADAGAST) project, African Monsoon Multidisciplinary Analyses (AMMA) project and the Department of Energy's Atmospheric Radiation Measurement (ARM) program, the ARM Mobile Facility was deployed during 2006 in this climatically sensitive region, thereby providing a unique opportunity to evaluate remotely sensed algorithms for deriving LST.

Cloud properties and their seasonal and diurnal variability from TOVS path-B

Eight years of cloud properties retrieved from Television Infrared Observation Satellite-N (TIROS-N) Observational Vertical Sounder (TOVS) observations aboard the NOAA polar orbiting satellites are presented. The relatively high spectral resolution of these instruments in the infrared allows especially reliable cirrus identification day and night. This dataset therefore provides complementary information to the International Satellite Cloud Climatology Project (ISCCP). According to this dataset, cirrus clouds cover about 27% of the earth and 45% of the Tropics, whereas ISCCP reports 19% and 25%, respectively. Both global datasets agree within 5% on the amount of single-layer low clouds, at 30%. From 1987 to 1995, global cloud amounts remained stable to within 2%. The seasonal cycle of cloud amount is in general stronger than its diurnal cycle and it is stronger than the one of effective cloud amount, the latter the relevant variable for radiative transfer. Maximum effective low cloud amount over ocean occurs in winter in SH subtropics in the early morning hours and in NH midlatitudes without diurnal cycle. Over land in winter the maximum is in the early afternoon, accompanied in the midlatitudes by thin cirrus. Over tropical land and in the other regions in summer, the maximum of mesoscale high opaque clouds occurs in the evening. Cirrus also increases during the afternoon and persists during night and early morning. The maximum of thin cirrus is in the early afternoon, then decreases slowly while cirrus and high opaque clouds increase. TOVS extends information of ISCCP during night, indicating that high cloudiness, increasing during the afternoon, persists longer during night in the Tropics and subtropics than in midlatitudes. A comparison of seasonal and diurnal cycle of high cloud amount between South America, Africa, and Indonesia during boreal winter has shown strong similarities between the two land regions, whereas the Indonesian islands show a seasonal and diurnal behavior strongly influenced by the surrounding ocean. Deeper precipitation systems over Africa than over South America do not seem to be directly reflected in the horizontal coverage and mesoscale effective emissivity of high clouds.

On the evolution of the solar wind between 1 and 5 AU at the time of the Cassini Jupiter flyby: Multispacecraft observations of interplanetary coronal mass ejections including the formation of a merged interaction region

The Cassini flyby of Jupiter occurred at a time near solar maximum. Consequently, the pre-Jupiter data set reveals clear and numerous transient perturbations to the Parker Spiral solar wind structure. Limited plasma data are available at Cassini for this period due to pointing restrictions imposed on the instrument. This renders the identification of the nature of such structures ambiguous, as determinations based on the magnetic field data alone are unreliable. However, a fortuitous alignment of the planets during this encounter allowed us to trace these structures back to those observed previously by the Wind spacecraft near the Earth. Of the phenomena that we are satisfactorily able to trace back to their manifestation at 1 AU, two are identified as being due to interplanetary coronal mass ejections. One event at Cassini is shown to be a merged interaction region, which is formed from the compression of a magnetic cloud by two anomalously fast solar wind streams. The flux-rope structure associated with this magnetic cloud is not as apparent at Cassini and has most likely been compressed and deformed. Confirmation of the validity of the ballistic projections used here is provided by results obtained from a one-dimensional magnetohydrodynamic projection of solar wind parameters measured upstream near the Earth. It is found that when the Earth and Cassini are within a few tens of degrees in heliospheric longitude, the results of this one-dimensional model predict the actual conditions measured at 5 AU to an impressive degree. Finally, the validity of the use of such one-dimensional projections in obtaining quasi-solar wind parameters at the outer planets is discussed.

Global climate and ocean circulation on an aquaplanet ocean-atmosphere general circulation model

A low resolution coupled ocean-atmosphere general circulation model OAGCM is used to study the characteristics of the large scale ocean circulation and its climatic impacts in a series of global coupled aquaplanet experiments. Three configurations, designed to produce fundamentally different ocean circulation regimes, are considered. The first has no obstruction to zonal flow, the second contains a low barrier that blocks zonal flow in the ocean at all latitudes, creating a single enclosed basin, whilst the third contains a gap in the barrier to allow circumglobal flow at high southern latitudes. Warm greenhouse climates with a global average air surface temperature of around 27C result in all cases. Equator to pole temperature gradients are shallower than that of a current climate simulation. Whilst changes in the land configuration cause regional changes in temperature, winds and rainfall, heat transports within the system are little affected. Inhibition of all ocean transport on the aquaplanet leads to a reduction in global mean surface temperature of 8C, along with a sharpening of the meridional temperature gradient. This results from a reduction in global atmospheric water vapour content and an increase in tropical albedo, both of which act to reduce global surface temperatures. Fitting a simple radiative model to the atmospheric characteristics of the OAGCM solutions suggests that a simpler atmosphere model, with radiative parameters chosen a priori based on the changing surface configuration, would have produced qualitatively different results. This implies that studies with reduced complexity atmospheres need to be guided by more complex OAGCM results on a case by case basis.

Vibrational intensities VII: ethane and ethane-d6

Absolute intensity measurements have been made on the fundamental vibrations of C2H6 and C2D6, using the extrapolation method of Wilson and Wells and using nitrogen at pressures up to 50 atmospheres to broaden the bands. The absorption coefficient was integrated against the logarithm of the frequency. Normal coordinates were calculated from the potential function of Hansen and Dennison, and were used to interpret the results in terms of quantities (∂p/∂Si) giving the change of dipole moment with respect to the symmetry coordinates Si. Consistency of data between the isotopes was used both to eliminate ambiguities in the interpretation, and as a criterion in separating overlapping pairs of absorption bands. The results have been interpreted in terms of bond effective moments.

Controls on boundary layer ventilation: Boundary layer processes and large-scale dynamics

Midlatitude cyclones are important contributors to boundary layer ventilation. However, it is uncertain how efficient such systems are at transporting pollutants out of the boundary layer, and variations between cyclones are unexplained. In this study 15 idealized baroclinic life cycles, with a passive tracer included, are simulated to identify the relative importance of two transport processes: horizontal divergence and convergence within the boundary layer and large-scale advection by the warm conveyor belt. Results show that the amount of ventilation is insensitive to surface drag over a realistic range of values. This indicates that although boundary layer processes are necessary for ventilation they do not control the magnitude of ventilation. A diagnostic for the mass flux out of the boundary layer has been developed to identify the synoptic-scale variables controlling the strength of ascent in the warm conveyor belt. A very high level of correlation (R-2 values exceeding 0.98) is found between the diagnostic and the actual mass flux computed from the simulations. This demonstrates that the large-scale dynamics control the amount of ventilation, and the efficiency of midlatitude cyclones to ventilate the boundary layer can be estimated using the new mass flux diagnostic. We conclude that meteorological analyses, such as ERA-40, are sufficient to quantify boundary layer ventilation by the large-scale dynamics.

Laboratory measurements of the water vapor continuum in the 1200–8000 cm−1 region between 293 K and 351 K

Laboratory Fourier transform spectroscopy of pure water vapor and water vapor mixed with air has been conducted between 1200 and 8000 cm−1 and at temperatures between 293 and 351 K with the purpose of detecting and characterizing the water vapor continuum. The spectral features of the continuum within the major water absorption bands are presented and compared where possible to those from previous experimental studies and to the commonly used MT_CKD and CKD models. It was observed that in the main, both models adequately capture the general spectral form of the continuum; however, there were a number of exceptions. Overall, there is no evidence to indicate that MT_CKD is an improvement upon the older CKD model in these spectral regions. There was generally good agreement between our results and those of other experimental investigators. The general mathematical forms of the self-continuum temperature dependence, given by both Roberts et al. (1976) and CKD/MT_CKD, fit well to the experimental continuum in these spectral regions. However, the range of temperatures over which we made measurements is not sufficient to discriminate between these two forms or to exclude the possibility of other forms of temperature dependence being more appropriate. At the same time, the actual parameters currently used in CKD/MT_CKD to describe the temperature dependence in many spectral regions cannot reproduce the observed strong spectral variation in the temperature dependence. It has not been possible to make definitive conclusions about the magnitude of the continuum absorption in the far wings of the absorption bands investigated here.

Integrated optics for nulling interferometry in the thermal infrared: progress and recent achievements

The search for Earth-like exoplanets, orbiting in the habitable zone of stars other than our Sun and showing biological activity, is one of the most exciting and challenging quests of the present time. Nulling interferometry from space, in the thermal infrared, appears as a promising candidate technique for the task of directly observing extra-solar planets. It has been studied for about 10 years by ESA and NASA in the framework of the Darwin and TPF-I missions respectively. Nevertheless, nulling interferometry in the thermal infrared remains a technological challenge at several levels. Among them, the development of the "modal filter" function is mandatory for the filtering of the wavefronts in adequacy with the objective of rejecting the central star flux to an efficiency of about 105. Modal filtering takes benefit of the capability of single-mode waveguides to transmit a single amplitude function, to eliminate virtually any perturbation of the interfering wavefronts, thus making very high rejection ratios possible. The modal filter may either be based on single-mode Integrated Optics (IO) and/or Fiber Optics. In this paper, we focus on IO, and more specifically on the progress of the on-going "Integrated Optics" activity of the European Space Agency.

Estimating relationships between air mass origin and chemical composition

Observations of a chemical at a point in the atmosphere typically show sudden transitions between episodes of high and low concentration. Often these are associated with a rapid change in the origin of air arriving at the site. Lagrangian chemical models riding along trajectories can reproduce such transitions, but small timing errors from trajectory phase errors dramatically reduce the correlation between modeled concentrations and observations. Here the origin averaging technique is introduced to obtain maps of average concentration as a function of air mass origin for the East Atlantic Summer Experiment 1996 (EASE96, a ground-based chemistry campaign). These maps are used to construct origin averaged time series which enable comparison between a chemistry model and observations with phase errors factored out. The amount of the observed signal explained by trajectory changes can be quantified, as can the systematic model errors as a function of air mass origin. The Cambridge Tropospheric Trajectory model of Chemistry and Transport (CiTTyCAT) can account for over 70% of the observed ozone signal variance during EASE96 when phase errors are side-stepped by origin averaging. The dramatic increase in correlation (from 23% without averaging) cannot be achieved by time averaging. The success of the model is attributed to the strong relationship between changes in ozone along trajectories and their origin and its ability to simulate those changes. The model performs less well for longer-lived chemical constituents because the initial conditions 5 days before arrival are insufficiently well known.

A simple relationship between volcanic sulfate aerosol optical depth and surface temperature change simulated in an atmosphere-ocean general circulation model.

In this study we quantify the relationship between the aerosol optical depth increase from a volcanic eruption and the severity of the subsequent surface temperature decrease. This investigation is made by simulating 10 different sizes of eruption in a global circulation model (GCM) by changing stratospheric sulfate aerosol optical depth at each time step. The sizes of the simulated eruptions range from Pinatubo‐sized up to the magnitude of supervolcanic eruptions around 100 times the size of Pinatubo. From these simulations we find that there is a smooth monotonic relationship between the global mean maximum aerosol optical depth anomaly and the global mean temperature anomaly and we derive a simple mathematical expression which fits this relationship well. We also construct similar relationships between global mean aerosol optical depth and the temperature anomaly at every individual model grid box to produce global maps of best‐fit coefficients and fit residuals. These maps are used with caution to find the eruption size at which a local temperature anomaly is clearly distinct from the local natural variability and to approximate the temperature anomalies which the model may simulate following a Tambora‐sized eruption. To our knowledge, this is the first study which quantifies the relationship between aerosol optical depth and resulting temperature anomalies in a simple way, using the wealth of data that is available from GCM simulations.

A comparison among four different retrieval methods for ice-cloud properties using data from CloudSat, CALIPSO, and MODIS

The A-Train constellation of satellites provides a new capability to measure vertical cloud profiles that leads to more detailed information on ice-cloud microphysical properties than has been possible up to now. A variational radar–lidar ice-cloud retrieval algorithm (VarCloud) takes advantage of the complementary nature of the CloudSat radar and Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) lidar to provide a seamless retrieval of ice water content, effective radius, and extinction coefficient from the thinnest cirrus (seen only by the lidar) to the thickest ice cloud (penetrated only by the radar). In this paper, several versions of the VarCloud retrieval are compared with the CloudSat standard ice-only retrieval of ice water content, two empirical formulas that derive ice water content from radar reflectivity and temperature, and retrievals of vertically integrated properties from the Moderate Resolution Imaging Spectroradiometer (MODIS) radiometer. The retrieved variables typically agree to within a factor of 2, on average, and most of the differences can be explained by the different microphysical assumptions. For example, the ice water content comparison illustrates the sensitivity of the retrievals to assumed ice particle shape. If ice particles are modeled as oblate spheroids rather than spheres for radar scattering then the retrieved ice water content is reduced by on average 50% in clouds with a reflectivity factor larger than 0 dBZ. VarCloud retrieves optical depths that are on average a factor-of-2 lower than those from MODIS, which can be explained by the different assumptions on particle mass and area; if VarCloud mimics the MODIS assumptions then better agreement is found in effective radius and optical depth is overestimated. MODIS predicts the mean vertically integrated ice water content to be around a factor-of-3 lower than that from VarCloud for the same retrievals, however, because the MODIS algorithm assumes that its retrieved effective radius (which is mostly representative of cloud top) is constant throughout the depth of the cloud. These comparisons highlight the need to refine microphysical assumptions in all retrieval algorithms and also for future studies to compare not only the mean values but also the full probability density function.