scholarly journals The effect of metallicity on the atmospheres of exoplanets with fully coupled 3D hydrodynamics, equilibrium chemistry, and radiative transfer

2018 ◽  
Vol 612 ◽  
pp. A105 ◽  
Author(s):  
B. Drummond ◽  
N. J. Mayne ◽  
I. Baraffe ◽  
P. Tremblin ◽  
J. Manners ◽  
...  

In this work, we have performed a series of simulations of the atmosphere of GJ 1214b assuming different metallicities using the Met Office Unified Model (UM). The UM is a general circulation model (GCM) that solves the deep, non-hydrostatic equations of motion and uses a flexible and accurate radiative transfer scheme, based on the two-stream and correlated-k approximations, to calculate the heating rates. In this work we consistently couple a well-tested Gibbs energy minimisation scheme to solve for the chemical equilibrium abundances locally in each grid cell for a general set of elemental abundances, further improving the flexibility and accuracy of the model. As the metallicity of the atmosphere is increased we find significant changes in the dynamical and thermal structure, with subsequent implications for the simulated phase curve. The trends that we find are qualitatively consistent with previous works, though with quantitative differences. We investigate in detail the effect of increasing the metallicity by splitting the mechanism into constituents, involving the mean molecular weight, the heat capacity and the opacities. We find the opacity effect to be the dominant mechanism in altering the circulation and thermal structure. This result highlights the importance of accurately computing the opacities and radiative transfer in 3D GCMs.

2014 ◽  
Vol 7 (3) ◽  
pp. 883-907 ◽  
Author(s):  
R. Fischer ◽  
S. Nowicki ◽  
M. Kelley ◽  
G. A. Schmidt

Abstract. The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere–ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.


Icarus ◽  
2019 ◽  
Vol 321 ◽  
pp. 232-250 ◽  
Author(s):  
Masaru Yamamoto ◽  
Kohei Ikeda ◽  
Masaaki Takahashi ◽  
Takeshi Horinouchi

2019 ◽  
Vol 12 (8) ◽  
pp. 3745-3758 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  

