Effects of grid spacing on high-frequency precipitation variance in coupled high-resolution global ocean-atmosphere models.

High-frequency precipitation variance is calculated in 12 different free-running (non-data-assimilative) coupled high resolution atmosphere-ocean model simulations, an assimilative coupled atmosphere-ocean weather forecast model, and an assimilative reanalysis. The results are compared with results from satellite estimates of precipitation and rain gauge observations. An analysis of irregular sub-daily fluctuations, which was applied by Covey et al. (Geophys Res Lett 45:12514-12522, 2018. 10.1029/2018GL078926) to satellite products and low-resolution climate models, is applied here to rain gauges and higher-resolution models. In contrast to lower-resolution climate simulations, which Covey et al. (2018) found to be lacking with respect to variance in irregular sub-daily fluctuations, the highest-resolution simulations examined here display an irregular sub-daily fluctuation variance that lies closer to that found in satellite products. Most of the simulations used here cannot be analyzed via the Covey et al. (2018) technique, because they do not output precipitation at sub-daily intervals. Thus the remainder of the paper focuses on frequency power spectral density of precipitation and on cumulative distribution functions over time scales (2-100 days) that are still relatively "high-frequency" in the context of climate modeling. Refined atmospheric or oceanic model grid spacing is generally found to increase high-frequency precipitation variance in simulations, approaching the values derived from observations. Mesoscale-eddy-rich ocean simulations significantly increase precipitation variance only when the atmosphere grid spacing is sufficiently fine (< 0.5°). Despite the improvements noted above, all of the simulations examined here suffer from the "drizzle effect", in which precipitation is not temporally intermittent to the extent found in observations.

GEOS-S2S Version 2: The GMAO High Resolution Coupled Model and Assimilation System for Seasonal Prediction.

The Global Modeling and Assimilation Office (GMAO) has recently released a new version of the Goddard Earth Observing System (GEOS) Sub-seasonal to Seasonal prediction (S2S) system, GEOS-S2S-2, that represents a substantial improvement in performance and infrastructure over the previous system. The system is described here in detail, and results are presented from forecasts, climate equillibrium simulations and data assimilation experiments. The climate or equillibrium state of the atmosphere and ocean showed a substantial reduction in bias relative to GEOS-S2S-1. The GEOS-S2S-2 coupled reanalysis also showed substantial improvements, attributed to the assimilation of along-track Absolute Dynamic Topography. The forecast skill on subseasonal scales showed a much-improved prediction of the Madden-Julian Oscillation in GEOS-S2S-2, and on a seasonal scale the tropical Pacific forecasts show substantial improvement in the east and comparable skill to GEOS-S2S-1 in the central Pacific. GEOS-S2S-2 anomaly correlations of both land surface temperature and precipitation were comparable to GEOS-S2S-1, and showed substantially reduced root mean square error of surface temperature. The remaining issues described here are being addressed in the development of GEOS-S2S Version 3, and with that system GMAO will continue its tradition of maintaining a state of the art seasonal prediction system for use in evaluating the impact on seasonal and decadal forecasts of assimilating newly available satellite observations, as well as to evaluate additional sources of predictability in the earth system through the expanded coupling of the earth system model and assimilation components.

The Roles of Climate Change and Climate Variability in the 2017 Atlantic Hurricane Season.

The 2017 Atlantic hurricane season was extremely active with six major hurricanes, the third most on record. The sea-surface temperatures (SSTs) over the eastern Main Development Region (EMDR), where many tropical cyclones (TCs) developed during active months of August/September, were ~0.96 °C above the 1901-2017 average (warmest on record): about ~0.42 °C from a long-term upward trend and the rest (~80%) attributed to the Atlantic Meridional Mode (AMM). The contribution to the SST from the North Atlantic Oscillation (NAO) over the EMDR was a weak warming, while that from El Niño-Southern Oscillation (ENSO) was negligible. Nevertheless, ENSO, the NAO, and the AMM all contributed to favorable wind shear conditions, while the AMM also produced enhanced atmospheric instability. Compared with the strong hurricane years of 2005/2010, the ocean heat content (OHC) during 2017 was larger across the tropics, with higher SST anomalies over the EMDR and Caribbean Sea. On the other hand, the dynamical/thermodynamical atmospheric conditions, while favorable for enhanced TC activity, were less prominent than in 2005/2010 across the tropics. The results suggest that unusually warm SST in the EMDR together with the long fetch of the resulting storms in the presence of record-breaking OHC may be key factors in driving the strong TC activity in 2017.

