(Table 1) Nd isotope ratios of ODP Leg 208 samples

We present new Nd isotope records from Walvis Ridge Ocean Drilling Program (ODP) sites 1262-1264 (southeastern Atlantic) spanning the past 24 Ma to investigate the Neogene evolution of Atlantic thermohaline circulation. The new data indicate that deepwater epsilon-Nd(t) values from ODP Site 1262 decrease from -11.0 at 10.6 Ma to -12.5 by 7.3 Ma. This decrease parallels the Nd isotope trends contained in Fe-Mn crust records from the northwestern Atlantic; however, the shift at ODP Site 1262 (4755 m water depth) occurred up to ~6 Ma earlier than the corresponding changes in crust records from the Atlantic and earlier than any Atlantic site shallower than 2700 m paleowater depth. Recent interpretations of the rapid decrease in Fe-Mn crust Nd epsilon-Nd(t)ss values invoke changes in weathering inputs to the Labrador Sea region rather than a fundamental change in deepwater convection in the Labrador Sea. However, the new evidence for significant depth stratification of the Nd isotope signal in the southeastern Atlantic between 10.6 and 7.3 Ma suggests that the onset of deepwater convection in the Labrador Sea may have played a role in the deepwater decrease in Nd isotopic composition. Climatic conditions during the middle to late Miocene likely favored an increase in the importance of glacially induced mechanical weathering, while at the same time promoting deep convection in the Labrador Sea.

Convective Cold Pools over the Atlas Mountains and Their Influence on the Saharan Heat Low

The West African Monsoon (WAM) and its representation in numerical models are strongly influenced by the Saharan Heat Low (SHL), a low-pressure system driven by radiative heating over the central Sahara and ventilated by the cold and moist inflow from adjacent oceans. It has recently been shown that a significant part of the southerly moisture flux into the SHL originates from convective cold pools over the Sahel. These density currents driven by evaporation of rain are largely absent in models with parameterized convection. This crucial issue has been hypothesized to contribute to the inability of many climate models to reproduce the variability of the WAM. Here, the role of convective cold pools approaching the SHL from the Atlas Mountains, which are a strong orographic trigger for deep convection in Northwest Africa, is analyzed. Knowledge about the frequency of these events, as well as their impact on large-scale dynamics, is required to understand their contribution to the variability of the SHL and to known model uncertainties. The first aspect is addressed through the development of an objective and automated method for the generation of multi-year climatologies not available before. The algorithm combines freely available standard surface observations with satellite microwave data. Representativeness of stations and influence of their spatial density are addressed by comparison to a satellite-only climatology. Applying this algorithm to data from automated weather stations and manned synoptic stations in and south of the Atlas Mountains reveals the frequent occurrence. On the order of 6 events per month are detected from May to September when the SHL is in its northernmost position. The events tend to cluster into several-days long convectively active periods, often with strong events on consecutive days. This study is the first to diagnose dynamical impacts of such periods on the SHL, based on simulations of two example cases using the Weather Research and Forecast (WRF) model at convection-permitting resolution. Sensitivity experiments with artificially removed cold pools as well as different resolutions and parameterizations are conducted. Results indicate increases in surface pressure of more than 1 hPa and significant moisture transports into the desert over several days. This moisture affects radiative heating and thus the energy balance of the SHL. Even though cold pool events north of the SHL are less frequent when compared to their Sahelian counterparts, it is shown that they gain importance due to their temporal clustering on synoptic timescale. Together with studies focusing on the Sahel, this work emphasizes the need for improved parameterization schemes for deep convection in order to produce more reliable climate projections for the WAM.

Vertical thermal structure of the atmosphere during active and weak phases of convection over the north Bay of Bengal: Observation and model results

Owing to the lack of atmospheric vertical profile data with sufficient accuracy and vertical resolution, the response of the deep atmosphere to passage of monsoon systems over the Bay of Bengal. had not been satisfactorily elucidated. Under the Indian Climate Research Programme, a special observational programme called 'Bay of Bengal Monsoon Experiment' (BOBMEX), was conducted during July-August 1999. The present study is based on the high-resolution radiosondes launched during BOBMEX in the north Bay. Clear changes in the vertical thermal structure of the atmosphere between active and weak phases of convection have been observed. The atmosphere cooled below 6 km height and became warmer between 6 and 13 km height. The warmest layer was located between 8 and 10 km height, and the coldest layer was found just below 5 km height. The largest fluctuations in the humidity field occurred in the mid-troposphere. The observed changes between active and weak phases of convection are compared with the results from an atmospheric general circulation model, which is similar to that used at the National Centre for Medium Range Weather Forecasting, New Delhi. The model is not able to capture realistically some important features of the temperature and humidity profiles in the lower troposphere and in the boundary layer during the active and weak spells.

