Search NASA⌕ Search

SEARCH · Search NASA

Results for “Convective parameterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

On the Vertical Distribution of Local and Remote Sources of Water for Precipitation

The vertical distribution of local and remote sources of water for precipitation and total column water over the United States are evaluated in a general circulation model simulation. The Goddard Earth Observing System (GEOS) general circulation model (GCM) includes passive constituent tracers to determine the geographical sources of the water in the column. Results show that the local percentage of precipitable water and local percentage of precipitation can be very different. The transport of water vapor from remote oceanic sources at mid and upper levels is important to the total water in the column over the central United States, while the access of locally evaporated water in convective precipitation processes is important to the local precipitation ratio. This result resembles the conceptual formulation of the convective parameterization. However, the formulations of simple models of precipitation recycling include the assumption that the ratio of the local water in the column is equal to the ratio of the local precipitation. The present results demonstrate the uncertainty in that assumption, as locally evaporated water is more concentrated near the surface.

Bosilovich, Michael G.↗

Using Satellite-derived Ice Concentration to Represent Antarctic Coastal Polynyas in Ocean Climate Models

The focus of this paper is on the representation of Antarctic coastal polynyas in global ice-ocean general circulation models (OGCMs), in particular their local, regional, and high-frequency behavior. This is verified with the aid of daily ice concentration derived from satellite passive microwave data using the NASATeam 2 (NT2) and the bootstrap (BS) algorithms. Large systematic regional and temporal discrepancies arise, some of which are related to the type of convection parameterization used in the model. An attempt is made to improve the fresh-water flux associated with melting and freezing in Antarctic coastal polynyas by ingesting (assimilating) satellite ice concentration where it comes to determining the thermodynamics of the open-water fraction of a model grid cell. Since the NT2 coastal open-water fraction (polynyas) tends to be less extensive than the simulated one in the decisive season and region, assimilating NT2 coastal ice concentration yields overall reduced net freezing rates, smaller formation rates of Antarctic Bottom Water, and a stronger southward flow of North Atlantic Deep Water across 30 S. Enhanced net freezing rates occur regionally when NT2 coastal ice concentration is assimilated, concomitant with a more realistic ice thickness distribution and accumulation of High-Salinity Shelf Water. Assimilating BS rather than NT2 coastal ice concentration, the differences to the non-assimilated simulation are generally smaller and of opposite sign. This suggests that the model reproduces coastal ice concentration in closer agreement with the BS data than with the NT2 data, while more realistic features emerge when NT2 data are assimilated.

Stoessel, Achim↗

Simulations of Hurricane Katrina (2005) with the 0.125 degree finite-volume General Circulation Model on the NASA Columbia Supercomputer

Hurricane Katrina was the sixth most intense hurricane in the Atlantic. Katrina's forecast poses major challenges, the most important of which is its rapid intensification. Hurricane intensity forecast with General Circulation Models (GCMs) is difficult because of their coarse resolution. In this article, six 5-day simulations with the ultra-high resolution finite-volume GCM are conducted on the NASA Columbia supercomputer to show the effects of increased resolution on the intensity predictions of Katrina. It is found that the 0.125 degree runs give comparable tracks to the 0.25 degree, but provide better intensity forecasts, bringing the center pressure much closer to observations with differences of only plus or minus 12 hPa. In the runs initialized at 1200 UTC 25 AUG, the 0.125 degree simulates a more realistic intensification rate and better near-eye wind distributions. Moreover, the first global 0.125 degree simulation without convection parameterization (CP) produces even better intensity evolution and near-eye winds than the control run with CP.

Shen, B.-W.↗

Hurricane Forecasts with a Global Mesoscale-resolving Model on the NASA Columbia Supercomputer Preliminary Simulations of Hurricane Katrina (2005)

It is known that the General Circulation Models (GCMs) have sufficient resolution to accurately simulate hurricane near-eye structure and intensity. To overcome this limitation, the mesoscale-resolving finite-element GCM (fvGCM) has been experimentally deployed on the NASA Columbia supercomputer, and its performance is evaluated choosing hurricane Katrina as an example in this study. On late August 2005 Katrina underwent two stages of rapid intensification and became the sixth most intense hurricane in the Atlantic. Six 5-day simulations of Katrina at both 0.25 deg and 0.125 deg show comparable track forecasts, but the 0,125 deg runs provide much better intensity forecasts, producing center pressure with errors of only +/- 12 hPa. The 0.125 deg simulates better near-eye wind distributions and a more realistic average intensification rate. A convection parameterization (CP) is one of the major limitations in a GCM, the 0.125 deg run with CP disabled produces very encouraging results.