Errors and improvements in the use of archived meteorological data for chemical transport modeling: an analysis using GEOS-Chem v11-01 driven by GEOS-5 meteorology.

Global simulations of atmospheric chemistry are commonly conducted with off-line chemical transport models (CTMs) driven by archived meteorological data from general circulation models (GCMs). The off-line approach has advantages of simplicity and expediency, but incurs errors due to temporal averaging in the meteorological archive and the inability to reproduce the GCM transport algorithms exactly. The CTM simulation is also often conducted at coarser grid resolution than the parent GCM. Here we investigate this cascade of CTM errors by using 222Rn-210Pb-7Be chemical tracer simulations offline in the GEOS-Chem CTM at rectilinear 0.25° ×0.3125° (≈25 km) and 2° ×2.5° (≈200 km) resolutions, and on-line in the parent GEOS-5 GCM at cubed-sphere c360 (≈25 km) and c48 (≈200 km) horizontal resolutions. The c360 GEOS-5 GCM meteorological archive, updated every 3 hours and remapped to 0.25° ×0.3125°, is the standard operational product generated by the NASA Global Modeling and Assimilation Office (GMAO) and used as input by GEOS-Chem. We find that the GEOS-Chem 222Rn simulation at native 0.25° ×0.3125° resolution is affected by vertical transport errors of up to 20% relative to the GEOS-5 c360 on-line simulation, in part due to loss of transient organized vertical motions in the GCM (resolved convection) that are temporally averaged out in the 3-hour meteorological archive. There is also significant error caused by operational remapping of the meteorological archive from cubed-sphere to rectilinear grid. Decreasing the GEOS-Chem resolution from 0.25°×0.3125° to 2°×2.5° induces further weakening of vertical transport as transient vertical motions are averaged out spatially as well as temporally. The resulting 222Rn concentrations simulated by the coarse-resolution GEOS-Chem are overestimated by up to 40% in surface air relative to the on-line c360 simulations, and underestimated by up to 40% in the upper troposphere, while the tropospheric lifetimes of 210Pb and 7Be against aerosol deposition are affected by 5-10%. The lost vertical transport in the coarse-resolution GEOS-Chem simulation can be partly restored by re-computing the convective mass fluxes at the appropriate resolution to replace the archived convective mass fluxes, and by correcting for bias 20 in spatial averaging of boundary layer mixing depths.

Assessing the Grell-Freitas Convection Parameterization in the NASA GEOS Modeling System.

We implemented and began to evaluate an alternative convection parameterization for the NASA Goddard Earth Observing System (GEOS) general circulation model (GCM). The proposed parameterization follows the mass flux approach with several closures, for equilibrium and nonequilibrium convection, and includes scale and aerosol aware functionalities. Recently, we extended the scheme to a trimodal spectral size distribution of allowed convective plumes to simulate the transition among shallow, congestus, and deep convection regimes. In addition, the inclusion of a new closure for nonequilibrium convection resulted in a substantial gain of realism in the model representation of the diurnal cycle of convection over the land. We demonstrated the scale-dependence functionality with a cascade of global-scale simulations from a nominal horizontal resolution of 50 km down to 6 km. The ability to realistically simulate the diurnal cycle of precipitation over various regions of the earth was verified against several remote sensing-derived intradiurnal precipitation estimates. We extended the model performance evaluation for weather-scale applications by bringing together some available operational short-range weather forecast models and global atmospheric reanalyses. Our results demonstrate that the GEOS GCM with the alternative convective parameterization has good properties and competitive skill in comparison with state-of-the-art observations and numerical simulations.

The impact of SST-forced and unforced teleconnections on 2015/16 El Niño winter precipitation over the western United States.

The factors impacting western U.S. winter precipitation during the 2015/16 El Niño are investigated using the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data, and simulations with the Goddard Earth Observing System version 5 (GEOS-5) atmospheric general circulation model forced with specified sea surface temperatures (SSTs). Results reveal that the simulated response to the tropical Pacific SST associated with the 2015/16 El Niño was to produce wetter than normal conditions over much of the west coast including California - a result at odds with the negative precipitation anomalies observed over much of the Southwestern U.S. It is shown that two factors acted to partly counter the canonical ENSO response in that region. First, a potentially predictable but modest response to the unusually strong and persistent warm SST in the northeastern Pacific decreased precipitation in the Southwestern U.S. by increasing sea level pressure, driving anticyclonic circulation and atmospheric descent, and reducing moisture transport into that region. Second, large-scale unforced (by SST) components of atmospheric variability (consisting of the leading modes of unpredictable intra-ensemble variability) resembling the positive phase of the North Atlantic Oscillation and Arctic Oscillation are found to be an important contributor to the drying over the western U.S. While a statistical reconstruction of the precipitation from our simulations that account for internal atmospheric variability does much to close the gap between the ensemble mean and observed precipitation in the Southwestern U.S., some differences remain, indicating that model error is also playing a role.