The diurnal cycle of convection and atmospheric tides in an aquaplanet

The diurnal cycle of tropical convection and its relationship to the atmospheric tides is investigated using an aquaplanet GCM. The diurnal and semidiurnal harmonics of precipitation are both found to contribute significantly to the total diurnal variability of precipitation in the model, which is broadly consistent with observations of the diurnal cycle of convection over the open ocean. The semidiurnal tide is found to be the dominant forcing for the semidiurnal harmonic of precipitation. In contrast the diurnal tide plays only a small role in forcing the diurnal harmonic of precipitation, which is dominated by the variations in shortwave and longwave heating. In both the diurnal and semidiurnal harmonics, the feedback onto the convection by the humidity tendencies due to the convection is found to be important in determining the phase of the harmonics. Further experiments show that the diurnal cycle of precipitation is sensitive to the choice of closure in the convection scheme. While the surface pressure signal of the simulated atmospheric tides in the model agree well with both theory and observations in their magnitude and phase, sensitivity experiments suggest that the role of the stratospheric ozone in forcing the semidiurnal tide is much reduced compared to theoretical predictions. Furthermore, the influence of the cloud radiative effects seems small. It is suggested that the radiative heating profile in the troposphere, associated primarily with the water vapor distribution, is more important than previously thought for driving the semidiurnal tide. However, this result may be sensitive to the vertical resolution and extent of the model.

Atmospheric factors governing banded orographic convection

Atmospheric factors Governing Banded Orographic Convection The three-dimensional structure of shallow orographic convection is investigated through simulations performed with a cloud-resolving numerical model. In moist flows that overcome a given topographic barrier to form statically unstable cap clouds, the organization of the convection depends on both the atmospheric structure and the mechanism by which the convection is initiated. Convection initiated by background thermal fluctuations embedded in the flow over a smooth mountain (without any small-scale topographic features) tends to be cellular and disorganized except that shear-parallel bands may form in flows with strong unidirectional vertical shear. The development of well-organized bands is favored when there is weak static instability inside the cloud and when the dry air surrounding the cloud is strongly stable. These bands move with the flow and distribute their cumulative precipitation evenly over the mountain upslope. Similar shear-parallel bands also develop in flows where convection is initiated by small-scale topographic noise superimposed onto the main mountain profile, but in this case stronger circulations are also triggered that create stationary rainbands parallel to the low-level flow. This second dominant mode, which is less sensitive to the atmospheric structure and the strength of forcing, is triggered by lee waves that form over small-scale topographic bumps near the upstream edge of the main orographic cloud. Due to their stationarity, these flow-parallel bands can produce locally heavy precipitation amounts.

Understanding and representing atmospheric convection across scales: recommendations from the meeting held at Dartington Hall, Devon, UK, 28–30 January 2013

As weather and climate models move toward higher resolution, there is growing excitement about potential future improvements in the understanding and prediction of atmospheric convection and its interaction with larger-scale phenomena. A meeting in January 2013 in Dartington, Devon was convened to address the best way to maximise these improvements, specifically in a UK context but with international relevance. Specific recommendations included increased convective-scale observations, high-resolution virtual laboratories, and a system of parameterization test beds with a range of complexities. The main recommendation was to facilitate the development of physically based convective parameterizations that are scale-aware, non-local, non-equilibrium, and stochastic.

Basic theory behind parameterizing atmospheric convection

