Journal of the Atmospheric Sciences
Latest Publications


TOTAL DOCUMENTS

13295
(FIVE YEARS 722)

H-INDEX

234
(FIVE YEARS 9)

Published By American Meteorological Society

1520-0469, 0022-4928

Abstract Intermittent transitions between turbulent and non-turbulent states are ubiquitous in the stable atmospheric surface layer (ASL). Data from two field experiments in Utqiagvik, Alaska, and from direct numerical simulations are used to probe these state transitions so as to (i) identify statistical metrics for the detection of intermittency, (ii) probe the physical origin of turbulent bursts, and (iii) quantify intermittency effects on overall fluxes and their representation in closure models. The analyses reveal three turbulence regimes, two of which correspond to weakly turbulent periods accompanied by intermittent behavior (regime 1: intermittent, regime 2: transitional), while the third is associated with a fully turbulent flow. Based on time series of the turbulence kinetic energy (TKE), two non-dimensional parameters are proposed to diagnostically categorize the ASL state into these regimes; the first characterizes the weakest turbulence state, while the second describes the range of turbulence variability. The origins of intermittent turbulence activity are then investigated based on the TKE budget over the identified bursts. While the quantitative results depend on the height, the analyses indicate that these bursts are predominantly advected by the mean flow, produced locally by mechanical shear, or lofted from lower levels by turbulent ejections. Finally, a new flux model is proposed using the vertical velocity variance in combination with different mixing length scales. The model provides improved representation (correlation coefficients with observations of 0.61 for momentum and 0.94 for sensible heat) compared to Monin–Obukhov similarity (correlation coefficients of 0.0047 for momentum and 0.49 for sensible heat), thus opening new pathways for improved parametrizations in coarse atmospheric models.


Abstract A dry-core idealized general circulation model with a stratospheric polar vortex in the northern hemisphere is run with a combination of simplified topography and imposed tropospheric temperature perturbations, each located in the northern hemisphere with a zonal wave number of one. The phase difference between the imposed temperature wave and the topography is varied to understand what effect this has on the occurrence of polar vortex displacements. Geometric moments are used to identify the centroid of the polar vortex for the purposes of classifying whether or not the polar vortex is displaced. Displacements of the polar vortex are a response to increased tropospheric wave activity. Compared to a model run with only topography, the likelihood of the polar vortex being displaced increases when the warm region is located west of the topography peak, and decreases when the cold region is west of the topography peak. This response from the polar vortex is due to the modulation of vertically propogating wave activity by the temperature forcing. When the southerly winds on the western side of the topographically forced anticyclone are collocated with warm or cold temperature forcing, the vertical wave activity flux in the troposphere becomes more positive or negative, respectively. This is in line with recent reanalysis studies which showed that anomalous warming west of the surface pressure high, in the climatological standing wave, precedes polar vortex disturbances.


Abstract In a recent study by Wang et al. (2021a) that introduced a dynamical efficiency to the intensification potential of a tropical cyclone (TC) system, a simplified energetically based dynamical system (EBDS) model was shown to be able to capture the intensity-dependence of TC potential intensification rate (PIR) in both idealized numerical simulations and observations. Although the EBDS model can capture the intensity-dependence of TC intensification as in observations, a detailed evaluation has not yet been done. This study provides an evaluation of the EBDS model in reproducing the intensity-dependent feature of the observed TC PIR based on the best-track data for TCs over the North Atlantic, central, eastern and western North Pacific during 1982–2019. Results show that the theoretical PIR estimated by the EBDS model can capture basic features of the observed PIR reasonably well. The TC PIR in the best-track data increases with increasing relative TC intensity (intensity normalized by its corresponding maximum potential intensity–MPI) and reaches a maximum at an intermediate relative intensity around 0.6, and then decreases with increasing relative intensity to zero as the TC approaches its MPI, as in idealized numerical simulations. Results also show that the PIR for a given relative intensity increases with the increasing MPI and thus increasing sea surface temperature, which is also consistent with the theoretical PIR implied by the EBDS model. In addition, future directions to include environmental effects and make the EBDS model applicable to predict intensity change of real TCs are also discussed.