Abstract. The ice thickness distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to discretize it remains an open question. Here, we use the ocean–sea ice general circulation model, Nucleus for European Modelling of the Ocean (NEMO) version 3.6 and Louvain-la-Neuve sea Ice Model (LIM) version 3 (NEMO3.6-LIM3), forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ∼30 % and ∼10 % of the reference values (run with five categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ∼4 and ∼2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean sea ice state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence. It is also found that the current default discretization of the NEMO3.6-LIM3 model is sufficient for large-scale present-day climate applications. In all cases, the role of the ITD discretization on the simulated mean sea ice state has to be appreciated relative to other influences (parameter uncertainty, forcing uncertainty, internal climate variability).


2008 ◽  
Vol 21 (22) ◽  
pp. 5797-5806 ◽  
Author(s):  
Paul A. O’Gorman ◽  
Tapio Schneider

Abstract As the climate changes, changes in static stability, meridional temperature gradients, and availability of moisture for latent heat release may exert competing effects on the energy of midlatitude transient eddies. This paper examines how the eddy kinetic energy in midlatitude baroclinic zones responds to changes in radiative forcing in simulations with an idealized moist general circulation model. In a series of simulations in which the optical thickness of the longwave absorber is varied over a wide range, the eddy kinetic energy has a maximum for a climate with mean temperature similar to that of present-day earth, with significantly smaller values both for warmer and for colder climates. In a series of simulations in which the meridional insolation gradient is varied, the eddy kinetic energy increases monotonically with insolation gradient. In both series of simulations, the eddy kinetic energy scales approximately linearly with the dry mean available potential energy averaged over the baroclinic zones. Changes in eddy kinetic energy can therefore be related to the changes in the atmospheric thermal structure that affect the mean available potential energy.


2018 ◽  
Vol 75 (1) ◽  
pp. 3-20 ◽  
Author(s):  
Nicholas J. Lutsko

An equatorial heat source mimicking the strong diabatic heating above the west Pacific is added to an idealized, dry general circulation model. For small (<0.5 K day−1) heating rates the responses closely match the expectations from linear Matsuno–Gill theory, though the amplitudes of the responses increase sublinearly. This “linear” regime breaks down for larger heating rates and it is found that this is because the stability of the tropical atmosphere increases. At the same time, the equatorial winds increasingly superrotate. This superrotation is driven by stationary eddy momentum fluxes by the waves excited by the heating and is damped by the vertical advection of low-momentum air by the mean flow and, at large heating rates, by the divergence of momentum by transient eddies. These dynamics are explored in additional experiments in which the equator-to-pole temperature gradient is varied. Very strong superrotation is produced when a large heating rate is applied to a setup with a relatively weak equator-to-pole temperature gradient, though there is no evidence that this is a case of “runaway” superrotation.


2009 ◽  
Vol 66 (11) ◽  
pp. 3277-3296 ◽  
Author(s):  
James J. Benedict ◽  
David A. Randall

Abstract The detailed dynamic and thermodynamic space–time structures of the Madden–Julian oscillation (MJO) as simulated by the superparameterized Community Atmosphere Model version 3.0 (SP-CAM) are analyzed. Superparameterization involves substituting conventional boundary layer, moist convection, and cloud parameterizations with a configuration of cloud-resolving models (CRMs) embedded in each general circulation model (GCM) grid cell. Unlike most GCMs that implement conventional parameterizations, the SP-CAM displays robust atmospheric variability on intraseasonal space and time (30–60 days) scales. The authors examine a 19-yr SP-CAM simulation based on the Atmospheric Model Intercomparison Project protocol, forced by prescribed sea surface temperatures. Overall, the space–time structures of MJO convective disturbances are very well represented in the SP-CAM. Compared to observations, the model produces a similar vertical progression of increased moisture, warmth, and heating from the boundary layer to the upper troposphere as deep convection matures. Additionally, important advective and convective processes in the SP-CAM compare favorably with those in observations. A deficiency of the SP-CAM is that simulated convective intensity organized on intraseasonal space–time scales is overestimated, particularly in the west Pacific. These simulated convective biases are likely due to several factors including unrealistic boundary layer interactions, a lack of weakening of the simulated disturbance over the Maritime Continent, and mean state differences.


2021 ◽  
Author(s):  
Alexei Belochitski ◽  
Vladimir Krasnopolsky

Abstract. The ability of Machine-Learning (ML) based model components to generalize to the previously unseen inputs, and the resulting stability of the models that use these components, has been receiving a lot of recent attention, especially when it comes to ML-based parameterizations. At the same time, ML-based emulators of existing parameterizations can be stable, accurate, and fast when used in the model they were specifically designed for. In this work we show that shallow-neural-network-based emulators of radiative transfer parameterizations developed almost a decade ago for a state-of-the-art GCM are robust with respect to the substantial structural and parametric change in the host model: when used in two seven month-long experiments with the new model, they not only remain stable, but generate realistic output. Aspects of neural network architecture and training set design potentially contributing to stability of ML-based model components are discussed.


1990 ◽  
Vol 13 ◽  
pp. 57-67
Author(s):  
Julio Buchmann ◽  
Lawrence E. Buja ◽  
Jan Paegle ◽  
Robert E. Dickinson

A series of real data experiments is performed with a general circulation model in order to ascertain the sensitivity of extended range rain forecasts over the Americas to the structure and magnitude of tropical heating anomalies. The emphasis is upon heat inputs over the tropical Atlantic which have shown particularly significant drying influences over North America in our prior simulations. The heating imposed in the prior experiments is compared to the condensation heating rates that naturally occur in the forecast model, and shown to be excessive by approximately a factor of two. Present experiments reduce the imposed anomaly by a factor of three, and also incorporate sea-surface temperature decreases over the eastern tropical Pacific Ocean. The new experimental results are in many ways consistent with our prior results. The dry North American response is statistically more significant than the South American response, and occurs at least as frequently in the different members of the experimental ensembles as in our prior experiments. The drying effect is accentuated by the presence of East Pacific cooling, but this does not appear to be the dominant influence. Over tropical South America, the Pacific and Atlantic modifications produce compensating influences, with the former dominating dominant, and allowing increased rainfall over the Amazon Basin.


Sign in / Sign up

Export Citation Format

Share Document