A length scale defining partially-resolved boundary-layer turbulence simulations

PublishedJournal ArticleThis is the author accepted manuscript. The final version is available from Springer via the DOI in this record.Numerical weather prediction (NWP) model forecasts at horizontal grid lengths in the range of 100 m to 1 km are now possible. Within this range of grid lengths, the convective boundary layer (CBL) is partially resolved and thus in the so-called 'grey zone'. For simulations in the grey zone, numerical dissipation sources from both the advection scheme and the subgrid model are likely to be significant. Until now, these effects have not been incorporated fully into our understanding of the grey zone. In order to quantify these effects, a dissipation length scale is defined based on the second moment of the turbulent kinetic energy (TKE) spectrum. An ensemble of simulations of a CBL are performed using a large-eddy model across the grey-zone resolutions and for a range of subgrid model, advection scheme and vertical grid configurations. The dissipation length scale distinguishes the effects of the different model configurations in the grey zone. In the middle of the boundary layer, the resolved TKE is strongly controlled by the numerical dissipation. This leads to a similarity law for the resolved TKE in the grey zone using the dissipation length scale. A new definition of the grey zone emerges where the inversion depth and dissipation length scale are the same size. This contrasts with the typical definition using the horizontal grid length. At the inversion, however, the variation of the dissipation length scale with grid length is less predictable, reflecting significant challenges for modelling entrainment in the grey zone. The dissipation length scale is thus a simple diagnostic to aid both NWP and large-eddy modellers in understanding the grey zone. © 2013 Springer Science+Business Media Dordrecht

Version

n/

The z < 1.2 optical luminosity function from a sample of ∼410,000 galaxies in Boötes

Using a sample of ~410,000 galaxies to a depth of IAB=24 over 8.26 deg2 in the Boötes field (~10 times larger than the z~1 luminosity function (LF) studies in the prior literature), we have accurately measured the evolving B-band LF of red galaxies at z&lt;1.2 and blue galaxies at z&lt;1.0 In addition to the large sample size, we utilize photometry that accounts for the varying angular sizes of galaxies, photometric redshifts verified with spectroscopy, and absolute magnitudes that should have very small random and systematic errors. Our results are consistent with the migration of galaxies from the blue cloud to the red sequence as they cease to form stars and with downsizing in which more massive and luminous blue galaxies cease star formation earlier than fainter less massive ones. Comparing the observed fading of red galaxies with that expected from passive evolution alone, we find that the stellar mass contained within the red galaxy population has increased by a factor of ~3.6 from z~1.1 to z~0.1 The bright end of the red galaxy LF fades with decreasing redshift, with the rate of fading increasing from ~0.2 mag per unit redshift at z = 1.0 to ~0.8 at z = 0.2. The overall decrease in luminosity implies that the stellar mass in individual highly luminous red galaxies increased by a factor of ~2.2 from z = 1.1 to z = 0.1

Diagnosing coherent structures in the convective boundary layer by optimizing their vertical turbulent scalar transfer

This is the final version. Available on open access from Springer via the DOI in this recordA new method is introduced to identify coherent structures in the convective boundary layer, based on optimizing the vertical scalar flux in a two-fluid representation of turbulent motions as simulated by a large-eddy simulation. The new approach partitions the joint frequency distribution (JFD) of the vertical velocity and a transported scalar into coherent structures (fluid 2) and their environment (fluid 1) by maximizing that part of the scalar flux resolved by the mean properties in fluid 2 and fluid 1. The proposed method does not rely on any a priori criteria for the partitioning of the flow nor any pre-assumptions about the shape of the JFD. Different flavours of the optimization approach are examined based on maximizing either the total (fluid 1 + fluid 2) or the fluid-2 resolved scalar flux, and on whether all possible partitions or only a subset are considered. These options can result in different derived area fractions for the coherent structures. The properties of coherent structures diagnosed by the optimization method are compared to the conditional sampling of a surface-emitted decaying tracer, in which coherent structures are defined as having tracer perturbation greater than some height-dependent threshold. Results show that the optimization method is able to smoothly define coherent thermal structures in both the horizontal and the vertical. Moreover, optimizing the turbulent transfer by the fluid-2 resolved flux produces very similar coherent structures to the tracer threshold method, especially in terms of their area fraction and updraft velocities. Nonetheless, further analysis of the partitioning of the JFD reveals that, even though the area fraction of coherent structures might be similar, their definition can occupy different quadrants of the JFD, implying the contribution of different physical mechanisms to the turbulent transfer in the boundary layer. Finally, the kinematic and thermodynamic characteristics of the coherent structures are examined based on their definition criteria.Natural Environment Research Council (NERC

Idealized large-eddy simulations of nocturnal low-level jets over subtropical desert regions and implications for dust-generating winds

Nocturnal low-level jets (LLJs) are maxima in the wind profile, which often form above the stable nocturnal boundary layer. Over the Sahara, the world's largest source of mineral dust, this phenomenon is of particular importance to the emission and transport of desert aerosol. We present the first ever detailed large-eddy simulations of dust-generating LLJs. Using sensitivity studies with the UK Met Office large-eddy model (LEM), two key controls of the nocturnal LLJ are investigated: surface roughness and the Coriolis force. Functional relationships derived from the LEM results help to identify optimal latitude-roughness configurations for a maximum LLJ enhancement. Ideal conditions are found in regions between 20 and 27°N with roughness lengths >0.0001 m providing long oscillation periods and large jet amplitudes. Typical LLJ enhancements reach up to 3.5 m s-1 for geostrophic winds of 10 m s-1. The findings are largely consistent with results from a theoretical LLJ model applied for comparison. The results demonstrate the importance of latitude and roughness in creating regional patterns of LLJ influence. Combining the functional relationships with high-resolution roughness data over northern Africa gives good agreement with the location of morning dust uplift in satellite observations. It is shown that shear-induced mixing plays an important role for the LLJ evolution and surface gustiness. With decreasing latitude the LLJ oscillation period is longer and, thus, shear-induced mixing is weaker, allowing a more stable nocturnal stratification to develop. This causes a later and more abrupt LLJ breakdown in the morning with stronger gusts, which can compensate for the slower LLJ evolution that leads to a weaker jet maximum. The findings presented here can serve as the first step towards a parametrization to improve the representation of the effects of nocturnal LLJs on dust emission in coarser-resolution models.European Research Counci