Climate models struggle to accurately represent the highly reflective boundary layer clouds overlying the remote and stormy Southern Ocean. We use in situ aircraft observations from the Southern Ocean Clouds, Radiation and Aerosol Transport Experimental Study (SOCRATES) to evaluate Southern Ocean clouds in a cloud‐resolving large‐eddy simulation (LES) and two coarse resolution global atmospheric models, the CESM Community Atmosphere Model (CAM6) and the GFDL Atmosphere Model (AM4), run in a nudged hindcast framework. We develop six case studies from SOCRATES data which span the range of observed cloud and boundary layer properties. For each case, the LES is run once forced purely using reanalysis data (fifth generation European Centre for Medium‐Range Weather Forecasts atmospheric reanalysis, “ERA5 based”) and once strongly nudged to an aircraft profile(“Obs based”). The ERA5‐based LES can be compared with the global models, which are also nudged to reanalysis data and are better for simulating cumulus. The Obs‐based LES closely matches an observed cloud profile and is useful for microphysical comparisons and sensitivity tests and simulating multilayer stratiform clouds. We use two‐moment Morrison microphysics in the LES and find that it simulates too few frozen particles in clouds occurring within the Hallett‐Mossop temperature range. We tweak the Hallett‐Mossop parameterization so that it activates within boundary layer clouds, and we achieve better agreement between observed and simulated microphysics. The nudged global climate models (GCMs) simulate liquid‐dominated mixed‐phase clouds in the stratiform cases but excessively glaciate cumulus clouds. Both GCMs struggle to represent two‐layer clouds, and CAM6 has low droplet concentrations in all cases and underpredicts stratiform cloud‐driven turbulence.
- Award ID(s):
- 1938108
- NSF-PAR ID:
- 10299814
- Date Published:
- Journal Name:
- Monthly Weather Review
- ISSN:
- 0027-0644
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
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
Abstract In an effort to better represent aerosol transport in mesoscale and global‐scale models, large eddy simulations (LES) from the National Center for Atmospheric Research (NCAR) Turbulence with Particles (NTLP) code are used to develop a Markov chain random walk model that predicts aerosol particle profiles in a cloud‐free marine atmospheric boundary layer (MABL). The evolution of vertical concentration profiles are simulated for a range of aerosol particle sizes and in a neutral and an unstable boundary layer. For the neutral boundary layer we find, based on the LES statistics and a specific model time step, that there exist significant correlation for particle positions, meaning that particles near the bottom of the boundary are more likely to remain near the bottom of the boundary layer than being abruptly transported to the top, and vice versa. For the unstable boundary layer, a similar time interval exhibits a weaker tendency for an aerosol particle to remain close to its current location compared to the neutral case due to the strong nonlocal convective motions. In the limit of a large time interval, particles have been mixed throughout the MABL and virtually no temporal correlation exists. We leverage this information to parameterize a Markov chain random walk model that accurately predicts the evolution of vertical concentration profiles. The new methodology has significant potential to be applied at the subgrid level for coarser‐scale weather and climate models, the utility of which is shown by comparison to airborne field data and global aerosol models.
-
Abstract. The observing system design of multidisciplinary fieldmeasurements involves a variety of considerations on logistics, safety, andscience objectives. Typically, this is done based on investigator intuitionand designs of prior field measurements. However, there is potential forconsiderable increases in efficiency, safety, and scientific success byintegrating numerical simulations in the design process. Here, we present anovel numerical simulation–environmental response function (NS–ERF)approach to observing system simulation experiments that aidssurface–atmosphere synthesis at the interface of mesoscale and microscalemeteorology. In a case study we demonstrate application of the NS–ERFapproach to optimize the Chequamegon Heterogeneous Ecosystem Energy-balanceStudy Enabled by a High-density Extensive Array of Detectors 2019(CHEESEHEAD19). During CHEESEHEAD19 pre-field simulation experiments, we considered theplacement of 20 eddy covariance flux towers, operations for 72 h oflow-altitude flux aircraft measurements, and integration of various remotesensing data products. A 2 h high-resolution large eddy simulationcreated a cloud-free virtual atmosphere for surface and meteorologicalconditions characteristic of the field campaign domain and period. Toexplore two specific design hypotheses we super-sampled this virtualatmosphere as observed by 13 different yet simultaneous observing systemdesigns consisting of virtual ground, airborne, and satellite observations.We then analyzed these virtual observations through ERFs to yield an optimalaircraft flight strategy for augmenting a stratified random flux towernetwork in combination with satellite retrievals. We demonstrate how the novel NS–ERF approach doubled CHEESEHEAD19'spotential to explore energy balance closure and spatial patterning scienceobjectives while substantially simplifying logistics. Owing to its modularextensibility, NS–ERF lends itself to optimizing observing system designs alsofor natural climate solutions, emission inventory validation, urban airquality, industry leak detection, and multi-species applications, among otheruse cases.more » « less
-
more » « less
Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios (
qc andqr , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulatedqc closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulatedqr properties, though, with bulk simulatedqr consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produceqr that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1175/MWR-D-20-0328.1
