- Award ID(s):
- 1712290
- PAR ID:
- 10158496
- Date Published:
- Journal Name:
- Journal of the Atmospheric Sciences
- ISSN:
- 0022-4928
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract We investigate the sensitivity of mesoscale atmospheric predictability to the slope of the background kinetic energy spectrum E by adding initial errors to simulations of idealized moist midlatitude cyclones at several wavenumbers k for which the slope of E (k) is significantly different. These different slopes arise from 1) differences in the E (k) generated by cyclones growing in two different moist baroclinically unstable environments, and 2) differences in the horizontal scale at which initial perturbations are added, with E (k) having steeper slopes at larger scales. When small-amplitude potential temperature perturbations are added, the error growth through the subsequent 36-h simulation is not sensitive to the slope of E (k) nor to the horizontal scale of the initial error. In all cases with small-amplitude perturbations, the error growth in physical space is dominated by moist convection along frontal boundaries. As such, the error field is localized in physical space and broad in wavenumber (spectral) space. In moist midlatitude cyclones, these broadly distributed errors in wavenumber space limit mesoscale predictability by growing up-amplitude rather than by cascading upscale to progressively longer wavelengths. In contrast, the error distribution in homogeneous turbulence is broad in physical space and localized in wavenumber space, and dimensional analysis can be used to estimate the error growth rate at a specific wavenumber k from E (k). Predictability estimates derived in this manner, and from the numerical solutions of idealized models of homogeneous turbulence, depend on whether the slope of E (k) is shallower or steeper than k^ −3 , which differs from the slope-insensitive behavior exhibited by moist midlatitude cyclones.more » « less
-
Abstract The predictability of passive scalar dispersion is of both theoretical interest and practical importance, for example for high‐resolution numerical weather prediction and air quality modeling. However, the implications for the numerical modeling of urban areas remain relatively unexplored. Using obstacle‐resolving large‐eddy simulations (LES), we conducted twin experiments, with and without a velocity perturbation, to investigate how the presence of urban roughness affects error growth in streamwise velocity ( u ) and passive scalar ( θ ) fields, as well as the differences between error evolutions in u and θ fields. The predictability limit is characterized using the signal‐to‐noise ratio (SNR) as a continuous metric to indicate when error reaches saturation. The presence of urban roughness decreases of the passive scalar by around 20% compared to cases without them. The error statistics of θ indicate that urban roughness‐induced flow structures and different scalar source locations affect the scalar dispersion and relative fluctuations, which subsequently dictate the evolution of the SNR. Analysis of the passive scalar error energy ( ϵ θ 2 ) budget indicates that the contributions from advective transport by the velocity and velocity error dominate. The error energy spectra of both u and θ exhibit a −5/3 slope in flat‐wall cases, but not in the presence of urban roughness, thereby highlighting the deviation from the assumption of locally isotropic turbulence. This study reveals that urban roughness can decrease the predictability of the passive scalar and destroy the similarity between the error statistics of the velocity and the passive scalar.more » « less
-
It is recognized that the atmosphere’s predictability is intrinsically limited by unobservably small uncertainties that are beyond our capability to eliminate. However, there have been discussions in recent years on whether forecast error grows upscale (small-scale error grows faster and transfers to progressively larger scales) or up-amplitude (grows at all scales at the same time) when unobservably small-amplitude initial uncertainties are imposed at the large scales and limit the intrinsic predictability. This study uses large-scale small-amplitude initial uncertainties of two different structures—one idealized, univariate, and isotropic, the other realistic, multivariate, and flow dependent—to examine the error growth characteristics in the intrinsic predictability regime associated with a record-breaking rainfall event that happened on 19–20 July 2021 in China. Results indicate upscale error growth characteristics regardless of the structure of the initial uncertainties: the errors at smaller scales grow fastest first; as the forecasts continue, the wavelengths of the fastest error growth gradually shift toward larger scales with reduced error growth rates. Therefore, error growth from smaller to larger scales was more important than the growth directly at the large scales of the initial errors. These upscale error growth characteristics also depend on the perturbed and examined quantities: if the examined quantity is perturbed, then its errors grow upscale; if there is no initial uncertainty in the examined quantity, then its errors grow at all scales at the same time, although its smaller-scale errors still grow faster for the first several hours, suggesting the existence of the upscale error growth.more » « less
Significance Statement This study compared the error growth characteristics associated with the atmosphere’s intrinsic predictability under two different structures of unobservably small-amplitude, large-scale initial uncertainties: one idealized, univariate, and isotropic, the other realistic, multivariate, and flow dependent. The characteristics of the errors growing upscale rather than up-amplitude regardless of the initial uncertainties’ structure are apparent. The large-scale errors do not grow if their initial amplitudes are much bigger than the small-scale errors. This study also examined how the error growth characteristics will change when the quantity that is used to describe the error growth is inconsistent with the quantity that contains uncertainty, suggesting the importance of including multivariate, covariant uncertainties of state variables in atmospheric predictability studies.
-
Abstract Analytical solutions for the three‐dimensional groundwater flow equation have been widely used to gain insight about subsurface flow structure and as an alternative to computationally expensive numerical models. Of particular interest are solutions that decompose prescribed hydraulic head boundaries (e.g., Dirichlet boundary condition) into a collection of harmonic functions. Previous studies estimate the frequencies and amplitudes of these harmonics with a least‐square approach where the amplitudes are fitted given a pre‐assigned set of frequencies. In these studies, an ad hoc and structured discretization of the frequency domain is typically used, excluding dominant frequencies while assigning importance to spurious frequencies, with significant consequences for estimating the fluxes and residence times. This study demonstrates the advantages of using a pre‐assigned frequency spectrum that targets the dominant frequencies based on rigorous statistical analysis with predefined significance levels. The new approach is tested for three hydrologic conceptualizations: (a) a synthetic periodic basin, (b) synthetic bedforms, and (c) a natural mountainous watershed. The performance of the frequency spectrum selection is compared with exact analytical or approximate numerical solutions. We found that the new approach better describes the fluxes and residence times for Dirichlet boundaries with well‐defined characteristics spatial scales (e.g., periodic basins and bedforms). For more complex scenarios, such as natural mountainous watersheds, both pre‐assigned frequency spectrums present similar performance. The spectral solutions presented here can play a central role in developing reduced‐complexity models for assessing regional water and solute fluxes within mountain watersheds and hyporheic zones.
-
Abstract In the atmospheric surface layer (ASL), a characteristic wavelength marking the limit between energy‐containing and inertial subrange scales can be defined from the vertical velocity spectrum. This wavelength is related to the integral length scale of turbulence, used in turbulence closure approaches for the ASL. The scaling laws describing the displacement of this wavelength with changes in atmospheric stability have eluded theoretical treatment and are considered here. Two derivations are proposed for mildly unstable to mildly stable ASL flows one that only makes use of normalizing constraints on the vertical velocity variance along with idealized spectral shapes featuring production to inertial subrange regimes, while another utilizes a co‐spectral budget with a return‐to‐isotropy closure. The expressions agree with field experiments and permit inference of the variations of the wavelength with atmospheric stability. This methodology offers a new perspective for numerical and theoretical modeling of ASL flows and for experimental design.