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).
Mesoscale climate models provide indispensable tools to understand land‐atmosphere interactions over urban regions. However, uncertainties in urban canopy parameters (UCPs) and parameterization schemes lead to degraded representation of the drag effect in complex built terrains. In particular, for the widely applied single‐layer urban canopy model (SLUCM) coupled with the Weather Research and Forecasting (WRF) model, near‐surface horizontal wind speed is known to be overestimated systematically. In this study, idealized large eddy simulations (LES) and WRF‐SLUCM simulations are conducted to study the separate effect of UCPs and aerodynamic parameterization on atmospheric boundary layer processes and rainfall variabilities in Phoenix, Arizona. For LES that explicitly resolves surface geometry, significant differences between three‐dimensional (3D) versus two‐dimensional (2D) representation of urban morphology are found in the surface layer and above. When surface drag is parameterized following SLUCM, surface morphologies have little impacts on the mean momentum transfer. WRF‐SLUCM simulation results, incorporated with 3D urban morphology data, indicate that simply refining the frontal area index will reduce the surface drag, which further amplifies the systematic positive bias of SLUCM in predicting horizontal wind speed. Replacing the drag parameterization in SLUCM by LES‐based aerodynamic parameters has evident impacts on near‐surface wind speed. The impact of urban roughness representation becomes the most evident during rainfall periods, due to the important role of surface drag in dictating moisture convergence. Our study underlines that apart from intensive efforts in obtaining detailed UCPs, it is also critical to enhance the urban momentum exchange parameterization schemes.
more » « less- NSF-PAR ID:
- 10446485
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Atmospheres
- Volume:
- 126
- Issue:
- 4
- ISSN:
- 2169-897X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
Urban canopy models (UCMs) in mesoscale numerical weather prediction models need evaluation to understand biases in urban environments under a range of conditions. The authors evaluate a new drag formula in the Weather Research and Forecasting (WRF) model’s multilayer UCM, the Building Effect Parameterization combined with the Building Energy Model (BEP+BEM), against both in-situ measurements in the urban environment as well as simulations with a simple bulk scheme and BEP+BEM using the old drag formula. The new drag formula varies with building packing density, while the old drag formula is constant. The case study is a strong cold frontal passage that occurred in Houston during the winter of 2017, causing high winds. It is found that both BEP+BEM simulations have lower peak wind speeds, consistent with near-surface measurements, while the bulk simulation has winds that are too strong. The constant-drag BEP+BEM simulation has a near-zero wind speed bias, while the new-drag simulation has a negative bias. Although the focus is on the impact of drag on the urban wind speeds, both BEP+BEM simulations have larger negative biases in the near-surface temperature than the bulk-scheme simulation. Reasons for the different performances are discussed.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.
-
null (Ed.)Abstract The simulated winds within the urban canopy of landfalling tropical cyclones are sensitive to the representation of the planetary-boundary and urban-canopy layers in numerical weather prediction models. In order to assess the sub-grid-scale parameterizations of these layers, mesoscale model simulations were executed and evaluated against near-surface observations as the outer wind field of Hurricane Irma (2017) interacted with the built-up region from downtown Miami northward to West Palm Beach. Four model simulations were examined, comprised of two different planetary boundary layer (PBL) parameterizations (a local closure scheme with turbulent kinetic energy prediction and a nonlocal closure scheme) and two different urban canopy models (UCMs) [a zeroth order bulk scheme and a multilayer Building Effect Parameterization (BEP) that mimics the three-dimensionality of buildings]. Overall, the simulated urban canopy winds were weakly sensitive to the PBL scheme and strongly sensitive to the UCM. The bulk simulations compared most favorably to an analyzed wind swath in the urban environment, while the BEP simulations had larger negative biases in the same region. There is uncertainty in magnitude of the urban environment biases due to the lack of many urban sheltered measurements in the wind swath analysis. Biases in the rural environment were similar among the bulk and BEP simulations. An improved comparison with the analyzed wind swath in the urban region was obtained by reducing the drag coefficient in BEP in one of the PBL schemes. The usefulness of BEP was demonstrated in its ability to predict realistic heterogeneous near-surface velocity patterns in urban regions.more » « less
-
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
