Modelling monsoons: understanding and predicting current and future behaviour

The global monsoon system is so varied and complex that understanding and predicting its diverse behaviour remains a challenge that will occupy modellers for many years to come. Despite the difficult task ahead, an improved monsoon modelling capability has been realized through the inclusion of more detailed physics of the climate system and higher resolution in our numerical models. Perhaps the most crucial improvement to date has been the development of coupled ocean-atmosphere models. From subseasonal to interdecadal time scales, only through the inclusion of air-sea interaction can the proper phasing and teleconnections of convection be attained with respect to sea surface temperature variations. Even then, the response to slow variations in remote forcings (e.g., El Niño—Southern Oscillation) does not result in a robust solution, as there are a host of competing modes of variability that must be represented, including those that appear to be chaotic. Understanding the links between monsoons and land surface processes is not as mature as that explored regarding air-sea interactions. A land surface forcing signal appears to dominate the onset of wet season rainfall over the North American monsoon region, though the relative role of ocean versus land forcing remains a topic of investigation in all the monsoon systems. Also, improved forecasts have been made during periods in which additional sounding observations are available for data assimilation. Thus, there is untapped predictability that can only be attained through the development of a more comprehensive observing system for all monsoon regions. Additionally, improved parameterizations - for example, of convection, cloud, radiation, and boundary layer schemes as well as land surface processes - are essential to realize the full potential of monsoon predictability. A more comprehensive assessment is needed of the impact of black carbon aerosols, which may modulate that of other anthropogenic greenhouse gases. Dynamical considerations require ever increased horizontal resolution (probably to 0.5 degree or higher) in order to resolve many monsoon features including, but not limited to, the Mei-Yu/Baiu sudden onset and withdrawal, low-level jet orientation and variability, and orographic forced rainfall. Under anthropogenic climate change many competing factors complicate making robust projections of monsoon changes. Absent aerosol effects, increased land-sea temperature contrast suggests strengthened monsoon circulation due to climate change. However, increased aerosol emissions will reflect more solar radiation back to space, which may temper or even reduce the strength of monsoon circulations compared to the present day. Precipitation may behave independently from the circulation under warming conditions in which an increased atmospheric moisture loading, based purely on thermodynamic considerations, could result in increased monsoon rainfall under climate change. The challenge to improve model parameterizations and include more complex processes and feedbacks pushes computing resources to their limit, thus requiring continuous upgrades of computational infrastructure to ensure progress in understanding and predicting current and future behaviour of monsoons.

The statistical properties of tropical ice clouds generated by the West-African and Australian monsoons from ground-based radar-lidar observations

In the present paper we characterize the statistical properties of non-precipitating tropical ice clouds (deep ice anvils resulting from deep convection and cirrus clouds) over Niamey, Niger, West Africa, and Darwin, northern Australia, using ground-based radar–lidar observations from the Atmospheric Radiation Measurement (ARM) programme. The ice cloud properties analysed in this paper are the frequency of ice cloud occurrence, cloud fraction, the morphological properties (cloud-top height, base height, and thickness), the microphysical and radiative properties (ice water content, visible extinction, effective radius, terminal fall speed, and concentration), and the internal cloud dynamics (in-cloud vertical air velocity). The main highlight of the paper is that it characterizes for the first time the probability density functions of the tropical ice cloud properties, their vertical variability and their diurnal variability at the same time. This is particularly important over West Africa, since the ARM deployment in Niamey provides the first vertically resolved observations of non-precipitating ice clouds in this crucial area in terms of redistribution of water and energy in the troposphere. The comparison between the two sites also provides an additional observational basis for the evaluation of the parametrization of clouds in large-scale models, which should be able to reproduce both the statistical properties at each site and the differences between the two sites. The frequency of ice cloud occurrence is found to be much larger over Darwin when compared to Niamey, and with a much larger diurnal variability, which is well correlated with the diurnal cycle of deep convective activity. The diurnal cycle of the ice cloud occurrence over Niamey is also much less correlated with that of deep convective activity than over Darwin, probably owing to the fact that Niamey is further away from the deep convective sources of the region. The frequency distributions of cloud fraction are strongly bimodal and broadly similar over the two sites, with a predominance of clouds characterized either by a very small cloud fraction (less than 0.3) or a very large cloud fraction (larger than 0.9). The ice clouds over Darwin are also much thicker (by 1 km or more statistically) and are characterized by a much larger diurnal variability than ice clouds over Niamey. Ice clouds over Niamey are also characterized by smaller particle sizes and fall speeds but in much larger concentrations, thereby carrying more ice water and producing more visible extinction than the ice clouds over Darwin. It is also found that there is a much larger occurrence of downward in-cloud air motions less than 1 m s−1 over Darwin, which together with the larger fall speeds retrieved over Darwin indicates that the life cycle of ice clouds is probably shorter over Darwin than over Niamey.