Shen, B.-W.↗

Clouds in GEOS-5

The GEOS-5 atmospheric model is being developed as a weather-and-climate capable model. It must perform well in assimilation mode as well as in weather and climate simulations and forecasts and in coupled chemistry-climate simulations. In developing GEOS-5, attention has focused on the representation of moist processes. The moist physics package uses a single phase prognostic condensate and a prognostic cloud fraction. Two separate cloud types are distinguished by their source: "anvil" cloud originates in detraining convection, and large-scale cloud originates in a PDF-based condensation calculation. Ice and liquid phases for each cloud type are considered. Once created, condensate and fraction from the anvil and statistical cloud types experience the same loss processes: evaporation of condensate and fraction, auto-conversion of liquid or mixed phase condensate, sedimentation of frozen condensate, and accretion of condensate by falling precipitation. The convective parameterization scheme is the Relaxed Arakawa-Schubert, or RAS, scheme. Satellite data are used to evaluate the performance of the moist physics packages and help in their tuning. In addition, analysis of and comparisons to cloud-resolving models such as the Goddard Cumulus Ensemble model are used to help improve the PDFs used in the moist physics. The presentation will show some of our evaluations including precipitation diagnostics.

Bacmeister, Julio↗

The Effect of Environmental Conditions on Tropical Deep Convective Systems Observed from the TRMM Satellite

This study uses measurements of radiation and cloud properties taken between January and August 1998 by three Tropical Rainfall Measuring Mission (TRMM) instruments, the Clouds and the Earth's Radiant Energy System (CERES) scanner, the TRMM Microwave Imager (TMI), and the Visible and InfraRed Scanner (VIRS), to evaluate the variations of tropical deep convective systems (DCS) with sea surface temperature (SST) and precipitation. This study finds that DCS precipitation efficiency increases with SST at a rate of approx. 2%/K. Despite increasing rainfall efficiency, the cloud areal coverage rises with SST at a rate of about 7%/K in the warm tropical seas. There, the boundary layer moisture supply for deep convection and the moisture transported to the upper troposphere for cirrus-anvil cloud formation increase by approx. 6.3%/K and approx. 4.0%/K, respectively. The changes in cloud formation efficiency, along with the increased transport of moisture available for cloud formation, likely contribute to the large rate of increasing DCS areal coverage. Although no direct observations are available, the increase of cloud formation efficiency with rising SST is deduced indirectly from measurements of changes in the ratio of DCS ice water path and boundary layer water vapor amount with SST. Besides the cloud areal coverage, DCS cluster effective sizes also increase with precipitation. Furthermore, other cloud properties, such as cloud total water and ice water paths, increase with SST. These changes in DCS properties will produce a negative radiative feedback for the earth's climate system due to strong reflection of shortwave radiation by the DCS. These results significantly differ from some previous hypothesized dehydration scenarios for warmer climates, and have great potential in testing current cloud-system resolving models and convective parameterizations of general circulation models.

Lin, Bing↗

Comparison of the Effects of RAS vs. Kain-Fritsch Convective Schemes on Katrina Forecasts with GEOS-5

Global forecasts were made with the 0.25-degree latitude version of GEOS-5, with the RAS scheme and with the Kain-Fritsch scheme. Examination was made of the Katrina (2005) hurricane simulation. Replacement of the RAS convective scheme with the K-F scheme results in a much more vigorous Katrina, closer to reality. Still, the result is not as vigorous as reality. In terms of wind maximum, the gap was closed by ~50%. The result seems to be due to the RAS scheme drying out the boundary layer, thus hampering the grid-scale secondary circulation and attending cyclone development. The RAS case never developed a full warm core, whereas the K-F case did. Not shown here: The K-F scheme also resulted in a more vigorous storm than when GEOS-5 is run with no convective parameterization. Also not shown: An experiment in which the RAS firing level was moved up by 3 model levels resulted in a stronger, warm-core storm, though not as strong as the K-F case. Effects on storm track were noticed, but not studied.