Last fall, a network of the European Cooperation in Science and Technology (COST), called “Basic Concepts for Convection Parameterization in Weather Forecast and Climate Models” (COST Action ES0905; see http://w3.cost.esf.org/index.php?id=205&action_number=ES0905), organized a 10-day training course on atmospheric convection and its parameterization. The aim of the workshop, held on the island of Brac, Croatia, was to help young scientists develop an in-depth understanding of the core theory underpinning convection parameterizations. The speakers also sought to impart an appreciation of the various approximations, compromises, and ansatz necessary to translate theory into operational practice for numerical models.

Including the sub-grid scale plume rise of vegetation fires in low resolution atmospheric transport models

We describe and begin to evaluate a parameterization to include the vertical transport of hot gases and particles emitted from biomass burning in low resolution atmospheric-chemistry transport models. This sub-grid transport mechanism is simulated by embedding a 1-D cloud-resolving model with appropriate lower boundary conditions in each column of the 3-D host model. Through assimilation of remote sensing fire products, we recognize which columns have fires. Using a land use dataset appropriate fire properties are selected. The host model provides the environmental conditions, allowing the plume rise to be simulated explicitly. The derived height of the plume is then used in the source emission field of the host model to determine the effective injection height, releasing the material emitted during the flaming phase at this height. Model results are compared with CO aircraft profiles from an Amazon basin field campaign and with satellite data, showing the huge impact that this mechanism has on model performance. We also show the relative role of each main vertical transport mechanisms, shallow and deep moist convection and the pyro-convection (dry or moist) induced by vegetation fires, on the distribution of biomass burning CO emissions in the troposphere.

Covariation of deep Southern Ocean oxygenation and atmospheric CO2 through the last ice age

No single mechanism can account for the full amplitude of past atmospheric carbon dioxide (CO2) concentration variability over glacial–interglacial cycles. A build-up of carbon in the deep ocean has been shown to have occurred during the Last Glacial Maximum. However, the mechanisms responsible for the release of the deeply sequestered carbon to the atmosphere at deglaciation, and the relative importance of deep ocean sequestration in regulating millennial-timescale variations in atmospheric CO2 concentration before the Last Glacial Maximum, have remained unclear. Here we present sedimentary redox-sensitive trace metal records from the Antarctic Zone of the Southern Ocean that provide a reconstruction of transient changes in deep ocean oxygenation and, by inference, respired carbon storage throughout the last glacial cycle. Our data suggest that respired carbon was removed from the abyssal Southern Ocean during the Northern Hemisphere cold phases of the deglaciation, when atmospheric CO2 concentration increased rapidly, reflecting—at least in part— a combination of dwindling iron fertilization by dust and enhanced deep ocean ventilation. Furthermore, our records show that the observed covariation between atmospheric CO2 concentration and abyssal Southern Ocean oxygenation was maintained throughout most of the past 80,000 years. This suggests that on millennial timescales deep ocean circulation and iron fertilization in the Southern Ocean played a consistent role in modifying atmospheric CO2 concentration.

Towards convection-resolving, global atmospheric simulations with the Model for Prediction Across Scales (MPAS) v3.1: an extreme scaling experiment, links to sourcecode of model and grides

The Model for Prediction Across Scales (MPAS) is a novel set of Earth system simulation components and consists of an atmospheric model, an ocean model and a land-ice model. Its distinct features are the use of unstructured Voronoi meshes and C-grid discretisation to address shortcomings of global models on regular grids and the use of limited area models nested in a forcing data set, with respect to parallel scalability, numerical accuracy and physical consistency. This concept allows one to include the feedback of regional land use information on weather and climate at local and global scales in a consistent way, which is impossible to achieve with traditional limited area modelling approaches. Here, we present an in-depth evaluation of MPAS with regards to technical aspects of performing model runs and scalability for three medium-size meshes on four different high-performance computing (HPC) sites with different architectures and compilers. We uncover model limitations and identify new aspects for the model optimisation that are introduced by the use of unstructured Voronoi meshes. We further demonstrate the model performance of MPAS in terms of its capability to reproduce the dynamics of the West African monsoon (WAM) and its associated precipitation in a pilot study. Constrained by available computational resources, we compare 11-month runs for two meshes with observations and a reference simulation from the Weather Research and Forecasting (WRF) model. We show that MPAS can reproduce the atmospheric dynamics on global and local scales in this experiment, but identify a precipitation excess for the West African region. Finally, we conduct extreme scaling tests on a global 3?km mesh with more than 65 million horizontal grid cells on up to half a million cores. We discuss necessary modifications of the model code to improve its parallel performance in general and specific to the HPC environment. We confirm good scaling (70?% parallel efficiency or better) of the MPAS model and provide numbers on the computational requirements for experiments with the 3?km mesh. In doing so, we show that global, convection-resolving atmospheric simulations with MPAS are within reach of current and next generations of high-end computing facilities.