skip to main content


Title: Coupling Warm Rain With an Eddy Diffusivity/Mass Flux Parameterization: 1. Model Description and Validation
Abstract

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),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
NSF-PAR ID:
10371063
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
14
Issue:
8
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The polarFregion 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‐latitudeFregion 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)Fregion 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
  2. Abstract

    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,https://doi.org/10.1175/1520-0477(1994)075〈1825:apftio〉2.0.co;2).

     
    more » « less
  3. Abstract

    Canopy‐snow unloading is the complex physical process of snow unloading from the canopy through meltwater drip, sublimation to the atmosphere, or solid snow unloading to the snowpack below. This process is difficult to parameterize due to limited observations. Time‐lapse photographs of snow in the canopy were characterized by citizen scientists to create a data set of snow interception observations at multiple locations across the western United States. This novel interception data set was used to evaluate three snow unloading parameterizations in the Structure for Unifying Multiple Modeling Alternatives (SUMMA) modular hydrologic modeling framework. SUMMA was modified to include a third snow unloading parameterization, termed Wind‐Temperature (Roesch et al., 2001,https://doi.org/10.1007/s003820100153), which includes wind‐dependent and temperature‐dependent unloading functions. It was compared to a meltwater drip unloading parameterization, termed Melt (Andreadis et al., 2009,https://doi.org/10.1029/2008wr007042), and a time‐dependent unloading parameterization, termed Exponential‐Decay (Hedstrom & Pomeroy, 1998,https://doi.org/10.1002/(SICI)1099-1085(199808/09)12:10/11<1611::AID-HYP684>3.0.CO;2-4). Wind‐Temperature performed well without calibration across sites, specifically in cold climates, where wind dominates unloading and rime accretion is low. At rime prone sites, Wind‐Temperature should be calibrated to account for longer interception events with less sensitivity to wind, otherwise Melt can be used without calibration. The absence of model physics in Exponential‐Decay requires local calibration that can only be transferred to sites with similar unloading patterns. The choice of unloading parameterization can result in 20% difference in SWE on the ground below the canopy and 10% difference in estimated average winter canopy albedo. These novel observations shed light on processes that are often overlooked in hydrology.

     
    more » « less
  4. Abstract

    Wave breaking induced bubbles contribute a significant part of air‐sea gas fluxes. Recent modeling of the sea state dependent CO2flux found that bubbles contribute up to ∼40% of the total CO2air‐sea fluxes (Reichl & Deike, 2020,https://doi.org/10.1029/2020gl087267). In this study, we implement the sea state dependent bubble gas transfer formulation of Deike and Melville (2018,https://doi.org/10.1029/2018gl078758) into a spectral wave model (WAVEWATCH III) incorporating the spectral modeling of the wave breaking distribution from Romero (2019,https://doi.org/10.1029/2019gl083408). We evaluate the accuracy of the sea state dependent gas transfer parameterization against available measurements of CO2gas transfer velocity from 9 data sets (11 research cruises, see Yang et al. (2022,https://doi.org/10.3389/fmars.2022.826421)). The sea state dependent parameterization for CO2gas transfer velocity is consistent with observations, while the traditional wind‐only parameterization used in most global models slightly underestimates the observations of gas transfer velocity. We produce a climatology of the sea state dependent gas transfer velocity using reanalysis wind and wave data spanning 1980–2017. The climatology shows that the enhanced gas transfer velocity occurs frequently in regions with developed sea states (with strong wave breaking and high significant wave height). The present study provides a general sea state dependent parameterization for gas transfer, which can be implemented in global coupled models.

     
    more » « less
  5. 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.1and 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/2021JD036035explores the impacts of decreasing resolution on responses in the thermosphere.

     
    more » « less