Miller, Timothy L.↗

Classification of Clouds and Deep Convection from GEOS-5 Using Satellite Observations

With the increased resolution of global atmospheric models and the push toward global cloud resolving models, the resemblance of model output to satellite observations has become strikingly similar. As we progress with our adaptation of the Goddard Earth Observing System Model, Version 5 (GEOS-5) as a high resolution cloud system resolving model, evaluation of cloud properties and deep convection require in-depth analysis beyond a visual comparison. Outgoing long-wave radiation (OLR) provides a sufficient comparison with infrared (IR) satellite imagery to isolate areas of deep convection. We have adopted a binning technique to generate a series of histograms for OLR which classify the presence and fraction of clear sky versus deep convection in the tropics that can be compared with a similar analyses of IR imagery from composite Geostationary Operational Environmental Satellite (GOES) observations. We will present initial results that have been used to evaluate the amount of deep convective parameterization required within the model as we move toward cloud system resolving resolutions of 10- to 1-km globally.

Putman, William↗

Tropical and Subtropical Cloud Transitions in Weather and Climate Prediction Models: The GCSS/WGNE Pacific Cross-Section Intercomparison (GPCI)

A model evaluation approach is proposed in which weather and climate prediction models are analyzed along a Pacific Ocean cross section, from the stratocumulus regions off the coast of California, across the shallow convection dominated trade winds, to the deep convection regions of the ITCZ the Global Energy and Water Cycle Experiment Cloud System Study/Working Group on Numerical Experimentation (GCSS/ WGNE) Pacific Cross-Section Intercomparison (GPCI). The main goal of GPCI is to evaluate and help understand and improve the representation of tropical and subtropical cloud processes in weather and climate prediction models. In this paper, a detailed analysis of cloud regime transitions along the cross section from the subtropics to the tropics for the season June July August of 1998 is presented. This GPCI study confirms many of the typical weather and climate prediction model problems in the representation of clouds: underestimation of clouds in the stratocumulus regime by most models with the corresponding consequences in terms of shortwave radiation biases; overestimation of clouds by the 40-yr ECMWF Re-Analysis (ERA-40) in the deep tropics (in particular) with the corresponding impact in the outgoing longwave radiation; large spread between the different models in terms of cloud cover, liquid water path and shortwave radiation; significant differences between the models in terms of vertical cross sections of cloud properties (in particular), vertical velocity, and relative humidity. An alternative analysis of cloud cover mean statistics is proposed where sharp gradients in cloud cover along the GPCI transect are taken into account. This analysis shows that the negative cloud bias of some models and ERA-40 in the stratocumulus regions [as compared to the first International Satellite Cloud Climatology Project (ISCCP)] is associated not only with lower values of cloud cover in these regimes, but also with a stratocumulus-to-cumulus transition that occurs too early along the trade wind Lagrangian trajectory. Histograms of cloud cover along the cross section differ significantly between models. Some models exhibit a quasi-bimodal structure with cloud cover being either very large (close to 100%) or very small, while other models show a more continuous transition. The ISCCP observations suggest that reality is in-between these two extreme examples. These different patterns reflect the diverse nature of the cloud, boundary layer, and convection parameterizations in the participating weather and climate prediction models.

Teixeira, J.↗

The Tethered Moon

