skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

Title: A New Theoretical Framework for Understanding Multiscale Atmospheric Predictability
Here we present a new theoretical framework that connects the error growth behavior in numerical weather prediction (NWP) with the atmospheric kinetic energy spectrum. Building on previous studies, our newly proposed framework applies to the canonical observed atmospheric spectrum that has a -3 slope at synoptic scales and a -5/3 slope at smaller scales. Based on this realistic hybrid energy spectrum, our new experiment using hybrid numerical models provides reasonable estimations for the finite predictable ranges at different scales. We further derive an analytical equation that helps understand the error growth behavior. Despite its simplicity, this new analytical error growth equation is capable of capturing the results of previous comprehensive theoretical and observational studies of atmospheric predictability. The success of this new theoretical framework highlights the combined effects of quasi-two-dimensional dynamics at synoptic-scales (-3 slope) and three-dimensional turbulence-like small-scale chaotic flows (-5/3 slope) in dictating the error growth. It is proposed that this new framework could serve as a guide for understanding and estimating the predictability limit in the real world.  more » « less
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of the Atmospheric Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. We analyse the pressure-driven flow of the Oldroyd-B fluid in slowly varying arbitrarily shaped, narrow channels and present a theoretical framework for calculating the relationship between the flow rate $q$ and pressure drop $\Delta p$ . We first identify the characteristic scales and dimensionless parameters governing the flow in the lubrication limit. Employing a perturbation expansion in powers of the Deborah number ( $De$ ), we provide analytical expressions for the velocity, stress and the $q$ – $\Delta p$ relation in the weakly viscoelastic limit up to $O(De^2)$ . Furthermore, we exploit the reciprocal theorem derived by Boyko $\&$ Stone ( Phys. Rev. Fluids , vol. 6, 2021, L081301) to obtain the $q$ – $\Delta p$ relation at the next order, $O(De^3)$ , using only the velocity and stress fields at the previous orders. We validate our analytical results with two-dimensional numerical simulations in the case of a hyperbolic, symmetric contracting channel and find excellent agreement. While the velocity remains approximately Newtonian in the weakly viscoelastic limit (i.e. the theorem of Tanner and Pipkin), we reveal that the pressure drop strongly depends on the viscoelastic effects and decreases with $De$ . We elucidate the relative importance of different terms in the momentum equation contributing to the pressure drop along the symmetry line and identify that a pressure drop reduction for narrow contracting geometries is primarily due to gradients in the viscoelastic shear stresses. We further show that, although for narrow geometries the viscoelastic axial stresses are negligible along the symmetry line, they are comparable or larger than shear stresses in the rest of the domain. 
    more » « less
  4. 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.

    more » « less
  5. null (Ed.)
    The marine-based West Antarctic Ice Sheet (WAIS) is currently locally retreating because of shifting wind-driven oceanic currents that transport warm waters toward the ice margin, resulting in ice shelf thinning and accelerated mass loss. Previous results from geologic drilling on Antarctica’s continental margins show significant variability in ice sheet extent during the late Neogene and Quaternary. Climate and ice sheet models indicate a fundamental role for oceanic heat in controlling ice sheet variability over at least the past 20 My. Although evidence for past ice sheet variability is available from ice-proximal marine settings, sedimentary sequences from the continental shelf and rise are required to evaluate the extent of past ice sheet variability and the associated forcings and feedbacks. International Ocean Discovery Program Expedition 374 drilled a latitudinal and depth transect of five sites from the outer continental shelf to rise in the central Ross Sea to resolve Neogene and Quaternary relationships between climatic and oceanic change and WAIS evolution. The Ross Sea was targeted because numerical ice sheet models indicate that this sector of Antarctica responds sensitively to changes in ocean heat flux. Expedition 374 was designed for optimal data-model integration to enable an improved understanding of Antarctic Ice Sheet (AIS) mass balance during warmer-than-present climates (e.g., the Pleistocene “super interglacials,” the mid-Pliocene, and the Miocene Climatic Optimum). The principal goals of Expedition 374 were to: 1. Evaluate the contribution of West Antarctica to far-field ice volume and sea level estimates; 2. Reconstruct ice-proximal oceanic and atmospheric temperatures to quantify past polar amplification; 3. Assess the role of oceanic forcing (e.g., temperature and sea level) on AIS variability; 4. Identify the sensitivity of the AIS to Earth’s orbital configuration under a variety of climate boundary conditions; and 5. Reconstruct Ross Sea paleobathymetry to examine relationships between seafloor geometry, ice sheet variability, and global climate. To achieve these objectives, postcruise studies will: 1. Use data and models to reconcile intervals of maximum Neogene and Quaternary ice advance and retreat with far-field records of eustatic sea level; 2. Reconstruct past changes in oceanic and atmospheric temperatures using a multiproxy approach; 3. Reconstruct Neogene and Quaternary sea ice margin fluctuations and correlate these records to existing inner continental shelf records; 4. Examine relationships among WAIS variability, Earth’s orbital configuration, oceanic temperature and circulation, and atmospheric pCO2; and 5. Constrain the timing of Ross Sea continental shelf overdeepening and assess its impact on Neogene and Quaternary ice dynamics. Expedition 374 departed from Lyttelton, New Zealand, in January 2018 and returned in March 2018. We recovered 1292.70 m of high-quality core from five sites spanning the early Miocene to late Quaternary. Three sites were cored on the continental shelf (Sites U1521, U1522, and U1523). At Site U1521, we cored a 650 m thick sequence of interbedded diamictite and diatom-rich mudstone penetrating seismic Ross Sea Unconformity 4 (RSU4). The depositional reconstructions of past glacial and open-marine conditions at this site will provide unprecedented insight into environmental change on the Antarctic continental shelf during the late early and middle Miocene. At Site U1522, we cored a discontinuous late Miocene to Pleistocene sequence of glacial and glaciomarine strata from the outer shelf with the primary objective of penetrating and dating RSU3, which is interpreted to reflect the first continental shelf–wide expansion of East and West Antarctic ice streams. Site U1523, located on the outer continental shelf, targeted a sediment drift beneath the westward-flowing Antarctic Slope Current (ASC) to test the hypothesis that changes in ASC vigor regulate ocean heat flux onto the continental shelf and thus ice sheet mass balance. We also cored two sites on the continental rise and slope. At Site U1524, we recovered a Plio–Pleistocene sedimentary sequence from the levee of the Hillary Canyon, one of the largest conduits of Antarctic Bottom Water from the continental shelf to the abyssal ocean. Site U1524 was designed to penetrate into middle Miocene and older strata, but coring was initially interrupted by drifting sea ice that forced us to abandon coring in Hole U1524A at 399.5 m drilling depth below seafloor (DSF). We moved to a nearby alternate site on the continental slope (Site U1525) to core a single hole designed to complement the record at Site U1524. We returned to Site U1524 after the sea ice cleared and cored Hole U1524C with the rotary core barrel system with the intention of reaching the target depth of 1000 m DSF. However, we were forced to terminate Hole U1524C at 441.9 m DSF because of a mechanical failure with the vessel that resulted in termination of all drilling operations and forced us to return to Lyttelton 16 days earlier than scheduled. The loss of 39% of our operational days significantly impacted our ability to achieve all Expedition 374 objectives. In particular, we were not able to recover continuous middle Miocene sequences from the continental rise designed to complement the discontinuous record from continental shelf Site U1521. The mechanical failure also meant we could not recover cores from proposed Site RSCR-19A, which was targeted to obtain a high-fidelity, continuous record of upper Neogene and Quaternary pelagic/hemipelagic sedimentation. Despite our failure to recover a continental shelf-to-rise Miocene transect, records from Sites U1522, U1524, and U1525 and legacy cores from the Antarctic Geological Drilling Project (ANDRILL) can be integrated to develop a shelf-to-rise Plio–Pleistocene transect. 
    more » « less