A new version of the stochastic multiplume Jet Propulsion Laboratory Eddy‐Diffusivity/Mass‐Flux (JPL‐EDMF) parameterization which consistently couples the simplified Khairoutdinov and Kogan (2000),
Dynamical cores used to study the circulation of the atmosphere employ various numerical methods ranging from finite‐volume, spectral element, global spectral, and hybrid methods. In this work, we explore the use of Flux‐Differencing Discontinuous Galerkin (FDDG) methods to simulate a fully compressible dry atmosphere at various resolutions. We show that the method offers a judicious compromise between high‐order accuracy and stability for large‐eddy simulations and simulations of the atmospheric general circulation. In particular, filters, divergence damping, diffusion, hyperdiffusion, or sponge‐layers are not required to ensure stability; only the numerical dissipation naturally afforded by FDDG is necessary. We apply the method to the simulation of dry convection in an atmospheric boundary layer and in a global atmospheric dynamical core in the standard benchmark of Held and Suarez (1994,
- NSF-PAR ID:
- 10409638
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 15
- Issue:
- 4
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract https://doi.org/10.1175/1520-0493(2000)128<0229:ANCPPI>2.0.CO;2 , warm phase cloud microphysical parameterization with the parameterization of cloud macrophysical and subgrid scale dynamical processes is described. The new parameterization combines the EDMF approach with an assumed shape of a joint probability density function of thermodynamic and kinematic variables which provide the basis for the computation of all parameterized processes. As far as we are aware this is the first attempt to consistently couple all of these parameterized processes in the EDMF framework. This paper is part one of a two paper series. Here, the JPL‐EDMF parameterization is described and benchmark simulations of precipitating stratocumulus and cumulus convection are performed in a single‐column‐model framework. The parameterization results compare favorably to the reference large‐eddy‐simulation results. In the second part (Smalley et al., 2022,https://doi.org/10.1029/2021MS002729 ) the JPL‐EDMF parameterization is validated for a wide range of observation‐based scenarios covering the continuous transition from subtropical stratocumulus to cumulus convection derived from global reanalysis, and parameterization uncertainties are studied in detail. -
more » « less
Abstract The polar
F region ionosphere frequently exhibits sporadic variability (e.g., Meek, 1949,https://doi.org/10.1029/JZ054i004p00339 ; Hill, 1963,https://doi.org/10.1175/1520‐0469(1963)020<0492:SEOLII>2.0.CO;2 ). Recent satellite data analysis (Noja et al., 2013,https://doi.org/10.1002/rds.20033 ; Chartier et al., 2018,https://doi.org/10.1002/2017JA024811 ) showed that the high‐latitudeF region ionosphere exhibits sporadic enhancements more frequently in January than in July in both the northern and southern hemispheres. The same pattern has been seen in statistics of the degradation and total loss of GPS service onboard low‐Earth orbit satellites (Xiong et al. 2018,https://doi.org/10.5194/angeo‐36‐679‐2018 ). Here, we confirm the existence of this annual pattern using ground GPS‐based images of TEC from the MIDAS algorithm. Images covering January and July 2014 confirm that the high‐latitude (>70 MLAT)F region exhibits a substantially larger range of values in January than in July in both the northern and southern hemispheres. The range of TEC values observed in the polar caps is 38–57 TECU (north‐south) in January versus 25–37 TECU in July. First‐principle modeling using SAMI3 reproduces this pattern, and indicates that it is caused by an asymmetry in plasma levels (30% higher in January than in July across both polar caps), as well as 17% longer O+plasma lifetimes in northern hemisphere winter, compared to southern hemisphere winter. -
more » « less
Abstract This paper describes an implementation of the combined hybrid‐parallel prediction (CHyPP) approach of Wikner et al. (2020),
https://doi.org/10.1063/5.0005541 on a low‐resolution atmospheric global circulation model (AGCM). The CHyPP approach combines a physics‐based numerical model of a dynamical system (e.g., the atmosphere) with a computationally efficient type of machine learning (ML) called reservoir computing to construct a hybrid model. This hybrid atmospheric model produces more accurate forecasts of most atmospheric state variables than the host AGCM for the first 7–8 forecast days, and for even longer times for the temperature and humidity near the earth's surface. It also produces more accurate forecasts than a model based only on ML, or a model that combines linear regression, rather than ML, with the AGCM. The potential of the CHyPP approach for climate research is demonstrated by a 10‐year long hybrid model simulation of the atmospheric general circulation, which shows that the hybrid model can simulate the general circulation with substantially smaller systematic errors and more realistic variability than the host AGCM. -
more » « less
Abstract Long‐term efforts have sought to extend global model resolution to smaller scales enabling more accurate descriptions of gravity wave (GW) sources and responses, given their major roles in coupling and variability throughout the atmosphere. Such studies reveal significant improvements accompanying increasing resolution, but no guidance on what is sufficient to approximate reality. We take the opposite approach, using a finite‐volume model solving the Navier‐Stokes equations exactly. The reference simulation addresses mountain wave (MW) generation and responses over the Southern Andes described using isotropic 500 m, central resolution by Fritts et al. (2021),
https://doi.org/10.1175/JAS-D-20-0207.1 and Lund et al. (2020),https://doi.org/10.1175/JAS-D-19-0356.1 . Reductions of horizontal resolution to 1 and 2 km result in (a) systematic increases in initial MW breaking altitudes, (b) weaker, larger‐scale generation of secondary GWs and acoustic waves accompanying these dynamics, and (c) significantly weaker and less extended responses in the mesosphere in latitude and longitude. Horizontal resolution of 4 km largely suppresses instabilities, but allows weak, sustained mean‐flow interactions. Responses for 8 km resolution are very weak and fail to capture any aspects of the high‐resolution responses. The chosen mean winds allow efficient MW penetration into the mesosphere and lower thermosphere, hence only exhibit strong pseudo‐momentum deposition and mean wind decelerations at higher altitudes. A companion paper by Fritts et al. (2022),https://doi.org/10.1029/2021JD036035 explores the impacts of decreasing resolution on responses in the thermosphere. -
more » « less
Abstract Analysis of 40 years of tide gauge data and reanalysis wind stresses from the Middle Atlantic Bight (MAB) indicate that along‐shelf wind stresses are a dominant driver of coastal dynamic sea level (sea level plus atmospheric pressure) variability at daily to yearly time scales. The sea‐level response to along‐shelf wind stress varies substantially along the coast and is accurately reproduced by a steady, barotropic, depth‐averaged model (Csanady, 1978,
https://doi.org/10.1175/1520‐0485(1978)008<0047:tatw>2.0.co;2 , Arrested Topographic Wave). The model indicates that the sea‐level response in the MAB depends primarily on the along‐shelf distribution of the along‐shelf wind stress, the Coriolis frequency, the bottom drag coefficient, and the cross‐shelf bottom slope. The along‐shelf wind stress varies along the MAB shelf due primarily to changes in the shelf orientation. The sea‐level response depends on both the local and upstream (in the sense of Kelvin wave propagation) along‐shelf wind stresses. Consequently, sea‐level variability at daily, monthly and yearly time scales along much of the central MAB coast is more strongly driven by upstream winds along the southern New England shelf than by local winds along the central MAB shelf. The residual coastal sea‐level variability, after removing the wind‐driven response and the trend, is roughly uniform along the MAB coast. The along‐coast average of the residual sea level at monthly and yearly time scales is caused by variations in shelf water densities primarily associated with the large annual cycle in water temperature and interannual variations in salinity.
