Search

This paper introduces the Flexible Global Ocean‐Atmosphere‐Land System Model: Grid‐Point Version 3 (FGOALS‐g3) and evaluates its basic performance based on some of its participation in the sixth phase of the Coupled Model Intercomparison Project (CMIP6) experiments. Our results show that many significant improvements have been achieved by FGOALS‐g3 in terms of climatological mean states, variabilities, and long‐term trends. For example, FGOALS‐g3 has a small (−0.015°C/100 yr) climate drift in 700‐yr preindustrial control (piControl) runs and smaller biases in climatological mean variables, such as the land/sea surface temperatures (SSTs) and seasonal soil moisture cycle, compared with its previous version FGOALS‐g2 during the historical period. The characteristics of climate variabilities, for example, Madden‐Julian oscillation (MJO) eastward/westward propagation ratios, spatial patterns of interannual variability of tropical SST anomalies, and relationship between the East Asian Summer Monsoon and El Niño–Southern Oscillation (ENSO), are well captured by FGOALS‐g3. In particular, the cooling trend of globally averaged surface temperature during 1940–1970, which is a challenge for most CMIP3 and CMIP5 models, is well reproduced by FGOALS‐g3 in historical runs. In addition to the external forcing factors recommended by CMIP6, anthropogenic groundwater forcing from 1965 to 2014 was incorporated into the FGOALS‐g3 historical runs.

The objective of this paper is to develop an optimized implicit‐explicit (IMEX) Runge‐Kutta scheme for atmospheric applications focusing on stability and accuracy. Following the common terminology, the proposed method is called IMEX‐SSP2(2,3,2), as it has second‐order accuracy and is composed of diagonally implicit two‐stage and explicit three‐stage parts. This scheme enjoys the Strong Stability Preserving (SSP) property for both parts. This new scheme is applied to nonhydrostatic compressible Boussinesq equations in two different arrangements, including (i) semiimplicit and (ii) Horizontally Explicit‐Vertically Implicit (HEVI) forms. The new scheme preserves the SSP property for larger regions of absolute monotonicity compared to the well‐studied scheme in the same class. In addition, numerical tests confirm that the IMEX‐SSP2(2,3,2) improves the maximum stable time step as well as the level of accuracy and computational cost compared to other schemes in the same class. It is demonstrated that the A‐stability property as well as satisfying “second‐stage order” and stiffly accurate conditions lead the proposed scheme to better performance than existing schemes for the applications examined herein.

Energetic particle precipitation (EPP) and ion chemistry affect the neutral composition of the polar middle atmosphere. For example, production of odd nitrogen and odd hydrogen during strong events can decrease ozone by tens of percent. However, the standard ion chemistry parameterization used in atmospheric models neglects the effects on some important species, such as nitric acid. We present WACCM-D, a variant of the Whole Atmosphere Community Climate Model, which includes a set of lower ionosphere (D-region) chemistry: 307 reactions of 20 positive ions and 21 negative ions. We consider realistic ionization scenarios and compare the WACCM-D results to those from the Sodankyla Ion and Neutral Chemistry (SIC), a state-of-the-art 1-D model of the D-region chemistry. We show that WACCM-D produces well the main characteristics of the D-region ionosphere, as well as the overall proportion of important ion groups, in agreement with SIC. Comparison of ion concentrations shows that the WACCM-D bias is typically within ±10% or less below 70 km. At 70–90 km, when strong altitude gradients in ionization rates and/or ion concentrations exist, the bias can be larger for some groups but is still within tens of percent. Based on the good agreement overall and the fact that part of the differences are caused by different model setups, WACCM-D provides a state-of-the-art global representation of D-region ion chemistry and is therefore expected to improve EPP modeling considerably. These improvements are demonstrated in a companion paper by Andersson et al.

