skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Parameterizing the Impact of Unresolved Temperature Variability on the Large‐Scale Density Field: Part 1. Theory.
Abstract Unresolved temperature and salinity fluctuations interact with a nonlinear seawater equation of state to produce significant errors in the ocean model evaluation of the large‐scale density field. It is shown that the impact of temperature fluctuations dominates the impact of salinity fluctuations and that the error in density is, to leading order, proportional to the product of a subgrid‐scale temperature variance and a second derivative of the equation of state. Two parameterizations are proposed to correct the large‐scale density field: one deterministic and one stochastic. Free parameters in both parameterizations are fit using fine‐resolution model data. Both parameterizations are computationally efficient as they require only one additional evaluation of a nonlinear equation at each grid cell. A companion paper will discuss the climate impacts of the parameterizations proposed here.  more » « less
Award ID(s):
1736708
PAR ID:
10363583
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Volume:
12
Issue:
12
ISSN:
1942-2466
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The Effective Field Theory of Large Scale Structure (EFTofLSS) has found tremendous success as a perturbative framework for the evolution of large scale structure, and it is now routinely used to compare theoretical predictions against cosmological observations. The model for the total matter field includes one nuisance parameter at 1-loop order, the effective sound speed, which can be extracted by matching the EFT to full N-body simulations. In this work we first leverage the Layzer-Irvine cosmic energy equation to show that the equation of state can be exactly computed with knowledge of the fully nonlinear power spectrum. When augmented with separate universe methods, we show one can estimate the effective sound speed. This estimate is in good agreement with simulation results, with errors at the few tens of percent level. We apply our method to investigate the cosmology dependence of the effective sound speed and to shed light on what cosmic structures shape its value. 
    more » « less
  2. Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed. 
    more » « less
  3. Abstract Improved modeling of permafrost active layer freeze‐thaw plays a crucial role in understanding the response of the Arctic ecosystem to the accelerating warming trend in the region over the past decades. However, modeling the dynamics of the active layer at diurnal time scale remains challenging using the traditional models of freeze‐thaw processes. In this study, a physically based analytical model is formulated to simulate the thaw depth of the active layer under changing boundary conditions of soil heat flux. Conservation of energy for the active layer leads to a nonlinear integral equation of the thaw depth using a temperature profile approximated from the analytical solution of the heat transfer equation forced by ground heat flux. Temporally variable ground heat flux is estimated using non‐gradient models when field observations are not available. Validation of the proposed model conducted against field data obtained from three Arctic forest and tundra sites demonstrates that the model is able to simulate both thaw depth and soil temperature profiles accurately. The model has the potential to estimate regional variability of the thaw depth for permafrost related applications. 
    more » « less
  4. Ocean mesoscale eddies are often poorly represented in climate models, and therefore, their effects on the large scale circulation must be parameterized. Traditional parameterizations, which represent the bulk effect of the unresolved eddies, can be improved with new subgrid models learned directly from data. Zanna and Bolton (ZB20) applied an equation‐discovery algorithm to reveal an interpretable expression parameterizing the subgrid momentum fluxes by mesoscale eddies through the components of the velocity‐gradient tensor. In this work, we implement the ZB20 parameterization into the primitive‐equation GFDL MOM6 ocean model and test it in two idealized configurations with significantly different dynamical regimes and topography. The original parameterization was found to generate excessive numerical noise near the grid scale. We propose two filtering approaches to avoid the numerical issues and additionally enhance the strength of large‐scale energy backscatter. The filtered ZB20 parameterizations led to improved climatological mean state and energy distributions, compared to the current state‐of‐the‐art energy backscatter parameterizations. The filtered ZB20 parameterizations are scale‐aware and, consequently, can be used with a single value of the non‐dimensional scaling coefficient for a range of resolutions. The successful application of the filtered ZB20 parameterizations to parameterize mesoscale eddies in two idealized configurations offers a promising opportunity to reduce long‐standing biases in global ocean simulations in future studies. 
    more » « less
  5. Abstract Parameterization of mesoscale eddies in coarse resolution ocean models is necessary to include the effect of eddies on the large‐scale oceanic circulation. We propose to use a multiple‐scale Quasi‐Geostrophic (MSQG) model to capture the eddy dynamics that develop in response to a prescribed large‐scale flow. The MSQG model consists in extending the traditional quasi geostrophic (QG) dynamics to include the effects of a variable Coriolis parameter and variable background stratification. Solutions to this MSQG equation are computed numerically and compared to a full primitive equation model. The large‐scale flow field permits baroclinically unstable QG waves to grow. These instabilities saturate due to non‐linearities and a filtering method is applied to remove large‐scale structures that develop due to the upscale cascade. The resulting eddy field represents a dynamically consistent response to the prescribed background flow, and can be used to rectify the large‐scale dynamics. Comparisons between Gent‐McWilliams eddy parameterization and the present solutions show large regions of agreement, while also indicating areas where the eddies feed back onto the large scale in a manner that the Gent‐McWilliams parameterization cannot capture. Also of interest is the time variability of the eddy feedback which can be used to build stochastic eddy parameterizations. 
    more » « less