A reasonable initial condition on Earth after the Moonforming impact is that it begins as a hot global magma ocean1,2. We therefore begin our study with the mantle as a liquid ocean with a surface temperature on the order of 3000- 4000 K at a time some 100-1000 years after the impact, by which point we can hope that early transients have settled down. A 2nd initial condition is a substantial atmosphere, 100-1000 bars of H2O and CO2, supplemented by smaller amounts of CO, H2, N2, various sulfur-containing gases, and a suite of geochemical volatiles evaporated from the magma. Third, we start the Moon with its current mass at the relevant Roche limit. The 4th initial condition is the angular momentum of the Earth-Moon system. Canonical models hold this constant, whilst some recent models begin with considerably more angular momentum than is present today. Here we present a ruthlessly simplified model of Earth's cooling magmasphere based on a full-featured atmosphere and including tidal heating by the newborn Moon. Thermal blanketing by H2O-CO2 atmospheres slows cooling of a magma ocean. Geochemical volatiles - chiefly S, Na, and Cl - raise the opacity of the magma ocean's atmosphere and slow cooling still more. We assume a uniform mantle with a single internal (potential) temperature and a global viscosity. The important "freezing point" is the sharp rheological transition between a fluid carrying suspended crystals and a solid matrix through which fluids percolate. Most tidal heating takes place at this "freezing point" in a gel that is both pliable and viscous. Parameterized convection links the cooling rate to the temperature and heat generation inside the Earth. Tidal heating is a major effect. Tidal dissipation in the magma ocean is described by viscosity. The Moon is entwined with Earth by the negative feedback between thermal blanketing and tidal heating that comes from the temperature-dependent viscosity of the magma ocean. Because of this feedback, the rate that the Moon's orbit evolves is limited by the modest radiative cooling rate of Earth's atmosphere, which in effect tethers the Moon to the Earth. Consequently the Moon's orbit evolves orders of magnitude more slowly than in conventional models. Slow orbital evolution promotes capture by orbital resonances that may have been important in the Earth-Moon system

Earth-Moon system↗

Diagnosing Warm Frontal Cloud Formation in a GCM: A Novel Approach Using Conditional Subsetting

This study analyzes characteristics of clouds and vertical motion across extratropical cyclone warm fronts in the NASA Goddard Institute for Space Studies general circulation model. The validity of the modeled clouds is assessed using a combination of satellite observations from CloudSat, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E), and the NASA Modern-Era Retrospective Analysis for Research and Applications (MERRA) reanalysis. The analysis focuses on developing cyclones, to test the model's ability to generate their initial structure. To begin, the extratropical cyclones and their warm fronts are objectively identified and cyclone-local fields are mapped into a vertical transect centered on the surface warm front. To further isolate specific physics, the cyclones are separated using conditional subsetting based on additional cyclone-local variables, and the differences between the subset means are analyzed. Conditional subsets are created based on 1) the transect clouds and 2) vertical motion; 3) the strength of the temperature gradient along the warm front, as well as the storm-local 4) wind speed and 5) precipitable water (PW). The analysis shows that the model does not generate enough frontal cloud, especially at low altitude. The subsetting results reveal that, compared to the observations, the model exhibits a decoupling between cloud formation at high and low altitudes across warm fronts and a weak sensitivity to moisture. These issues are caused in part by the parameterized convection and assumptions in the stratiform cloud scheme that are valid in the subtropics. On the other hand, the model generates proper covariability of low-altitude vertical motion and cloud at the warm front and a joint dependence of cloudiness on wind and PW.

cyclones↗

Assessment of NASA GISS CMIP5 and Post-CMIP5 Simulated Clouds and TOA Radiation Budgets Using Satellite Observations. Part I: Cloud Fraction and Properties

Although many improvements have been made in phase 5 of the Coupled Model Intercomparison Project (CMIP5), clouds remain a significant source of uncertainty in general circulation models (GCMs) because their structural and optical properties are strongly dependent upon interactions between aerosol/cloud microphysics and dynamics that are unresolved in such models. Recent changes to the planetary boundary layer (PBL) turbulence and moist convection parameterizations in the NASA GISS Model E2 atmospheric GCM(post-CMIP5, hereafter P5) have improved cloud simulations significantly compared to its CMIP5 (hereafter C5) predecessor. A study has been performed to evaluate these changes between the P5 and C5 versions of the GCM, both of which used prescribed sea surface temperatures. P5 and C5 simulated cloud fraction (CF), liquid water path (LWP), ice water path (IWP), cloud water path (CWP), precipitable water vapor (PWV), and relative humidity (RH) have been compared to multiple satellite observations including the Clouds and the Earth's Radiant Energy System-Moderate Resolution Imaging Spectroradiometer (CERES-MODIS, hereafter CM), CloudSat- Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO; hereafter CC), Atmospheric Infrared Sounder (AIRS), and Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E). Although some improvements are observed in the P5 simulation on a global scale, large improvements have been found over the southern midlatitudes (SMLs), where correlations increased and both bias and root-mean-square error (RMSE) significantly decreased, in relation to the previous C5 simulation, when compared to observations. Changes to the PBL scheme have resulted in improved total column CFs, particularly over the SMLs where marine boundary layer (MBL) CFs have increased by nearly 20% relative to the previous C5 simulation. Globally, the P5 simulated CWPs are 25 gm22 lower than the previous C5 results. The P5 version of the GCM simulates PWV and RH higher than its C5 counterpart and agrees well with the AMSR-E and AIRS observations. The moister atmospheric conditions simulated by P5 are consistent with the CF comparison and provide a strong support for the increase in MBL clouds over the SMLs. Over the tropics, the P5 version of the GCM simulated total column CFs and CWPs are slightly lower than the previous C5 results, primarily as a result of the shallower tropical boundary layer in P5 relative to C5 in regions outside the marine stratocumulus decks.

