A laboratory-scale experiment and a numerical simulation of unusual spiral plumes in a High-Prandtl-number fluid, links supplementary material

We experimentally and numerically investigated the generation of plumes from a local heat source (LHS) and studied the interaction of these plumes with cellular convective motion (CCM) in a rectangular cavity filled with silicon oil at a Prandtl number (Pr) of approximately two thousand. The LHS is generated using a 0.2-W green laser beam. A roll-type CCM is generated by vertically heating one side of the cavity. The CCM may lead to the formation of an unusual spiral convective plume that resembles a vertical Archimedes spiral. A similar plume is obtained in a direct numerical simulation. We discuss the physical mechanism for the formation of a spiral plume and the application of the results to mantle convection problems. We also estimate the Reynolds (Re) and Rayleigh (Ra) numbers and apply self-similarity theory to convection in the Earth's mantle. Spiral plumes can be used to interpret mantle tomography results over the last decade.

Results of calculated surface air temperature anomaly of COSMOS Pliocene experiment 2 in NetCDF format

In this manuscript we describe the experimental procedure employed at the Alfred Wegener Institute in Germany in the preparation of the simulations for the Pliocene Model Intercomparison Project (PlioMIP). We present a description of the utilized Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) and document the procedures that we applied to transfer the Pliocene Research, Interpretation and Synoptic Mapping (PRISM) Project mid-Pliocene reconstruction into model forcing fields. The model setup and spin-up procedure are described for both the paleo- and preindustrial (PI) time slices of PlioMIP experiments 1 and 2, and general results that depict the performance of our model setup for mid-Pliocene conditions are presented. The mid-Pliocene, as simulated with our COSMOS setup and PRISM boundary conditions, is both warmer and wetter in the global mean than the PI. The globally averaged annual mean surface air temperature in the mid-Pliocene standalone atmosphere (fully coupled atmosphere-ocean) simulation is 17.35 °C (17.82 °C), which implies a warming of 2.23 °C (3.40 °C) relative to the respective PI control simulation.

Results of analysis and model simulation of the bed profiles

Bedforms such as dunes and ripples are ubiquitous in rivers and coastal seas, and commonly described as triangular shapes from which height and length are calculated to estimate hydrodynamic and sediment dynamic parameters. Natural bedforms, however, present a far more complicated morphology; the difference between natural bedform shape and the often assumed triangular shape is usually neglected, and how this may affect the flow is unknown. This study investigates the shapes of natural bedforms and how they influence flow and shear stress, based on four datasets extracted from earlier studies on two rivers (the Rio Paraná in Argentina, and the Lower Rhine in The Netherlands). The most commonly occurring morphological elements are a sinusoidal stoss side made of one segment and a lee side made of two segments, a gently sloping upper lee side and a relatively steep (6 to 21°) slip face. A non-hydrostatic numerical model, set up using Delft3D, served to simulate the flow over fixed bedforms with various morphologies derived from the identified morphological elements. Both shear stress and turbulence increase with increasing slip face angle and are only marginally affected by the dimensions and positions of the upper and lower lee side. The average slip face angle determined from the bed profiles is 14°, over which there is no permanent flow separation. Shear stress and turbulence above natural bedforms are higher than above a flat bed but much lower than over the often assumed 30° lee side angle.

Holocene experiment with coupled atmosphere-ocean-model ECHAM5/MPI-OM

Changes in the Earth's orbit lead to changes in the seasonal and meridional distribution of insolation. We quantify the influence of orbitally induced changes on the seasonal temperature cycle in a transient simulation of the last 6000 years - from the mid-Holocene to today - using a coupled atmosphere-ocean general circulation model (ECHAM5/MPI-OM) including a land surface model (JSBACH). The seasonal temperature cycle responds directly to the insolation changes almost everywhere. In the Northern Hemisphere, its amplitude decreases according to an increase in winter insolation and a decrease in summer insolation. In the Southern Hemisphere, the opposite is true. Over the Arctic Ocean, decreasing summer insolation leads to an increase in sea-ice cover. The insulating effect of sea ice between the ocean and the atmosphere leads to decreasing heat flux and favors more "continental" conditions over the Arctic Ocean in winter, resulting in strongly decreasing temperatures. Consequently, there are two competing effects: the direct response to insolation changes and a sea-ice insulation effect. The sea-ice insulation effect is stronger, and thus an increase in the amplitude of the seasonal temperature cycle over the Arctic Ocean occurs. This increase is strongest over the Barents Shelf and influences the temperature response over northern Europe. We compare our modeled seasonal temperatures over Europe to paleo reconstructions. We find better agreements in winter temperatures than in summer temperatures and better agreements in northern Europe than in southern Europe, since the model does not reproduce the southern European Holocene summer cooling inferred from the paleo reconstructions. The temperature reconstructions for northern Europe support the notion of the influence of the sea-ice insulation effect on the evolution of the seasonal temperature cycle.

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.