Abstract The properties of diurnal variability in tropical cyclones (TCs) and the mechanisms behind them remain an intriguing aspect of TC research. This study provides a comprehensive analysis of diurnal variability in two simulations of TCs to explore these mechanisms. One simulation is a well known Hurricane Nature Run, which is a realistic simulation of a TC produced using the Weather Research and Forecasting model (WRF). The other simulation is a realistic simulation produced using WRF of Hurricane Florence (2018) using hourly ERA5 reanalysis data as input. Empirical orthogonal functions and Fourier filtering are used to analyze diurnal variability in the TCs. In both simulations a diurnal squall forms at sunrise in the inner core and propagates radially outwards and intensifies until midday. At midday the upper-level outflow strengthens, surface inflow weakens, and the cirrus canopy reaches its maximum height and radial extent. At sunset and overnight, the surface inflow is stronger, and convection inside the RMW peaks. Therefore, two diurnal cycles of convection exist in the TCs with different phases of maxima: eyewall convection at sunset and at night, and rainband convection in the early morning. This study finds that the diurnal pulse in the cirrus canopy is not advectively-driven, nor can it be attributed to weaker inertial stability at night; rather, the results indicate direct solar heating as a mechanism for cirrus canopy lifting and enhanced daytime outflow. These results show a strong diurnal modulation of tropical cyclone structure, and are consistent with other recent observational and modeling studies of the TC diurnal cycle.


Abstract Understanding ice development in Cumulus Congestus (CuCg) clouds, which are ubiquitous globally, is critical for improving our knowledge of cloud physics, cloud resolution and climate prediction models. Results presented here are representative of data collected in 1,008 penetrations of moderate to strong updrafts in CuCg clouds by five research aircraft in six geographic locations. The results show that CuCg with warm (> ∼20°C) cloud base temperatures, such as in tropical marine environments, experience a strong collision-coalescence process. Development of coalescence is also correlated with drop effective radius > ∼12 to 14 µm in diameter. Increasing the cloud-base drop concentration with diameters from 15 to 35 µm and decreasing the drop concentration < 15 µm appears to enhance coalescence. While the boundary-layer aerosol population is not a determinate factor in development of coalescence in tropical marine environments, its impact on coalescence is not yet fully determined. Some supercooled large drops generated via coalescence fracture when freezing, producing a secondary ice process (SIP) with production of copious small ice particles that naturally seed the cloud. The SIP produces an avalanche effect, freezing the majority of supercooled liquid water before fresh updrafts reach the −16°C level. Conversely, CuCg with cloud base temperatures ≤ ∼8°C develop significant concentrations of ice particles at colder temperatures, so that small supercooled water drops are lofted to higher elevations before freezing. Recirculation of ice in downdrafts at the edges of updrafts appears to be the primary mechanism for development of precipitation in CuCg with colder cloud base temperatures.


Abstract Four state-of-the-science numerical weather prediction (NWP) models were used to perform mountain wave- (MW) resolving hind-casts over the Drake Passage of a 10-day period in 2010 with numerous observed MW cases. The Integrated Forecast System (IFS) and the Icosahedral Nonhydrostatic (ICON) model were run at Δx ≈ 9 and 13 km globally. TheWeather Research and Forecasting (WRF) model and the Met Office Unified Model (UM) were both configured with a Δx = 3 km regional domain. All domains had tops near 1 Pa (z ≈ 80 km). These deep domains allowed quantitative validation against Atmospheric InfraRed Sounder (AIRS) observations, accounting for observation time, viewing geometry, and radiative transfer. All models reproduced observed middle-atmosphere MWs with remarkable skill. Increased horizontal resolution improved validations. Still, all models underrepresented observed MW amplitudes, even after accounting for model effective resolution and instrument noise, suggesting even at Δx ≈ 3 km resolution, small-scale MWs are under-resolved and/or over-diffused. MWdrag parameterizations are still necessary in NWP models at current operational resolutions of Δx ≈ 10 km. Upper GW sponge layers in the operationally configured models significantly, artificially reduced MW amplitudes in the upper stratosphere and mesosphere. In the IFS, parameterized GW drags partly compensated this deficiency, but still, total drags were ≈ 6 time smaller than that resolved at Δx ≈ 3 km. Meridionally propagating MWs significantly enhance zonal drag over the Drake Passage. Interestingly, drag associated with meridional fluxes of zonal momentum (i.e. ) were important; not accounting for these terms results in a drag in the wrong direction at and below the polar night jet.


