Abstract

Numerical weather prediction (NWP) models are now capable of operating at horizontal resolutions in the 100-m to 1-km range, a grid spacing similar in scale to that of the turbulent eddies present in the atmospheric convective boundary layer (CBL). Known as the ‘grey zone’ of turbulence, this regime is characterized by significant contributions from both the resolved and subgrid components to represent the dominant motions of the system. This study examines how the initiation of resolved turbulence – a concept commonly referred to as ‘spin-up’– can be delayed during the evolution of a simulated CBL in the grey zone. We identify the importance of imposed pseudo-random perturbations of potential temperature (theta ) for the development of the resolved fields showing that without such perturbations, resolved turbulence does not become established at all. When the perturbations are organized, spin-up can develop more rapidly, and we find that the earliest spin-up times can be achieved by applying an idealized profile of variance to derive the theta perturbation values. The perturbation structures are shown to be most effective when applied at intervals following the mixed-layer time scale, t_{*}, rather than perturbing only at the initial time. We also propose a modification to the three-dimensional Smagorinsky turbulence closure, in which the Smagorinsky constant is replaced by a scale-dependent coefficient. Both the approaches of: (1) applying structured theta perturbations, and (2) using a dynamically-evolving Smagorinsky coefficient are shown to encourage faster spin-up independently of each other, but the best results clearly emerge when the two methods are applied concurrently.

Highlights

  • Since its inception, numerical weather prediction (NWP) has always been governed by a tight balance; on the one hand using the most up-to-date science in representing atmospheric processes, and on the other, the best use of available computational resources

  • The difference between omitting and including random perturbations at time t + 0 is immediately apparent, and it is clear from Fig. 3a that no resolved motion develops during the 9-h period in the unperturbed case

  • To address the question of whether or not this resolved turbulence should be allowed in the grey zone at all, our results show that allowing such resolved turbulence provides much of the necessary non-local transport to complement the 3D static Smagorinsky scheme, even in the absence of the countergradient correction term that is characteristic of 1D planetary boundary-layer (PBL) schemes

Read more

Summary

Introduction

Numerical weather prediction (NWP) has always been governed by a tight balance; on the one hand using the most up-to-date science in representing atmospheric processes, and on the other, the best use of available computational resources. Ever since the earliest days of operational NWP in the 1950s, parametrization has been key in maintaining this balance, such that subgrid-scale physical processes could be represented and their influence be passed back to the grid. The grey zone of turbulence (or terra incognita, Wyngaard 2004) refers to resolutions at which turbulent eddies in the atmospheric boundary layer are partially resolved and partially parametrized, a regime that is emerging in the highest resolution mesoscale models. Eddies that are represented within the grey zone often appear to evolve in a grid-dependent way (Ching et al 2014; Zhou et al 2014), and since neither one-dimensional (1D) planetary boundary-layer (PBL) schemes (e.g. Lock et al 2000) nor three-dimensional (3D) large-eddy simulation (LES, Lilly 1967) formulations are designed to work within this regime, a necessity is growing to find new ways to represent these processes

Methods
Results
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call