The effect of doubled CO2 and model basic state biases on the monsoon-ENSO system. II: Changing ENSO regimes

Integrations of a fully-coupled climate model with and without flux adjustments in the equatorial oceans are performed under 2×CO2 conditions to explore in more detail the impact of increased greenhouse gas forcing on the monsoon-ENSO system. When flux adjustments are used to correct some systematic model biases, ENSO behaviour in the modelled future climate features distinct irregular and periodic (biennial) regimes. Comparison with the observed record yields some consistency with ENSO modes primarily based on air-sea interaction and those dependent on basinwide ocean wave dynamics. Simple theory is also used to draw analogies between the regimes and irregular (stochastically forced) and self-excited oscillations respectively. Periodic behaviour is also found in the Asian-Australian monsoon system, part of an overall biennial tendency of the model under these conditions related to strong monsoon forcing and increased coupling between the Indian and Pacific Oceans. The tropospheric biennial oscillation (TBO) thus serves as a useful descriptor for the coupled monsoon-ENSO system in this case. The presence of obvious regime changes in the monsoon-ENSO system on interdecadal timescales, when using flux adjustments, suggests there may be greater uncertainty in projections of future climate, although further modelling studies are required to confirm the realism and cause of such changes.

Evaluation of the Met Office global forecast model using Geostationary Earth Radiation Budget (GERB) data

Simulations of the top-of-atmosphere radiative-energy budget from the Met Office global numerical weather-prediction model are evaluated using new data from the Geostationary Earth Radiation Budget (GERB) instrument on board the Meteosat-8 satellite. Systematic discrepancies between the model simulations and GERB measurements greater than 20 Wm-2 in outgoing long-wave radiation (OLR) and greater than 60 Wm-2 in reflected short-wave radiation (RSR) are identified over the period April-September 2006 using 12 UTC data. Convective cloud over equatorial Africa is spatially less organized and less reflective than in the GERB data. This bias depends strongly on convective-cloud cover, which is highly sensitive to changes in the model convective parametrization. Underestimates in model OLR over the Gulf of Guinea coincide with unrealistic southerly cloud outflow from convective centres to the north. Large overestimates in model RSR over the subtropical ocean, greater than 50 Wm-2 at 12 UTC, are explained by unrealistic radiative properties of low-level cloud relating to overestimation of cloud liquid water compared with independent satellite measurements. The results of this analysis contribute to the development and improvement of parametrizations in the global forecast model.

Origins of dry air in the Tropics and Subtropics

The humidity in the dry regions of the tropical and subtropical troposphere has a major impact on the ability of the atmosphere to radiate heat to space. The water vapour content in these regions is determined by their ``origins'', here defined as the last condensation event following air masses. Trajectory simulations are used to investigate such origins using ERA40 data for January 1993. It is shown that 96% of air parcels experience condensation within 24 days and most of the remaining 4% originate in the stratosphere. Dry air masses are shown to experience a net pressure increase since last condensation which is uniform with latitude, while the median time taken for descent is 5 days into the subtropics but exceeds 16 days into the equatorial lower troposphere. The associated rate of decrease in potential temperature is consistent with radiative cooling. The relationship between the drier regions in the tropics and subtropics and the geographical localization of their origin is investigated. Four transport processes are identified to explain these relationships.

The role of the basic state in the ENSO-monsoon relationship and implications for predictability

The impact of systematic model errors on a coupled simulation of the Asian Summer monsoon and its interannual variability is studied. Although the mean monsoon climate is reasonably well captured, systematic errors in the equatorial Pacific mean that the monsoon-ENSO teleconnection is rather poorly represented in the GCM. A system of ocean-surface heat flux adjustments is implemented in the tropical Pacific and Indian Oceans in order to reduce the systematic biases. In this version of the GCM, the monsoon-ENSO teleconnection is better simulated, particularly the lag-lead relationships in which weak monsoons precede the peak of El Nino. In part this is related to changes in the characteristics of El Nino, which has a more realistic evolution in its developing phase. A stronger ENSO amplitude in the new model version also feeds back to further strengthen the teleconnection. These results have important implications for the use of coupled models for seasonal prediction of systems such as the monsoon, and suggest that some form of flux correction may have significant benefits where model systematic error compromises important teleconnections and modes of interannual variability.

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.

An overview of Earth's global electric circuit and atmospheric conductivity