Abstract Large-eddy simulations are conducted to investigate and physically interpret the impacts of heterogeneous, low terrain on deep-convection initiation (CI). The simulations are based on a case of shallow-to-deep convective transition over the Amazon River basin, and use idealized terrains with varying levels of ruggedness. The terrain is designed by specifying its power-spectral shape in wavenumber space, inverting to physical space assuming random phases for all wave modes, and scaling the terrain to have a peak height of 200 m. For the case in question, these modest terrain fields expedite CI by up to 2-3 h, largely due to the impacts of the terrain on the size of, and subcloud support for, incipient cumuli. Terrain-induced circulations enhance subcloud kinetic energy on the mesoscale, which is realized as wider and longer-lived subcloud circulations. When the updraft branches of these circulations breach the level of free convection, they initiate wider and more persistent cumuli that subsequently undergo less entrainment-induced cloud dilution and detrainment-induced mass loss. As a result, the clouds become more vigorous and penetrate deeper into the troposphere. Larger-scale terrains are more effective than smaller-scale terrains in promoting CI because they induce larger enhancements in both the width and the persistence of subcloud updrafts.


Author(s):  
Joshua G. Gebauer ◽  
David B. Parsons

Abstract An analytical model is presented for the generation of a Blackadar-like nocturnal low-level jet in a broad baroclinic zone. The flow is forced from below (flat ground) by a surface buoyancy gradient and from above (free atmosphere) by a constant pressure gradient force. Diurnally-varying mixing coefficients are specified to increase abruptly at sunrise and decrease abruptly at sunset. With attention restricted to a surface buoyancy that varies linearly with a horizontal coordinate, the Boussinesq-approximated equations of motion, thermal energy, and mass conservation reduce to a system of one-dimensional equations that can be solved analytically. Sensitivity tests with southerly jets suggest that (i) stronger jets are associated with larger decreases of the eddy viscosity at sunset (as in Blackadar theory), (ii) the nighttime surface buoyancy gradient has little impact on jet strength, and (iii) for pure baroclinic forcing (no free-atmosphere geostrophic wind), the nighttime eddy diffusivity has little impact on jet strength, but the daytime eddy diffusivity is very important and has a larger impact than the daytime eddy viscosity. The model was applied to a jet that developed in fair weather conditions over the Great Plains from southern Texas to northern South Dakota on 1 May 2020. The ECMWF Reanalysis v5 (ERA5) for the afternoon prior to jet formation showed that a broad north-south-oriented baroclinic zone covered much of the region. The peak model-predicted winds were in good agreement with ERA5 winds and lidar data from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) central facility in north-central Oklahoma.


Abstract We consider the closure problem of representing the higher order moments (HOMs) in terms of lower-order moments, a central feature in turbulence modelling based on the Reynolds-Averaged Navier-Stokes (RANS) approach. Our focus is on models suited for the description of asymmetric, non-local and semi-organized turbulence in the dry atmospheric convective boundary layer (CBL). We establish a multivariate probability density function (PDF) describing populations of plumes which are embedded in a sea of weaker randomly spaced eddies, and apply an assumed Delta-PDF approximation. The main content of this approach consists of capturing the bulk properties of the PDF. We solve the closure problem analytically for all relevant higher order moments (HOMs) involving velocity components and temperature and establish a hierarchy of new non-Gaussian turbulence closure models of different content and complexity ranging from analytical to semi-analytical. All HOMs in the hierarchy have a universal and simple functional form. They refine the widely used Millionshchikov closure hypothesis and generalize the famous quadratic skewness-kurtosis relationship to higher-order. We examine the performance of the new closures by comparison with measurement, LES and DNS data and derive empirical constants for semi-analytical models, which are best for practical applications. We show that the new models have a good skill in predicting the HOMs for atmospheric CBL. Our closures can be implemented in second-, third- and fourth-order RANS turbulence closure models of bi-, tri-and four-variate levels of complexity. Finally, several possible generalizations of our approach are discussed.


Abstract Convective self-aggregation refers to a phenomenon in which random convection can self-organize into large-scale clusters over an ocean surface with uniform temperature in cloud-resolving models. Previous literature studies convective aggregation primarily by analyzing vertically integrated (VI) moist static energy (MSE) variance. That is the global MSE variance, including both the local MSE variance at a given altitude and the covariance of MSE anomalies between different altitudes. Here we present a vertically resolved (VR) MSE framework that focuses on the local MSE variance to study convective self-aggregation. Using a cloud-resolving simulation, we show that the development of self-aggregation is associated with an increase of local MSE variance, and that the diabatic and adiabatic generation of the MSE variance is mainly dominated by the boundary layer (BL, the lowest 2 km). The results agree with recent numerical simulation results and the available potential energy analyses showing that the BL plays a key role in the development of self-aggregation. Additionally, we find that the lower free troposphere (2 - 4 km) also generates significant MSE variance in the first 15 days. We further present a detailed comparison between the global and local MSE variance frameworks in their mathematical formulation and diagnostic results, highlighting their differences.


Sign in / Sign up

Export Citation Format

Share Document