Mid-Holocene monsoons: a multi-model analysis of the inter-hemispheric differences in the responses to orbital forcing and ocean feedbacks

The response of monsoon circulation in the northern and southern hemisphere to 6 ka orbital forcing has been examined in 17 atmospheric general circulation models and 11 coupled ocean–atmosphere general circulation models. The atmospheric response to increased summer insolation at 6 ka in the northern subtropics strengthens the northern-hemisphere summer monsoons and leads to increased monsoonal precipitation in western North America, northern Africa and China; ocean feedbacks amplify this response and lead to further increase in monsoon precipitation in these three regions. The atmospheric response to reduced summer insolation at 6 ka in the southern subtropics weakens the southern-hemisphere summer monsoons and leads to decreased monsoonal precipitation in northern South America, southern Africa and northern Australia; ocean feedbacks weaken this response so that the decrease in rainfall is smaller than might otherwise be expected. The role of the ocean in monsoonal circulation in other regions is more complex. There is no discernable impact of orbital forcing in the monsoon region of North America in the atmosphere-only simulations but a strong increase in precipitation in the ocean–atmosphere simulations. In contrast, there is a strong atmospheric response to orbital forcing over northern India but ocean feedback reduces the strength of the change in the monsoon although it still remains stronger than today. Although there are differences in magnitude and exact location of regional precipitation changes from model to model, the same basic mechanisms are involved in the oceanic modulation of the response to orbital forcing and this gives rise to a robust ensemble response for each of the monsoon systems. Comparison of simulated and reconstructed changes in regional climate suggest that the coupled ocean–atmosphere simulations produce more realistic changes in the northern-hemisphere monsoons than atmosphere-only simulations, though they underestimate the observed changes in precipitation in all regions. Evaluation of the southern-hemisphere monsoons is limited by lack of quantitative reconstructions, but suggest that model skill in simulating these monsoons is limited.

Global monsoons in the mid-Holocene and oceanic feedback

The response of the six major summer monsoon systems (the North American monsoon, the northern Africa monsoon, the Asia monsoon, the northern Australasian monsoon, the South America monsoon and the southern Africa monsoon) to mid-Holocene orbital forcing has been investigated using a coupled ocean–atmosphere general circulation model (FOAM), with the focus on the distinct roles of the direct insolation forcing and oceanic feedback. The simulation result is also found to compare well with the NCAR CSM. The direct effects of the change in insolation produce an enhancement of the Northern Hemisphere monsoons and a reduction of the Southern Hemisphere monsoons. Ocean feedbacks produce a further enhancement of the northern Africa monsoon and the North American monsoon. However, ocean feedbacks appear to weaken the Asia monsoon, although the overall effect (direct insolation forcing plus ocean feedback) remains a strengthened monsoon. The impact of ocean feedbacks on the South American and southern African monsoons is relatively small, and therefore these regions, especially the South America, experienced a reduced monsoon regime compared to present. However, there is a strong ocean feedback on the northern Australian monsoon that negates the direct effects of orbital changes and results in a strengthening of austral summer monsoon precipitation in this region. A new synthesis is made for mid-Holocene paleoenvironmental records and is compared with the model simulations. Overall, model simulations produce changes in regional climates that are generally consistent with paleoenvironmental observations.

The mean evolution and variability of the Asian summer monsoon: Comparison of ECMWF and NCEP-NCAR reanalyses

