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.

The role of the observed tropical convection in the generation of frost events in the southern cone of South America

Based on previous results obtained from observations and linear wave theory analysis, the hypothesis that large-scale patterns can generate extreme cold events in southeast South America through the propagation of remotely excited Rossby waves was already suggested. This work will confirm these findings and extend their analysis through a series of numerical experiments using a primitive equation model where waves are excited by a thermal forcing situated in positions chosen according to observed convection anomalies over the equatorial region. The basic state used for these experiments is a composite of austral winters with maximum and minimum frequency of occurrence of generalized frosts that can affect a large area known as the Wet Pampas located in the central and eastern part of Argentina. The results suggest that stationary Rossby waves may be one important mechanism linking anomalous tropical convection with the extreme cold events in the Wet Pampas. The combination of tropical convection and a specific basic state can generate the right environment to guide the Rossby waves trigged by the tropical forcing towards South America. Depending on the phase of the waves entering the South American continent, they can favour the advection of anomalous wind at low levels from the south carrying cold and dry air over the whole southern extreme of the continent, producing a generalized frost in the Wet Pampa region. On the other hand, when a basic state based on the composites of minimum frosts is used, an anomalous anticyclone over the southern part of the continent generates a circulation with a south-southeast wind which brings maritime air and therefore humidity over the Wet Pampas region, creating negative temperature anomalies only over the northeastern part of the region. Under these conditions even if frosts occur they would not be generalized, as observed for the other basic state with maximum frequency of occurrence of generalized frosts.

Robust relations between CCN and the vertical evolution of cloud drop size distribution in deep convective clouds

In-situ measurements in convective clouds (up to the freezing level) over the Amazon basin show that smoke from deforestation fires prevents clouds from precipitating until they acquire a vertical development of at least 4 km, compared to only 1-2 km in clean clouds. The average cloud depth required for the onset of warm rain increased by similar to 350 m for each additional 100 cloud condensation nuclei per cm(3) at a super-saturation of 0.5% (CCN0.5%). In polluted clouds, the diameter of modal liquid water content grows much slower with cloud depth (at least by a factor of similar to 2), due to the large number of droplets that compete for available water and to the suppressed coalescence processes. Contrary to what other studies have suggested, we did not observe this effect to reach saturation at 3000 or more accumulation mode particles per cm(3). The CCN0.5% concentration was found to be a very good predictor for the cloud depth required for the onset of warm precipitation and other microphysical factors, leaving only a secondary role for the updraft velocities in determining the cloud drop size distributions. The effective radius of the cloud droplets (r(e)) was found to be a quite robust parameter for a given environment and cloud depth, showing only a small effect of partial droplet evaporation from the cloud's mixing with its drier environment. This supports one of the basic assumptions of satellite analysis of cloud microphysical processes: the ability to look at different cloud top heights in the same region and regard their r(e) as if they had been measured inside one well developed cloud. The dependence of r(e) on the adiabatic fraction decreased higher in the clouds, especially for cleaner conditions, and disappeared at r(e)>=similar to 10 mu m. We propose that droplet coalescence, which is at its peak when warm rain is formed in the cloud at r(e)=similar to 10 mu m, continues to be significant during the cloud's mixing with the entrained air, cancelling out the decrease in r(e) due to evaporation.

Anomalous compressibility of ferropericlase throughout the iron spin cross-over

