Forest canopies play a critical role in affecting momentum and scalar transfer. Although there have been recent advances in numerical simulations of turbulent flows and scalar transfer across plant canopies and the atmosphere interface, few models have incorporated all important physical and physiological processes in subcanopy layers. Here we describe and evaluate an advanced multiple‐layer canopy module (MCANOPY), which is developed based largely on the Community Land Model version 4.5 and then coupled with the Weather Research and Forecasting model with large‐eddy simulations (WRF‐LES). The MCANOPY includes a suite of subcanopy processes, including radiation transfer, photosynthesis, canopy layer energy balance, momentum drag, and heat, water vapor, and CO2exchange between canopy layers and the canopy atmosphere. Numerical schemes for heat and water transport in soil, ground surface energy balance, and soil respiration are also included. Both the stand‐alone MCANOPY and the coupled system (the WRF‐LES‐MCANOPY) are evaluated against data measured in the Canopy Horizontal Array Turbulence Study field experiment. The MCANOPY performs reasonably well in reproducing vertical profiles of mean and turbulent flows as well as second‐order statistical quantities including heat and scalar fluxes within the canopy under unstable stability conditions. The coupled WRF‐LES‐MCANOPY captures major features of canopy edge flows under both neutral and unstable conditions. Limitations of the MCANOPY are discussed for our further work. Our results suggest that our model can be a promising modeling system for a variety of applications to study canopy flows and scalar transport (e.g., CO2).
- Award ID(s):
- 1821885
- NSF-PAR ID:
- 10353440
- Date Published:
- Journal Name:
- Journal of the Atmospheric Sciences
- ISSN:
- 0022-4928
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Sensitivity of large eddy simulations of tropical cyclone to sub-grid scale mixing parameterization.The surface wind structure and vertical turbulent transport processes in the eyewall of hurricane Isabel (2003) are investigated using six large-eddy simulations (LESs) with different horizontal grid spacing and three-dimensional (3D) sub-grid scale (SGS) turbulent mixing models and a convection permitting simulation that uses a coarser grid spacing and one-dimensional vertical turbulent mixing scheme. The mean radius-height distribution of storm tangential wind and radial flow, vertical velocity structure, and turbulent kinetic energy and momentum fluxes in the boundary layer generated by LESs are consistent with those derived from historical dropsonde composites, Doppler radar, and aircraft measurements. Unlike the convection permitting simulation that produces storm wind fields lacking small-scale disturbances, all LESs are able to produce sub-kilometer and kilometer scale eddy circulations in the eyewall. The inter-LES differences generally reduce with the decrease of model grid spacing. At 100-m horizontal grid spacing, the vertical momentum fluxes induced by the model-resolved eddies and the associated eddy exchange coefficients in the eyewall simulated by the LESs with different 3D SGS mixing schemes are fairly consistent. Although with uncertainties, the decomposition in terms of eddy scales suggests that sub-kilometer eddies are mainly responsible for the vertical turbulent transport within the boundary layer (~1 km depth following the conventional definition) whereas eddies greater than 1 km become the dominant contributors to the vertical momentum transport above the boundary layer in the eyewall. The strong dependence of vertical turbulent transport on eddy scales suggests that the vertical turbulent mixing parameterization in mesoscale simulations of tropical cyclones is ultimately a scale-sensitive problem.more » « less
-
null (Ed.)Abstract Surface friction contributes to tornado formation and maintenance by enhancing the convergence of angular momentum. The traditional lower boundary condition in atmospheric models typically assumes an instant equilibrium between the unresolved stress and the resolved shear. This assumption ignores the physics that turbulent motions are generated and dissipated at finite rates—in effect, turbulence has a memory through its lifetime. In this work, a modified lower boundary condition is proposed to account for the effect of turbulence memory. Specifically, when an air parcel moves along a curved trajectory, a normal surface-shear-stress component arises owing to turbulence memory. In the accompanying large-eddy simulation (LES) of idealized tornadoes, the normal surface-shear-stress component is a source of additional dynamic instability, which provides an extra pathway for the development of turbulent motions. The influence of turbulence memory on the intensity of quasi-steady-state tornadoes remains negligible as long as assumptions employed by the modified lower boundary condition hold over a relatively large fraction of the flow region of interest. However, tornadoes in a transient state may be especially sensitive to turbulence memory. Significance Statement Friction between the wind and the ground can influence atmospheric phenomena in important ways. For example, surface friction can be a significant source of rotation in some thunderstorms, and it can also help to intensify rotation when rotation is already present. Unfortunately, the representation of friction’s effects in atmospheric simulations is especially error-prone in phenomena characterized by rapid temporal evolution or strong spatial variations. Our work explores a new framework for representing friction to include the effect of the so-called turbulence memory. The approach is tested in idealized tornado simulations, but it may be applied to a wide range of atmospheric vortices.more » « less
-
more » « less
Abstract Populated urban areas along many coastal regions are vulnerable to landfalling tropical cyclones (TCs). To the detriment of surface parameterizations in mesoscale models, the complexities of turbulence at high TC wind speeds in urban canopies are presently poorly understood. Thus, this study explores the impacts of urban morphology on TC-strength winds and boundary layer turbulence in landfalling TCs. To better quantify how urban structures interact with TC winds, large-eddy simulations (LESs) are conducted with the Cloud Model 1 (CM1). This implementation of CM1 includes immersed boundary conditions (IBCs) to represent buildings and eddy recycling to maintain realistic turbulent flow perturbations. Within the IBCs, an idealized coastal city with varying scales is introduced. TC winds impinge perpendicularly to the urbanized coastline. Numerical experiments show that buildings generate distinct, intricate flow patterns that vary significantly as the city structure is varied. Urban IBCs produce much stronger turbulent kinetic energy than is produced by conventional surface parameterizations. Strong effective eddy viscosity due to resolved eddy mixing is displayed in the wake of buildings within the urban canopy, while deep and enhanced effective eddy viscosity is present downstream. Such effects are not seen in a comparison LES using a simple surface parameterization with high roughness values. Wind tunneling effects in streamwise canyons enhance pedestrian-level winds well beyond what is possible without buildings. In the arena of regional mesoscale modeling, this type of LES framework with IBCs can be used to improve parameters in surface and boundary layer schemes to more accurately represent the drag coefficient and the eddy viscosity in landfalling TC boundary layers.
Significance Statement This is among the first large-eddy simulation model studies to examine the impacts of tropical cyclone–like winds around explicitly resolved buildings. This work is a step forward in bridging the gap between engineering studies that use computational fluid dynamics models or laboratory experiments for flow through cities and mesoscale model simulations of landfalling tropical cyclones that use surface parameterizations specialized for urban land use.
-
more » « less
This study utilizes a large-eddy simulation (LES) approach to systematically assess the directional variability of wave-driven Langmuir turbulence (LT) in the ocean surface boundary layer (OSBL) under tropical cyclones (TCs). The Stokes drift vector, which drives LT through the Craik–Leibovich vortex force, is obtained through spectral wave simulations. LT’s direction is identified by horizontally elongated turbulent structures and objectively determined from horizontal autocorrelations of vertical velocities. In spite of a TC’s complex forcing with great wind and wave misalignments, this study finds that LT is approximately aligned with the wind. This is because the Reynolds stress and the depth-averaged Lagrangian shear (Eulerian plus Stokes drift shear) that are key in determining the LT intensity (determined by normalized depth-averaged vertical velocity variances) and direction are also approximately aligned with the wind relatively close to the surface. A scaling analysis of the momentum budget suggests that the Reynolds stress is approximately constant over a near-surface layer with predominant production of turbulent kinetic energy by Stokes drift shear, which is confirmed from the LES results. In this layer, Stokes drift shear, which dominates the Lagrangian shear, is aligned with the wind because of relatively short, wind-driven waves. On the contrary, Stokes drift exhibits considerable amount of misalignments with the wind. This wind–wave misalignment reduces LT intensity, consistent with a simple turbulent kinetic energy model. Our analysis shows that both the Reynolds stress and LT are aligned with the wind for different reasons: the former is dictated by the momentum budget, while the latter is controlled by wind-forced waves.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/JAS-D-21-0033.1