Direct estimation of the global distribution of vertical velocity within cirrus clouds.

Cirrus clouds determine the radiative balance of the upper troposphere and the transport of water vapor across the tropopause. The representation of vertical wind velocity, W, in atmospheric models constitutes the largest source of uncertainty in the calculation of the cirrus formation rate. Using global atmospheric simulations with a spatial resolution of 7 km we obtain for the first time a direct estimate of the distribution of W at the scale relevant for cirrus formation, validated against long-term observations at two different ground sites. The standard deviation in W, σ w , varies widely over the globe with the highest values resulting from orographic uplift and convection, and the lowest occurring in the Arctic. Globally about 90% of the simulated σ w values are below 0.1 m s-1 and about one in 104 cloud formation events occur in environments with σ w > 0.8 m s-1. Combining our estimate with reanalysis products and an advanced cloud formation scheme results in lower homogeneous ice nucleation frequency than previously reported, and a decreasing average ice crystal concentration with decreasing temperature. These features are in agreement with observations and suggest that the correct parameterization of σ w is critical to simulate realistic cirrus properties.

Chemical Mechanisms and Their Applications in the Goddard Earth Observing System (GEOS) Earth System Model.

NASA's Goddard Earth Observing System (GEOS) Earth System Model (ESM) is a modular, general circulation model (GCM), and data assimilation system (DAS) that is used to simulate and study the coupled dynamics, physics, chemistry, and biology of our planet. GEOS is developed by the Global Modeling and Assimilation Office (GMAO) at NASA Goddard Space Flight Center. It generates near-real-time analyzed data products, reanalyses, and weather and seasonal forecasts to support research targeted to understanding interactions among Earth System processes. For chemistry, our efforts are focused on ozone and its influence on the state of the atmosphere and oceans, and on trace gas data assimilation and global forecasting at mesoscale discretization. Several chemistry and aerosol modules are coupled to the GCM, which enables GEOS to address topics pertinent to NASA's Earth Science Mission. This paper describes the atmospheric chemistry components of GEOS and provides an overview of its Earth System Modeling Framework (ESMF)-based software infrastructure, which promotes a rich spectrum of feedbacks that influence circulation and climate, and impact human and ecosystem health. We detail how GEOS allows model users to select chemical mechanisms and emission scenarios at run time, establish the extent to which the aerosol and chemical components communicate, and decide whether either or both influence the radiative transfer calculations. A variety of resolutions facilitates research on spatial and temporal scales relevant to problems ranging from hourly changes in air quality to trace gas trends in a changing climate. Samples of recent GEOS chemistry applications are provided.

The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2).

The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) is the latest atmospheric reanalysis of the modern satellite era produced by NASA's Global Modeling and Assimilation Office (GMAO). MERRA-2 assimilates observation types not available to its predecessor, MERRA, and includes updates to the Goddard Earth Observing System (GEOS) model and analysis scheme so as to provide a viable ongoing climate analysis beyond MERRA's terminus. While addressing known limitations of MERRA, MERRA-2 is also intended to be a development milestone for a future integrated Earth system analysis (IESA) currently under development at GMAO. This paper provides an overview of the MERRA-2 system and various performance metrics. Among the advances in MERRA-2 relevant to IESA are the assimilation of aerosol observations, several improvements to the representation of the stratosphere including ozone, and improved representations of cryospheric processes. Other improvements in the quality of MERRA-2 compared with MERRA include the reduction of some spurious trends and jumps related to changes in the observing system, and reduced biases and imbalances in aspects of the water cycle. Remaining deficiencies are also identified. Production of MERRA-2 began in June 2014 in four processing streams, and converged to a single near-real time stream in mid 2015. MERRA-2 products are accessible online through the NASA Goddard Earth Sciences Data Information Services Center (GES DISC).

Practice and philosophy of climate model tuning across six U.S. modeling centers.

Model calibration (or "tuning") is a necessary part of developing and testing coupled ocean-atmosphere climate models regardless of their main scientific purpose. There is an increasing recognition that this process needs to become more transparent for both users of climate model output and other developers. Knowing how and why climate models are tuned and which targets are used is essential to avoiding possible misattributions of skillful predictions to data accommodation and vice versa. This paper describes the approach and practice of model tuning for the six major U.S. climate modeling centers. While details differ among groups in terms of scientific missions, tuning targets and tunable parameters, there is a core commonality of approaches. However, practices differ significantly on some key aspects, in particular, in the use of initialized forecast analyses as a tool, the explicit use of the historical transient record, and the use of the present day radiative imbalance vs. the implied balance in the pre-industrial as a target.