General circulation models↗

Application of the CERES Flux-by-Cloud Type Simulator to GCM Output

The CERES Flux By CloudType data product produces CERES top-of-atmosphere (TOA) fluxes by region and cloud type. Here, the cloud types are defined by cloud optical depth (t) and cloud top pressure (pc), with bins similar to those used by ISCCP (International Satellite Cloud Climatology Project). This data product has the potential to be a powerful tool for the evaluation of the clouds produced by climate models by helping to identify which physical parameterizations have problems (e.g., boundary-layer parameterizations, convective clouds, processes that affect surface albedo). Also, when the flux-by-cloud type and frequency of cloud types are simultaneously used to evaluate a model, the results can determine whether an unrealistically large or small occurrence of a given cloud type has an important radiative impact for a given region. A simulator of the flux-by-cloud type product has been applied to three-hourly data from the year 2008 from the UK Met Office HadGEM2-A model using the Langley Fu-Lour radiative transfer model to obtain TOA SW and LW fluxes.

Eitzen, Zachary↗

Large-Scale Transport Responses to Tropospheric Circulation Changes Using GEOS-5

The mean age since air was last at the Northern Hemisphere midlatitude surface is a fundamental property of tropospheric transport. Recent comparisons among chemistry climate models, however, reveal that there are large differences in the mean age among models and that these differences are most likely related to differences in tropical (parameterized) convection. Here we use aquaplanet simulations of the Goddard Earth Observing System Model Version 5 (GEOS-5) to explore the sensitivity of the mean age to changes in the tropical circulation. Tropical circulation changes are forced by prescribed localized off-equatorial warm sea surface temperature anomalies that (qualitatively) reproduce the convection and circulation differences among the comprehensive models. Idealized chemical species subject to prescribed OH loss are also integrated in parallel in order to illustrate the impact of tropical transport changes on interhemispheric constituent transport.

Orbe, Clara↗

Overview of Large-scale Tropospheric Transport in the Chemistry Climate Model Initiative (CCMI) Simulations

The transport of chemicals is a major uncertainty in the modeling of tropospheric composition. Here we compare the large-scale tropospheric transport properties among different models in the Chemistry Climate Modeling Initiative (CCMI) with a focus on transport defined with respect to the Northern Hemisphere (NH) midlatitude surface. Among simulations of the recent past (1980-2010) we show that there are substantial differences in their global-scale tropospheric transport properties. For example, the mean transit time since southern hemisphere air last contacted the NH midlatitude surface differs by ~30-40% among simulations. We show that these differences are most likely associated with differences in parameterized convection over the oceans, such that the spread in transport among simulations constrained with analysis fields is as large as the spread among free-running simulations.

large-scale tropospheric transport↗

Large-Scale Transport into the Arctic: The Roles of the Midlatitude Jet and the Hadley Cell