The behavior of the Asian summer monsoon is documented and compared using the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA) and the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) Reanalysis. In terms of seasonal mean climatologies the results suggest that, in several respects, the ERA is superior to the NCEP-NCAR Reanalysis. The overall better simulation of the precipitation and hence the diabatic heating field over the monsoon domain in ERA means that the analyzed circulation is probably nearer reality. In terms of interannual variability, inconsistencies in the definition of weak and strong monsoon years based on typical monsoon indices such as All-India Rainfall (AIR) anomalies and the large-scale wind shear based dynamical monsoon index (DMI) still exist. Two dominant modes of interannual variability have been identified that together explain nearly 50% of the variance. Individually, they have many features in common with the composite flow patterns associated with weak and strong monsoons, when defined in terms of regional AIR anomalies and the large-scale DMI. The reanalyses also show a common dominant mode of intraseasonal variability that describes the latitudinal displacement of the tropical convergence zone from its oceanic-to-continental regime and essentially captures the low-frequency active/break cycles of the monsoon. The relationship between interannual and intraseasonal variability has been investigated by considering the probability density function (PDF) of the principal component of the dominant intraseasonal mode. Based on the DMI, there is an indication that in years with a weaker monsoon circulation, the PDF is skewed toward negative values (i,e., break conditions). Similarly, the PDFs for El Nino and La Nina years suggest that El Nino predisposes the system to more break spells, although the sample size may limit the statistical significance of the results.

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.

MJO Simulation Diagnostics

The Madden–Julian oscillation (MJO) interacts with and influences a wide range of weather and climate phenomena (e.g., monsoons, ENSO, tropical storms, midlatitude weather), and represents an important, and as yet unexploited, source of predictability at the subseasonal time scale. Despite the important role of the MJO in climate and weather systems, current global circulation models (GCMs) exhibit considerable shortcomings in representing this phenomenon. These shortcomings have been documented in a number of multimodel comparison studies over the last decade. However, diagnosis of model performance has been challenging, and model progress has been difficult to track, because of the lack of a coherent and standardized set of MJO diagnostics. One of the chief objectives of the U.S. Climate Variability and Predictability (CLIVAR) MJO Working Group is the development of observation-based diagnostics for objectively evaluating global model simulations of the MJO in a consistent framework. Motivation for this activity is reviewed, and the intent and justification for a set of diagnostics is provided, along with specification for their calculation, and illustrations of their application. The diagnostics range from relatively simple analyses of variance and correlation to more sophisticated space–time spectral and empirical orthogonal function analyses. These diagnostic techniques are used to detect MJO signals, to construct composite life cycles, to identify associations of MJO activity with the mean state, and to describe interannual variability of the MJO.

Using idealized snow forcing to test teleconnections with the Indian summer monsoon in the Hadley Centre GCM

Anomalous heavy snow during winter or spring has long been regarded as a possible precursor of deficient Indian monsoon rainfall during the subsequent summer. However previous work in this field is inconclusive, in terms of the mechanism that communicates snow anomalies to the monsoon summer, and even the region from which snow has the most impact. In this study we explore these issues in coupled and atmosphere-only versions of the Hadley Centre model. A 1050-year control integration of the HadCM3 coupled model, which well represents the seasonal cycle of snow cover over the Eurasian continent, is analysed and shows evidence for weakened monsoons being preceded by strong snow forcing (in the absence of ENSO) over either the Himalaya/Tibetan Plateau or north/west Eurasia regions. However, empirical orthogonal function (EOF) analysis of springtime interannual variability in snow depth shows the leading mode to have opposite signs between these two regions, suggesting that competing mechanisms may be possible. To determine the dominant region, ensemble integrations are carried out using HadAM3, the atmospheric component of HadCM3, and a variety of anomalous snow forcing initial conditions obtained from the control integration of the coupled model. Forcings are applied during spring in separate experiments over the Himalaya/Tibetan Plateau and north/west Eurasia regions, in conjunction with climatological SSTs in order to avoid the direct effects of ENSO. With the aid of idealized forcing conditions in sensitivity tests, we demonstrate that forcing from the Himalaya region is dominant in this model via a Blanford-type mechanism involving reduced surface sensible heat and longwave fluxes, reduced heating of the troposphere over the Tibetan Plateau and consequently a reduced meridional tropospheric temperature gradient which weakens the monsoon during early summer. Snow albedo is shown to be key to the mechanism, explaining around 50% of the perturbation in sensible heating over the Tibetan Plateau, and accounting for the majority of cooling through the troposphere.

