ALBERT

Description: We analyse the wind and boundary layer properties of turbulent Rayleigh-Bénard convection in a cylindrical container with aspect ratio one for Prandtl number Pr = 0.786 and Rayleigh numbers (Ra) up to 10 9 by means of highly resolved direct numerical simulations. We identify time periods in which the orientation of the large-scale circulation (LSC) is nearly constant in order to perform a statistical analysis of the LSC. The analysis is then reduced to two dimensions by considering only the plane of the LSC. Within this plane the LSC is treated as a wind with thermal and viscous boundary layers developing close to the horizontal plates. Special focus is on the spatial development of the wind magnitude and the boundary layer thicknesses along the bottom plate. A method for the local analysis of the instantaneous boundary layer thicknesses is introduced which shows a dramatically changing wind magnitude along the wind path. Furthermore a linear increase of the viscous and thermal boundary layer thickness along the wind direction is observed for all Ra considered while their ratio is spatially constant but depends weakly on Ra. A possible explanation is a strong spatial variation of the wind magnitude and fluctuations in the boundary layer region. © 2012 Cambridge University Press.

Description: To approximate the velocity and temperature within the boundary layers in turbulent thermal convection at moderate Rayleigh numbers, we consider the Falkner-Skan ansatz, which is a generalization of the Prandtl-Blasius one to a non-zero-pressure-gradient case. This ansatz takes into account the influence of the angle of attack β of the large-scale circulation of a fluid inside a convection cell against the heated/cooled horizontal plate. With respect to turbulent Rayleigh-Bénard convection, we derive several theoretical estimates, among them the limiting cases of the temperature profiles for all angles β, for infinite and for infinitesimal Prandtl numbers Pr. Dependences on Pr and β of the ratio of the thermal to viscous boundary layers are obtained from the numerical solutions of the boundary layers equations. For particular cases of β, accurate approximations are developed as functions on Pr. The theoretical results are corroborated by our direct numerical simulations for Pr= 0. 786 (air) and Pr= 4. 38 (water). The angle of attack β is estimated based on the information on the locations within the plane of the large-scale circulation where the time-averaged wall shear stress vanishes. For the fluids considered it is found that β 0. 7 π and the theoretical predictions based on the Falkner-Skan approximation for this β leads to better agreement with the DNS results, compared with the Prandtl-Blasius approximation for β = π. © 2013 Cambridge University Press.

Description: Direct numerical simulations (DNS) of turbulent thermal convection in a box-shaped domain with regular surface roughness at the heated bottom and cooled top surfaces are conducted for Prandtl number Pr = 0.786 and Rayleigh numbers Ra between 106 and 108. The surface roughness is introduced by four parallelepiped equidistantly distributed obstacles attached to the bottom plate, and four obstacles located symmetrically at the top plate. By varying Ra and the height and width of the obstacles, we investigate the influence of the regular wall roughness on the turbulent heat transport, measured by the Nusselt number Nu. For fixed Ra, the change in the value of Nu is determined not only by the covering area of the surface, i.e. the obstacle height, but also by the distance between the obstacles. The heat flux enhancement is found to be largest for wide cavities between the obstacles which can be 'washed out' by the flow. This is also manifested in an empirical relation, which is based on the DNS data. We further discuss theoretical limiting cases for very wide and very narrow obstacles and combine them into a simple model for the heat flux enhancement due to the wall roughness, without introducing any free parameters. This model predicts well the general trends and the order of magnitude of the heat flux enhancement obtained in the DNS. In the Nu versus Ra scaling, the obstacles work in two ways: for smaller Ra an increase of the scaling exponent compared to the smooth case is found, which is connected to the heat flux entering the cavities from below. For larger Ra the scaling exponent saturates to the one for smooth plates, which can be understood as a full washing-out of the cavities. The latter is also investigated by considering the strength of the mean secondary flow in the cavities and its relation to the wind (i.e. the large-scale circulation), that develops in the core part of the domain. Generally, an increase in the roughness height leads to stronger flows both in the cavities and in the bulk region, while an increase in the width of the obstacles strengthens only the large-scale circulation of the fluid and weakens the secondary flows. An increase of the Rayleigh number always leads to stronger flows, both in the cavities and in the bulk. © 2014 Cambridge University Press.

Description: The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Description: The large variety of atmospheric circulation systems affecting the eastern Asian climate is reflected by the complex Asian vegetation distribution. Particularly in the transition zones of these circulation systems, vegetation is supposed to be very sensitive to climate change. Since proxy records are scarce, hitherto a mechanistic understanding of the past spatio-temporal climate–vegetation relationship is lacking. To assess the Holocene vegetation change and to obtain an ensemble of potential mid-Holocene biome distributions for eastern Asia, we forced the diagnostic biome model BIOME4 with climate anomalies of different transient Holocene climate simulations performed in coupled atmosphere–ocean(–vegetation) models. The simulated biome changes are compared with pollen-based biome records for different key regions. In all simulations, substantial biome shifts during the last 6000 years are confined to the high northern latitudes and the monsoon–westerly wind transition zone, but the temporal evolution and amplitude of change strongly depend on the climate forcing. Large parts of the southern tundra are replaced by taiga during the mid-Holocene due to a warmer growing season and the boreal treeline in northern Asia is shifted northward by approx. 4° in the ensemble mean, ranging from 1.5 to 6° in the individual simulations, respectively. This simulated treeline shift is in agreement with pollen-based reconstructions from northern Siberia. The desert fraction in the transition zone is reduced by 21 % during the mid-Holocene compared to pre-industrial due to enhanced precipitation. The desert–steppe margin is shifted westward by 5° (1–9° in the individual simulations). The forest biomes are expanded north-westward by 2°, ranging from 0 to 4° in the single simulations. These results corroborate pollen-based reconstructions indicating an extended forest area in north-central China during the mid-Holocene. According to the model, the forest-to-non-forest and steppe-to-desert changes in the climate transition zones are spatially not uniform and not linear since the mid-Holocene.

Description: Climate reconstructions are means to extract the signal from uncertain paleo-observations, so-called proxies. It is essential to evaluate these reconstructions to understand and quantify their uncertainties. Similarly, comparing climate simulations and proxies requires approaches to bridge the temporal and spatial differences between both and to address their specific uncertainties. One way to achieve these two goals is so-called pseudoproxies. These are surrogate proxy records within the virtual reality of a climate simulation. They in turn depend on an understanding of the uncertainties of the real proxies including the noise characteristics disturbing the original environmental signal. Common pseudoproxy approaches so far concentrate on data with high temporal resolution over the last approximately 2000 years. Here we provide a simple but flexible noise model for potentially low-resolution sedimentary climate proxies for temperature on millennial timescales, the code for calculating a set of pseudoproxies from a simulation, and one example of pseudoproxies. The noise model considers the influence of other environmental variables, a dependence on the climate state, a bias due to changing seasonality, modifications of the archive (for example bioturbation), potential sampling variability, and a measurement error. Model, code, and data allow us to develop new ways of comparing simulation data with proxies on long timescales. Code and data are available at https://doi.org/10.17605/OSF.IO/ZBEHX (Bothe et al., 2018).