Transport from the Northern Hemisphere (NH) midlatitudes to the Arctic plays a crucial role in determining the abundance of trace gases and aerosols that are important to Arctic climate via impacts on radiation and chemistry. Here we examine this transport using an idealized tracer with a fixed lifetime and predominantly midlatitude land-based sources in models participating in the Chemistry Climate Model Initiative (CCMI). We show that there is a 25%-45% difference in the Arctic concentrations of this tracer among the models. This spread is correlated with the spread in the location of the Pacific jet, as well as the spread in the location of the Hadley Cell (HC) edge, which varies consistently with jet latitude. Our results suggest that it is likely that the HC-related zonal-mean meridional transport rather than the jet-related eddy mixing is the major contributor to the inter-model spread in the transport of land-based tracers into the Arctic. Specifically, in models with a more northern jet, the HC generally extends further north and the tracer source region is mostly covered by surface southward flow associated with the lower branch of the HC, resulting in less efficient transport poleward to the Arctic. During boreal summer, there are poleward biases in jet location in free-running models, and these models likely underestimate the rate of transport into the Arctic. Models using specified dynamics do not have biases in the jet location, but do have biases in the surface meridional flow, which may result in differences in transport into the Arctic. In addition to the land-based tracer, the midlatitude-to-Arctic transport is further examined by another idealized tracer with zonally uniform sources. With equal sources from both land and ocean, the inter-model spread of this zonally uniform tracer is more related to variations in parameterized convection over oceans rather than variations in HC extent, particularly during boreal winter. This suggests that transport of land-based and oceanic tracers or aerosols towards the Arctic differs in pathways and therefore their corresponding inter-model variabilities result from different physical processes.

chemistry climate models↗

GISS Model E2.2: A Climate Model Optimized for the Middle Atmosphere. Part 2: Validation of Large-Scale Transport and Evaluation of Climate Response

Here we examine the large‐scale transport characteristics of the new ``Middle Atmosphere" NASA Goddard Institute for Space Studies climate model (E2.2). First we evaluate the stratospheric transport circulation in historical atmosphere‐only simulations integrated with interactive trace gas and aerosol chemistry. Compared to lower vertical resolution model versions, E2.2 exhibits improved tropical ascent and older stratospheric mean ages that are more consistent with observed values. In the troposphere, poleward transport to the Arctic and interhemispheric mean ages in E2.2 are comparable to models participating in the Chemistry Climate Modeling Initiative. In addition to validating E2.2 we also assess its ``transport sensitivity" using the coupled atmosphere‐ocean abrupt 4xCO2 and transient 1%CO2 simulations submitted to the Coupled Model Intercomparison Project, Phase 6, along with a 2xCO2 simulation used to evaluate the linearity of the transport circulation's response to increased CO2. We show that decreases (increases) in a stratospheric mean age (idealized surface loss) tracer scale linearly with increased lower stratospheric upwelling, which also increases linearly with warming tropical sea surface temperatures (SSTs). Abrupt 2xCO2 and 4xCO2 experiments constrained with (fixed) pre‐industrial SSTs are also used to quantify the relative importance of rapid adjustments versus SST feedbacks to the transport circulation responses in the model. Finally, sensitivity experiments are presented to illustrate the impact of changes in the convective parameterization on stratospheric transport.

GISS Model E2.2↗

Dynamical and Trace Gas Responses of the Quasi-Biennial Oscillation to Increased CO2

Long-term projections of the Quasi-Biennial Oscillation (QBO) remain uncertain. Using the high-top \Middle Atmosphere" (MA) version of the NASA Goddard Institute for Space Studies (GISS) ModelE (E2.2-AP), we show that increased CO2 reduces the QBO period and weakens its amplitude, consistent with results from recent multimodel studies. The amplitude response is asymmetric, with the easterly phase weakening more than the westerly phase. Results from both coupled atmosphere{ocean and climatologically fixed sea surface temperature simulations link the period response to changes in parameterized convective gravity wave fluxes, whereas the amplitude changes are more closely tied to changes in tropical upwelling. Comparisons of simulations with and without ozone-dynamical feedbacks reveal that ozone feedbacks increase the QBO amplitude and period in pre-industrial control simulations, and reduce the response of the amplitude to increased CO2. The response of stratospheric tracers to projected changes in the QBO is also investigated. Similarity between results using fully interactive (trace gas and aerosol) versus linearized ozone chemical mechanisms suggests that the QBO response to CO2 does not depend on additional chemical feedbacks beyond those associated with ozone, which implies that projected QBO changes may be credibly represented using simpli fied stratospheric ozone schemes.

GISS Model E2.2-AP↗