Solar influences on climate

Understanding the influence of solar variability on the Earth’s climate requires knowledge of solar variability, solar-terrestrial interactions and the mechanisms determining the response of the Earth’s climate system. We provide a summary of our current understanding in each of these three areas. Observations and mechanisms for the Sun's variability are described, including solar irradiance variations on both decadal and centennial timescales and their relation to galactic cosmic rays. Corresponding observations of variations of the Earth’s climate on associated timescales are described, including variations in ozone, temperatures, winds, clouds, precipitation and regional modes of variability such as the monsoons and the North Atlantic Oscillation. A discussion of the available solar and climate proxies is provided. Mechanisms proposed to explain these climate observations are described, including the effects of variations in solar irradiance and of charged particles. Finally, the contribution of solar variations to recent observations of global climate change are discussed.

Dependence of Indian monsoon rainfall on moisture fluxes across the Arabian Sea and the impact of coupled model sea surface temperature biases

The Arabian Sea is an important moisture source for Indian monsoon rainfall. The skill of climate models in simulating the monsoon and its variability varies widely, while Arabian Sea cold sea surface temperature (SST) biases are common in coupled models and may therefore influence the monsoon and its sensitivity to climate change. We examine the relationship between monsoon rainfall, moisture fluxes and Arabian Sea SST in observations and climate model simulations. Observational analysis shows strong monsoons depend on moisture fluxes across the Arabian Sea, however detecting consistent signals with contemporaneous summer SST anomalies is complicated in the observed system by air/sea coupling and large-scale induced variability such as the El Niño-Southern Oscillation feeding back onto the monsoon through development of the Somali Jet. Comparison of HadGEM3 coupled and atmosphere-only configurations suggests coupled model cold SST biases significantly reduce monsoon rainfall. Idealised atmosphere-only experiments show that the weakened monsoon can be mainly attributed to systematic Arabian Sea cold SST biases during summer and their impact on the monsoon-moisture relationship. The impact of large cold SST biases on atmospheric moisture content over the Arabian Sea, and also the subsequent reduced latent heat release over India, dominates over any enhancement in the land-sea temperature gradient and results in changes to the mean state. We hypothesize that a cold base state will result in underestimation of the impact of larger projected Arabian Sea SST changes in future climate, suggesting that Arabian Sea biases should be a clear target for model development.

Aspects of the large scale tropical atmospheric circulation

An account is given of a number of recent studies with idealised models whose aim is to further understanding of the large-scale tropical atmospheric circulation. Initial-value integrations with a model with imposed heating are used to discuss aspects of the Asian summer monsoon, including constraints on cross-equatorial flow into the monsoon. The summer descent in the Mediterranean region and on the eastern sides of the summer subtropical anticyclones are seen to be associated with the monsoons to their east. An aqua-planet GCM is used to investigate the relationship between simple SST distributions and tropical convection and circulation. The existence of strong equatorial convection and Hadley cells is found to depend sensitively on the curvature of the meridional profile in SST. Zonally confined SST maxima produce convective maxima centred to the west and suppression of convection elsewhere. Strong equatorial zonal flow changes are found in some experiments and three mechanisms for producing these are investigated in a model with imposed heating. 1.

The effect of Arabian Sea optical properties on SST biases and the South Asian summer monsoon in a coupled GCM

This study examines the effect of seasonally varying chlorophyll on the climate of the Arabian Sea and South Asian monsoon. The effect of such seasonality on the radiative properties of the upper ocean is often a missing process in coupled general circulation models and its large amplitude in the region makes it a pertinent choice for study to determine any impact on systematic biases in the mean and seasonality of the Arabian Sea. In this study we examine the effects of incorporating a seasonal cycle in chlorophyll due to phytoplankton blooms in the UK Met Office coupled atmosphere-ocean GCM HadCM3. This is achieved by performing experiments in which the optical properties of water in the Arabian Sea - a key signal of the semi-annual cycle of phytoplankton blooms in the region - are calculated from a chlorophyll climatology derived from Sea-viewing Wide Field-of-View Sensor (SeaWiFS) data. The SeaWiFS chlorophyll is prescribed in annual mean and seasonally-varying experiments. In response to the chlorophyll bloom in late spring, biases in mixed layer depth are reduced by up to 50% and the surface is warmed, leading to increases in monsoon rainfall during the onset period. However when the monsoons are fully established in boreal winter and summer and there are strong surface winds and a deep mixed layer, biases in the mixed layer depth are reduced but the surface undergoes cooling. The seasonality of the response of SST to chlorophyll is found to depend on the relative depth of the mixed layer to that of the anomalous penetration depth of solar fluxes. Thus the inclusion of the effects of chlorophyll on radiative properties of the upper ocean acts to reduce biases in mixed layer depth and increase seasonality in SST.