Large-scale Controls on Atlantic Tropical Cyclone Activity on Seasonal Time Scales.

Interannual variations in seasonal tropical cyclone (TC) activity (e.g., genesis frequency and location, track pattern, and landfall) over the Atlantic are explored by employing observationally-constrained simulations with the NASA Goddard Earth Observing System version (GEOS-5) atmospheric general circulation model. The climate modes investigated are El Niño-Southern Oscillation (ENSO), the North Atlantic Oscillation (NAO), and the Atlantic Meridional Mode (AMM). The results show that the NAO and AMM can strongly modify and even oppose the well-known ENSO impacts, like in 2005, when a strong positive AMM (associated with warm SSTs and a negative SLP anomaly over the western tropical Atlantic), led to a very active TC season with enhanced TC genesis over the Caribbean Sea and a number of landfalls over North America, under a neutral ENSO condition. On the other end, the weak TC activity during 2013 (characterized by weak negative Niño index) appears caused by a NAO-induced positive SLP anomaly with enhanced vertical wind shear over the tropical North Atlantic. During 2010, the combined impact of the three modes produced positive SST anomalies across the entire low- latitudinal Atlantic and a weaker subtropical high, leading to more early recurvers and thus fewer landfalls despite enhanced TC genesis. The study provides evidence that TC number and track are very sensitive to the relative phases and intensities of these three modes, and not just to ENSO alone. Examination of seasonal predictability reveals that predictive skill of the three modes is limited over tropics to sub-tropics, with the AMM having the highest predictability over the North Atlantic, followed by ENSO and NAO.

Structure and Dynamics of the Quasi-Biennial Oscillation in MERRA-2.

The structure, dynamics, and ozone signal of the Quasi-Biennial Oscillation produced by the 35-year NASA MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications) reanalysis are examined based on monthly mean output. Along with the analysis of the QBO in assimilation winds and ozone, the QBO forcings created by assimilated observations, dynamics, parameterized gravity wave drag, and ozone chemistry parameterization are examined and compared with the original MERRA system. Results show that the MERRA-2 reanalysis produces a realistic QBO in the zonal winds, mean meridional circulation, and ozone over the 1980-2015 time period. In particular, the MERRA-2 zonal winds show improved representation of the QBO 50 hPa westerly phase amplitude at Singapore when compared to MERRA. The use of limb ozone observations creates improved vertical structure and realistic downward propagation of the ozone QBO signal during times when the MLS ozone limb observations are available (October 2004 to present). The increased equatorial GWD in MERRA-2 has reduced the zonal wind data analysis contribution compared to MERRA so that the QBO mean meridional circulation can be expected to be more physically forced and therefore more physically consistent. This can be important for applications in which MERRA-2 winds are used to drive transport experiments.

Frequency and Impact of Summertime Stratospheric Intrusions over Maryland during DISCOVER-AQ (2011): New Evidence from NASA's GEOS-5 Simulations.

Aircraft observations and ozonesonde profiles collected on July 14 and 27, 2011, during the Maryland month-long DISCOVER-AQ campaign, indicate the presence of stratospheric air just above the planetary boundary layer (PBL). This raises the question of whether summer stratospheric intrusions (SIs) elevate surface ozone levels and to what degree they influence background ozone levels and contribute to ozone production. We used idealized stratospheric air tracers, along with observations, to determine the frequency and extent of SIs in Maryland during July 2011. On 4 of 14 flight days, SIs were detected in layers that the aircraft encountered above the PBL from the coincidence of enhanced ozone, moderate CO, and low moisture. Satellite observations of lower tropospheric humidity confirmed the occurrence of synoptic scale influence of SIs as do simulations with the GEOS-5 Atmospheric General Circulation Model. The evolution of GEOS-5 stratospheric air tracers agree with the timing and location of observed stratospheric influence and indicate that more than 50% of air in SI layers above the PBL had resided in the stratosphere within the previous 14 days. Despite having a strong influence in the lower free troposphere, these events did not significantly affect surface ozone, which remained low on intrusion days. The model indicates similar frequencies of stratospheric influence during all summers from 2009-2013. GEOS-5 results suggest that, over Maryland, the strong inversion capping the summer PBL limits downward mixing of stratospheric air during much of the day, helping to preserve low surface ozone associated with frontal passages that precede SIs.