The thermoelastic properties of ferropericlase Mg(1-x)Fe(x)O (x = 0.1875) throughout the iron high-to-low spin cross-over have been investigated by first principles at Earth`s lower mantle conditions. This cross-over has important consequences for elasticity such as an anomalous bulk modulus (K(S)) reduction. At room temperature the anomaly is somewhat sharp in pressure but broadens with increasing temperature. Along a typical geotherm it occurs across most of the lower mantle with a more significant K(S) reduction at approximate to 1,400-1,600 km depth. This anomaly might also cause a reduction in the effective activation energy for diffusion creep and lead to a viscosity minimum in the mid-lower mantle, in apparent agreement with results from inversion of data related with mantle convection and postglacial rebound.

The role of mobile belts for the longevity of deep cratonic lithosphere: The crumple zone model

Many Archean cratons are surrounded by Proterozoic mobile belts that have experienced episodes of tectonic re-activation over their lifetimes. This suggests that mobile belt lithosphere may be associated with long lived, inherited weakness. It is proposed that the proximity of this weakness can increase the longevity of deep Archean lithosphere by buffering Archean cratons from mantle derived stresses. The physical plausibility of this idea is explored through numerical simulations of mantle convection that include continents and allow for material rheologies that model the combined brittle and ductile behavior of the lithosphere. Within the simulations, the longevity of deep cratonic lithosphere does increase if it is buffered by mobile belts that can fail at relatively low stress levels.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Host-symbiont interactions in the deep-sea vent mussel Bathymodiolus azoricus : a molecular approach

Tese de Doutoramento, Ciências do Mar, especialidade de Biologia Marinha, 19 de Dezembro de 2015, Universidade dos Açores.

Development of procedures for the simulation of atmospheric flows over complex terrain, using OpenFOAM

The main purpose of this work was the development of procedures for the simulation of atmospheric ows over complex terrain, using OpenFOAM. For this aim, tools and procedures were developed apart from this code for the preprocessing and data extraction, which were thereafter applied in the simulation of a real case. For the generation of the computational domain, a systematic method able to translate the terrain elevation model to a native OpenFOAM format (blockMeshDict) was developed. The outcome was a structured mesh, in which the user has the ability to de ne the number of control volumes and its dimensions. With this procedure, the di culties of case set up and the high computation computational e ort reported in literature associated to the use of snappyHexMesh, the OpenFOAM resource explored until then for the accomplishment of this task, were considered to be overwhelmed. Developed procedures for the generation of boundary conditions allowed for the automatic creation of idealized inlet vertical pro les, de nition of wall functions boundary conditions and the calculation of internal eld rst guesses for the iterative solution process, having as input experimental data supplied by the user. The applicability of the generated boundary conditions was limited to the simulation of turbulent, steady-state, incompressible and neutrally strati ed atmospheric ows, always recurring to RaNS (Reynolds-averaged Navier-Stokes) models. For the modelling of terrain roughness, the developed procedure allowed to the user the de nition of idealized conditions, like an uniform aerodynamic roughness length or making its value variable as a function of topography characteristic values, or the using of real site data, and it was complemented by the development of techniques for the visual inspection of generated roughness maps. The absence and the non inclusion of a forest canopy model limited the applicability of this procedure to low aerodynamic roughness lengths. The developed tools and procedures were then applied in the simulation of a neutrally strati ed atmospheric ow over the Askervein hill. In the performed simulations was evaluated the solution sensibility to di erent convection schemes, mesh dimensions, ground roughness and formulations of the k - ε and k - ω models. When compared to experimental data, calculated values showed a good agreement of speed-up in hill top and lee side, with a relative error of less than 10% at a height of 10 m above ground level. Turbulent kinetic energy was considered to be well simulated in the hill windward and hill top, and grossly predicted in the lee side, where a zone of ow separation was also identi ed. Despite the need of more work to evaluate the importance of the downstream recirculation zone in the quality of gathered results, the agreement between the calculated and experimental values and the OpenFOAM sensibility to the tested parameters were considered to be generally in line with the simulations presented in the reviewed bibliographic sources.

Effects of convection on the removal of the multiplicity of stable states in one-dimensional vertical models

Here I develop a model of a radiative-convective atmosphere with both radiative and convective schemes highly simplified. The atmospheric absorption of radiation at selective wavelengths makes use of constant mass absorption coefficients in finite width spectral bands. The convective regime is introduced by using a prescribed lapse rate in the troposphere. The main novelty of the radiative-convective model developed here is that it is solved without using any angular approximation for the radiation field. The solution obtained in the purely radiation mode (i. e. with convection ignored) leads to multiple equilibria of stable states, being very similar to some results recently found in simple models of planetary atmospheres. However, the introduction of convective processes removes the multiple equilibria of stable states. This shows the importance of taking convective processes into account even for qualitative analyses of planetary atmosphere