[1] In this paper, we discuss the importance of the surface boundary condition (fixed versus interactive surface temperature) for the long time scale of approach to Radiative-Convective Equilibrium (RCE). Using a simple linearized two-variable model for surface-atmosphere interaction, we derive an analytic expression for sC ,a long climate relaxation time scale that remains well defined and much longer than either mixing time scale of Tompkins and Craig (1998b), even in the limit that the heat capacity of the surface vanishes. We show that the size of sC is an intrinsic property of the coupling between the atmosphere and surface, and not a result of the thermal inertia of the surface alone. When the surface heat capacity is low, sC can be several times longer than expected, due to the effects of moisture on the effective heat capacity of the atmosphere. We also show that the theoretical expression for sC is a good predictor of best fit exponential relaxation time scales in a single-column model with full physics, across a range of surface temperatures and surface heat capacities.

[1] The ECHAM6 atmospheric general circulation model is the atmosphere component of the Max Planck Institute Earth System Model (MPI-ESM) that is used in the Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations. As ECHAM6 has its uppermost layer centered at 0.01 hPa in the upper mesosphere, these simulations offer the opportunity to study the middle atmosphere climate change and its relation to the troposphere on the basis of a very comprehensive set of state-of-the-art model simulations. The goals of this paper are (a) to introduce those new features of ECHAM6 particularly relevant for the middle atmosphere, including external forcing data, and (b) to evaluate the simulated middle atmosphere and describe the simulated response to natural and anthropogenic forcings. New features in ECHAM6 with respect to ECHAM5 include a new short-wave radiation scheme, the option to vary spectral irradiance independent of total solar irradiance, and a latitude-dependent gravity-wave source strength. The description of external forcing data focuses on solar irradiance and ozone. Stratospheric temperature trends simulated with the MPI-ESM for the last decades of the 20th century agree well with observations. The future projections depend strongly on the scenario. Under the high emission scenario RCP8.5, simulated temperatures are locally lower by more than 20 K than preindustrial values. Many of the simulated patterns of the responses to natural forcings as provided by solar variability, volcanic aerosols, and El Nino–Southern Oscillation, largely agree with the observations.

[1] In the paper ‘‘Atmospheric dynamics of Earth-like tidally locked aquaplanets’’ by Timothy M. Merlis and Tapio Schneider (Journal of Advances in Modeling Earth Systems, 2, 13, doi:10.3894/JAMES.2010.2.13, 2010), the shortwave optical depth ts~0:22 is missing from the expression (2.2) for the solar radiative flux. The second paragraph of section 2.1 should read: [2] Absorption of solar radiation in the atmosphere is represented by attenuation of the downward flux of radiation with an optical depth that varies quadratically with pressure, so that S p ð Þ~STOA exp {ts p

[1] This paper examines the ability of recent versions of the Geophysical Fluid Dynamics Laboratory Operational Hurricane Forecast Model (GHM) to reproduce the observed relationship between hurricane intensity and hurricane-induced Sea Surface Temperature (SST) cooling. The analysis was performed by taking a Lagrangian composite of all hurricanes in the North Atlantic from 1998–2009 in observations and 2005–2009 for the GHM. A marked improvement in the intensity-SST relationship for the GHM compared to observations was found between the years 2005 and 2006–2009 due to the introduction of warm-core eddies, a representation of the loop current, and changes to the drag coefficient parameterization for bulk turbulent flux computation. A Conceptual Hurricane Intensity Model illustrates the essential steady-state characteristics of the intensity-SST relationship and is explained by two coupled equations for the atmosphere and ocean. The conceptual model qualitatively matches observations and the 2006–2009 period in the GHM, and presents supporting evidence for the conclusion that weaker upper oceanic thermal stratification in the Gulf of Mexico, caused by the introduction of the loop current and warm core eddies, is crucial to explaining the observed SST-intensity pattern. The diagnostics proposed by the conceptual model offer an independent set of metrics for comparing operational hurricane forecast models to observations.