The Earth’s global atmospheric electric circuit depends on the upper and lower atmospheric boundaries formed by the ionosphere and the planetary surface. Thunderstorms and electrified rain clouds drive a DC current (∼1 kA) around the circuit, with the current carried by molecular cluster ions; lightning phenomena drive the AC global circuit. The Earth’s near-surface conductivity ranges from 10−7 S m−1 (for poorly conducting rocks) to 10−2 S m−1 (for clay or wet limestone), with a mean value of 3.2 S m−1 for the ocean. Air conductivity inside a thundercloud, and in fair weather regions, depends on location (especially geomagnetic latitude), aerosol pollution and height, and varies from ∼10−14 S m−1 just above the surface to 10−7 S m−1 in the ionosphere at ∼80 km altitude. Ionospheric conductivity is a tensor quantity due to the geomagnetic field, and is determined by parameters such as electron density and electron–neutral particle collision frequency. In the current source regions, point discharge (coronal) currents play an important role below electrified clouds; the solar wind-magnetosphere dynamo and the unipolar dynamo due to the terrestrial rotating dipole moment also apply atmospheric potential differences. Detailed measurements made near the Earth’s surface show that Ohm’s law relates the vertical electric field and current density to air conductivity. Stratospheric balloon measurements launched from Antarctica confirm that the downward current density is ∼1 pA m−2 under fair weather conditions. Fortuitously, a Solar Energetic Particle (SEP) event arrived at Earth during one such balloon flight, changing the observed atmospheric conductivity and electric fields markedly. Recent modelling considers lightning discharge effects on the ionosphere’s electric potential (∼+250 kV with respect to the Earth’s surface) and hence on the fair weather potential gradient (typically ∼130 V m−1 close to the Earth’s surface. We conclude that cloud-to-ground (CG) lightning discharges make only a small contribution to the ionospheric potential, and that sprites (namely, upward lightning above energetic thunderstorms) only affect the global circuit in a miniscule way. We also investigate the effects of mesoscale convective systems on the global circuit.

Water vapor variability in the tropics and its links to dynamics and precipitation

The distribution and variability of water vapor and its links with radiative cooling and latent heating via precipitation are crucial to understanding feedbacks and processes operating within the climate system. Column-integrated water vapor (CWV) and additional variables from the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year reanalysis (ERA40) are utilized to quantify the spatial and temporal variability in tropical water vapor over the period 1979–2001. The moisture variability is partitioned between dynamical and thermodynamic influences and compared with variations in precipitation provided by the Climate Prediction Center Merged Analysis of Precipitation (CMAP) and the Global Precipitation Climatology Project (GPCP). The spatial distribution of CWV is strongly determined by thermodynamic constraints. Spatial variability in CWV is dominated by changes in the large-scale dynamics, in particular associated with the El Niño–Southern Oscillation (ENSO). Trends in CWV are also dominated by dynamics rather than thermodynamics over the period considered. However, increases in CWV associated with changes in temperature are significant over the equatorial east Pacific when analyzing interannual variability and over the north and northwest Pacific when analyzing trends. Significant positive trends in CWV tend to predominate over the oceans while negative trends in CWV are found over equatorial Africa and Brazil. Links between changes in CWV and vertical motion fields are identified over these regions and also the equatorial Atlantic. However, trends in precipitation are generally incoherent and show little association with the CWV trends. This may in part reflect the inadequacies of the precipitation data sets and reanalysis products when analyzing decadal variability. Though the dynamic component of CWV is a major factor in determining precipitation variability in the tropics, in some regions/seasons the thermodynamic component cancels its effect on precipitation variability.

Vertical current flow through extensive layer clouds

The global atmospheric electrical circuit sustains a vertical current density between the ionosphere and the Earth's surface, the existence of which is well-established from measurements made in fair-weather conditions. In overcast, but non-thunderstorm, non-precipitating conditions, the current travels through the cloud present, despite cloud layers having low electrical conductivity. For extensive layer clouds, this leads to space charge at the upper and lower cloud boundaries. Using a combination of atmospheric electricity and solar radiation measurements at three UK sites, vertical current measurements have been categorised into clear, broken, and overcast cloud conditions. This approach shows that the vertical “fair weather” current is maintained despite the presence of cloud. In fully overcast conditions with thick cloud, the vertical current is reduced compared to thin cloud overcast conditions, associated with the cloud's resistance contributions. Contribution of cloud to the columnar resistance depends both on cloud thickness, and the cloud's height.

Observed atmospheric electricity effect on clouds

The atmosphere's fair weather electric field is a permanent feature, arising from the combination of distant thunderstorms, Earth's conducting surface, a charged ionosphere and cosmic ray ionization. Despite its ubiquity, no fair weather electricity effect on clouds has been hitherto demonstrated. Here we report surface measurements of radiation emitted and scattered by extensive thin continental cloud, which, after ~2 min delay, shows changes closely following the fair weather electric field. For typical fluctuations in the fair weather electric field, changes of about 10% are subsequently induced in the diffuse short-wave radiation. These observations are consistent with enhanced production of large cloud droplets from charging at layer cloud edges.

Impact of Resolution on the Tropical Pacific Circulation in a Matrix of Coupled Models

Results are presented from a matrix of coupled model integrations, using atmosphere resolutions of 135 and 90 km, and ocean resolutions of 1° and 1/3°, to study the impact of resolution on simulated climate. The mean state of the tropical Pacific is found to be improved in the models with a higher ocean resolution. Such an improved mean state arises from the development of tropical instability waves, which are poorly resolved at low resolution; these waves reduce the equatorial cold tongue bias. The improved ocean state also allows for a better simulation of the atmospheric Walker circulation. Several sensitivity studies have been performed to further understand the processes involved in the different component models. Significantly decreasing the horizontal momentum dissipation in the coupled model with the lower-resolution ocean has benefits for the mean tropical Pacific climate, but decreases model stability. Increasing the momentum dissipation in the coupled model with the higher-resolution ocean degrades the simulation toward that of the lower-resolution ocean. These results suggest that enhanced ocean model resolution can have important benefits for the climatology of both the atmosphere and ocean components of the coupled model, and that some of these benefits may be achievable at lower ocean resolution, if the model formulation allows.

Stratospheric temperature and radiative forcing response to 11-year solar cycle changes in irradiance and ozone

The 11-yr solar cycle temperature response to spectrally resolved solar irradiance changes and associated ozone changes is calculated using a fixed dynamical heating (FDH) model. Imposed ozone changes are from satellite observations, in contrast to some earlier studies. A maximum of 1.6 K is found in the equatorial upper stratosphere and a secondary maximum of 0.4 K in the equatorial lower stratosphere, forming a double peak in the vertical. The upper maximum is primarily due to the irradiance changes while the lower maximum is due to the imposed ozone changes. The results compare well with analyses using the 40-yr ECMWF Re-Analysis (ERA-40) and NCEP/NCAR datasets. The equatorial lower stratospheric structure is reproduced even though, by definition, the FDH calculations exclude dynamically driven temperature changes, suggesting an important role for an indirect dynamical effect through ozone redistribution. The results also suggest that differences between the Stratospheric Sounding Unit (SSU)/Microwave Sounding Unit (MSU) and ERA-40 estimates of the solar cycle signal can be explained by the poor vertical resolution of the SSU/MSU measurements. The adjusted radiative forcing of climate change is also investigated. The forcing due to irradiance changes was 0.14 W m−2, which is only 78% of the value obtained by employing the standard method of simple scaling of the total solar irradiance (TSI) change. The difference arises because much of the change in TSI is at wavelengths where ozone absorbs strongly. The forcing due to the ozone change was only 0.004 W m−2 owing to strong compensation between negative shortwave and positive longwave forcings.

Why the western Pacific subtropical high has extended westward since the late 1970s

The western Pacific subtropical high (WPSH) is closely related to Asian climate. Previous examination of changes in the WPSH found a westward extension since the late 1970s, which has contributed to the inter-decadal transition of East Asian climate. The reason for the westward extension is unknown, however. The present study suggests that this significant change of WPSH is partly due to the atmosphere's response to the observed Indian Ocean-western Pacific (IWP) warming. Coordinated by a European Union's Sixth Framework Programme, Understanding the Dynamics of the Coupled Climate System (DYNAMITE), five AGCMs were forced by identical idealized sea surface temperature patterns representative of the IWP warming and cooling. The results of these numerical experiments suggest that the negative heating in the central and eastern tropical Pacific and increased convective heating in the equatorial Indian Ocean/ Maritime Continent associated with IWP warming are in favor of the westward extension of WPSH. The SST changes in IWP influences the Walker circulation, with a subsequent reduction of convections in the tropical central and eastern Pacific, which then forces an ENSO/Gill-type response that modulates the WPSH. The monsoon diabatic heating mechanism proposed by Rodwell and Hoskins plays a secondary reinforcing role in the westward extension of WPSH. The low-level equatorial flank of WPSH is interpreted as a Kelvin response to monsoon condensational heating, while the intensified poleward flow along the western flank of WPSH is in accord with Sverdrup vorticity balance. The IWP warming has led to an expansion of the South Asian high in the upper troposphere, as seen in the reanalysis.

Climate impacts of recent multidecadal changes in Atlantic Ocean Sea Surface Temperature: A multimodel comparison

During the twentieth century sea surface temperatures in the Atlantic Ocean exhibited prominent multidecadal variations. The source of such variations has yet to be rigorously established—but the question of their impact on climate can be investigated. Here we report on a set of multimodel experiments to examine the impact of patterns of warming in the North Atlantic, and cooling in the South Atlantic, derived from observations, that is characteristic of the positive phase of the Atlantic Multidecadal Oscillation (AMO). The experiments were carried out with six atmospheric General Circulation Models (including two versions of one model), and a major goal was to assess the extent to which key climate impacts are consistent between the different models. The major climate impacts are found over North and South America, with the strongest impacts over land found over the United States and northern parts of South America. These responses appear to be driven by a combination of an off-equatorial Gill response to diabatic heating over the Caribbean due to increased rainfall within the region and a Northward shift in the Inter Tropical Convergence Zone (ITCZ) due to the anomalous cross-equatorial SST gradient. The majority of the models show warmer US land temperatures and reduced Mean Sea Level Pressure during summer (JJA) in response to a warmer North Atlantic and a cooler South Atlantic, in line with observations. However the majority of models show no significant impact on US rainfall during summer. Over northern South America, all models show reduced rainfall in southern hemisphere winter (JJA), whilst in Summer (DJF) there is a generally an increase in rainfall. However, there is a large spread amongst the models in the magnitude of the rainfall anomalies over land. Away from the Americas, there are no consistent significant modelled responses. In particular there are no significant changes in the North Atlantic Oscillation (NAO) over the North Atlantic and Europe in Winter (DJF). Additionally, the observed Sahel drying signal in African rainfall is not seen in the modelled responses. Suggesting that, in contrast to some studies, the Atlantic Multidecadal Oscillation was not the primary driver of recent reductions in Sahel rainfall.

Coupled simulations of the Indian monsoon intraseasonal oscillation using a fine-resolution mixed-layer ocean model

The intraseasonal variability (ISV) of the Indian summer monsoon is dominated by a 30–50 day oscillation between “active” and “break” events of enhanced and reduced rainfall over the subcontinent, respectively. These organized convective events form in the equatorial Indian Ocean and propagate north to India. Atmosphere–ocean coupled processes are thought to play a key role the intensity and propagation of these events. A high-resolution, coupled atmosphere–mixed-layer-oceanmodel is assembled: HadKPP. HadKPP comprises the Hadley Centre Atmospheric Model (HadAM3) and the K Profile Parameterization (KPP) mixed-layer ocean model. Following studies that upper-ocean vertical resolution and sub-diurnal coupling frequencies improve the simulation of ISV in SSTs, KPP is run at 1 m vertical resolution near the surface; the atmosphere and ocean are coupled every three hours. HadKPP accurately simulates the 30–50 day ISV in rainfall and SSTs over India and the Bay of Bengal, respectively, but suffers from low ISV on the equator. This is due to the HadAM3 convection scheme producing limited ISV in surface fluxes. HadKPP demonstrates little of the observed northward propagation of intraseasonal events, producing instead a standing oscillation. The lack of equatorial ISV in convection in HadAM3 constrains the ability of KPP to produce equatorial SST anomalies, which further weakens the ISV of convection. It is concluded that while atmosphere–ocean interactions are undoubtedly essential to an accurate simulation of ISV, they are not a panacea for model deficiencies. In regions where the atmospheric forcing is adequate, such as the Bay of Bengal, KPP produces SST anomalies that are comparable to the Tropical Rainfall Measuring Mission Microwave Imager (TMI) SST analyses in both their magnitude and their timing with respect to rainfall anomalies over India. HadKPP also displays a much-improved phase relationship between rainfall and SSTs over a HadAM3 ensemble forced by observed SSTs, when both are compared to observations. Coupling to mixed-layer models such as KPP has the potential to improve operational predictions of ISV, particularly when the persistence time of SST anomalies is shorter than the forecast lead time.

The role of eddies in driving the tropospheric response to stratospheric heating perturbations

A simplified general circulation model has been used to investigate the chain of causality whereby changes in tropospheric circulation and temperature are produced in response to stratospheric heating perturbations. Spinup ensemble experiments have been performed to examine the evolution of the tropospheric circulation in response to such perturbations. The primary aim of these experiments is to investigate the possible mechanisms whereby a tropospheric response to changing solar activity over the 11-yr solar cycle could be produced in response to heating of the equatorial lower stratosphere. This study therefore focuses on a stratospheric heating perturbation in which the heating is largest in the tropics. For comparison, experiments are also performed in which the stratosphere is heated uniformly at all latitudes and in which it is heated preferentially in the polar region. Thus, the mechanisms discussed have a wider relevance for the impact of stratospheric perturbations on the troposphere. The results demonstrate the importance of changing eddy momentum fluxes in driving the tropospheric response. This is confirmed by the lack of a similar response in a zonally symmetric model with fixed eddy forcing. Furthermore, it is apparent that feedback between the tropospheric eddy fluxes and tropospheric circulation changes is required to produce the full model response. The quasigeostrophic index of refraction is used to diagnose the cause of the changes in eddy behavior. It is demonstrated that the latitudinal extent of stratospheric heating is important in determining the direction of displacement of the tropospheric jet and storm track.