Predictability of the western North Pacific summer climate demonstrated by the coupled models of ENSEMBLES

The Asian monsoon system, including the western North Pacific (WNP), East Asian, and Indian monsoons, dominates the climate of the Asia-Indian Ocean-Pacific region, and plays a significant role in the global hydrological and energy cycles. The prediction of monsoons and associated climate features is a major challenge in seasonal time scale climate forecast. In this study, a comprehensive assessment of the interannual predictability of the WNP summer climate has been performed using the 1-month lead retrospective forecasts (hindcasts) of five state-of-the-art coupled models from ENSEMBLES for the period of 1960–2005. Spatial distribution of the temporal correlation coefficients shows that the interannual variation of precipitation is well predicted around the Maritime Continent and east of the Philippines. The high skills for the lower-tropospheric circulation and sea surface temperature (SST) spread over almost the whole WNP. These results indicate that the models in general successfully predict the interannual variation of the WNP summer climate. Two typical indices, the WNP summer precipitation index and the WNP lower-tropospheric circulation index (WNPMI), have been used to quantify the forecast skill. The correlation coefficient between five models’ multi-model ensemble (MME) mean prediction and observations for the WNP summer precipitation index reaches 0.66 during 1979–2005 while it is 0.68 for the WNPMI during 1960–2005. The WNPMI-regressed anomalies of lower-tropospheric winds, SSTs and precipitation are similar between observations and MME. Further analysis suggests that prediction reliability of the WNP summer climate mainly arises from the atmosphere–ocean interaction over the tropical Indian and the tropical Pacific Ocean, implying that continuing improvement in the representation of the air–sea interaction over these regions in CGCMs is a key for long-lead seasonal forecast over the WNP and East Asia. On the other hand, the prediction of the WNP summer climate anomalies exhibits a remarkable spread resulted from uncertainty in initial conditions. The summer anomalies related to the prediction spread, including the lower-tropospheric circulation, SST and precipitation anomalies, show a Pacific-Japan or East Asia-Pacific pattern in the meridional direction over the WNP. Our further investigations suggest that the WNPMI prediction spread arises mainly from the internal dynamics in air–sea interaction over the WNP and Indian Ocean, since the local relationships among the anomalous SST, circulation, and precipitation associated with the spread are similar to those associated with the interannual variation of the WNPMI in both observations and MME. However, the magnitudes of these anomalies related to the spread are weaker, ranging from one third to a half of those anomalies associated with the interannual variation of the WNPMI in MME over the tropical Indian Ocean and subtropical WNP. These results further support that the improvement in the representation of the air–sea interaction over the tropical Indian Ocean and subtropical WNP in CGCMs is a key for reducing the prediction spread and for improving the long-lead seasonal forecast over the WNP and East Asia.

Isomap nonlinear dimensionality reduction and bimodality of Asian monsoon convection

It is known that the empirical orthogonal function method is unable to detect possible nonlinear structure in climate data. Here, isometric feature mapping (Isomap), as a tool for nonlinear dimensionality reduction, is applied to 1958–2001 ERA-40 sea-level pressure anomalies to study nonlinearity of the Asian summer monsoon intraseasonal variability. Using the leading two Isomap time series, the probability density function is shown to be bimodal. A two-dimensional bivariate Gaussian mixture model is then applied to identify the monsoon phases, the obtained regimes representing enhanced and suppressed phases, respectively. The relationship with the large-scale seasonal mean monsoon indicates that the frequency of monsoon regime occurrence is significantly perturbed in agreement with conceptual ideas, with preference for enhanced convection on intraseasonal time scales during large-scale strong monsoons. Trend analysis suggests a shift in concentration of monsoon convection, with less emphasis on South Asia and more on the East China Sea.